(* Generated by JWS Online *)

(* This is an experimental feature of JWS Online. Please report any mistakes.*)

(* Note that the following notable SBML entities or features are not supported in notebook outputyet: *)
    (* Events *)
    (* Constraints *)
    (* Units and UnitDefinitions *)
    (* AlgebraicRules *)
    (* conversionFactors *)

variables = {
B\[LetterSpace]alan[t], B\[LetterSpace]ffa[t], B\[LetterSpace]gluc[t], B\[LetterSpace]glucgn[t], B\[LetterSpace]ins[t], B\[LetterSpace]ket[t], B\[LetterSpace]lac[t], Cvar[t], F\[LetterSpace]TG[t], F\[LetterSpace]acyl[t], F\[LetterSpace]g6p[t], GPa[t], GSa[t], M\[LetterSpace]g6p[t], M\[LetterSpace]glycgn[t], M\[LetterSpace]ket[t], M\[LetterSpace]lac[t], M\[LetterSpace]pyr[t], PKa[t], PP1[t], PP1\[LetterSpace]GPa[t], R2C2[t], R2\[LetterSpace]C\[LetterSpace]cAMP2[t], R2\[LetterSpace]cAMP4[t], aK[t], acet\[LetterSpace]c[t], acet\[LetterSpace]m[t], alan[t], cAMP[t], citrate[t], g6p[t], gluc[t], glutamate[t], glycgn[t], ket[t], lac[t], malate[t], malonyl[t], oa\[LetterSpace]c[t], oa\[LetterSpace]m[t], palm[t], palmCoA[t], pep[t], pyr[t]
};

initialValues = {
B\[LetterSpace]alan[0] == 0.1, B\[LetterSpace]ffa[0] == 0.3, B\[LetterSpace]gluc[0] == 4, B\[LetterSpace]glucgn[0] == 0, B\[LetterSpace]ins[0] == 0, B\[LetterSpace]ket[0] == 0.08, B\[LetterSpace]lac[0] == 0.2, Cvar[0] == 0, F\[LetterSpace]TG[0] == 0, F\[LetterSpace]acyl[0] == 0, F\[LetterSpace]g6p[0] == 0, GPa[0] == 0.035, GSa[0] == 0.0015, M\[LetterSpace]g6p[0] == 8, M\[LetterSpace]glycgn[0] == 100, M\[LetterSpace]ket[0] == 0.1, M\[LetterSpace]lac[0] == 0.2, M\[LetterSpace]pyr[0] == 0.2, PKa[0] == 0, PP1[0] == 0.00025, PP1\[LetterSpace]GPa[0] == 0, R2C2[0] == 0.00025, R2\[LetterSpace]C\[LetterSpace]cAMP2[0] == 0, R2\[LetterSpace]cAMP4[0] == 0, aK[0] == 1, acet\[LetterSpace]c[0] == 6, acet\[LetterSpace]m[0] == 7, alan[0] == 0.1, cAMP[0] == 0, citrate[0] == 0.5, g6p[0] == 0.2, gluc[0] == 5, glutamate[0] == 0.1, glycgn[0] == 10, ket[0] == 0.08, lac[0] == 1, malate[0] == 3, malonyl[0] == 0.3, oa\[LetterSpace]c[0] == 0.0004, oa\[LetterSpace]m[0] == 0.0002, palm[0] == 0.3, palmCoA[0] == 0.3, pep[0] == 0.001, pyr[0] == 3
};

