// This model generated automatically from SBML

// unit definitions
import nsrunit;
unit conversion off;

// SBML property definitions
property sbmlRole=string;
property sbmlName=string;
property sbmlCompartment=string;

// SBML reactions
// v1: alpha => p 
// v2: Ste2 => Ste2a 
// v3: Ste2a => Ste2 
// v4: Ste2a => p 
// v5: Ste2 => p 
// v6: Gabc => GaGTP Gbc 
// v7: GaGTP => GaGDP 
// v8: GaGTP => GaGDP 
// v9: GaGDP Gbc => Gabc 
// v10: Gbc complexC => complexD 
// v11: complexD => Gbc complexC 
// v12: Ste5 Ste11 => complexA 
// v13: complexA => Ste5 Ste11 
// v14: Ste7 Fus3 => complexB 
// v15: complexB => Ste7 Fus3 
// v16: complexA complexB => complexC 
// v17: complexC => Ste11 Ste5 Ste7 Fus3 
// v18: complexD Ste20 => complexE 
// v19: complexE => complexD Ste20 
// v20: complexE 
// v21: complexE => Gbc Ste5 Ste11 Ste7 Fus3 Ste20 
// v22: complexF => complexG 
// v23: complexF => Gbc Ste5 Ste11 Ste7 Fus3 Ste20 
// v24: complexG => complexH 
// v25: complexG => Gbc Ste5 Ste11 Ste7 Fus3 Ste20 
// v26: complexH => complexI 
// v27: complexH => Gbc Ste5 Ste11 Ste7 Fus3 Ste20 
// v28: complexI => complexL Fus3PP 
// v29: complexL Fus3 => complexK 
// v30: complexK => Fus3 complexL 
// v31: complexK => complexI 
// v32: complexL => Gbc Ste5 Ste11 Ste7 Ste20 
// v33: Fus3PP => Fus3 
// v34: Ste12 Fus3PP => Ste12a 
// v35: Ste12a => Ste12 Fus3PP 
// v36: Bar1 => Bar1a 
// v37: Bar1a => Bar1 
// v38: Bar1a => Bar1aex 
// v39: Far1 => Far1PP 
// v40: Far1PP => Far1 
// v41: Far1 => Far1U 
// v42: Far1PP Gbc => complexM 
// v43: complexM => Gbc Far1PP 
// v44: complexN => Cdc28 Far1PP 
// v45: Cdc28 Far1PP => complexN 
// v46: p => Sst2 
// v47: Sst2 => p 

