(* 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 = {
Bar1[t], Bar1a[t], Bar1aex[t], Cdc28[t], Far1[t], Far1PP[t], Far1U[t], Fus3[t], Fus3PP[t], GaGDP[t], GaGTP[t], Gabc[t], Gbc[t], Sst2[t], Ste11[t], Ste12[t], Ste12a[t], Ste2[t], Ste20[t], Ste2a[t], Ste5[t], Ste7[t], alpha[t], complexA[t], complexB[t], complexC[t], complexD[t], complexE[t], complexF[t], complexG[t], complexH[t], complexI[t], complexK[t], complexL[t], complexM[t], complexN[t]
};

initialValues = {
Bar1[0] == 200.0, Bar1a[0] == 0.0, Bar1aex[0] == 0.0, Cdc28[0] == 300.0, Far1[0] == 500.0, Far1PP[0] == 0.0, Far1U[0] == 0.0, Fus3[0] == 686.399701640513, Fus3PP[0] == 0.0, GaGDP[0] == 0.0, GaGTP[0] == 0.0, Gabc[0] == 1666.6666667, Gbc[0] == 0.0, Sst2[0] == 0.0, Ste11[0] == 158.33176608789, Ste12[0] == 200.0, Ste12a[0] == 0.0, Ste2[0] == 1666.6666667, Ste20[0] == 1000.0, Ste2a[0] == 0.0, Ste5[0] == 158.33176608789, Ste7[0] == 36.3997016405141, alpha[0] == 100.0, complexA[0] == 105.943298120207, complexB[0] == 77.8753625675829, complexC[0] == 235.724935791903, complexD[0] == 0.0, complexE[0] == 0.0, complexF[0] == 0.0, complexG[0] == 0.0, complexH[0] == 0.0, complexI[0] == 0.0, complexK[0] == 0.0, complexL[0] == 0.0, complexM[0] == 0.0, complexN[0] == 0.0
};

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]40, v\[LetterSpace]41, v\[LetterSpace]42, v\[LetterSpace]43, v\[LetterSpace]44, v\[LetterSpace]45, v\[LetterSpace]46, v\[LetterSpace]47, v\[LetterSpace]5, v\[LetterSpace]6, v\[LetterSpace]7, v\[LetterSpace]8, v\[LetterSpace]9
};

rateEquations = {
v\[LetterSpace]1 -> k1*alpha[t]*Bar1aex[t], v\[LetterSpace]10 -> k10*complexC[t]*Gbc[t], v\[LetterSpace]11 -> k11*complexD[t], v\[LetterSpace]12 -> k12*Ste11[t]*Ste5[t], v\[LetterSpace]13 -> k13*complexA[t], v\[LetterSpace]14 -> k14*Fus3[t]*Ste7[t], v\[LetterSpace]15 -> k15*complexB[t], v\[LetterSpace]16 -> k16*complexA[t]*complexB[t], v\[LetterSpace]17 -> k17*complexC[t], v\[LetterSpace]18 -> k18*complexD[t]*Ste20[t], v\[LetterSpace]19 -> k19*complexE[t], v\[LetterSpace]2 -> k2*alpha[t]*Ste2[t], v\[LetterSpace]20 -> k20*complexE[t], v\[LetterSpace]21 -> k21*complexE[t], v\[LetterSpace]22 -> k22*complexF[t], v\[LetterSpace]23 -> k23*complexF[t], v\[LetterSpace]24 -> k24*complexG[t], v\[LetterSpace]25 -> k25*complexG[t], v\[LetterSpace]26 -> k26*complexH[t], v\[LetterSpace]27 -> k27*complexH[t], v\[LetterSpace]28 -> k28*complexI[t], v\[LetterSpace]29 -> k29*complexL[t]*Fus3[t], v\[LetterSpace]3 -> k3*Ste2a[t], v\[LetterSpace]30 -> k30*complexK[t], v\[LetterSpace]31 -> k31*complexK[t], v\[LetterSpace]32 -> k32*complexL[t], v\[LetterSpace]33 -> k33*Fus3PP[t], v\[LetterSpace]34 -> k34*Fus3PP[t]*Ste12[t], v\[LetterSpace]35 -> k35*Ste12a[t], v\[LetterSpace]36 -> k36*Bar1[t]*Ste12a[t], v\[LetterSpace]37 -> k37*Bar1a[t], v\[LetterSpace]38 -> k38*Bar1a[t], v\[LetterSpace]39 -> (k39*Far1[t]*Fus3PP[t]^2)/(10000 + Fus3PP[t]^2), v\[LetterSpace]4 -> k4*Ste2a[t], v\[LetterSpace]40 -> k40*Far1PP[t], v\[LetterSpace]41 -> k41*Cdc28[t]*Far1[t], v\[LetterSpace]42 -> k42*Far1PP[t]*Gbc[t], v\[LetterSpace]43 -> k43*complexM[t], v\[LetterSpace]44 -> k44*complexN[t], v\[LetterSpace]45 -> k45*Cdc28[t]*Far1PP[t], v\[LetterSpace]46 -> (k46*Fus3PP[t]^2)/(16 + Fus3PP[t]^2), v\[LetterSpace]47 -> k47*Sst2[t], v\[LetterSpace]5 -> k5*Ste2[t], v\[LetterSpace]6 -> k6*Gabc[t]*Ste2a[t], v\[LetterSpace]7 -> k7*GaGTP[t], v\[LetterSpace]8 -> k8*GaGTP[t]*Sst2[t], v\[LetterSpace]9 -> k9*GaGDP[t]*Gbc[t]
};

