(* 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 = {
APCP[t], BCKI[t], CKI[t], Cdc20A[t], Cdc20in[t], Cdh1[t], CycA[t], CycB[t], CycE[t], Mass[t], TriA[t], TriE[t], pB[t], pBCKI[t]
};

initialValues = {
APCP[0] == 0.0759140253067017, BCKI[0] == 0.679449200630188, CKI[0] == 0.02882070094347, Cdc20A[0] == 0.357272386550903, Cdc20in[0] == 0.770238757133484, Cdh1[0] == 0.718939363956451, CycA[0] == 0.0945030823349953, CycB[0] == 0.190358594059944, CycE[0] == 0.209202438592911, Mass[0] == 1.33826780319214, TriA[0] == 0.349222421646118, TriE[0] == 0.0, pB[0] == 0.0123442625626922, pBCKI[0] == 0.0479593835771084
};

rates = {
v\[LetterSpace]1, v\[LetterSpace]10, v\[LetterSpace]11, v\[LetterSpace]12, v\[LetterSpace]13, v\[LetterSpace]14, v\[LetterSpace]15, v\[LetterSpace]16, v\[LetterSpace]17, v\[LetterSpace]18, v\[LetterSpace]19, v\[LetterSpace]2, v\[LetterSpace]20, v\[LetterSpace]21, v\[LetterSpace]22, v\[LetterSpace]23, v\[LetterSpace]24, v\[LetterSpace]25, v\[LetterSpace]26, v\[LetterSpace]27, v\[LetterSpace]28, v\[LetterSpace]29, v\[LetterSpace]3, v\[LetterSpace]30, v\[LetterSpace]31, v\[LetterSpace]32, v\[LetterSpace]33, v\[LetterSpace]34, v\[LetterSpace]35, v\[LetterSpace]36, v\[LetterSpace]37, v\[LetterSpace]38, v\[LetterSpace]39, v\[LetterSpace]4, v\[LetterSpace]5, v\[LetterSpace]6, v\[LetterSpace]7, v\[LetterSpace]8, v\[LetterSpace]9
};

rateEquations = {
v\[LetterSpace]1 -> mu*Mass[t]*(1 - Mass[t]/MaxMass), v\[LetterSpace]10 -> kdissb*pBCKI[t], v\[LetterSpace]11 -> V25*pBCKI[t], v\[LetterSpace]12 -> Vwee*BCKI[t], v\[LetterSpace]13 -> Vdb*BCKI[t], v\[LetterSpace]14 -> Vdi*BCKI[t], v\[LetterSpace]15 -> Vdb*pBCKI[t], v\[LetterSpace]16 -> Vdi*pBCKI[t], v\[LetterSpace]17 -> Vsi, v\[LetterSpace]18 -> Vdi*CKI[t], v\[LetterSpace]19 -> kassa*CKI[t]*CycA[t], v\[LetterSpace]2 -> Vsb, v\[LetterSpace]20 -> kdissa*TriA[t], v\[LetterSpace]21 -> Vdi*TriA[t], v\[LetterSpace]22 -> Vda*TriA[t], v\[LetterSpace]23 -> kasse*CKI[t]*CycE[t], v\[LetterSpace]24 -> kdisse*TriE[t], v\[LetterSpace]25 -> Vdi*TriE[t], v\[LetterSpace]26 -> Vde*TriE[t], v\[LetterSpace]27 -> Vsa, v\[LetterSpace]28 -> Vda*CycA[t], v\[LetterSpace]29 -> Vse, v\[LetterSpace]3 -> Vdb*CycB[t], v\[LetterSpace]30 -> Vde*CycE[t], v\[LetterSpace]31 -> Michaelis[CycB[t], Jaie, kaie, 1 - APCP[t]], v\[LetterSpace]32 -> Michaelis[1, Jiie, kiie, APCP[t]], v\[LetterSpace]33 -> ks20p + (ks20pp*CycB[t]^n20)/(J20^n20 + CycB[t]^n20), v\[LetterSpace]34 -> kd20*Cdc20in[t], v\[LetterSpace]35 -> Michaelis[APCP[t], Ja20, ka20, Cdc20in[t]], v\[LetterSpace]36 -> Michaelis[1, Ji20, ki20, Cdc20A[t]], v\[LetterSpace]37 -> kd20*Cdc20A[t], v\[LetterSpace]38 -> Michaelis[Vah1, Jah1, 1, 1 - Cdh1[t]], v\[LetterSpace]39 -> Michaelis[Vih1, Jih1, 1, Cdh1[t]], v\[LetterSpace]4 -> V25*pB[t], v\[LetterSpace]5 -> Vwee*CycB[t], v\[LetterSpace]6 -> kassb*CKI[t]*CycB[t], v\[LetterSpace]7 -> kdissb*BCKI[t], v\[LetterSpace]8 -> Vdb*pB[t], v\[LetterSpace]9 -> kassb*CKI[t]*pB[t]
};

parameters = {
CycD0 -> 0.108, J20 -> 10.0, Ja20 -> 1.0, Ja25 -> 0.1, Jafb -> 0.1, Jafi -> 1.0, Jah1 -> 0.03, Jaie -> 0.1, Jatf -> 0.01, Jawee -> 0.05, Ji20 -> 1.0, Ji25 -> 0.1, Jifb -> 0.1, Jifi -> 2.0, Jih1 -> 0.03, Jiie -> 0.1, Jitf -> 0.01, Jiwee -> 0.05, KEZ -> 0.2, MaxMass -> 10000.0, k14di -> 12.0, k25p -> 0.01, k25pp -> 5.0, ka20 -> 1.0, ka25p -> 0.0, ka25pp -> 1.0, kafb -> 1.0, kafi -> 6.0, kah1p -> 0.02, kah1pp -> 0.8, kaie -> 0.1, kassa -> 50.0, kassb -> 60.0, kasse -> 0.0, katfapp -> 1.5, katfdpp -> 3.0, katfepp -> 0.38, katfp -> 0.0, kaweep -> 0.3, kaweepp -> 0.0, kd20 -> 0.05, kdap -> 0.01, kdapp -> 0.16, kdappp -> 0.0, kdbcpp -> 0.15, kdbhpp -> 0.4, kdbp -> 0.003, kdeapp -> 0.0, kdebpp -> 0.0, kdeepp -> 0.0, kdep -> 0.12, kdiapp -> 0.1, kdibpp -> 0.8, kdidpp -> 0.1, kdiepp -> 0.12, kdip -> 0.002, kdissa -> 0.06, kdissb -> 0.05, kdisse -> 0.0, ki20 -> 0.16, ki25p -> 0.3, ki25pp -> 0.0, kifb -> 0.15, kifibpp -> 0.05, kifip -> 0.008, kih1app -> 0.35, kih1bpp -> 0.1, kih1dpp -> 0.005, kih1epp -> 0.06, kih1p -> 0.001, kiie -> 0.15, kitfapp -> 0.0, kitfbpp -> 8.0, kitfp -> 0.75, kiweep -> 0.0, kiweepp -> 1.0, ks20p -> 0.001, ks20pp -> 1.0, ksap -> 0.0015, ksapp -> 0.01, ksbp -> 0.004, ksbpp -> 0.04, ksep -> 0.0, ksepp -> 0.15, ksip -> 0.018, ksipp -> 0.18, kweep -> 0.02, kweepp -> 0.2, mu -> 0.005776, n20 -> 1.0, sink -> 1.0, source -> 1.0, default\[LetterSpace]compartment -> 1.0
};

assignments = {
Michaelis[M1_,J1_,k1_,S1_] -> (k1*M1*S1)/(J1 + S1), GK[a_,b_,c_,d_] -> (2*a*d)/(-a + b + b*c + a*d + Sqrt[-4*a*(-a + b)*d + (-a + b + b*c + a*d)^2]), Cdc25 -> GK[ka25p + ka25pp*CycB[t], ki25p + Cdc14*ki25pp, Ja25, Ji25], Vatf -> CycD*katfdpp + katfp + katfapp*CycA[t] + katfepp*CycE[t], Vah1 -> kah1p + Cdc14*kah1pp, Vsi -> ksip + ksipp*TFI, Cdc14 -> Cdc20A[t], CycD -> CycD0*Mass[t], TFI -> GK[Cdc14*kafi, kifip + kifibpp*CycB[t], Jafi, Jifi], Vdi -> (CycD*kdidpp + kdip + kdiapp*CycA[t] + kdibpp*CycB[t] + kdiepp*CycE[t])/(1 + Cdc14*k14di), Vwee -> kweep + kweepp*Wee1, Vse -> (ksep + ksepp*TFE)*Mass[t], Vda -> kdap + (kdapp + kdappp)*Cdc20A[t] + kdappp*Cdc20in[t], V25 -> k25p + Cdc25*k25pp, Vsa -> (ksap + ksapp*TFE)*Mass[t], TFE -> GK[Vatf, Vitf, Jatf, Jitf], TFB -> GK[kafb*CycB[t], kifb, Jafb, Jifb], Vsb -> (ksbp + ksbpp*TFB)*Mass[t], Vih1 -> CycD*kih1dpp + kih1p + kih1app*CycA[t] + kih1bpp*CycB[t] + kih1epp*CycE[t], Vitf -> kitfp + kitfapp*CycA[t] + kitfbpp*CycB[t], Vdb -> kdbp + kdbcpp*Cdc20A[t] + kdbhpp*Cdh1[t], Vde -> kdep + kdeapp*CycA[t] + kdebpp*CycB[t] + kdeepp*CycE[t], Wee1 -> GK[kaweep + Cdc14*kaweepp, kiweep + kiweepp*CycB[t], Jawee, Jiwee]
};

events = {

};

speciesAnnotations = {

};

reactionAnnotations = {

};

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

(* Time evolution *)
odes = {
APCP'[t] == 1.0*v\[LetterSpace]31 -1.0*v\[LetterSpace]32, BCKI'[t] == 1.0*v\[LetterSpace]11 +1.0*v\[LetterSpace]6 -1.0*v\[LetterSpace]13 -1.0*v\[LetterSpace]12 -1.0*v\[LetterSpace]7 -1.0*v\[LetterSpace]14, CKI'[t] == 1.0*v\[LetterSpace]22 +1.0*v\[LetterSpace]13 +1.0*v\[LetterSpace]15 +1.0*v\[LetterSpace]20 +1.0*v\[LetterSpace]7 +1.0*v\[LetterSpace]24 +1.0*v\[LetterSpace]10 +1.0*v\[LetterSpace]17 +1.0*v\[LetterSpace]26 -1.0*v\[LetterSpace]18 -1.0*v\[LetterSpace]19 -1.0*v\[LetterSpace]9 -1.0*v\[LetterSpace]23 -1.0*v\[LetterSpace]6, Cdc20A'[t] == 1.0*v\[LetterSpace]35 -1.0*v\[LetterSpace]36 -1.0*v\[LetterSpace]37, Cdc20in'[t] == 1.0*v\[LetterSpace]36 +1.0*v\[LetterSpace]33 -1.0*v\[LetterSpace]34 -1.0*v\[LetterSpace]35, Cdh1'[t] == 1.0*v\[LetterSpace]38 -1.0*v\[LetterSpace]39, CycA'[t] == 1.0*v\[LetterSpace]20 +1.0*v\[LetterSpace]27 +1.0*v\[LetterSpace]21 -1.0*v\[LetterSpace]28 -1.0*v\[LetterSpace]19, CycB'[t] == 1.0*v\[LetterSpace]2 +1.0*v\[LetterSpace]4 +1.0*v\[LetterSpace]7 +1.0*v\[LetterSpace]14 -1.0*v\[LetterSpace]5 -1.0*v\[LetterSpace]3 -1.0*v\[LetterSpace]6, CycE'[t] == 1.0*v\[LetterSpace]25 +1.0*v\[LetterSpace]29 +1.0*v\[LetterSpace]24 -1.0*v\[LetterSpace]30 -1.0*v\[LetterSpace]23, Mass'[t] == 1.0*v\[LetterSpace]1 , TriA'[t] == 1.0*v\[LetterSpace]19 -1.0*v\[LetterSpace]22 -1.0*v\[LetterSpace]20 -1.0*v\[LetterSpace]21, TriE'[t] == 1.0*v\[LetterSpace]23 -1.0*v\[LetterSpace]25 -1.0*v\[LetterSpace]24 -1.0*v\[LetterSpace]26, pB'[t] == 1.0*v\[LetterSpace]16 +1.0*v\[LetterSpace]5 +1.0*v\[LetterSpace]10 -1.0*v\[LetterSpace]8 -1.0*v\[LetterSpace]4 -1.0*v\[LetterSpace]9, pBCKI'[t] == 1.0*v\[LetterSpace]12 +1.0*v\[LetterSpace]9 -1.0*v\[LetterSpace]11 -1.0*v\[LetterSpace]16 -1.0*v\[LetterSpace]15 -1.0*v\[LetterSpace]10
};
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]}]