rates = {
Jg10, Jg11, Jg3, Jg4, Jg5, Jg6, Jg7, Jg8, Jg9, v\[LetterSpace]F1, v\[LetterSpace]F3, v\[LetterSpace]F4, v\[LetterSpace]F5, v\[LetterSpace]Ggluc, v\[LetterSpace]Igluc, v\[LetterSpace]L10, v\[LetterSpace]L11f, v\[LetterSpace]L11r, v\[LetterSpace]L12, v\[LetterSpace]L13, v\[LetterSpace]L14, v\[LetterSpace]L15, v\[LetterSpace]L16, v\[LetterSpace]L17, v\[LetterSpace]L18, v\[LetterSpace]L19, v\[LetterSpace]L1f, v\[LetterSpace]L1r, v\[LetterSpace]L20, v\[LetterSpace]L21f, v\[LetterSpace]L21r, v\[LetterSpace]L22, v\[LetterSpace]L2f, v\[LetterSpace]L2r, v\[LetterSpace]L3f, v\[LetterSpace]L3r, v\[LetterSpace]L4, v\[LetterSpace]L5f, v\[LetterSpace]L5r, v\[LetterSpace]L6, v\[LetterSpace]L7, v\[LetterSpace]L8, v\[LetterSpace]L9, v\[LetterSpace]S1f, v\[LetterSpace]S1r, v\[LetterSpace]S2, v\[LetterSpace]S3f, v\[LetterSpace]S3r, v\[LetterSpace]S4f, v\[LetterSpace]S4r, v\[LetterSpace]c1cAMP, v\[LetterSpace]c2cAMP, v\[LetterSpace]d\[LetterSpace]Balan, v\[LetterSpace]d\[LetterSpace]Bffa, v\[LetterSpace]d\[LetterSpace]Bgluc, v\[LetterSpace]d\[LetterSpace]Bglucgn, v\[LetterSpace]d\[LetterSpace]Bins, v\[LetterSpace]d\[LetterSpace]Bket, v\[LetterSpace]d\[LetterSpace]Blac, v\[LetterSpace]feed, v\[LetterSpace]glucgn, v\[LetterSpace]ins, v\[LetterSpace]s\[LetterSpace]dket, v\[LetterSpace]tF1, v\[LetterSpace]tF3, v\[LetterSpace]tL1, v\[LetterSpace]tL2, v\[LetterSpace]tL3, v\[LetterSpace]tL5, v\[LetterSpace]tL6, v\[LetterSpace]tS1, v\[LetterSpace]tS2, v\[LetterSpace]tS3, v\[LetterSpace]tS4
};

