// This model generated automatically from SBML
// WARNING: Renamed variable source to source_ due to JSim namespace restrictions

// unit definitions
import nsrunit;
unit conversion off;

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

// SBML reactions
// v1: source => Mass 
// v2: source => CycB 
// v3: CycB => sink 
// v4: pB => CycB 
// v5: CycB => pB 
// v6: CycB CKI => BCKI 
// v7: BCKI => CycB CKI 
// v8: pB => sink 
// v9: pB CKI => pBCKI 
// v10: pBCKI => pB CKI 
// v11: pBCKI => BCKI 
// v12: BCKI => pBCKI 
// v13: BCKI => CKI 
// v14: BCKI => CycB 
// v15: pBCKI => CKI 
// v16: pBCKI => pB 
// v17: source => CKI 
// v18: CKI => sink 
// v19: CKI CycA => TriA 
// v20: TriA => CKI CycA 
// v21: TriA => CycA 
// v22: TriA => CKI 
// v23: CKI CycE => TriE 
// v24: TriE => CKI CycE 
// v25: TriE => CycE 
// v26: TriE => CKI 
// v27: source => CycA 
// v28: CycA => sink 
// v29: source => CycE 
// v30: CycE => sink 
// v31: source => APCP 
// v32: APCP => sink 
// v33: source => Cdc20in 
// v34: Cdc20in => sink 
// v35: Cdc20in => Cdc20A 
// v36: Cdc20A => Cdc20in 
// v37: Cdc20A => sink 
// v38: source => Cdh1 
// v39: Cdh1 => sink 

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

  // variable definitions
  real compartment = 1 L;
  real CycD(time);
  real Cdc14(time);
  real Wee1(time);
  real Vwee(time);
  real Cdc25(time);
  real V25(time);
  real TFB(time);
  real Vatf(time);
  real Vitf(time);
  real TFE(time);
  real TFI(time);
  real Vsb(time);
  real Vsa(time);
  real Vse(time);
  real Vah1(time);
  real Vih1(time);
  real Vdb(time);
  real Vda(time);
  real Vde(time);
  real Vsi(time);
  real Vdi(time);
  real CycD0 = .108;
  real ksep = 0;
  real ksepp = .15;
  real kdep = .12;
  real kdeapp = 0;
  real kdebpp = 0;
  real kdeepp = 0;
  real kasse = 0;
  real kdisse = 0;
  real katfp = 0;
  real katfapp = 1.5;
  real katfdpp = 3;
  real katfepp = .38;
  real kitfp = .75;
  real kitfapp = 0;
  real kitfbpp = 8;
  real Jatf = .01;
  real Jitf = .01;
  real ksbp = .004;
  real ksbpp = .04;
  real kdbp = .003;
  real kdbhpp = .4;
  real kdbcpp = .15;
  real kassb = 60;
  real kdissb = .05;
  real kweep = .02;
  real kweepp = .2;
  real k25p = .01;
  real k25pp = 5;
  real kafb = 1;
  real kifb = .15;
  real Jafb = .1;
  real Jifb = .1;
  real ksap = .0015;
  real ksapp = .01;
  real kdap = .01;
  real kdapp = .16;
  real kdappp = 0;
  real kassa = 50;
  real kdissa = .06;
  real ksip = .018;
  real ksipp = .18;
  real kdip = .002;
  real kdiapp = .1;
  real kdibpp = .8;
  real kdidpp = .1;
  real kdiepp = .12;
  real k14di = 12;
  real kafi = 6;
  real kifip = .008;
  real kifibpp = .05;
  real Jafi = 1;
  real Jifi = 2;
  real kaie = .1;
  real kiie = .15;
  real Jaie = .1;
  real Jiie = .1;
  real ks20p = .001;
  real ks20pp = 1;
  real n20 = 1;
  real J20 = 10;
  real kd20 = .05;
  real ka20 = 1;
  real ki20 = .16;
  real Ja20 = 1;
  real Ji20 = 1;
  real kah1p = .02;
  real kah1pp = .8;
  real kih1p = .001;
  real kih1app = .35;
  real kih1bpp = .1;
  real kih1dpp = .005;
  real kih1epp = .06;
  real Jah1 = .03;
  real Jih1 = .03;
  real kaweep = .3;
  real kaweepp = 0;
  real kiweep = 0;
  real kiweepp = 1;
  real Jawee = .05;
  real Jiwee = .05;
  real ka25p = 0;
  real ka25pp = 1;
  real ki25p = .3;
  real ki25pp = 0;
  real Ja25 = .1;
  real Ji25 = .1;
  real KEZ = .2;
  real mu = .005776;
  real MaxMass = 1E4;
  real source_ = 1 M;
  real Mass(time) M;
  real CycB(time) M;
  real sink = 1 M;
  real pB(time) M;
  real CKI(time) M;
  real BCKI(time) M;
  real pBCKI(time) M;
  real CycA(time) M;
  real TriA(time) M;
  real CycE(time) M;
  real TriE(time) M;
  real APCP(time) M;
  real Cdc20in(time) M;
  real Cdc20A(time) M;
  real Cdh1(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;
  private real Michaelis_call0(time);
  real v32(time) katal;
  private real Michaelis_call1(time);
  real v33(time) katal;
  real v34(time) katal;
  real v35(time) katal;
  private real Michaelis_call2(time);
  real v36(time) katal;
  private real Michaelis_call3(time);
  real v37(time) katal;
  real v38(time) katal;
  private real Michaelis_call4(time);
  real v39(time) katal;
  private real Michaelis_call5(time);
  private real GK_call0(time);
  private real GK_call1(time);
  private real GK_call2(time);
  private real GK_call3(time);
  private real GK_call4(time);

  // equations
  CycD = CycD0*Mass;
  Cdc14 = Cdc20A;
  Wee1 = GK_call0;
  Vwee = kweep+kweepp*Wee1;
  Cdc25 = GK_call1;
  V25 = k25p+k25pp*Cdc25;
  TFB = GK_call2;
  Vatf = katfp+katfapp*CycA+katfepp*CycE+katfdpp*CycD;
  Vitf = kitfp+kitfapp*CycA+kitfbpp*CycB;
  TFE = GK_call3;
  TFI = GK_call4;
  Vsb = (ksbp+ksbpp*TFB)*Mass;
  Vsa = (ksap+ksapp*TFE)*Mass;
  Vse = (ksep+ksepp*TFE)*Mass;
  Vah1 = kah1p+kah1pp*Cdc14;
  Vih1 = kih1p+kih1app*CycA+kih1bpp*CycB+kih1epp*CycE+kih1dpp*CycD;
  Vdb = kdbp+kdbhpp*Cdh1+kdbcpp*Cdc20A;
  Vda = kdap+(kdapp+kdappp)*Cdc20A+kdappp*Cdc20in;
  Vde = kdep+kdeepp*CycE+kdeapp*CycA+kdebpp*CycB;
  Vsi = ksip+ksipp*TFI;
  Vdi = (kdip+kdiapp*CycA+kdibpp*CycB+kdiepp*CycE+kdidpp*CycD)/(1+k14di*Cdc14);
  when (time=time.min) Mass = 1.3382678;
  Mass:time = (v1)/compartment;
  when (time=time.min) CycB = .19035859;
  CycB:time = (v2 + -1*v3 + v4 + -1*v5 + -1*v6 + v7 + v14)/compartment;
  when (time=time.min) pB = .01234426;
  pB:time = (-1*v4 + v5 + -1*v8 + -1*v9 + v10 + v16)/compartment;
  when (time=time.min) CKI = .0288207;
  CKI:time = (-1*v6 + v7 + -1*v9 + v10 + v13 + v15 + v17 + -1*v18 + -1*v19 + v20 + v22 + -1*v23 + v24 + v26)/compartment;
  when (time=time.min) BCKI = .6794492;
  BCKI:time = (v6 + -1*v7 + v11 + -1*v12 + -1*v13 + -1*v14)/compartment;
  when (time=time.min) pBCKI = .04795938;
  pBCKI:time = (v9 + -1*v10 + -1*v11 + v12 + -1*v15 + -1*v16)/compartment;
  when (time=time.min) CycA = .09450308;
  CycA:time = (-1*v19 + v20 + v21 + v27 + -1*v28)/compartment;
  when (time=time.min) TriA = .34922242;
  TriA:time = (v19 + -1*v20 + -1*v21 + -1*v22)/compartment;
  when (time=time.min) CycE = .20920244;
  CycE:time = (-1*v23 + v24 + v25 + v29 + -1*v30)/compartment;
  when (time=time.min) TriE = 0;
  TriE:time = (v23 + -1*v24 + -1*v25 + -1*v26)/compartment;
  when (time=time.min) APCP = .07591403;
  APCP:time = (v31 + -1*v32)/compartment;
  when (time=time.min) Cdc20in = .77023876;
  Cdc20in:time = (v33 + -1*v34 + -1*v35 + v36)/compartment;
  when (time=time.min) Cdc20A = .35727239;
  Cdc20A:time = (v35 + -1*v36 + -1*v37)/compartment;
  when (time=time.min) Cdh1 = .71893936;
  Cdh1:time = (v38 + -1*v39)/compartment;
  v1 = mu*Mass*(1-Mass/MaxMass);
  v2 = Vsb;
  v3 = Vdb*CycB;
  v4 = V25*pB;
  v5 = Vwee*CycB;
  v6 = kassb*CycB*CKI;
  v7 = kdissb*BCKI;
  v8 = Vdb*pB;
  v9 = kassb*pB*CKI;
  v10 = kdissb*pBCKI;
  v11 = V25*pBCKI;
  v12 = Vwee*BCKI;
  v13 = Vdb*BCKI;
  v14 = Vdi*BCKI;
  v15 = Vdb*pBCKI;
  v16 = Vdi*pBCKI;
  v17 = Vsi;
  v18 = Vdi*CKI;
  v19 = kassa*CKI*CycA;
  v20 = kdissa*TriA;
  v21 = Vdi*TriA;
  v22 = Vda*TriA;
  v23 = kasse*CKI*CycE;
  v24 = kdisse*TriE;
  v25 = Vdi*TriE;
  v26 = Vde*TriE;
  v27 = Vsa;
  v28 = Vda*CycA;
  v29 = Vse;
  v30 = Vde*CycE;
  v31 = Michaelis_call0;
  Michaelis_call0 = kaie*(1-APCP)*CycB/(Jaie+(1-APCP));
  v32 = Michaelis_call1;
  Michaelis_call1 = kiie*APCP*1/(Jiie+APCP);
  v33 = ks20p+ks20pp*CycB^n20/(J20^n20+CycB^n20);
  v34 = kd20*Cdc20in;
  v35 = Michaelis_call2;
  Michaelis_call2 = ka20*Cdc20in*APCP/(Ja20+Cdc20in);
  v36 = Michaelis_call3;
  Michaelis_call3 = ki20*Cdc20A*1/(Ji20+Cdc20A);
  v37 = kd20*Cdc20A;
  v38 = Michaelis_call4;
  Michaelis_call4 = 1*(1-Cdh1)*Vah1/(Jah1+(1-Cdh1));
  v39 = Michaelis_call5;
  Michaelis_call5 = 1*Cdh1*Vih1/(Jih1+Cdh1);
  GK_call0 = 2*(kaweep+kaweepp*Cdc14)*Jiwee/(kiweep+kiweepp*CycB-(kaweep+kaweepp*Cdc14)+(kiweep+kiweepp*CycB)*Jawee+(kaweep+kaweepp*Cdc14)*Jiwee+((kiweep+kiweepp*CycB-(kaweep+kaweepp*Cdc14)+(kiweep+kiweepp*CycB)*Jawee+(kaweep+kaweepp*Cdc14)*Jiwee)^2-4*(kaweep+kaweepp*Cdc14)*Jiwee*(kiweep+kiweepp*CycB-(kaweep+kaweepp*Cdc14)))^(1/2));
  GK_call1 = 2*(ka25p+ka25pp*CycB)*Ji25/(ki25p+ki25pp*Cdc14-(ka25p+ka25pp*CycB)+(ki25p+ki25pp*Cdc14)*Ja25+(ka25p+ka25pp*CycB)*Ji25+((ki25p+ki25pp*Cdc14-(ka25p+ka25pp*CycB)+(ki25p+ki25pp*Cdc14)*Ja25+(ka25p+ka25pp*CycB)*Ji25)^2-4*(ka25p+ka25pp*CycB)*Ji25*(ki25p+ki25pp*Cdc14-(ka25p+ka25pp*CycB)))^(1/2));
  GK_call2 = 2*(kafb*CycB)*Jifb/(kifb-kafb*CycB+kifb*Jafb+kafb*CycB*Jifb+((kifb-kafb*CycB+kifb*Jafb+kafb*CycB*Jifb)^2-4*(kafb*CycB)*Jifb*(kifb-kafb*CycB))^(1/2));
  GK_call3 = 2*Vatf*Jitf/(Vitf-Vatf+Vitf*Jatf+Vatf*Jitf+((Vitf-Vatf+Vitf*Jatf+Vatf*Jitf)^2-4*Vatf*Jitf*(Vitf-Vatf))^(1/2));
  GK_call4 = 2*(kafi*Cdc14)*Jifi/(kifip+kifibpp*CycB-kafi*Cdc14+(kifip+kifibpp*CycB)*Jafi+kafi*Cdc14*Jifi+((kifip+kifibpp*CycB-kafi*Cdc14+(kifip+kifibpp*CycB)*Jafi+kafi*Cdc14*Jifi)^2-4*(kafi*Cdc14)*Jifi*(kifip+kifibpp*CycB-kafi*Cdc14))^(1/2));

  // Used function definitions
  // Function definition Michaelis(M1,J1,k1,S1)=k1*S1*M1/(J1+S1);
  // Function definition GK(a,b,c,d)=2*a*d/(b-a+b*c+a*d+((b-a+b*c+a*d)^2-4*a*d*(b-a))^(1/2));

  // variable properties
  compartment.sbmlRole="compartment";
  CycD.sbmlRole="parameter";
  Cdc14.sbmlRole="parameter";
  Wee1.sbmlRole="parameter";
  Vwee.sbmlRole="parameter";
  Cdc25.sbmlRole="parameter";
  V25.sbmlRole="parameter";
  TFB.sbmlRole="parameter";
  Vatf.sbmlRole="parameter";
  Vitf.sbmlRole="parameter";
  TFE.sbmlRole="parameter";
  TFI.sbmlRole="parameter";
  Vsb.sbmlRole="parameter";
  Vsa.sbmlRole="parameter";
  Vse.sbmlRole="parameter";
  Vah1.sbmlRole="parameter";
  Vih1.sbmlRole="parameter";
  Vdb.sbmlRole="parameter";
  Vda.sbmlRole="parameter";
  Vde.sbmlRole="parameter";
  Vsi.sbmlRole="parameter";
  Vdi.sbmlRole="parameter";
  CycD0.sbmlRole="parameter";
  ksep.sbmlRole="parameter";
  ksepp.sbmlRole="parameter";
  kdep.sbmlRole="parameter";
  kdeapp.sbmlRole="parameter";
  kdebpp.sbmlRole="parameter";
  kdeepp.sbmlRole="parameter";
  kasse.sbmlRole="parameter";
  kdisse.sbmlRole="parameter";
  katfp.sbmlRole="parameter";
  katfapp.sbmlRole="parameter";
  katfdpp.sbmlRole="parameter";
  katfepp.sbmlRole="parameter";
  kitfp.sbmlRole="parameter";
  kitfapp.sbmlRole="parameter";
  kitfbpp.sbmlRole="parameter";
  Jatf.sbmlRole="parameter";
  Jitf.sbmlRole="parameter";
  ksbp.sbmlRole="parameter";
  ksbpp.sbmlRole="parameter";
  kdbp.sbmlRole="parameter";
  kdbhpp.sbmlRole="parameter";
  kdbcpp.sbmlRole="parameter";
  kassb.sbmlRole="parameter";
  kdissb.sbmlRole="parameter";
  kweep.sbmlRole="parameter";
  kweepp.sbmlRole="parameter";
  k25p.sbmlRole="parameter";
  k25pp.sbmlRole="parameter";
  kafb.sbmlRole="parameter";
  kifb.sbmlRole="parameter";
  Jafb.sbmlRole="parameter";
  Jifb.sbmlRole="parameter";
  ksap.sbmlRole="parameter";
  ksapp.sbmlRole="parameter";
  kdap.sbmlRole="parameter";
  kdapp.sbmlRole="parameter";
  kdappp.sbmlRole="parameter";
  kassa.sbmlRole="parameter";
  kdissa.sbmlRole="parameter";
  ksip.sbmlRole="parameter";
  ksipp.sbmlRole="parameter";
  kdip.sbmlRole="parameter";
  kdiapp.sbmlRole="parameter";
  kdibpp.sbmlRole="parameter";
  kdidpp.sbmlRole="parameter";
  kdiepp.sbmlRole="parameter";
  k14di.sbmlRole="parameter";
  kafi.sbmlRole="parameter";
  kifip.sbmlRole="parameter";
  kifibpp.sbmlRole="parameter";
  Jafi.sbmlRole="parameter";
  Jifi.sbmlRole="parameter";
  kaie.sbmlRole="parameter";
  kiie.sbmlRole="parameter";
  Jaie.sbmlRole="parameter";
  Jiie.sbmlRole="parameter";
  ks20p.sbmlRole="parameter";
  ks20pp.sbmlRole="parameter";
  n20.sbmlRole="parameter";
  J20.sbmlRole="parameter";
  kd20.sbmlRole="parameter";
  ka20.sbmlRole="parameter";
  ki20.sbmlRole="parameter";
  Ja20.sbmlRole="parameter";
  Ji20.sbmlRole="parameter";
  kah1p.sbmlRole="parameter";
  kah1pp.sbmlRole="parameter";
  kih1p.sbmlRole="parameter";
  kih1app.sbmlRole="parameter";
  kih1bpp.sbmlRole="parameter";
  kih1dpp.sbmlRole="parameter";
  kih1epp.sbmlRole="parameter";
  Jah1.sbmlRole="parameter";
  Jih1.sbmlRole="parameter";
  kaweep.sbmlRole="parameter";
  kaweepp.sbmlRole="parameter";
  kiweep.sbmlRole="parameter";
  kiweepp.sbmlRole="parameter";
  Jawee.sbmlRole="parameter";
  Jiwee.sbmlRole="parameter";
  ka25p.sbmlRole="parameter";
  ka25pp.sbmlRole="parameter";
  ki25p.sbmlRole="parameter";
  ki25pp.sbmlRole="parameter";
  Ja25.sbmlRole="parameter";
  Ji25.sbmlRole="parameter";
  KEZ.sbmlRole="parameter";
  mu.sbmlRole="parameter";
  MaxMass.sbmlRole="parameter";
  source_.sbmlRole="species";
  source_.sbmlCompartment="compartment";
  Mass.sbmlRole="species";
  Mass.sbmlCompartment="compartment";
  CycB.sbmlRole="species";
  CycB.sbmlCompartment="compartment";
  sink.sbmlRole="species";
  sink.sbmlCompartment="compartment";
  pB.sbmlRole="species";
  pB.sbmlCompartment="compartment";
  CKI.sbmlRole="species";
  CKI.sbmlCompartment="compartment";
  BCKI.sbmlRole="species";
  BCKI.sbmlCompartment="compartment";
  pBCKI.sbmlRole="species";
  pBCKI.sbmlCompartment="compartment";
  CycA.sbmlRole="species";
  CycA.sbmlCompartment="compartment";
  TriA.sbmlRole="species";
  TriA.sbmlCompartment="compartment";
  CycE.sbmlRole="species";
  CycE.sbmlCompartment="compartment";
  TriE.sbmlRole="species";
  TriE.sbmlCompartment="compartment";
  APCP.sbmlRole="species";
  APCP.sbmlCompartment="compartment";
  Cdc20in.sbmlRole="species";
  Cdc20in.sbmlCompartment="compartment";
  Cdc20A.sbmlRole="species";
  Cdc20A.sbmlCompartment="compartment";
  Cdh1.sbmlRole="species";
  Cdh1.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";
}