parameters = {
k1 -> 0.03, k10 -> 0.1, k11 -> 5.0, k12 -> 1.0, k13 -> 3.0, k14 -> 1.0, k15 -> 3.0, k16 -> 3.0, k17 -> 100.0, k18 -> 5.0, k19 -> 1.0, k2 -> 0.0012, k20 -> 10.0, k21 -> 5.0, k22 -> 47.0, k23 -> 5.0, k24 -> 345.0, k25 -> 5.0, k26 -> 50.0, k27 -> 5.0, k28 -> 140.0, k29 -> 10.0, k3 -> 0.6, k30 -> 1.0, k31 -> 250.0, k32 -> 5.0, k33 -> 50.0, k34 -> 18.0, k35 -> 10.0, k36 -> 0.1, k37 -> 0.1, k38 -> 0.01, k39 -> 18.0, k4 -> 0.24, k40 -> 1.0, k41 -> 0.02, k42 -> 0.1, k43 -> 0.01, k44 -> 0.01, k45 -> 0.1, k46 -> 200.0, k47 -> 1.0, k5 -> 0.024, k6 -> 0.0036, k7 -> 0.24, k8 -> 0.033, k9 -> 2000.0, p -> 0.0, default\[LetterSpace]compartment -> 1.0
};

assignments = {

};

events = {

};

speciesAnnotations = {

};

reactionAnnotations = {

};

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