rateEquations = {
Jg10 -> kgc1*cAMP[t]^2*R2C2[t] - k\[LetterSpace]gc1*Cvar[t]*R2\[LetterSpace]C\[LetterSpace]cAMP2[t], Jg11 -> -(k\[LetterSpace]gc2*Cvar[t]*R2\[LetterSpace]cAMP4[t]) + kgc2*cAMP[t]^2*R2\[LetterSpace]C\[LetterSpace]cAMP2[t], Jg3 -> (kg3*Cvar[t]*(kt - PKa[t]))/(kmg3 + kt - PKa[t]), Jg4 -> (kg4*PKa[t]*(PP1[t] + PP1\[LetterSpace]GPa[t]))/(kmg4 + PKa[t]), Jg5 -> (kg5*(pt - GPa[t])*PKa[t])/(kmg5s + pt - GPa[t]), Jg6 -> (kg6*GPa[t]*(PP1[t] + PP1\[LetterSpace]GPa[t]))/(kmg6s + GPa[t]), Jg7 -> (kg7*GSa[t]*(Cvar[t] + PKa[t]))/(kmg7s + GSa[t]), Jg8 -> (kg8*(st - GSa[t])*PP1[t])/(kmg8s + st - GSa[t]), Jg9 -> ka*GPa[t]*PP1[t] - k\[LetterSpace]a*PP1\[LetterSpace]GPa[t], v\[LetterSpace]F1 -> k\[LetterSpace]f1*(1 + B\[LetterSpace]ins[t]^ep14/(k\[LetterSpace]Dins^ep14 + B\[LetterSpace]ins[t]^ep14))*F\[LetterSpace]g6p[t], v\[LetterSpace]F3 -> k\[LetterSpace]f3*(1 + B\[LetterSpace]ins[t]^ep15/(k\[LetterSpace]Dins^ep15 + B\[LetterSpace]ins[t]^ep15))*F\[LetterSpace]acyl[t], v\[LetterSpace]F4 -> (k\[LetterSpace]Dins^en12*k\[LetterSpace]f4*F\[LetterSpace]TG[t])/(k\[LetterSpace]Dins^en12 + B\[LetterSpace]ins[t]^en12), v\[LetterSpace]F5 -> F\[LetterSpace]ffa^3*k\[LetterSpace]f5, v\[LetterSpace]Ggluc -> (k1glucgn*B\[LetterSpace]gluc[t]^ng)/(k\[LetterSpace]mGlgn^ng + B\[LetterSpace]gluc[t]^ng), v\[LetterSpace]Igluc -> (k1ins*B\[LetterSpace]gluc[t]^ni)/(k\[LetterSpace]mIns^ni + B\[LetterSpace]gluc[t]^ni), v\[LetterSpace]L10 -> k\[LetterSpace]L10*aK[t], v\[LetterSpace]L11f -> k\[LetterSpace]L11f*(1 + B\[LetterSpace]glucgn[t]^ep5/(k\[LetterSpace]Dglucgn^ep5 + B\[LetterSpace]glucgn[t]^ep5))*oa\[LetterSpace]m[t], v\[LetterSpace]L11r -> k\[LetterSpace]L11r*malate[t], v\[LetterSpace]L12 -> (k\[LetterSpace]Dins^en4*k\[LetterSpace]L12*(1 + B\[LetterSpace]glucgn[t]^ep6/(k\[LetterSpace]Dglucgn^ep6 + B\[LetterSpace]glucgn[t]^ep6))*malate[t])/(k\[LetterSpace]Dins^en4 + B\[LetterSpace]ins[t]^en4), v\[LetterSpace]L13 -> (k\[LetterSpace]Dins^en5*k\[LetterSpace]L13*(1 + B\[LetterSpace]glucgn[t]^ep7/(k\[LetterSpace]Dglucgn^ep7 + B\[LetterSpace]glucgn[t]^ep7))*oa\[LetterSpace]c[t])/(k\[LetterSpace]Dins^en5 + B\[LetterSpace]ins[t]^en5), v\[LetterSpace]L14 -> (k\[LetterSpace]Dglucgn^en7*k\[LetterSpace]L14*(1 + B\[LetterSpace]ins[t]^ep10/(k\[LetterSpace]Dins^ep10 + B\[LetterSpace]ins[t]^ep10))*citrate[t])/(k\[LetterSpace]Dglucgn^en7 + B\[LetterSpace]glucgn[t]^en7), v\[LetterSpace]L15 -> (k\[LetterSpace]i5*k\[LetterSpace]L15*acet\[LetterSpace]c[t]*(1 + (p1*citrate[t])/(k\[LetterSpace]p1 + citrate[t])))/(k\[LetterSpace]i5 + palmCoA[t]), v\[LetterSpace]L16 -> k\[LetterSpace]L16*acet\[LetterSpace]c[t]*malonyl[t]^7, v\[LetterSpace]L17 -> (k\[LetterSpace]i1*k\[LetterSpace]L17*palm[t])/(k\[LetterSpace]i1 + malonyl[t]), v\[LetterSpace]L18 -> (k\[LetterSpace]i2*k\[LetterSpace]L18*palmCoA[t])/(k\[LetterSpace]i2 + malonyl[t]), v\[LetterSpace]L19 -> (k\[LetterSpace]Dins^en8*k\[LetterSpace]L19*acet\[LetterSpace]m[t]^2*(1 + cAMP[t]^ep11/(k\[LetterSpace]dcAMP^ep11 + cAMP[t]^ep11)))/(k\[LetterSpace]Dins^en8 + B\[LetterSpace]ins[t]^en8), v\[LetterSpace]L1f -> (k\[LetterSpace]L1f*(1 + B\[LetterSpace]ins[t]^ep1/(k\[LetterSpace]Dins^ep1 + B\[LetterSpace]ins[t]^ep1))*gluc[t])/(k\[LetterSpace]mL1f + gluc[t]), v\[LetterSpace]L1r -> (k\[LetterSpace]L1r*(1 + B\[LetterSpace]glucgn[t]^ep9/(k\[LetterSpace]Dglucgn^ep9 + B\[LetterSpace]glucgn[t]^ep9))*g6p[t])/(k\[LetterSpace]mL1r + g6p[t]), v\[LetterSpace]L20 -> (k\[LetterSpace]Dins2^en3*k\[LetterSpace]L20*alan[t])/(k\[LetterSpace]Dins2^en3 + B\[LetterSpace]ins[t]^en3), v\[LetterSpace]L21f -> (k\[LetterSpace]Dins^en3*k\[LetterSpace]L21f*glutamate[t]*pyr[t])/((k\[LetterSpace]Dins^en3 + B\[LetterSpace]ins[t]^en3)*(k\[LetterSpace]mL21g + glutamate[t])*(k\[LetterSpace]mL21p + pyr[t])), v\[LetterSpace]L21r -> (k\[LetterSpace]Dins^en3*k\[LetterSpace]L21r*acet\[LetterSpace]m[t]*aK[t]*alan[t])/((k\[LetterSpace]mL21k + aK[t])*(k\[LetterSpace]mL21a + alan[t])*(k\[LetterSpace]Dins^en3 + B\[LetterSpace]ins[t]^en3)), v\[LetterSpace]L22 -> k\[LetterSpace]L22*glutamate[t], v\[LetterSpace]L2f -> (k\[LetterSpace]L2f*g6p[t]*GSa[t])/(k\[LetterSpace]mL2f + g6p[t]), v\[LetterSpace]L2r -> (k\[LetterSpace]L2r*glycgn[t]*GPa[t])/(k\[LetterSpace]mL2r + glycgn[t]), v\[LetterSpace]L3f -> (k\[LetterSpace]Dglucgn^en1*k\[LetterSpace]L3f*(1 + B\[LetterSpace]ins[t]^ep2/(k\[LetterSpace]Dins^ep2 + B\[LetterSpace]ins[t]^ep2))*g6p[t])/((k\[LetterSpace]Dglucgn^en1 + B\[LetterSpace]glucgn[t]^en1)*(k\[LetterSpace]mL3f + g6p[t])), v\[LetterSpace]L3r -> (k\[LetterSpace]Dins2^en6*k\[LetterSpace]L3r*(1 + B\[LetterSpace]glucgn[t]^ep8/(k\[LetterSpace]Dglucgn^ep8 + B\[LetterSpace]glucgn[t]^ep8))*pep[t])/((k\[LetterSpace]Dins2^en6 + B\[LetterSpace]ins[t]^en6)*(k\[LetterSpace]mL3r + pep[t])), v\[LetterSpace]L4 -> (k\[LetterSpace]Dglucgn^en2*k\[LetterSpace]i13*k\[LetterSpace]L4*(1 + B\[LetterSpace]ins[t]^ep3/(k\[LetterSpace]Dins^ep3 + B\[LetterSpace]ins[t]^ep3))*pep[t])/((k\[LetterSpace]i13 + alan[t])*(k\[LetterSpace]Dglucgn^en2 + B\[LetterSpace]glucgn[t]^en2)*(k\[LetterSpace]mL4 + pep[t])), v\[LetterSpace]L5f -> (k\[LetterSpace]L5f*pyr[t])/(k\[LetterSpace]mL5f + pyr[t]), v\[LetterSpace]L5r -> (k\[LetterSpace]L5r*lac[t])/(k\[LetterSpace]mL5r + lac[t]), v\[LetterSpace]L6 -> (k\[LetterSpace]L6*(1 + (p2*acet\[LetterSpace]m[t])/(k\[LetterSpace]p2 + acet\[LetterSpace]m[t]))*(1 + B\[LetterSpace]glucgn[t]^ep4/(k\[LetterSpace]Dglucgn^ep4 + B\[LetterSpace]glucgn[t]^ep4))*pyr[t])/(k\[LetterSpace]mL6 + pyr[t]), v\[LetterSpace]L7 -> (k\[LetterSpace]i8*k\[LetterSpace]L7*pyr[t])/((k\[LetterSpace]i8 + acet\[LetterSpace]m[t])*(k\[LetterSpace]mL7 + pyr[t])), v\[LetterSpace]L8 -> (k\[LetterSpace]i4*k\[LetterSpace]L8*acet\[LetterSpace]m[t]*oa\[LetterSpace]m[t])/(k\[LetterSpace]i4 + palmCoA[t]), v\[LetterSpace]L9 -> k\[LetterSpace]L9*citrate[t], v\[LetterSpace]S1f -> k\[LetterSpace]s1f*M\[LetterSpace]g6p[t], v\[LetterSpace]S1r -> k\[LetterSpace]s1r*M\[LetterSpace]glycgn[t], v\[LetterSpace]S2 -> k\[LetterSpace]s2*M\[LetterSpace]g6p[t], v\[LetterSpace]S3f -> k\[LetterSpace]s3f*M\[LetterSpace]pyr[t], v\[LetterSpace]S3r -> k\[LetterSpace]s3r*M\[LetterSpace]lac[t], v\[LetterSpace]S4f -> k\[LetterSpace]s4f*M\[LetterSpace]pyr[t], v\[LetterSpace]S4r -> k\[LetterSpace]s4r*M\[LetterSpace]alan*M\[LetterSpace]ket[t], v\[LetterSpace]c1cAMP -> (kc1*B\[LetterSpace]glucgn[t]^ng)/(kcm1^ng + B\[LetterSpace]glucgn[t]^ng), v\[LetterSpace]c2cAMP -> (kc2*B\[LetterSpace]ins[t]^ni*cAMP[t])/(kcm2^ni + B\[LetterSpace]ins[t]^ni), v\[LetterSpace]d\[LetterSpace]Balan -> kd\[LetterSpace]Balan*B\[LetterSpace]alan[t], v\[LetterSpace]d\[LetterSpace]Bffa -> kd\[LetterSpace]Bffa*B\[LetterSpace]ffa[t], v\[LetterSpace]d\[LetterSpace]Bgluc -> kd\[LetterSpace]Bgluc*B\[LetterSpace]gluc[t], v\[LetterSpace]d\[LetterSpace]Bglucgn -> kd\[LetterSpace]Bglucgn*B\[LetterSpace]glucgn[t], v\[LetterSpace]d\[LetterSpace]Bins -> kd\[LetterSpace]Bins*B\[LetterSpace]ins[t], v\[LetterSpace]d\[LetterSpace]Bket -> kd\[LetterSpace]Bket*B\[LetterSpace]ket[t], v\[LetterSpace]d\[LetterSpace]Blac -> kd\[LetterSpace]Blac*B\[LetterSpace]lac[t], v\[LetterSpace]feed -> Piecewise[{{kf, t < 20000}, {kresume, t > tresume}}, (1.*kf)/E^(0.0006157521601035995*kf^2*(-20000 + t)^2)], v\[LetterSpace]glucgn -> k\[LetterSpace]glucgn, v\[LetterSpace]ins -> k\[LetterSpace]ins, v\[LetterSpace]s\[LetterSpace]dket -> ks\[LetterSpace]dket*M\[LetterSpace]ket[t], v\[LetterSpace]tF1 -> k\[LetterSpace]tF1*B\[LetterSpace]gluc[t]*(1 + B\[LetterSpace]ins[t]^ep12/(k\[LetterSpace]Dins^ep12 + B\[LetterSpace]ins[t]^ep12)), v\[LetterSpace]tF3 -> (F\[LetterSpace]ffa^8*k\[LetterSpace]Dins2^en9*k\[LetterSpace]tF3)/(k\[LetterSpace]Dins2^en9 + B\[LetterSpace]ins[t]^en9), v\[LetterSpace]tL1 -> k\[LetterSpace]tL1*(B\[LetterSpace]gluc[t] - gluc[t]), v\[LetterSpace]tL2 -> k\[LetterSpace]tL2*B\[LetterSpace]lac[t], v\[LetterSpace]tL3 -> k\[LetterSpace]tL3*ket[t], v\[LetterSpace]tL5 -> k\[LetterSpace]tL5*B\[LetterSpace]ffa[t], v\[LetterSpace]tL6 -> (k\[LetterSpace]Dins2^en10*k\[LetterSpace]tL6*B\[LetterSpace]alan[t])/(k\[LetterSpace]Dins2^en10 + B\[LetterSpace]ins[t]^en10), v\[LetterSpace]tS1 -> k\[LetterSpace]tS1*B\[LetterSpace]gluc[t]*(1 + B\[LetterSpace]ins[t]^ep13/(k\[LetterSpace]Dins^ep13 + B\[LetterSpace]ins[t]^ep13)), v\[LetterSpace]tS2 -> k\[LetterSpace]tS2*B\[LetterSpace]ket[t], v\[LetterSpace]tS3 -> k\[LetterSpace]tS3*M\[LetterSpace]lac[t], v\[LetterSpace]tS4 -> (k\[LetterSpace]Dins2^en11*k\[LetterSpace]tS4*M\[LetterSpace]alan)/(k\[LetterSpace]Dins2^en11 + B\[LetterSpace]ins[t]^en11)
};

