(* 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 = {
IE\[LetterSpace]p[t], UbE\[LetterSpace]star[t], cdc25\[LetterSpace]p[t], cyclin[t], dimer[t], dimer\[LetterSpace]p[t], p\[LetterSpace]dimer[t], p\[LetterSpace]dimer\[LetterSpace]p[t], wee1\[LetterSpace]p[t]
};

initialValues = {
IE\[LetterSpace]p[0] == 0.0, UbE\[LetterSpace]star[0] == 0.0, cdc25\[LetterSpace]p[0] == 0.0, cyclin[0] == 100.0, dimer[0] == 0.0, dimer\[LetterSpace]p[0] == 0.0, p\[LetterSpace]dimer[0] == 0.0, p\[LetterSpace]dimer\[LetterSpace]p[0] == 0.0, wee1\[LetterSpace]p[0] == 0.0
};

rates = {
R1, R10, R11, R12, R13, R14, R15, R17, R18, R19, R2, R20, R21, R22, R23, R24, R3, R4, R5, R6, R7, R8, R9
};

rateEquations = {
R1 -> k1AA, R10 -> k2*p\[LetterSpace]dimer[t], R11 -> kinh*p\[LetterSpace]dimer\[LetterSpace]p[t], R12 -> kwee*dimer\[LetterSpace]p[t], R13 -> k25*p\[LetterSpace]dimer\[LetterSpace]p[t], R14 -> k2*p\[LetterSpace]dimer\[LetterSpace]p[t], R15 -> k2*dimer\[LetterSpace]p[t], R17 -> (ka*(total\[LetterSpace]cdc25 - cdc25\[LetterSpace]p[t])*dimer\[LetterSpace]p[t])/(K\[LetterSpace]a + total\[LetterSpace]cdc25 - cdc25\[LetterSpace]p[t]), R18 -> (kbPPase*cdc25\[LetterSpace]p[t])/(K\[LetterSpace]b + cdc25\[LetterSpace]p[t]), R19 -> (ke*dimer\[LetterSpace]p[t]*(total\[LetterSpace]wee1 - wee1\[LetterSpace]p[t]))/(K\[LetterSpace]e + total\[LetterSpace]wee1 - wee1\[LetterSpace]p[t]), R2 -> k2*cyclin[t], R20 -> (kfPPase*wee1\[LetterSpace]p[t])/(K\[LetterSpace]f + wee1\[LetterSpace]p[t]), R21 -> (kg*dimer\[LetterSpace]p[t]*(total\[LetterSpace]IE - IE\[LetterSpace]p[t]))/(K\[LetterSpace]g + total\[LetterSpace]IE - IE\[LetterSpace]p[t]), R22 -> (khPPAse*IE\[LetterSpace]p[t])/(K\[LetterSpace]h + IE\[LetterSpace]p[t]), R23 -> (kc*IE\[LetterSpace]p[t]*(total\[LetterSpace]UbE - UbE\[LetterSpace]star[t]))/(K\[LetterSpace]c + total\[LetterSpace]UbE - UbE\[LetterSpace]star[t]), R24 -> (kd\[LetterSpace]anti\[LetterSpace]IE*UbE\[LetterSpace]star[t])/(K\[LetterSpace]d + UbE\[LetterSpace]star[t]), R3 -> cdc2*k3*cyclin[t], R4 -> kinh*dimer\[LetterSpace]p[t], R5 -> kwee*dimer[t], R6 -> kcak*dimer[t], R7 -> k2*dimer[t], R8 -> k25*p\[LetterSpace]dimer[t], R9 -> kcak*p\[LetterSpace]dimer[t]
};

parameters = {
K\[LetterSpace]a -> 0.1, K\[LetterSpace]b -> 0.1, K\[LetterSpace]c -> 0.01, K\[LetterSpace]d -> 0.01, K\[LetterSpace]e -> 0.3, K\[LetterSpace]f -> 0.3, K\[LetterSpace]g -> 0.01, K\[LetterSpace]h -> 0.01, V25\[LetterSpace]double\[LetterSpace]prime -> 2.0, V25\[LetterSpace]prime -> 0.1, V2\[LetterSpace]double\[LetterSpace]prime -> 1.0, V2\[LetterSpace]prime -> 0.015, Vwee\[LetterSpace]double\[LetterSpace]prime -> 1.0, Vwee\[LetterSpace]prime -> 0.1, k1AA -> 1.0, k3 -> 0.01, ka -> 0.01, kbPPase -> 0.125, kc -> 0.1, kcak -> 0.25, kd\[LetterSpace]anti\[LetterSpace]IE -> 0.095, ke -> 0.0133, kfPPase -> 0.1, kg -> 0.0065, khPPAse -> 0.087, kinh -> 0.025, total\[LetterSpace]IE -> 1.0, total\[LetterSpace]UbE -> 1.0, total\[LetterSpace]cdc2 -> 100.0, total\[LetterSpace]cdc25 -> 1.0, total\[LetterSpace]wee1 -> 1.0, cytoplasm -> 1.0
};

assignments = {
wee1 -> total\[LetterSpace]wee1 - wee1\[LetterSpace]p[t], cdc25 -> total\[LetterSpace]cdc25 - cdc25\[LetterSpace]p[t], cdc2 -> total\[LetterSpace]cdc2 - dimer[t] - dimer\[LetterSpace]p[t] - p\[LetterSpace]dimer[t] - p\[LetterSpace]dimer\[LetterSpace]p[t], UbE -> total\[LetterSpace]UbE - UbE\[LetterSpace]star[t], IE -> total\[LetterSpace]IE - IE\[LetterSpace]p[t], k2 -> V2\[LetterSpace]prime*(total\[LetterSpace]UbE - UbE\[LetterSpace]star[t]) + V2\[LetterSpace]double\[LetterSpace]prime*UbE\[LetterSpace]star[t], kwee -> Vwee\[LetterSpace]double\[LetterSpace]prime*(total\[LetterSpace]wee1 - wee1\[LetterSpace]p[t]) + Vwee\[LetterSpace]prime*wee1\[LetterSpace]p[t], Y15P -> p\[LetterSpace]dimer[t] + p\[LetterSpace]dimer\[LetterSpace]p[t], total\[LetterSpace]cyclin -> cyclin[t] + dimer[t] + dimer\[LetterSpace]p[t] + p\[LetterSpace]dimer[t] + p\[LetterSpace]dimer\[LetterSpace]p[t], k25 -> V25\[LetterSpace]prime*(total\[LetterSpace]cdc25 - cdc25\[LetterSpace]p[t]) + V25\[LetterSpace]double\[LetterSpace]prime*cdc25\[LetterSpace]p[t]
};

events = {

};

speciesAnnotations = {
cyclin[t]->"http://identifiers.org/interpro/IPR015454"
};

reactionAnnotations = {
R19->"http://identifiers.org/ec-code/2.7.11.22", R19->"http://identifiers.org/go/GO:0044257", R21->"http://identifiers.org/ec-code/2.7.11.22"
};

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

(* Time evolution *)
odes = {
IE\[LetterSpace]p'[t] == 1.0*R21 -1.0*R22, UbE\[LetterSpace]star'[t] == 1.0*R23 -1.0*R24, cdc25\[LetterSpace]p'[t] == 1.0*R17 -1.0*R18, cyclin'[t] == 1.0*R1 -1.0*R2 -1.0*R3, dimer'[t] == 1.0*R3 +1.0*R4 +1.0*R8 -1.0*R5 -1.0*R6 -1.0*R7, dimer\[LetterSpace]p'[t] == 1.0*R13 +1.0*R6 -1.0*R4 -1.0*R12 -1.0*R15, p\[LetterSpace]dimer'[t] == 1.0*R5 +1.0*R11 -1.0*R8 -1.0*R9 -1.0*R10, p\[LetterSpace]dimer\[LetterSpace]p'[t] == 1.0*R9 +1.0*R12 -1.0*R11 -1.0*R13 -1.0*R14, wee1\[LetterSpace]p'[t] == 1.0*R19 -1.0*R20
};
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]}]