(* 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 = {
X1[t], X3[t], X5[t], X7[t], X6[t], X4[t], X2[t]
};

initialValues = {
X1[0] == 0.03, X3[0] == 0.1, X5[0] == 1.0, X7[0] == 0.05, X6[0] == 0.02, X4[0] == 0.7, X2[0] == 1.0
};

rates = {

};

rateEquations = {

};

parameters = {
heat\[LetterSpace]shock -> 0.0, X0 -> 1.0, X1 -> 0.03, X2 -> 1.0, X3 -> 0.1, X4 -> 0.7, X5 -> 1.0, X6 -> 0.02, X7 -> 0.05, cell -> 1.0, external -> 1.0
};

assignments = {
flux\[LetterSpace]X5\[LetterSpace]in -> 11.06*X14^0.6*X15f^0.4*X2[t]^0.04*X3[t]^0.32*X4[t]^0.16, flux\[LetterSpace]X3\[LetterSpace]out -> (76.434*X12r^0.718*X13^0.18*X15f^0.103*X3[t]^0.593)/X2[t]^0.412, flux\[LetterSpace]X2\[LetterSpace]out -> (30.12*X10^0.5111*X11^0.0667*X12f^0.411*X17^0.0111*X2[t]^0.575*X4[t]^0.00333)/(X1[t]^0.00333*X3[t]^0.17), flux\[LetterSpace]X1\[LetterSpace]in -> (31.912*X0^0.968*X19^0.0323*X8^0.968*X7[t]^0.00968)/X2[t]^0.194, flux\[LetterSpace]X1\[LetterSpace]out -> (89.935*X9*X1[t]^0.75)/X6[t]^0.4, flux\[LetterSpace]X2\[LetterSpace]in -> (142.72*X12r^0.311*X9^0.689*X1[t]^0.517*X3[t]^0.183)/(X2[t]^0.179*X6[t]^0.276), flux\[LetterSpace]X3\[LetterSpace]in -> (7.8819*X12f^0.949*X15r^0.0513*X2[t]^0.394*X5[t]^0.0128)/(X3[t]^0.392*X4[t]^0.01), flux\[LetterSpace]X4\[LetterSpace]in -> 11.07*X13*X3[t]^0.5, flux\[LetterSpace]X4\[LetterSpace]out -> (3.4556*X14^0.857*X17^0.143*X2[t]^0.214*X4[t]^0.386)/X1[t]^0.0429, flux\[LetterSpace]X5\[LetterSpace]out -> (4.929*X15r^0.2*X16^0.8*X5[t]^0.25)/(X2[t]^0.04*X4[t]^0.04), flux\[LetterSpace]X6\[LetterSpace]in -> (0.19424*X17*X2[t]^0.3*X4[t]^0.3)/X1[t]^0.3, flux\[LetterSpace]X6\[LetterSpace]out -> 1.0939*X18*X6[t]^0.2, flux\[LetterSpace]X7\[LetterSpace]in -> 1.0939*X18*X6[t]^0.2, flux\[LetterSpace]X7\[LetterSpace]out -> 1.2288*X19*X7[t]^0.3, X8 -> Piecewise[{{8, heat\[LetterSpace]shock == 1}}, 1], X9 -> 1, X10 -> 1, X11 -> 1, X12f -> 1, X12r -> Piecewise[{{16, heat\[LetterSpace]shock == 1}}, 1], X13 -> Piecewise[{{16, heat\[LetterSpace]shock == 1}}, 1], X14 -> 1, X15f -> 1, X15r -> Piecewise[{{50, heat\[LetterSpace]shock == 1}}, 1], X16 -> Piecewise[{{16, heat\[LetterSpace]shock == 1}}, 1], X17 -> 1, X18 -> Piecewise[{{18, heat\[LetterSpace]shock == 1}}, 1], X19 -> 1
};

events = {

};

speciesAnnotations = {
X0[t]->"http://identifiers.org/chebi/CHEBI:17925", X1[t]->"http://identifiers.org/chebi/CHEBI:17925", X2[t]->"http://identifiers.org/chebi/CHEBI:17665", X3[t]->"http://identifiers.org/chebi/CHEBI:16077", X4[t]->"http://identifiers.org/chebi/CHEBI:18066", X5[t]->"http://identifiers.org/chebi/CHEBI:28087", X6[t]->"http://identifiers.org/chebi/CHEBI:18283", X7[t]->"http://identifiers.org/chebi/CHEBI:16551"
};

reactionAnnotations = {

};

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

(* Time evolution *)
odes = {
X1'[t] == flux\[LetterSpace]X1\[LetterSpace]in - flux\[LetterSpace]X1\[LetterSpace]out, X3'[t] == flux\[LetterSpace]X3\[LetterSpace]in - flux\[LetterSpace]X3\[LetterSpace]out, X5'[t] == flux\[LetterSpace]X5\[LetterSpace]in - flux\[LetterSpace]X5\[LetterSpace]out, X7'[t] == flux\[LetterSpace]X7\[LetterSpace]in - flux\[LetterSpace]X7\[LetterSpace]out, X6'[t] == flux\[LetterSpace]X6\[LetterSpace]in - flux\[LetterSpace]X6\[LetterSpace]out, X4'[t] == flux\[LetterSpace]X4\[LetterSpace]in - flux\[LetterSpace]X4\[LetterSpace]out, X2'[t] == flux\[LetterSpace]X2\[LetterSpace]in - flux\[LetterSpace]X2\[LetterSpace]out
};
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]}]