parameters = {
Gmax -> 5*^-08, Gmin -> 3*^-08, Imax -> 1.3*^-06, Imin -> 7*^-07, PP1t -> 0.00025, c0 -> 5.0, capkt -> 0.00025, cmax -> 0.0032, cmin -> 2*^-06, en1 -> 10.0, en10 -> 10.0, en11 -> 10.0, en12 -> 10.0, en2 -> 10.0, en3 -> 10.0, en4 -> 10.0, en5 -> 10.0, en6 -> 10.0, en7 -> 10.0, en8 -> 10.0, en9 -> 20.0, ep1 -> 10.0, ep10 -> 10.0, ep11 -> 10.0, ep12 -> 10.0, ep13 -> 10.0, ep14 -> 10.0, ep15 -> 10.0, ep2 -> 10.0, ep3 -> 10.0, ep4 -> 10.0, ep5 -> 10.0, ep6 -> 10.0, ep7 -> 10.0, ep8 -> 10.0, ep9 -> 10.0, index -> 5.0, k11 -> 4.3*^-05, k22 -> 0.0007, k\[LetterSpace]Dglucgn -> 4*^-08, k\[LetterSpace]Dins -> 1*^-06, k\[LetterSpace]Dins2 -> 7.5*^-07, k\[LetterSpace]L10 -> 0.1, k\[LetterSpace]L11f -> 0.6, k\[LetterSpace]L11r -> 0.01, k\[LetterSpace]L12 -> 0.6, k\[LetterSpace]L13 -> 0.5, k\[LetterSpace]L14 -> 0.01, k\[LetterSpace]L15 -> 0.01, k\[LetterSpace]L16 -> 0.01, k\[LetterSpace]L17 -> 0.01, k\[LetterSpace]L18 -> 0.01, k\[LetterSpace]L19 -> 0.01, k\[LetterSpace]L1f -> 3.0, k\[LetterSpace]L1r -> 4.0, k\[LetterSpace]L20 -> 0.2, k\[LetterSpace]L21f -> 0.001, k\[LetterSpace]L21r -> 0.2, k\[LetterSpace]L22 -> 0.01, k\[LetterSpace]L2f -> 200.0, k\[LetterSpace]L2r -> 20.0, k\[LetterSpace]L3f -> 0.1, k\[LetterSpace]L3r -> 0.3, k\[LetterSpace]L4 -> 2.0, k\[LetterSpace]L5f -> 0.0, k\[LetterSpace]L5r -> 0.1, k\[LetterSpace]L6 -> 1.0, k\[LetterSpace]L7 -> 1.0, k\[LetterSpace]L8 -> 0.1, k\[LetterSpace]L9 -> 0.1, k\[LetterSpace]a -> 60.0, k\[LetterSpace]dcAMP -> 3.16227*^-06, k\[LetterSpace]f1 -> 0.1, k\[LetterSpace]f3 -> 0.2, k\[LetterSpace]f4 -> 0.1, k\[LetterSpace]f5 -> 0.1, k\[LetterSpace]gc1 -> 60000.0, k\[LetterSpace]gc2 -> 60000.0, k\[LetterSpace]i1 -> 0.1, k\[LetterSpace]i13 -> 2.0, k\[LetterSpace]i2 -> 1.0, k\[LetterSpace]i4 -> 3.0, k\[LetterSpace]i5 -> 1.0, k\[LetterSpace]i8 -> 2.0, k\[LetterSpace]mGlgn -> 8.0, k\[LetterSpace]mIns -> 8.0, k\[LetterSpace]mL1f -> 7.7, k\[LetterSpace]mL1r -> 1.3, k\[LetterSpace]mL21a -> 21.0, k\[LetterSpace]mL21g -> 4.3, k\[LetterSpace]mL21k -> 0.22, k\[LetterSpace]mL21p -> 0.4, k\[LetterSpace]mL2f -> 0.57, k\[LetterSpace]mL2r -> 1.4, k\[LetterSpace]mL3f -> 0.01, k\[LetterSpace]mL3r -> 0.0034, k\[LetterSpace]mL4 -> 0.18, k\[LetterSpace]mL5f -> 0.03, k\[LetterSpace]mL5r -> 0.8, k\[LetterSpace]mL6 -> 0.22, k\[LetterSpace]mL7 -> 0.0204, k\[LetterSpace]p1 -> 0.5, k\[LetterSpace]p2 -> 2.0, k\[LetterSpace]s1f -> 0.0, k\[LetterSpace]s1r -> 0.0, k\[LetterSpace]s2 -> 0.02, k\[LetterSpace]s3f -> 0.01, k\[LetterSpace]s3r -> 0.01, k\[LetterSpace]s4f -> 0.07, k\[LetterSpace]s4r -> 0.0, k\[LetterSpace]tF1 -> 0.01, k\[LetterSpace]tF3 -> 0.008, k\[LetterSpace]tL1 -> 100.0, k\[LetterSpace]tL2 -> 0.1, k\[LetterSpace]tL3 -> 0.1, k\[LetterSpace]tL5 -> 0.1, k\[LetterSpace]tL6 -> 1.0, k\[LetterSpace]tS1 -> 0.01, k\[LetterSpace]tS2 -> 1.0, k\[LetterSpace]tS3 -> 0.01, k\[LetterSpace]tS4 -> 3.0, kc1 -> 1.0, kcm1 -> 4*^-08, kcm2 -> 1*^-06, kd\[LetterSpace]Balan -> 0.015, kd\[LetterSpace]Bffa -> 0.015, kd\[LetterSpace]Bgluc -> 0.015, kd\[LetterSpace]Bglucgn -> 0.015, kd\[LetterSpace]Bins -> 0.015, kd\[LetterSpace]Bket -> 0.015, kd\[LetterSpace]Blac -> 0.015, kf -> 0.5, kg2 -> 0.5, kg3 -> 1200.0, kg4 -> 300.0, kg5 -> 1200.0, kg6 -> 300.0, kg7 -> 1200.0, kg8 -> 300.0, kgi -> 10.0, kmg3 -> 0.0004, kmg4 -> 0.0011, kmg5 -> 0.01, kmg6 -> 0.005, kmg7 -> 0.015, kmg8 -> 0.00012, kresume -> 0.6, ks\[LetterSpace]dket -> 0.01, kt -> 0.0025, ng -> 10.0, ni -> 10.0, p1 -> 1.0, p2 -> 1.0, pt -> 0.07, s1 -> 100.0, s2 -> 0.001, st -> 0.003, tresume -> 20250.0, F\[LetterSpace]ffa -> 2, M\[LetterSpace]alan -> 0.4, PKb -> 0.0025, blood -> 1.0, fat -> 1.0, glycogen\[LetterSpace]regulation -> 1.0, liver -> 1.0, muscle -> 1.0
};