(* Time evolution *)
odes = {
Bar1'[t] == 1.0*v\[LetterSpace]37 -1.0*v\[LetterSpace]36, Bar1a'[t] == 1.0*v\[LetterSpace]36 -1.0*v\[LetterSpace]38 -1.0*v\[LetterSpace]37, Bar1aex'[t] == 1.0*v\[LetterSpace]38 , Cdc28'[t] == 1.0*v\[LetterSpace]44 -1.0*v\[LetterSpace]45, Far1'[t] == 1.0*v\[LetterSpace]40 -1.0*v\[LetterSpace]39 -1.0*v\[LetterSpace]41, Far1PP'[t] == 1.0*v\[LetterSpace]39 +1.0*v\[LetterSpace]44 +1.0*v\[LetterSpace]43 -1.0*v\[LetterSpace]40 -1.0*v\[LetterSpace]42 -1.0*v\[LetterSpace]45, Far1U'[t] == 1.0*v\[LetterSpace]41 , Fus3'[t] == 1.0*v\[LetterSpace]23 +1.0*v\[LetterSpace]17 +1.0*v\[LetterSpace]15 +1.0*v\[LetterSpace]33 +1.0*v\[LetterSpace]21 +1.0*v\[LetterSpace]30 +1.0*v\[LetterSpace]25 +1.0*v\[LetterSpace]27 -1.0*v\[LetterSpace]29 -1.0*v\[LetterSpace]14, Fus3PP'[t] == 1.0*v\[LetterSpace]35 +1.0*v\[LetterSpace]28 -1.0*v\[LetterSpace]33 -1.0*v\[LetterSpace]34, GaGDP'[t] == 1.0*v\[LetterSpace]7 +1.0*v\[LetterSpace]8 -1.0*v\[LetterSpace]9, GaGTP'[t] == 1.0*v\[LetterSpace]6 -1.0*v\[LetterSpace]7 -1.0*v\[LetterSpace]8, Gabc'[t] == 1.0*v\[LetterSpace]9 -1.0*v\[LetterSpace]6, Gbc'[t] == 1.0*v\[LetterSpace]23 +1.0*v\[LetterSpace]32 +1.0*v\[LetterSpace]11 +1.0*v\[LetterSpace]6 +1.0*v\[LetterSpace]43 +1.0*v\[LetterSpace]21 +1.0*v\[LetterSpace]25 +1.0*v\[LetterSpace]27 -1.0*v\[LetterSpace]10 -1.0*v\[LetterSpace]9 -1.0*v\[LetterSpace]42, Sst2'[t] == 1.0*v\[LetterSpace]46 -1.0*v\[LetterSpace]47, Ste11'[t] == 1.0*v\[LetterSpace]23 +1.0*v\[LetterSpace]17 +1.0*v\[LetterSpace]32 +1.0*v\[LetterSpace]13 +1.0*v\[LetterSpace]21 +1.0*v\[LetterSpace]25 +1.0*v\[LetterSpace]27 -1.0*v\[LetterSpace]12, Ste12'[t] == 1.0*v\[LetterSpace]35 -1.0*v\[LetterSpace]34, Ste12a'[t] == 1.0*v\[LetterSpace]34 -1.0*v\[LetterSpace]35, Ste2'[t] == 1.0*v\[LetterSpace]3 -1.0*v\[LetterSpace]5 -1.0*v\[LetterSpace]2, Ste20'[t] == 1.0*v\[LetterSpace]23 +1.0*v\[LetterSpace]32 +1.0*v\[LetterSpace]21 +1.0*v\[LetterSpace]19 +1.0*v\[LetterSpace]25 +1.0*v\[LetterSpace]27 -1.0*v\[LetterSpace]18, Ste2a'[t] == 1.0*v\[LetterSpace]2 -1.0*v\[LetterSpace]3 -1.0*v\[LetterSpace]4, Ste5'[t] == 1.0*v\[LetterSpace]23 +1.0*v\[LetterSpace]17 +1.0*v\[LetterSpace]32 +1.0*v\[LetterSpace]13 +1.0*v\[LetterSpace]21 +1.0*v\[LetterSpace]25 +1.0*v\[LetterSpace]27 -1.0*v\[LetterSpace]12, Ste7'[t] == 1.0*v\[LetterSpace]23 +1.0*v\[LetterSpace]17 +1.0*v\[LetterSpace]32 +1.0*v\[LetterSpace]15 +1.0*v\[LetterSpace]21 +1.0*v\[LetterSpace]25 +1.0*v\[LetterSpace]27 -1.0*v\[LetterSpace]14, alpha'[t] ==  -1.0*v\[LetterSpace]1, complexA'[t] == 1.0*v\[LetterSpace]12 -1.0*v\[LetterSpace]16 -1.0*v\[LetterSpace]13, complexB'[t] == 1.0*v\[LetterSpace]14 -1.0*v\[LetterSpace]16 -1.0*v\[LetterSpace]15, complexC'[t] == 1.0*v\[LetterSpace]16 +1.0*v\[LetterSpace]11 -1.0*v\[LetterSpace]17 -1.0*v\[LetterSpace]10, complexD'[t] == 1.0*v\[LetterSpace]10 +1.0*v\[LetterSpace]19 -1.0*v\[LetterSpace]11 -1.0*v\[LetterSpace]18, complexE'[t] == 1.0*v\[LetterSpace]18 -1.0*v\[LetterSpace]20 -1.0*v\[LetterSpace]21 -1.0*v\[LetterSpace]19, complexF'[t] == 1.0*v\[LetterSpace]20 -1.0*v\[LetterSpace]23 -1.0*v\[LetterSpace]22, complexG'[t] == 1.0*v\[LetterSpace]22 -1.0*v\[LetterSpace]24 -1.0*v\[LetterSpace]25, complexH'[t] == 1.0*v\[LetterSpace]24 -1.0*v\[LetterSpace]26 -1.0*v\[LetterSpace]27, complexI'[t] == 1.0*v\[LetterSpace]31 +1.0*v\[LetterSpace]26 -1.0*v\[LetterSpace]28, complexK'[t] == 1.0*v\[LetterSpace]29 -1.0*v\[LetterSpace]31 -1.0*v\[LetterSpace]30, complexL'[t] == 1.0*v\[LetterSpace]28 +1.0*v\[LetterSpace]30 -1.0*v\[LetterSpace]29 -1.0*v\[LetterSpace]32, complexM'[t] == 1.0*v\[LetterSpace]42 -1.0*v\[LetterSpace]43, complexN'[t] == 1.0*v\[LetterSpace]45 -1.0*v\[LetterSpace]44
};
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]}]