(* 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 = {
s10[t], s11[t], s12[t], s13[t], s14[t], s15[t], s16[t], s17[t], s18[t], s19[t], s2[t], s21[t], s23[t], s25[t], s26[t], s27[t], s29[t], s3[t], s31[t], s33[t], s35[t], s36[t], s37[t], s38[t], s4[t], s5[t], s6[t], s7[t], s8[t], s9[t]
};

initialValues = {
s10[0] == 17.7476118652, s11[0] == 0.0, s12[0] == 1.0, s13[0] == 0.759332005, s14[0] == 0.0, s15[0] == 0.0, s16[0] == 0.0, s17[0] == 0.0, s18[0] == 0.0, s19[0] == 0.0, s2[0] == 1.0, s21[0] == 0.0, s23[0] == 37.81914621, s25[0] == 2.507671281, s26[0] == 852.608418, s27[0] == 3.930695895, s29[0] == 0.012942682, s3[0] == 14.0774258421, s31[0] == 0.004015031, s33[0] == 0.000988831, s35[0] == 0.043352951, s36[0] == 0.0, s37[0] == 0.0, s38[0] == 0.0, s4[0] == 104.07239819, s5[0] == 108.094519859, s6[0] == 1.79487179487, s7[0] == 25.1885369533, s8[0] == 68.8788335846, s9[0] == 1.0
};

rates = {
re30, re31, re32, re33, re34, re35, re36, re37, re38, re39, re40, re41, re42, re43, re44, re45, re46, re47, re48, re49, re50, re51, re52, re53, re54, re55, re56, re57, re58, re59
};

rateEquations = {
re30 -> (K\[LetterSpace]PFKL\[LetterSpace]akg*K\[LetterSpace]PFKL\[LetterSpace]cit*K\[LetterSpace]PFKL\[LetterSpace]icit*K\[LetterSpace]PFKL\[LetterSpace]mal*K\[LetterSpace]PFKL\[LetterSpace]pep*s22*Vf\[LetterSpace]PFKL*s13[t]*s9[t])/((K\[LetterSpace]PFKL\[LetterSpace]f6p + s22)*(K\[LetterSpace]PFKL\[LetterSpace]cit + s10[t])*(K\[LetterSpace]PFKL\[LetterSpace]PHOS\[LetterSpace]S775 + s13[t])*(K\[LetterSpace]PFKL\[LetterSpace]pep + s5[t])*(K\[LetterSpace]PFKL\[LetterSpace]icit + s6[t])*(K\[LetterSpace]PFKL\[LetterSpace]akg + s7[t])*(K\[LetterSpace]PFKL\[LetterSpace]mal + s8[t])*(K\[LetterSpace]PFKL\[LetterSpace]f26bp + s9[t])), re31 -> (K\[LetterSpace]FBPase\[LetterSpace]f26bp*Vf\[LetterSpace]FBPase*s10[t]*s4[t])/((K\[LetterSpace]FBPase\[LetterSpace]cit + s10[t])*(K\[LetterSpace]FBPase\[LetterSpace]f16bp + s4[t])*(K\[LetterSpace]FBPase\[LetterSpace]f26bp + s9[t])), re32 -> k\[LetterSpace]ALDO*s4[t], re33 -> -k\[LetterSpace]mal, re34 -> -k\[LetterSpace]akg, re35 -> -k\[LetterSpace]pep, re36 -> -k\[LetterSpace]icit, re37 -> -k\[LetterSpace]f26bp, re38 -> -k\[LetterSpace]cit, re39 -> -k\[LetterSpace]f6p, re40 -> re40\[LetterSpace]k1*(s24 - s23[t]), re41 -> insulin*re41\[LetterSpace]k1*s23[t] - re41\[LetterSpace]k2*s25[t], re42 -> re42\[LetterSpace]k1*s23[t]*s31[t], re43 -> insulin*re43\[LetterSpace]k1*s26[t] - re43\[LetterSpace]k2*s27[t], re44 -> re44\[LetterSpace]k1*s25[t]*s31[t], re45 -> re45\[LetterSpace]k1*s27[t], re46 -> re46\[LetterSpace]k1*s26[t], re47 -> re47\[LetterSpace]k1*s27[t], re48 -> re48\[LetterSpace]k1*s25[t]*(s28 - s29[t]), re49 -> re49\[LetterSpace]k1*s29[t], re50 -> re50\[LetterSpace]k1*s26[t], re51 -> re51\[LetterSpace]k1*s25[t], re52 -> re52\[LetterSpace]k1*s29[t]*(s30 - s31[t]), re53 -> re53\[LetterSpace]k1*s31[t], re54 -> re54\[LetterSpace]k1*s31[t]*(s32 - s33[t]), re55 -> re55\[LetterSpace]k1*s33[t]*s35[t], re56 -> re56\[LetterSpace]k1*s31[t]*(s34 - s35[t]), re57 -> re57\[LetterSpace]k1*s35[t], re58 -> re58\[LetterSpace]k1*(s1 - s13[t])*s33[t], re59 -> re59\[LetterSpace]k1*s13[t]
};

parameters = {
K\[LetterSpace]FBPase\[LetterSpace]cit -> 0.0211646, K\[LetterSpace]FBPase\[LetterSpace]f16bp -> 0.104089638, K\[LetterSpace]FBPase\[LetterSpace]f26bp -> 17.51744342, K\[LetterSpace]PFKL\[LetterSpace]PHOS\[LetterSpace]S775 -> 6.283705757, K\[LetterSpace]PFKL\[LetterSpace]akg -> 24661.01154, K\[LetterSpace]PFKL\[LetterSpace]cit -> 41.30426029, K\[LetterSpace]PFKL\[LetterSpace]f26bp -> 1.282443082, K\[LetterSpace]PFKL\[LetterSpace]f6p -> 0.014114844, K\[LetterSpace]PFKL\[LetterSpace]icit -> 1784.508205, K\[LetterSpace]PFKL\[LetterSpace]mal -> 9.544729149, K\[LetterSpace]PFKL\[LetterSpace]pep -> 0.633518366, Vf\[LetterSpace]FBPase -> 9.932861302, Vf\[LetterSpace]PFKL -> 695063.7194, insulin -> 1.0, k\[LetterSpace]ALDO -> 0.008187906, k\[LetterSpace]akg -> -3.544494721, k\[LetterSpace]cit -> -0.351935646, k\[LetterSpace]f26bp -> -0.083430336, k\[LetterSpace]f6p -> -0.930115636, k\[LetterSpace]icit -> -0.038210156, k\[LetterSpace]mal -> 1.005530417, k\[LetterSpace]pep -> 43.99195576, s1 -> 1.235570941, s24 -> 46.2222520609226, s28 -> 4.33812187331668, s30 -> 0.0959163701057347, s32 -> 2.7769943784158, s34 -> 14.9913345914433, re40\[LetterSpace]k1 -> 0.0477985900779305, re41\[LetterSpace]k1 -> 7.78160761103111, re41\[LetterSpace]k2 -> 1.61147523779118, re42\[LetterSpace]k1 -> 9.93311225447353, re43\[LetterSpace]k1 -> 0.363030286526969, re43\[LetterSpace]k2 -> 0.408127912886804, re44\[LetterSpace]k1 -> 3.88248960751442*^-05, re45\[LetterSpace]k1 -> 0.28442597446931, re46\[LetterSpace]k1 -> 1.00000000282413*^-05, re47\[LetterSpace]k1 -> 0.0948960328385619, re48\[LetterSpace]k1 -> 0.00919578911309774, re49\[LetterSpace]k1 -> 7.70618517548016, re50\[LetterSpace]k1 -> 1.00277786609339*^-05, re51\[LetterSpace]k1 -> 0.00792717614041237, re52\[LetterSpace]k1 -> 0.419682384304397, re53\[LetterSpace]k1 -> 0.124330492920416, re54\[LetterSpace]k1 -> 0.00752464611370891, re55\[LetterSpace]k1 -> 1.95497593092361, re56\[LetterSpace]k1 -> 0.00105342379833469, re57\[LetterSpace]k1 -> 0.00145811601430322, re58\[LetterSpace]k1 -> 26.8316707654711, re59\[LetterSpace]k1 -> 0.0166525139097609, default -> 1.0
};

assignments = {
s22 -> s3[t]
};

events = {

};

speciesAnnotations = {

};

reactionAnnotations = {

};

units = {
{"time" -> "", "metabolite" -> "", "extent" -> ""}
};

(* Time evolution *)
odes = {
s10'[t] ==  -1.0*re38, s11'[t] == 1.0*re32 , s12'[t] == 0.0 , s13'[t] == 1.0*re58 -1.0*re59, s14'[t] == 1.0*re33 , s15'[t] == 1.0*re34 , s16'[t] == 1.0*re35 , s17'[t] == 1.0*re36 , s18'[t] == 1.0*re37 , s19'[t] == 1.0*re38 , s2'[t] == 0.0 , s21'[t] == 1.0*re39 , s23'[t] == 1.0*re40 +1.0*re46 -1.0*re41 -1.0*re42, s25'[t] == 1.0*re41 +1.0*re45 -1.0*re44 -1.0*re51, s26'[t] == 1.0*re42 -1.0*re43 -1.0*re46 -1.0*re50, s27'[t] == 1.0*re43 +1.0*re44 -1.0*re45 -1.0*re47, s29'[t] == 1.0*re48 -1.0*re49, s3'[t] ==  -1.0*re39, s31'[t] == 1.0*re52 -1.0*re53, s33'[t] == 1.0*re54 -1.0*re55, s35'[t] == 1.0*re56 -1.0*re57, s36'[t] == 1.0*re47 , s37'[t] == 1.0*re50 , s38'[t] == 1.0*re51 , s4'[t] == 1.0*re30 -1.0*re31 -1.0*re32, s5'[t] ==  -1.0*re35, s6'[t] ==  -1.0*re36, s7'[t] ==  -1.0*re34, s8'[t] ==  -1.0*re33, s9'[t] ==  -1.0*re37
};
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]}]