assignments = {
k1ins -> (Imax - Imin)*kd\[LetterSpace]Bins, kgc1 -> k\[LetterSpace]gc1/k11, kgc2 -> k\[LetterSpace]gc2/k22, kmg5s -> kmg5*(1 + (s1*g6p[t])/kg2), kmg6s -> kmg6/(1 + (s2*gluc[t])/kgi), kmg7s -> kmg7*(1 + (s1*g6p[t])/kg2), kmg8s -> kmg8/(1 + (s1*g6p[t])/kg2), kc2 -> 316227.7660168379*kc1, ka -> k\[LetterSpace]a/kd, k\[LetterSpace]glucgn -> Gmax*kd\[LetterSpace]Bglucgn, k1glucgn -> (Gmax - Gmin)*kd\[LetterSpace]Bglucgn, kd -> cmin + (c0^index*(cmax - cmin))/(c0^index + glycgn[t]^index), k\[LetterSpace]ins -> Imin*kd\[LetterSpace]Bins, GPaPlot -> 14.285714285714285*GPa[t], GSaPlot -> 333.3333333333333*GSa[t]
};

events = {

};

speciesAnnotations = {

};

reactionAnnotations = {

};

units = {
{"time" -> "min", "metabolite" -> "mmol/L", "extent" -> "mM"}
};

