(* 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 = {
IR[t], IRS1[t], IRS1p[t], IRS1ser307[t], IRS1ser307p[t], IRi[t], IRp[t], S6[t], S6K[t], S6KDN[t], S6Kp[t], S6p[t], mTOR[t], mTORp[t], mTORpS6KDN[t]
};

initialValues = {
IR[0] == 0.9181, IRS1[0] == 0.9999, IRS1p[0] == 0.0001, IRS1ser307[0] == 0.4614, IRS1ser307p[0] == 0.5386, IRi[0] == 0.0818, IRp[0] == 0.0001, S6[0] == 0.7765, S6K[0] == 0.9999, S6KDN[0] == 0.0, S6Kp[0] == 0.0001, S6p[0] == 0.2235, mTOR[0] == 0.9999, mTORp[0] == 0.0001, mTORpS6KDN[0] == 0.0
};

rates = {
v1, v1m, v2, v2m, v3, v3m, v4, v4m, v5, v5T, v5m, v6, v6m, v7, v7m, vM5T
};

rateEquations = {
v1 -> IRbasal*IR[t] + Ins*k1*IR[t], v1m -> k1m*IRp[t], v2 -> k2*IRp[t]*IRS1[t], v2m -> k2m*IRS1p[t], v3 -> k3*IRS1ser307[t]*S6Kp[t], v3m -> k3m*IRS1ser307p[t], v4 -> k4*IRS1p[t]*mTOR[t], v4m -> k4m*mTORp[t], v5 -> k5*mTORp[t]*S6K[t], v5T -> I*mTORp[t]*S6KDN[t], v5m -> k5m*S6Kp[t], v6 -> k6*S6[t]*S6Kp[t], v6m -> k6m*S6p[t], v7 -> k7*IRi[t], v7m -> k7m*IRp[t], vM5T -> k5tm*mTORpS6KDN[t]
};

parameters = {
I -> 0.0006, IRbasal -> 0.0635, Ins -> 100.0, k1 -> 6.0982, k1m -> 366.2795, k2 -> 0.4397, k2m -> 0.5174, k3 -> 20660.1506, k3m -> 1.9552, k4 -> 0.085, k4m -> 0.082, k5 -> 0.1965, k5m -> 0.1968, k5tm -> 1.7709, k6 -> 35.9042, k6m -> 0.0138, k7 -> 0.1499, k7m -> 97.5553, default -> 1.0
};

assignments = {
fig6c\[LetterSpace]IRS1ser307p -> 2.60412436101051*IRS1ser307p[t], fig6g\[LetterSpace]IRS1ser307p -> 283.31617806680646*(-0.538627018751669 + IRS1ser307p[t]), fig6f\[LetterSpace]S6p -> 223.46246244978784*(-0.000106747023438893 + S6p[t]), fig6e\[LetterSpace]IRS1p -> 83773.99869569197*(-0.000106747023438893 + IRS1p[t]), fig6d\[LetterSpace]S6p -> 5.0036269004157*S6p[t], fig6a\[LetterSpace]IRS1p -> 28765.5323430673*IRS1p[t], fig3b\[LetterSpace]S6p -> 140.688371115354*S6p[t], fig3b\[LetterSpace]IRS1ser307p -> 110.795803138294*IRS1ser307p[t], fig6b\[LetterSpace]S6Kp -> 4455.78636783405*S6Kp[t]
};

events = {

};

speciesAnnotations = {

};

reactionAnnotations = {

};

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

(* Time evolution *)
odes = {
IR'[t] == 1.0*v1m +1.0*v7 -1.0*v1, IRS1'[t] == 1.0*v2m -1.0*v2, IRS1p'[t] == 1.0*v2 -1.0*v2m, IRS1ser307'[t] == 1.0*v3m -1.0*v3, IRS1ser307p'[t] == 1.0*v3 -1.0*v3m, IRi'[t] == 1.0*v7m -1.0*v7, IRp'[t] == 1.0*v1 -1.0*v1m -1.0*v7m, S6'[t] == 1.0*v6m -1.0*v6, S6K'[t] == 1.0*v5m -1.0*v5, S6KDN'[t] == 1.0*vM5T -1.0*v5T, S6Kp'[t] == 1.0*v5 -1.0*v5m, S6p'[t] == 1.0*v6 -1.0*v6m, mTOR'[t] == 1.0*v4m -1.0*v4, mTORp'[t] == 1.0*v4 +1.0*vM5T -1.0*v4m -1.0*v5T, mTORpS6KDN'[t] == 1.0*v5T -1.0*vM5T
};
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]}]