(* 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 = {
Cdc15[t], Cdc20[t], Cdh1[t], Clb2[t], MEN[t], Net1[t], Net1Cdc14[t], Polo[t], PoloT[t], Tem1[t], securinT[t], securinseparase[t], separaseT[t]
};

initialValues = {
Cdc15[0] == 0.0257, Cdc20[0] == 0.0, Cdh1[0] == 3.8*^-05, Clb2[0] == 0.997, MEN[0] == 8*^-06, Net1[0] == 0.098, Net1Cdc14[0] == 0.485, Polo[0] == 0.945, PoloT[0] == 0.99, Tem1[0] == 0.00389, securinT[0] == 0.6, securinseparase[0] == 0.248, separaseT[0] == 0.25
};

rates = {
securinseparase\[LetterSpace]securindegradation, securinseparase\[LetterSpace]separase\[LetterSpace]degradation, v1, v10, v11, v13, v14, v15, v16, v17, v18, v19, v2, v20, v21, v22, v23, v24, v25, v26, v27, v28, v29, v3, v30, v4, v5, v6, v7, v8, v9
};

rateEquations = {
securinseparase\[LetterSpace]securindegradation -> (kdsecurin + kadsecurin*Cdc20[t])*securinseparase[t], securinseparase\[LetterSpace]separase\[LetterSpace]degradation -> kdseparase*securinseparase[t], v1 -> ksclb2, v10 -> kdseparase*separaseT[t], v11 -> lasecurin*securin*separase - ldsecurin*securinseparase[t], v13 -> (Net1P*(Cdc14*kad + kd*PP2A))/(Jnet + Net1P), v14 -> ((Cdk*kp + kap*MEN[t])*Net1[t])/(Jnet + Net1[t] + Net1Cdc14[t]), v15 -> Cdc14*lanet*Net1[t], v16 -> ldnet*Net1Cdc14[t], v17 -> ((Cdk*kp + kap*MEN[t])*Net1Cdc14[t])/(Jnet + Net1[t] + Net1Cdc14[t]), v18 -> kspolo, v19 -> (kdpolo + kadpolo*Cdh1[t])*PoloT[t], v2 -> (kdclb2 + kadclb2*Cdc20[t] + kaadclb2*Cdh1[t])*Clb2[t], v20 -> ((Cdk*kaapolo + kapolo)*(-Polo[t] + PoloT[t]))/(Jpolo - Polo[t] + PoloT[t]), v21 -> (kipolo*Polo[t])/(Jpolo + Polo[t]), v22 -> (kdpolo + kadpolo*Cdh1[t])*Polo[t], v23 -> ((katem + kaatem*Polo[t])*(Tem1T - Tem1[t]))/(Jtem1 + Tem1T - Tem1[t]), v24 -> ((kitem + kaitem*PP2A)*Tem1[t])/(Jtem1 + Tem1[t]), v25 -> ((Cdc14*kaacdc15 + kacdc15)*(Cdc15T - Cdc15[t]))/(Cdc15T + Jcdc15 - Cdc15[t]), v26 -> ((Cdk*kaicdc15 + kicdc15)*Cdc15[t])/(Jcdc15 + Cdc15[t]), v27 -> lamen*(Cdc15[t] - MEN[t])*(-MEN[t] + Tem1[t]), v28 -> ldmen*MEN[t], v29 -> ((kitem + kaitem*PP2A)*MEN[t])/(Jtem1 + Tem1[t]), v3 -> kscdc20, v30 -> ((Cdk*kaicdc15 + kicdc15)*MEN[t])/(Jcdc15 + Cdc15[t]), v4 -> Cdc20[t]*(kdcdc20 + kadcdc20*Cdh1[t]), v5 -> ((Cdc14*kadcdh + kdcdh)*(Cdh1T - Cdh1[t]))/(Cdh1T + Jcdh - Cdh1[t]), v6 -> (Cdk*kapcdh*Cdh1[t])/(Jcdh + Cdh1[t]), v7 -> kssecurin, v8 -> (kdsecurin + kadsecurin*Cdc20[t])*securinT[t], v9 -> ksseparase
};

parameters = {
Cdc14T -> 0.5, Cdc15T -> 1.0, Cdh1T -> 1.0, Inh -> 0.0, Jcdc15 -> 0.2, Jcdh -> 0.0015, Jnet -> 0.2, Jpolo -> 0.25, Jtem1 -> 0.005, Net1T -> 1.0, PP2AT -> 1.0, Tem1T -> 1.0, kaacdc15 -> 0.5, kaadclb2 -> 2.0, kaapolo -> 0.5, kaatem -> 0.5, kacdc15 -> 0.02, kad -> 0.1, kadcdc20 -> 2.0, kadcdh -> 1.0, kadclb2 -> 0.2, kadpolo -> 0.25, kadsecurin -> 2.0, kaicdc15 -> 0.2, kaitem -> 1.0, kap -> 2.0, kapcdh -> 1.0, kapolo -> 0.0, katem -> 0.0, kd -> 0.45, kdcdc20 -> 0.05, kdcdh -> 0.01, kdclb2 -> 0.03, kdpolo -> 0.01, kdsecurin -> 0.05, kdseparase -> 0.004, ki -> 20.0, kicdc15 -> 0.0, kipolo -> 0.1, kitem -> 0.1, kp -> 0.4, kpp -> 0.1, kscdc20 -> 0.015, ksclb2 -> 0.03, kspolo -> 0.01, kssecurin -> 0.03, ksseparase -> 0.001, lamen -> 10.0, lanet -> 500.0, lasecurin -> 500.0, ldmen -> 0.1, ldnet -> 1.0, ldsecurin -> 1.0, AA -> 1.0, degr -> 1.0, compartment -> 1.0
};

assignments = {
PP2A -> (PP2AT*(1 + ki*kpp*separase))/(1 + ki*separase), Cdc15\[LetterSpace]i -> Cdc15T - Cdc15[t], Cdc14 -> Cdc14T - Net1Cdc14[t], Net1P -> Net1T - Net1[t] - Net1Cdc14[t], Cdk -> Clb2[t]/(1 + Inh), Polo\[LetterSpace]i -> -Polo[t] + PoloT[t], securin -> -securinseparase[t] + securinT[t], Tem1\[LetterSpace]i -> Tem1T - Tem1[t], Cdh1\[LetterSpace]i -> Cdh1T - Cdh1[t], separase -> -securinseparase[t] + separaseT[t]
};

events = {

};

speciesAnnotations = {

};

reactionAnnotations = {

};

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

(* Time evolution *)
odes = {
Cdc15'[t] == 1.0*v25 -1.0*v26, Cdc20'[t] == 1.0*v3 -1.0*v4, Cdh1'[t] == 1.0*v5 -1.0*v6, Clb2'[t] == 1.0*v1 -1.0*v2, MEN'[t] == 1.0*v27 -1.0*v28 -1.0*v29 -1.0*v30, Net1'[t] == 1.0*v13 +1.0*v16 -1.0*v14 -1.0*v15, Net1Cdc14'[t] == 1.0*v15 -1.0*v16 -1.0*v17, Polo'[t] == 1.0*v20 -1.0*v21 -1.0*v22, PoloT'[t] == 1.0*v18 -1.0*v19, Tem1'[t] == 1.0*v23 -1.0*v24, securinT'[t] == 1.0*v7 -1.0*v8, securinseparase'[t] == 1.0*v11 -1.0*securinseparase\[LetterSpace]securindegradation -1.0*securinseparase\[LetterSpace]separase\[LetterSpace]degradation, separaseT'[t] == 1.0*v9 -1.0*v10
};
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]}]