(* Time evolution *)
odes = {
B\[LetterSpace]alan'[t] == 1.0*v\[LetterSpace]tS4 -1.0*v\[LetterSpace]d\[LetterSpace]Balan -1.0*v\[LetterSpace]tL6, B\[LetterSpace]ffa'[t] == 1.0*v\[LetterSpace]tF3 -1.0*v\[LetterSpace]d\[LetterSpace]Bffa -1.0*v\[LetterSpace]tL5, B\[LetterSpace]gluc'[t] == 1.0*v\[LetterSpace]feed -1.0*v\[LetterSpace]d\[LetterSpace]Bgluc -1.0*v\[LetterSpace]tF1 -1.0*v\[LetterSpace]tL1 -1.0*v\[LetterSpace]tS1, B\[LetterSpace]glucgn'[t] == 1.0*v\[LetterSpace]glucgn -1.0*v\[LetterSpace]Ggluc -1.0*v\[LetterSpace]d\[LetterSpace]Bglucgn, B\[LetterSpace]ins'[t] == 1.0*v\[LetterSpace]Igluc +1.0*v\[LetterSpace]ins -1.0*v\[LetterSpace]d\[LetterSpace]Bins, B\[LetterSpace]ket'[t] == 1.0*v\[LetterSpace]tL3 -1.0*v\[LetterSpace]d\[LetterSpace]Bket -1.0*v\[LetterSpace]tS2, B\[LetterSpace]lac'[t] == 1.0*v\[LetterSpace]tS3 -1.0*v\[LetterSpace]d\[LetterSpace]Blac -1.0*v\[LetterSpace]tL2, Cvar'[t] == 1.0*Jg10 +1.0*Jg11 , F\[LetterSpace]TG'[t] == 1.0*v\[LetterSpace]F3 -1.0*v\[LetterSpace]F4, F\[LetterSpace]acyl'[t] == 1.0*v\[LetterSpace]F1 +1.0*v\[LetterSpace]F5 -1.0*v\[LetterSpace]F3, F\[LetterSpace]g6p'[t] == 1.0*v\[LetterSpace]tF1 -1.0*v\[LetterSpace]F1, GPa'[t] == 1.0*Jg5 -1.0*Jg6 -1.0*Jg9, GSa'[t] == 1.0*Jg8 -1.0*Jg7, M\[LetterSpace]g6p'[t] == 1.0*v\[LetterSpace]S1r +1.0*v\[LetterSpace]tS1 -1.0*v\[LetterSpace]S1f -1.0*v\[LetterSpace]S2, M\[LetterSpace]glycgn'[t] == 1.0*v\[LetterSpace]S1f -1.0*v\[LetterSpace]S1r, M\[LetterSpace]ket'[t] == 1.0*v\[LetterSpace]tS2 -1.0*v\[LetterSpace]s\[LetterSpace]dket, M\[LetterSpace]lac'[t] == 1.0*v\[LetterSpace]S3f -1.0*v\[LetterSpace]S3r -1.0*v\[LetterSpace]tS3, M\[LetterSpace]pyr'[t] == 2.0*v\[LetterSpace]S2 +1.0*v\[LetterSpace]S3r +1.0*v\[LetterSpace]S4r -1.0*v\[LetterSpace]S3f -1.0*v\[LetterSpace]S4f, PKa'[t] == 1.0*Jg3 -1.0*Jg4, PP1'[t] ==  -1.0*Jg9, PP1\[LetterSpace]GPa'[t] == 1.0*Jg9 , R2C2'[t] ==  -1.0*Jg10, R2\[LetterSpace]C\[LetterSpace]cAMP2'[t] == 1.0*Jg10 -1.0*Jg11, R2\[LetterSpace]cAMP4'[t] == 1.0*Jg11 , aK'[t] == 1.0*v\[LetterSpace]L21f +1.0*v\[LetterSpace]L9 -1.0*v\[LetterSpace]L10 -1.0*v\[LetterSpace]L21r, acet\[LetterSpace]c'[t] == 1.0*v\[LetterSpace]L14 -1.0*v\[LetterSpace]L15 -1.0*v\[LetterSpace]L16, acet\[LetterSpace]m'[t] == 8.0*v\[LetterSpace]L18 +1.0*v\[LetterSpace]L7 -2.0*v\[LetterSpace]L19 -1.0*v\[LetterSpace]L8, alan'[t] == 1.0*v\[LetterSpace]L21f +1.0*v\[LetterSpace]tL6 -1.0*v\[LetterSpace]L20 -1.0*v\[LetterSpace]L21r, cAMP'[t] == 1.0*v\[LetterSpace]c1cAMP -2.0*Jg10 -2.0*Jg11 -1.0*v\[LetterSpace]c2cAMP, citrate'[t] == 1.0*v\[LetterSpace]L8 -1.0*v\[LetterSpace]L14 -1.0*v\[LetterSpace]L9, g6p'[t] == 1.0*v\[LetterSpace]L1f +1.0*v\[LetterSpace]L2r +1.0*v\[LetterSpace]L3r -1.0*v\[LetterSpace]L1r -1.0*v\[LetterSpace]L2f -1.0*v\[LetterSpace]L3f, gluc'[t] == 1.0*v\[LetterSpace]L1r +1.0*v\[LetterSpace]tL1 -1.0*v\[LetterSpace]L1f, glutamate'[t] == 1.0*v\[LetterSpace]L21r -1.0*v\[LetterSpace]L21f -1.0*v\[LetterSpace]L22, glycgn'[t] == 1.0*v\[LetterSpace]L2f -1.0*v\[LetterSpace]L2r, ket'[t] == 1.0*v\[LetterSpace]L19 -1.0*v\[LetterSpace]tL3, lac'[t] == 1.0*v\[LetterSpace]L5f +1.0*v\[LetterSpace]tL2 -1.0*v\[LetterSpace]L5r, malate'[t] == 1.0*v\[LetterSpace]L10 +1.0*v\[LetterSpace]L11f -1.0*v\[LetterSpace]L11r -1.0*v\[LetterSpace]L12, malonyl'[t] == 1.0*v\[LetterSpace]L15 -7.0*v\[LetterSpace]L16, oa\[LetterSpace]c'[t] == 1.0*v\[LetterSpace]L12 +1.0*v\[LetterSpace]L14 -1.0*v\[LetterSpace]L13, oa\[LetterSpace]m'[t] == 1.0*v\[LetterSpace]L11r +1.0*v\[LetterSpace]L6 -1.0*v\[LetterSpace]L11f -1.0*v\[LetterSpace]L8, palm'[t] == 1.0*v\[LetterSpace]L16 +1.0*v\[LetterSpace]tL5 -1.0*v\[LetterSpace]L17, palmCoA'[t] == 1.0*v\[LetterSpace]L17 -1.0*v\[LetterSpace]L18, pep'[t] == 1.0*v\[LetterSpace]L13 +2.0*v\[LetterSpace]L3f -2.0*v\[LetterSpace]L3r -1.0*v\[LetterSpace]L4, pyr'[t] == 1.0*v\[LetterSpace]L20 +1.0*v\[LetterSpace]L21r +1.0*v\[LetterSpace]L4 +1.0*v\[LetterSpace]L5r -1.0*v\[LetterSpace]L21f -1.0*v\[LetterSpace]L5f -1.0*v\[LetterSpace]L6 -1.0*v\[LetterSpace]L7
};
timeCourse = NDSolve[Join[odes, initialValues]//.rateEquations//.assignments//.parameters, variables, {t, 0, 100}];

(* Steady-state solution initialized with result of time evolution *)
findRootEquations = odes /.D[_[t],t]->0;
findRootVariables = Partition[Flatten[{#, #/.timeCourse/.t->100} &/@variables],2];
steadyStateVariables = FindRoot[findRootEquations//.rateEquations//.assignments//.parameters, findRootVariables, MaxIterations->100]
fluxes = #//.assignments//.parameters/.steadyStateVariables&/@rateEquations

(* Plot the time evolution of the variables *)
plotTable=Table[Plot[variables[[i]]/.parameters/.timeCourse,{t,0,100},PlotLegends->variables[[i]],PlotRange->Full],{i,Length[variables]}]