(* 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 = {
FZYa[t], IEa\[LetterSpace]1[t], MPFc[t], MPFn[t], N[t], StgPc[t], StgPn[t], Stgc[t], Stgm[t], Stgn[t], Wee1Pn[t], Wee1c[t], Wee1n[t], Xm[t], Xp[t], preMPFc[t], preMPFn[t]
};

initialValues = {
FZYa[0] == 0.0, IEa\[LetterSpace]1[0] == 0.0, MPFc[0] == 1.0, MPFn[0] == 0.0, N[0] == 1.0, StgPc[0] == 0.8, StgPn[0] == 0.0, Stgc[0] == 0.0, Stgm[0] == 1.0, Stgn[0] == 0.0, Wee1Pn[0] == 0.0, Wee1c[0] == 0.0, Wee1n[0] == 0.0, Xm[0] == 0.0, Xp[0] == 0.0, preMPFc[0] == 0.0, preMPFn[0] == 0.0
};

rates = {
Nuclei\[LetterSpace]1, R\[LetterSpace]1, R\[LetterSpace]10, R\[LetterSpace]12, R\[LetterSpace]13, R\[LetterSpace]14, R\[LetterSpace]15, R\[LetterSpace]16, R\[LetterSpace]19, R\[LetterSpace]2, R\[LetterSpace]20, R\[LetterSpace]21, R\[LetterSpace]22, R\[LetterSpace]23, R\[LetterSpace]24, R\[LetterSpace]25, R\[LetterSpace]26, R\[LetterSpace]27, R\[LetterSpace]28, R\[LetterSpace]29, R\[LetterSpace]3, R\[LetterSpace]30, R\[LetterSpace]33, R\[LetterSpace]34, R\[LetterSpace]37, R\[LetterSpace]38, R\[LetterSpace]39, R\[LetterSpace]40, R\[LetterSpace]41, R\[LetterSpace]42, R\[LetterSpace]43, R\[LetterSpace]44, R\[LetterSpace]45, R\[LetterSpace]46, R\[LetterSpace]6, R\[LetterSpace]7, R\[LetterSpace]8, R\[LetterSpace]9, \[LetterSpace]16, \[LetterSpace]18, \[LetterSpace]182\[LetterSpace]1, \[LetterSpace]20, \[LetterSpace]22, \[LetterSpace]24, \[LetterSpace]26, \[LetterSpace]28, \[LetterSpace]30, \[LetterSpace]32, \[LetterSpace]50, \[LetterSpace]51, exportofMPFfromnucleus\[LetterSpace]1, exportofpreMPFnfromnucleus\[LetterSpace]1, importofMPFintocytoplasm\[LetterSpace]1, importofpreMPFintocytoplaslm\[LetterSpace]1
};

rateEquations = {
Nuclei\[LetterSpace]1 -> 0, R\[LetterSpace]1 -> cytoplasm*ksc, R\[LetterSpace]10 -> nuclei*((kdmp*Stgm[t])/(Jm + Stgm[t]) + kdm*Stgm[t]*Xp[t]), R\[LetterSpace]12 -> cytoplasm*ksstg*Stgm[t], R\[LetterSpace]13 -> (cytoplasm*(kastgp + kastg*MPFc[t])*Stgc[t])/(Jastg + Stgc[t]), R\[LetterSpace]14 -> (cytoplasm*kistg*StgPc[t])/(Jistg + StgPc[t]), R\[LetterSpace]15 -> cytoplasm*kdstg*Stgc[t], R\[LetterSpace]16 -> cytoplasm*kdstg*StgPc[t], R\[LetterSpace]19 -> (cytoplasm*E\[LetterSpace]1*kt*MPFc[t]*N[t])/(1 - E\[LetterSpace]1*N[t]), R\[LetterSpace]2 -> cytoplasm*preMPFc[t]*(kstgp + kstg*StgPc[t]), R\[LetterSpace]20 -> (cytoplasm*E\[LetterSpace]1*kt*N[t]*preMPFc[t])/(1 - E\[LetterSpace]1*N[t]), R\[LetterSpace]21 -> koutw\[LetterSpace]1*nuclei*Wee1Pn[t], R\[LetterSpace]22 -> (cytoplasm*E\[LetterSpace]1*kinw\[LetterSpace]1*Wee1Pc*N[t])/(1 - E\[LetterSpace]1*N[t]), R\[LetterSpace]23 -> koutw\[LetterSpace]1*nuclei*Wee1n[t], R\[LetterSpace]24 -> (cytoplasm*E\[LetterSpace]1*kinw\[LetterSpace]1*N[t]*Wee1c[t])/(1 - E\[LetterSpace]1*N[t]), R\[LetterSpace]25 -> kouts\[LetterSpace]1*nuclei*StgPn[t], R\[LetterSpace]26 -> (cytoplasm*E\[LetterSpace]1*kins\[LetterSpace]1*N[t]*StgPc[t])/(1 - E\[LetterSpace]1*N[t]), R\[LetterSpace]27 -> kouts\[LetterSpace]1*nuclei*Stgn[t], R\[LetterSpace]28 -> (cytoplasm*E\[LetterSpace]1*kins\[LetterSpace]1*N[t]*Stgc[t])/(1 - E\[LetterSpace]1*N[t]), R\[LetterSpace]29 -> nuclei*preMPFn[t]*(kstgp + kstg*StgPn[t]), R\[LetterSpace]3 -> cytoplasm*MPFc[t]*(kweep + kwee*Wee1c[t]), R\[LetterSpace]30 -> nuclei*MPFn[t]*(kweep + kwee*Wee1n[t]), R\[LetterSpace]33 -> nuclei*(kdnp + kdn*FZYa[t])*preMPFn[t], R\[LetterSpace]34 -> nuclei*(kdnp + kdn*FZYa[t])*MPFn[t], R\[LetterSpace]37 -> (kiie*nuclei*IEa\[LetterSpace]1[t])/(Jiie + IEa\[LetterSpace]1[t]), R\[LetterSpace]38 -> (kaie*nuclei*(1 - IEa\[LetterSpace]1[t])*MPFn[t])/(1 + Jaie - IEa\[LetterSpace]1[t]), R\[LetterSpace]39 -> (kafzy*nuclei*(1 - FZYa[t])*IEa\[LetterSpace]1[t])/(1 + Jafzy - FZYa[t]), R\[LetterSpace]40 -> (kifzy*nuclei*FZYa[t])/(Jifzy + FZYa[t]), R\[LetterSpace]41 -> (nuclei*(kiweep + kiwee*MPFn[t])*Wee1n[t])/(Jiwee + Wee1n[t]), R\[LetterSpace]42 -> (kawee*nuclei*Wee1Pn[t])/(Jawee + Wee1Pn[t]), R\[LetterSpace]43 -> (nuclei*(kastgp + kastg*MPFn[t])*Stgn[t])/(Jastg + Stgn[t]), R\[LetterSpace]44 -> (kistg*nuclei*StgPn[t])/(Jistg + StgPn[t]), R\[LetterSpace]45 -> kdstg*nuclei*Stgn[t], R\[LetterSpace]46 -> kdstg*nuclei*StgPn[t], R\[LetterSpace]6 -> cytoplasm*kdc*preMPFc[t], R\[LetterSpace]7 -> cytoplasm*kdc*MPFc[t], R\[LetterSpace]8 -> (cytoplasm*(kiweep + kiwee*MPFc[t])*Wee1c[t])/(Jiwee + Wee1c[t]), R\[LetterSpace]9 -> (cytoplasm*kawee*Wee1Pc)/(Jawee + Wee1Pc), \[LetterSpace]16 -> cytoplasm*kt*MPFc[t], \[LetterSpace]18 -> cytoplasm*kt*preMPFc[t], \[LetterSpace]182\[LetterSpace]1 -> (E\[LetterSpace]1*koutw\[LetterSpace]1*nuclei*N[t]*Wee1Pn[t])/(1 - E\[LetterSpace]1*N[t]), \[LetterSpace]20 -> cytoplasm*kinw\[LetterSpace]1*Wee1Pc, \[LetterSpace]22 -> (E\[LetterSpace]1*koutw\[LetterSpace]1*nuclei*N[t]*Wee1n[t])/(1 - E\[LetterSpace]1*N[t]), \[LetterSpace]24 -> cytoplasm*kinw\[LetterSpace]1*Wee1c[t], \[LetterSpace]26 -> (E\[LetterSpace]1*kouts\[LetterSpace]1*nuclei*N[t]*StgPn[t])/(1 - E\[LetterSpace]1*N[t]), \[LetterSpace]28 -> cytoplasm*kins\[LetterSpace]1*StgPc[t], \[LetterSpace]30 -> (E\[LetterSpace]1*kouts\[LetterSpace]1*nuclei*N[t]*Stgn[t])/(1 - E\[LetterSpace]1*N[t]), \[LetterSpace]32 -> cytoplasm*kins\[LetterSpace]1*Stgc[t], \[LetterSpace]50 -> ksxm\[LetterSpace]1*nuclei*N[t], \[LetterSpace]51 -> cytoplasm*ksxp\[LetterSpace]1*Xm[t], exportofMPFfromnucleus\[LetterSpace]1 -> kout\[LetterSpace]1*nuclei*MPFn[t], exportofpreMPFnfromnucleus\[LetterSpace]1 -> kout\[LetterSpace]1*nuclei*preMPFn[t], importofMPFintocytoplasm\[LetterSpace]1 -> (E\[LetterSpace]1*kout\[LetterSpace]1*nuclei*MPFn[t]*N[t])/(1 - E\[LetterSpace]1*N[t]), importofpreMPFintocytoplaslm\[LetterSpace]1 -> (E\[LetterSpace]1*kout\[LetterSpace]1*nuclei*N[t]*preMPFn[t])/(1 - E\[LetterSpace]1*N[t])
};

parameters = {
E\[LetterSpace]1 -> 7*^-05, Jafzy -> 0.01, Jaie -> 0.01, Jastg -> 0.05, Jawee -> 0.05, Jifzy -> 0.01, Jiie -> 0.01, Jistg -> 0.05, Jiwee -> 0.05, Jm -> 0.05, StringT -> 0.8, Wee1T -> 0.8, factor\[LetterSpace]1 -> 1.95, kafzy -> 1.0, kaie -> 1.0, kastg -> 1.0, kastgp -> 0.0, kawee -> 0.3, kdc -> 0.01, kdm -> 0.2, kdmp -> 0.002, kdn -> 1.5, kdnp -> 0.01, kdstg -> 0.0, kez\[LetterSpace]1 -> 0.5, kifzy -> 0.2, kiie -> 0.4, kins\[LetterSpace]1 -> 0.08, kinw\[LetterSpace]1 -> 0.04, kistg -> 0.3, kiwee -> 1.0, kiweep -> 0.01, kout\[LetterSpace]1 -> 0.0, kouts\[LetterSpace]1 -> 0.02, koutw\[LetterSpace]1 -> 0.01, ksc -> 0.01, ksstg -> 0.0, kstg -> 2.0, kstgp -> 0.2, ksxm\[LetterSpace]1 -> 0.0005, ksxp\[LetterSpace]1 -> 0.001, kt -> 0.15, kwee -> 1.0, kweep -> 0.005, cytoplasm -> 1.0, nuclei -> 1.0
};

assignments = {
Wee1Pc -> -Wee1c[t] + (Wee1T - E\[LetterSpace]1*N[t]*(Wee1n[t] + Wee1Pn[t]))/(1 - E\[LetterSpace]1*N[t]), StgPT -> (1 - E\[LetterSpace]1*N[t])*StgPc[t] + E\[LetterSpace]1*N[t]*StgPn[t], CycBT -> (1 - E\[LetterSpace]1*N[t])*(MPFc[t] + preMPFc[t]) + E\[LetterSpace]1*N[t]*(MPFn[t] + preMPFn[t])
};

events = {

};

speciesAnnotations = {
FZYa[t]->"http://identifiers.org/uniprot/Q24044", FZYa[t]->"http://identifiers.org/interpro/IPR000002", StgPc[t]->"http://identifiers.org/uniprot/P20483", StgPn[t]->"http://identifiers.org/uniprot/P20483", StgPn[t]->"http://identifiers.org/interpro/IPR000751", Stgc[t]->"http://identifiers.org/uniprot/P20483", Stgn[t]->"http://identifiers.org/uniprot/P20483", Stgn[t]->"http://identifiers.org/interpro/IPR000751", Wee1Pc[t]->"http://identifiers.org/uniprot/P54350", Wee1Pc[t]->"http://identifiers.org/uniprot/P30291", Wee1Pn[t]->"http://identifiers.org/uniprot/P54350", Wee1Pn[t]->"http://identifiers.org/uniprot/P30291", Wee1c[t]->"http://identifiers.org/uniprot/P54350", Wee1c[t]->"http://identifiers.org/uniprot/P30291", Wee1n[t]->"http://identifiers.org/uniprot/P54350", Wee1n[t]->"http://identifiers.org/uniprot/P30291"
};

reactionAnnotations = {

};

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

(* Time evolution *)
odes = {
FZYa'[t] == 1.0*R\[LetterSpace]39 -1.0*R\[LetterSpace]40, IEa\[LetterSpace]1'[t] == 1.0*R\[LetterSpace]38 -1.0*R\[LetterSpace]37, MPFc'[t] == 1.0*R\[LetterSpace]1 +1.0*R\[LetterSpace]2 +1.0*importofMPFintocytoplasm\[LetterSpace]1 -1.0*R\[LetterSpace]3 -1.0*R\[LetterSpace]7 -1.0*R\[LetterSpace]19, MPFn'[t] == 1.0*\[LetterSpace]16 +1.0*R\[LetterSpace]29 -1.0*exportofMPFfromnucleus\[LetterSpace]1 -1.0*R\[LetterSpace]30 -1.0*R\[LetterSpace]34, N'[t] == 1.0*Nuclei\[LetterSpace]1 , StgPc'[t] == 1.0*R\[LetterSpace]13 +1.0*\[LetterSpace]26 -1.0*R\[LetterSpace]14 -1.0*R\[LetterSpace]16 -1.0*R\[LetterSpace]26, StgPn'[t] == 1.0*\[LetterSpace]28 +1.0*R\[LetterSpace]43 -1.0*R\[LetterSpace]25 -1.0*R\[LetterSpace]44 -1.0*R\[LetterSpace]46, Stgc'[t] == 1.0*R\[LetterSpace]12 +1.0*R\[LetterSpace]14 +1.0*\[LetterSpace]30 -1.0*R\[LetterSpace]13 -1.0*R\[LetterSpace]15 -1.0*R\[LetterSpace]28, Stgm'[t] ==  -1.0*R\[LetterSpace]10, Stgn'[t] == 1.0*\[LetterSpace]32 +1.0*R\[LetterSpace]44 -1.0*R\[LetterSpace]27 -1.0*R\[LetterSpace]43 -1.0*R\[LetterSpace]45, Wee1Pn'[t] == 1.0*\[LetterSpace]20 +1.0*R\[LetterSpace]41 -1.0*R\[LetterSpace]21 -1.0*R\[LetterSpace]42, Wee1c'[t] == 1.0*R\[LetterSpace]9 +1.0*\[LetterSpace]22 -1.0*R\[LetterSpace]8 -1.0*R\[LetterSpace]24, Wee1n'[t] == 1.0*\[LetterSpace]24 +1.0*R\[LetterSpace]42 -1.0*R\[LetterSpace]23 -1.0*R\[LetterSpace]41, Xm'[t] == 1.0*\[LetterSpace]50 , Xp'[t] == 1.0*\[LetterSpace]51 , preMPFc'[t] == 1.0*R\[LetterSpace]3 +1.0*importofpreMPFintocytoplaslm\[LetterSpace]1 -1.0*R\[LetterSpace]2 -1.0*R\[LetterSpace]6 -1.0*R\[LetterSpace]20, preMPFn'[t] == 1.0*\[LetterSpace]18 +1.0*R\[LetterSpace]30 -1.0*exportofpreMPFnfromnucleus\[LetterSpace]1 -1.0*R\[LetterSpace]29 -1.0*R\[LetterSpace]33
};
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]}]