math main {
  realDomain time second;
  time.min=0;
  extern time.max;
  extern time.delta;

  // variable definitions
  real compartment = 1 L;
  real k1 = .03;
  real k2 = .0012;
  real k3 = .6;
  real k4 = .24;
  real k5 = .024;
  real k6 = .0036;
  real k7 = .24;
  real k8 = .033;
  real k9 = 2E3;
  real k10 = .1;
  real k11 = 5;
  real k12 = 1;
  real k13 = 3;
  real k14 = 1;
  real k15 = 3;
  real k16 = 3;
  real k17 = 100;
  real k18 = 5;
  real k19 = 1;
  real k20 = 10;
  real k21 = 5;
  real k22 = 47;
  real k23 = 5;
  real k24 = 345;
  real k25 = 5;
  real k26 = 50;
  real k27 = 5;
  real k28 = 140;
  real k29 = 10;
  real k30 = 1;
  real k31 = 250;
  real k32 = 5;
  real k33 = 50;
  real k34 = 18;
  real k35 = 10;
  real k36 = .1;
  real k37 = .1;
  real k38 = .01;
  real k39 = 18;
  real k40 = 1;
  real k41 = .02;
  real k42 = .1;
  real k43 = .01;
  real k44 = .01;
  real k45 = .1;
  real k46 = 200;
  real k47 = 1;
  real alpha(time) M;
  real p = 0 M;
  real Bar1aex(time) M;
  real Ste2(time) M;
  real Ste2a(time) M;
  real Gabc(time) M;
  real GaGTP(time) M;
  real Gbc(time) M;
  real GaGDP(time) M;
  real Sst2(time) M;
  real complexC(time) M;
  real complexD(time) M;
  real Ste5(time) M;
  real Ste11(time) M;
  real complexA(time) M;
  real Ste7(time) M;
  real Fus3(time) M;
  real complexB(time) M;
  real Ste20(time) M;
  real complexE(time) M;
  real complexF(time) M;
  real complexG(time) M;
  real complexH(time) M;
  real complexI(time) M;
  real complexL(time) M;
  real Fus3PP(time) M;
  real complexK(time) M;
  real Ste12(time) M;
  real Ste12a(time) M;
  real Bar1(time) M;
  real Bar1a(time) M;
  real Far1(time) M;
  real Far1PP(time) M;
  real Far1U(time) M;
  real Cdc28(time) M;
  real complexM(time) M;
  real complexN(time) M;
  real v1(time) katal;
  real v2(time) katal;
  real v3(time) katal;
  real v4(time) katal;
  real v5(time) katal;
  real v6(time) katal;
  real v7(time) katal;
  real v8(time) katal;
  real v9(time) katal;
  real v10(time) katal;
  real v11(time) katal;
  real v12(time) katal;
  real v13(time) katal;
  real v14(time) katal;
  real v15(time) katal;
  real v16(time) katal;
  real v17(time) katal;
  real v18(time) katal;
  real v19(time) katal;
  real v20(time) katal;
  real v21(time) katal;
  real v22(time) katal;
  real v23(time) katal;
  real v24(time) katal;
  real v25(time) katal;
  real v26(time) katal;
  real v27(time) katal;
  real v28(time) katal;
  real v29(time) katal;
  real v30(time) katal;
  real v31(time) katal;
  real v32(time) katal;
  real v33(time) katal;
  real v34(time) katal;
  real v35(time) katal;
  real v36(time) katal;
  real v37(time) katal;
  real v38(time) katal;
  real v39(time) katal;
  real v40(time) katal;
  real v41(time) katal;
  real v42(time) katal;
  real v43(time) katal;
  real v44(time) katal;
  real v45(time) katal;
  real v46(time) katal;
  real v47(time) katal;

  // equations
  when (time=time.min) alpha = 100;
  alpha:time = (-1*v1)/compartment;
  when (time=time.min) Bar1aex = 0;
  Bar1aex:time = (v38)/compartment;
  when (time=time.min) Ste2 = 1.6666667E3;
  Ste2:time = (-1*v2 + v3 + -1*v5)/compartment;
  when (time=time.min) Ste2a = 0;
  Ste2a:time = (v2 + -1*v3 + -1*v4)/compartment;
  when (time=time.min) Gabc = 1.6666667E3;
  Gabc:time = (-1*v6 + v9)/compartment;
  when (time=time.min) GaGTP = 0;
  GaGTP:time = (v6 + -1*v7 + -1*v8)/compartment;
  when (time=time.min) Gbc = 0;
  Gbc:time = (v6 + -1*v9 + -1*v10 + v11 + v21 + v23 + v25 + v27 + v32 + -1*v42 + v43)/compartment;
  when (time=time.min) GaGDP = 0;
  GaGDP:time = (v7 + v8 + -1*v9)/compartment;
  when (time=time.min) Sst2 = 0;
  Sst2:time = (v46 + -1*v47)/compartment;
  when (time=time.min) complexC = 2.3572494E2;
  complexC:time = (-1*v10 + v11 + v16 + -1*v17)/compartment;
  when (time=time.min) complexD = 0;
  complexD:time = (v10 + -1*v11 + -1*v18 + v19)/compartment;
  when (time=time.min) Ste5 = 1.5833177E2;
  Ste5:time = (-1*v12 + v13 + v17 + v21 + v23 + v25 + v27 + v32)/compartment;
  when (time=time.min) Ste11 = 1.5833177E2;
  Ste11:time = (-1*v12 + v13 + v17 + v21 + v23 + v25 + v27 + v32)/compartment;
  when (time=time.min) complexA = 1.059433E2;
  complexA:time = (v12 + -1*v13 + -1*v16)/compartment;
  when (time=time.min) Ste7 = 36.39970164;
  Ste7:time = (-1*v14 + v15 + v17 + v21 + v23 + v25 + v27 + v32)/compartment;
  when (time=time.min) Fus3 = 6.863997E2;
  Fus3:time = (-1*v14 + v15 + v17 + v21 + v23 + v25 + v27 + -1*v29 + v30 + v33)/compartment;
  when (time=time.min) complexB = 77.87536257;
  complexB:time = (v14 + -1*v15 + -1*v16)/compartment;
  when (time=time.min) Ste20 = 1E3;
  Ste20:time = (-1*v18 + v19 + v21 + v23 + v25 + v27 + v32)/compartment;
  when (time=time.min) complexE = 0;
  complexE:time = (v18 + -1*v19 + -1*v20 + -1*v21)/compartment;
  when (time=time.min) complexF = 0;
  complexF:time = (-1*v22 + -1*v23)/compartment;
  when (time=time.min) complexG = 0;
  complexG:time = (v22 + -1*v24 + -1*v25)/compartment;
  when (time=time.min) complexH = 0;
  complexH:time = (v24 + -1*v26 + -1*v27)/compartment;
  when (time=time.min) complexI = 0;
  complexI:time = (v26 + -1*v28 + v31)/compartment;
  when (time=time.min) complexL = 0;
  complexL:time = (v28 + -1*v29 + v30 + -1*v32)/compartment;
  when (time=time.min) Fus3PP = 0;
  Fus3PP:time = (v28 + -1*v33 + -1*v34 + v35)/compartment;
  when (time=time.min) complexK = 0;
  complexK:time = (v29 + -1*v30 + -1*v31)/compartment;
  when (time=time.min) Ste12 = 200;
  Ste12:time = (-1*v34 + v35)/compartment;
  when (time=time.min) Ste12a = 0;
  Ste12a:time = (v34 + -1*v35)/compartment;
  when (time=time.min) Bar1 = 200;
  Bar1:time = (-1*v36 + v37)/compartment;
  when (time=time.min) Bar1a = 0;
  Bar1a:time = (v36 + -1*v37 + -1*v38)/compartment;
  when (time=time.min) Far1 = 500;
  Far1:time = (-1*v39 + v40 + -1*v41)/compartment;
  when (time=time.min) Far1PP = 0;
  Far1PP:time = (v39 + -1*v40 + -1*v42 + v43 + v44 + -1*v45)/compartment;
  when (time=time.min) Far1U = 0;
  Far1U:time = (v41)/compartment;
  when (time=time.min) Cdc28 = 300;
  Cdc28:time = (v44 + -1*v45)/compartment;
  when (time=time.min) complexM = 0;
  complexM:time = (v42 + -1*v43)/compartment;
  when (time=time.min) complexN = 0;
  complexN:time = (-1*v44 + v45)/compartment;
  v1 = alpha*Bar1aex*k1;
  v2 = Ste2*alpha*k2;
  v3 = Ste2a*k3;
  v4 = Ste2a*k4;
  v5 = Ste2*k5;
  v6 = Ste2a*Gabc*k6;
  v7 = GaGTP*k7;
  v8 = GaGTP*Sst2*k8;
  v9 = GaGDP*Gbc*k9;
  v10 = Gbc*complexC*k10;
  v11 = complexD*k11;
  v12 = Ste5*Ste11*k12;
  v13 = complexA*k13;
  v14 = Ste7*Fus3*k14;
  v15 = complexB*k15;
  v16 = complexA*complexB*k16;
  v17 = complexC*k17;
  v18 = complexD*Ste20*k18;
  v19 = complexE*k19;
  v20 = complexE*k20;
  v21 = complexE*k21;
  v22 = complexF*k22;
  v23 = complexF*k23;
  v24 = complexG*k24;
  v25 = complexG*k25;
  v26 = complexH*k26;
  v27 = complexH*k27;
  v28 = complexI*k28;
  v29 = complexL*Fus3*k29;
  v30 = complexK*k30;
  v31 = complexK*k31;
  v32 = complexL*k32;
  v33 = Fus3PP*k33;
  v34 = Ste12*Fus3PP*k34;
  v35 = Ste12a*k35;
  v36 = Ste12a*Bar1*k36;
  v37 = Bar1a*k37;
  v38 = Bar1a*k38;
  v39 = Far1*Fus3PP*Fus3PP/(100*100+Fus3PP*Fus3PP)*k39;
  v40 = Far1PP*k40;
  v41 = Far1*Cdc28*k41;
  v42 = Gbc*Far1PP*k42;
  v43 = complexM*k43;
  v44 = complexN*k44;
  v45 = Far1PP*Cdc28*k45;
  v46 = Fus3PP^2/(4^2+Fus3PP^2)*k46;
  v47 = Sst2*k47;

  // variable properties
  compartment.sbmlRole="compartment";
  k1.sbmlRole="parameter";
  k2.sbmlRole="parameter";
  k3.sbmlRole="parameter";
  k4.sbmlRole="parameter";
  k5.sbmlRole="parameter";
  k6.sbmlRole="parameter";
  k7.sbmlRole="parameter";
  k8.sbmlRole="parameter";
  k9.sbmlRole="parameter";
  k10.sbmlRole="parameter";
  k11.sbmlRole="parameter";
  k12.sbmlRole="parameter";
  k13.sbmlRole="parameter";
  k14.sbmlRole="parameter";
  k15.sbmlRole="parameter";
  k16.sbmlRole="parameter";
  k17.sbmlRole="parameter";
  k18.sbmlRole="parameter";
  k19.sbmlRole="parameter";
  k20.sbmlRole="parameter";
  k21.sbmlRole="parameter";
  k22.sbmlRole="parameter";
  k23.sbmlRole="parameter";
  k24.sbmlRole="parameter";
  k25.sbmlRole="parameter";
  k26.sbmlRole="parameter";
  k27.sbmlRole="parameter";
  k28.sbmlRole="parameter";
  k29.sbmlRole="parameter";
  k30.sbmlRole="parameter";
  k31.sbmlRole="parameter";
  k32.sbmlRole="parameter";
  k33.sbmlRole="parameter";
  k34.sbmlRole="parameter";
  k35.sbmlRole="parameter";
  k36.sbmlRole="parameter";
  k37.sbmlRole="parameter";
  k38.sbmlRole="parameter";
  k39.sbmlRole="parameter";
  k40.sbmlRole="parameter";
  k41.sbmlRole="parameter";
  k42.sbmlRole="parameter";
  k43.sbmlRole="parameter";
  k44.sbmlRole="parameter";
  k45.sbmlRole="parameter";
  k46.sbmlRole="parameter";
  k47.sbmlRole="parameter";
  alpha.sbmlRole="species";
  alpha.sbmlCompartment="compartment";
  p.sbmlRole="species";
  p.sbmlCompartment="compartment";
  Bar1aex.sbmlRole="species";
  Bar1aex.sbmlCompartment="compartment";
  Ste2.sbmlRole="species";
  Ste2.sbmlCompartment="compartment";
  Ste2a.sbmlRole="species";
  Ste2a.sbmlCompartment="compartment";
  Gabc.sbmlRole="species";
  Gabc.sbmlCompartment="compartment";
  GaGTP.sbmlRole="species";
  GaGTP.sbmlCompartment="compartment";
  Gbc.sbmlRole="species";
  Gbc.sbmlCompartment="compartment";
  GaGDP.sbmlRole="species";
  GaGDP.sbmlCompartment="compartment";
  Sst2.sbmlRole="species";
  Sst2.sbmlCompartment="compartment";
  complexC.sbmlRole="species";
  complexC.sbmlCompartment="compartment";
  complexD.sbmlRole="species";
  complexD.sbmlCompartment="compartment";
  Ste5.sbmlRole="species";
  Ste5.sbmlCompartment="compartment";
  Ste11.sbmlRole="species";
  Ste11.sbmlCompartment="compartment";
  complexA.sbmlRole="species";
  complexA.sbmlCompartment="compartment";
  Ste7.sbmlRole="species";
  Ste7.sbmlCompartment="compartment";
  Fus3.sbmlRole="species";
  Fus3.sbmlCompartment="compartment";
  complexB.sbmlRole="species";
  complexB.sbmlCompartment="compartment";
  Ste20.sbmlRole="species";
  Ste20.sbmlCompartment="compartment";
  complexE.sbmlRole="species";
  complexE.sbmlCompartment="compartment";
  complexF.sbmlRole="species";
  complexF.sbmlCompartment="compartment";
  complexG.sbmlRole="species";
  complexG.sbmlCompartment="compartment";
  complexH.sbmlRole="species";
  complexH.sbmlCompartment="compartment";
  complexI.sbmlRole="species";
  complexI.sbmlCompartment="compartment";
  complexL.sbmlRole="species";
  complexL.sbmlCompartment="compartment";
  Fus3PP.sbmlRole="species";
  Fus3PP.sbmlCompartment="compartment";
  complexK.sbmlRole="species";
  complexK.sbmlCompartment="compartment";
  Ste12.sbmlRole="species";
  Ste12.sbmlCompartment="compartment";
  Ste12a.sbmlRole="species";
  Ste12a.sbmlCompartment="compartment";
  Bar1.sbmlRole="species";
  Bar1.sbmlCompartment="compartment";
  Bar1a.sbmlRole="species";
  Bar1a.sbmlCompartment="compartment";
  Far1.sbmlRole="species";
  Far1.sbmlCompartment="compartment";
  Far1PP.sbmlRole="species";
  Far1PP.sbmlCompartment="compartment";
  Far1U.sbmlRole="species";
  Far1U.sbmlCompartment="compartment";
  Cdc28.sbmlRole="species";
  Cdc28.sbmlCompartment="compartment";
  complexM.sbmlRole="species";
  complexM.sbmlCompartment="compartment";
  complexN.sbmlRole="species";
  complexN.sbmlCompartment="compartment";
  v1.sbmlRole="rate";
  v2.sbmlRole="rate";
  v3.sbmlRole="rate";
  v4.sbmlRole="rate";
  v5.sbmlRole="rate";
  v6.sbmlRole="rate";
  v7.sbmlRole="rate";
  v8.sbmlRole="rate";
  v9.sbmlRole="rate";
  v10.sbmlRole="rate";
  v11.sbmlRole="rate";
  v12.sbmlRole="rate";
  v13.sbmlRole="rate";
  v14.sbmlRole="rate";
  v15.sbmlRole="rate";
  v16.sbmlRole="rate";
  v17.sbmlRole="rate";
  v18.sbmlRole="rate";
  v19.sbmlRole="rate";
  v20.sbmlRole="rate";
  v21.sbmlRole="rate";
  v22.sbmlRole="rate";
  v23.sbmlRole="rate";
  v24.sbmlRole="rate";
  v25.sbmlRole="rate";
  v26.sbmlRole="rate";
  v27.sbmlRole="rate";
  v28.sbmlRole="rate";
  v29.sbmlRole="rate";
  v30.sbmlRole="rate";
  v31.sbmlRole="rate";
  v32.sbmlRole="rate";
  v33.sbmlRole="rate";
  v34.sbmlRole="rate";
  v35.sbmlRole="rate";
  v36.sbmlRole="rate";
  v37.sbmlRole="rate";
  v38.sbmlRole="rate";
  v39.sbmlRole="rate";
  v40.sbmlRole="rate";
  v41.sbmlRole="rate";
  v42.sbmlRole="rate";
  v43.sbmlRole="rate";
  v44.sbmlRole="rate";
  v45.sbmlRole="rate";
  v46.sbmlRole="rate";
  v47.sbmlRole="rate";
}