// This model generated automatically from SBML
// WARNING: Renamed variable Jspn to J_spn 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: s => Cln2 
// v2: Cln2 => p 
// v3: s => Clb2T 
// v4: Clb2T => p 
// v5: s => Clb5T 
// v6: Clb5T => p 
// v7: s => Sic1T 
// v8: Sic1T => p 
// v9: s => Clb2Sic1 
// v10: Clb2Sic1 => p 
// v11: s => Clb5Sic1 
// v12: Clb5Sic1 => p 
// v13: s => Cdc20T 
// v14: Cdc20T => p 
// v15: s => Cdc20 
// v16: Cdc20 => p 
// v17: s => Hct1 
// v18: Hct1 => p 
// v19: s => mass 
// v20: s => ORI 
// v21: ORI => p 
// v22: s => BUD 
// v23: BUD => p 
// v24: s => SPN 
// v25: SPN => p 

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

  // variable definitions
  real compartment = 1 L;
  real Vdb2(time);
  real Vdb5(time);
  real Bck2(time);
  real Clb2(time);
  real Clb5(time);
  real Vd2c1(time);
  real Vasbf(time);
  real Vit1(time);
  real Sic1(time);
  real SBF(time);
  real MBF(time);
  real Mcm1(time);
  real Swi5(time);
  real ksn2a = 0;
  real ksb2a = .002;
  real kdb2a = .01;
  real ksb5a = .006;
  real ksc1a = .02;
  real kasb2 = 50;
  real kasb5 = 50;
  real ks20a = .005;
  real ka20 = 1;
  real kat1a = .04;
  real ksori = 2;
  real kdori = .06;
  real kdbud = .06;
  real kdspn = .06;
  real kisbfa = .5;
  real kimcm = .15;
  real ksn2aa = .05;
  real ksb2aa = .05;
  real kdb2aa = 2;
  real ksb5aa = .02;
  real ksc1aa = .1;
  real kdib2 = .05;
  real kdib5 = .05;
  real ks20aa = .06;
  real ki20a = .1;
  real kat1aa = 2;
  real ksbud = .3;
  real kisbfaa = 6;
  real u = .005776;
  real kdn2 = .1;
  real kdb2aaa = .05;
  real kdb5a = .1;
  real kd1c1 = .01;
  real kd20 = .08;
  real ki20aa = 10;
  real kit1a = 0;
  real ksspn = .08;
  real kasbf = 1;
  real kamcm = 1;
  real kaswi = 1;
  real kiswia = .3;
  real kdb5aa = .25;
  real kd2c1 = .3;
  real kit1aa = .64;
  real kiswiaa = 20;
  real Cln3max = .02;
  real J_spn = .2;
  real Jasbf = .01;
  real Jisbf = .01;
  real Bck20 = .0027;
  real Jd2c1 = .05;
  real Jamcm = 1;
  real Jimcm = 1;
  real Hct1T = 1;
  real Jaswi = .1;
  real Jiswi = .1;
  real Jat1 = .05;
  real Jit1 = .05;
  real ec1n3 = 20;
  real eit1n2 = 1;
  real eorib2 = .4;
  real ec1k2 = 2;
  real eit1b2 = 1;
  real ebudb5 = 1;
  real ec1b2 = .067;
  real eit1b5 = .5;
  real esbfn3 = 75;
  real ec1b5 = 1;
  real esbfb5 = .5;
  real f = .433;
  real Jn3 = 6;
  real Dn3 = 1;
  real Vi20 = 10;
  real s = 1 M;
  real Cln2(time) M;
  real p = 1 M;
  real Clb2T(time) M;
  real Clb5T(time) M;
  real Sic1T(time) M;
  real Clb2Sic1(time) M;
  real Clb5Sic1(time) M;
  real Cdc20T(time) M;
  real Cdc20(time) M;
  real Hct1(time) M;
  real mass(time) M;
  real ORI(time) M;
  real BUD(time) M;
  real SPN(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;
  extern real v16 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;
  private real G_call0(time);
  private real G_call1(time);
  private real G_call2(time);
  private real G_call3(time);

  // equations
  Vdb2 = kdb2a*(Hct1T-Hct1)+kdb2aa*Hct1+kdb2aaa*Cdc20;
  Vdb5 = kdb5a+kdb5aa*Cdc20;
  Bck2 = Bck20*mass;
  Clb2 = Clb2T-Clb2Sic1;
  Clb5 = Clb5T-Clb5Sic1;
  Vd2c1 = kd2c1*(ec1n3*(Cln3max*Dn3*mass/(Jn3+Dn3*mass))+ec1k2*Bck2+Cln2+ec1b5*Clb5+ec1b2*Clb2);
  Vasbf = kasbf*(Cln2+esbfn3*(Cln3max*Dn3*mass/(Jn3+Dn3*mass)+Bck2)+esbfb5*Clb5);
  Vit1 = kit1a+kit1aa*(Cln3max*Dn3*mass/(Jn3+Dn3*mass)+eit1n2*Cln2+eit1b5*Clb5+eit1b2*Clb2);
  Sic1 = Sic1T-(Clb2Sic1+Clb5Sic1);
  SBF = G_call0;
  MBF = G_call1;
  Mcm1 = G_call2;
  Swi5 = G_call3;
  when (time=time.min) Cln2 = .01322672;
  Cln2:time = (v1 + -1*v2)/compartment;
  when (time=time.min) Clb2T = 1.3987145;
  Clb2T:time = (v3 + -1*v4)/compartment;
  when (time=time.min) Clb5T = .07965319;
  Clb5T:time = (v5 + -1*v6)/compartment;
  when (time=time.min) Sic1T = .01595117;
  Sic1T:time = (v7 + -1*v8)/compartment;
  when (time=time.min) Clb2Sic1 = .01485261;
  Clb2Sic1:time = (v9 + -1*v10)/compartment;
  when (time=time.min) Clb5Sic1 = 8.0991225E-4;
  Clb5Sic1:time = (v11 + -1*v12)/compartment;
  when (time=time.min) Cdc20T = .83375003;
  Cdc20T:time = (v13 + -1*v14)/compartment;
  when (time=time.min) Cdc20 = .07507414;
  Cdc20:time = (v15 + -1*v16)/compartment;
  when (time=time.min) Hct1 = .01220496;
  Hct1:time = (v17 + -1*v18)/compartment;
  when (time=time.min) mass = 1.5222296;
  mass:time = (v19)/compartment;
  when (time=time.min) ORI = 15.6743543;
  ORI:time = (v20 + -1*v21)/compartment;
  when (time=time.min) BUD = 1.0339811;
  BUD:time = (v22 + -1*v23)/compartment;
  when (time=time.min) SPN = .98317565;
  SPN:time = (v24 + -1*v25)/compartment;
  v1 = (ksn2a+ksn2aa*SBF)*mass;
  v2 = kdn2*Cln2;
  v3 = (ksb2a+ksb2aa*Mcm1)*mass;
  v4 = Vdb2*Clb2T;
  v5 = (ksb5a+ksb5aa*MBF)*mass;
  v6 = Vdb5*Clb5T;
  v7 = ksc1a+ksc1aa*Swi5;
  v8 = (kd1c1+Vd2c1/(Jd2c1+Sic1T))*Sic1T;
  v9 = kasb2*Clb2*Sic1;
  v10 = (kdib2+Vdb2+kd1c1+Vd2c1/(Jd2c1+Sic1T))*Clb2Sic1;
  v11 = kasb5*Clb5*Sic1;
  v12 = (kdib5+Vdb5+kd1c1+Vd2c1/(Jd2c1+Sic1T))*Clb5Sic1;
  v13 = ks20a+ks20aa*Clb2;
  v14 = kd20*Cdc20T;
  v15 = ka20*(Cdc20T-Cdc20);
  v17 = (kat1a+kat1aa*Cdc20)*(Hct1T-Hct1)/(Jat1+Hct1T-Hct1);
  v18 = Vit1*Hct1/(Jit1+Hct1);
  v19 = u*mass;
  v20 = ksori*(Clb5+eorib2*Clb2);
  v21 = kdori*ORI;
  v22 = ksbud*(Cln2+Cln3max*Dn3*mass/(Jn3+Dn3*mass)+ebudb5*Clb5);
  v23 = kdbud*BUD;
  v24 = ksspn*Clb2/(J_spn+Clb2);
  v25 = kdspn*SPN;
  G_call0 = 2*Vasbf*Jisbf/(kisbfa+kisbfaa*Clb2-Vasbf+(kisbfa+kisbfaa*Clb2)*Jasbf+Vasbf*Jisbf+((kisbfa+kisbfaa*Clb2-Vasbf+(kisbfa+kisbfaa*Clb2)*Jasbf+Vasbf*Jisbf)^2-4*Vasbf*Jisbf*(kisbfa+kisbfaa*Clb2-Vasbf))^(1/2));
  G_call1 = 2*Vasbf*Jisbf/(kisbfa+kisbfaa*Clb2-Vasbf+(kisbfa+kisbfaa*Clb2)*Jasbf+Vasbf*Jisbf+((kisbfa+kisbfaa*Clb2-Vasbf+(kisbfa+kisbfaa*Clb2)*Jasbf+Vasbf*Jisbf)^2-4*Vasbf*Jisbf*(kisbfa+kisbfaa*Clb2-Vasbf))^(1/2));
  G_call2 = 2*(kamcm*Clb2)*Jimcm/(kimcm-kamcm*Clb2+kimcm*Jamcm+kamcm*Clb2*Jimcm+((kimcm-kamcm*Clb2+kimcm*Jamcm+kamcm*Clb2*Jimcm)^2-4*(kamcm*Clb2)*Jimcm*(kimcm-kamcm*Clb2))^(1/2));
  G_call3 = 2*(kaswi*Cdc20)*Jiswi/(kiswia+kiswiaa*Clb2-kaswi*Cdc20+(kiswia+kiswiaa*Clb2)*Jaswi+kaswi*Cdc20*Jiswi+((kiswia+kiswiaa*Clb2-kaswi*Cdc20+(kiswia+kiswiaa*Clb2)*Jaswi+kaswi*Cdc20*Jiswi)^2-4*(kaswi*Cdc20)*Jiswi*(kiswia+kiswiaa*Clb2-kaswi*Cdc20))^(1/2));

  // Used function definitions
  // Function definition G(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";
  Vdb2.sbmlRole="parameter";
  Vdb5.sbmlRole="parameter";
  Bck2.sbmlRole="parameter";
  Clb2.sbmlRole="parameter";
  Clb5.sbmlRole="parameter";
  Vd2c1.sbmlRole="parameter";
  Vasbf.sbmlRole="parameter";
  Vit1.sbmlRole="parameter";
  Sic1.sbmlRole="parameter";
  SBF.sbmlRole="parameter";
  MBF.sbmlRole="parameter";
  Mcm1.sbmlRole="parameter";
  Swi5.sbmlRole="parameter";
  ksn2a.sbmlRole="parameter";
  ksb2a.sbmlRole="parameter";
  kdb2a.sbmlRole="parameter";
  ksb5a.sbmlRole="parameter";
  ksc1a.sbmlRole="parameter";
  kasb2.sbmlRole="parameter";
  kasb5.sbmlRole="parameter";
  ks20a.sbmlRole="parameter";
  ka20.sbmlRole="parameter";
  kat1a.sbmlRole="parameter";
  ksori.sbmlRole="parameter";
  kdori.sbmlRole="parameter";
  kdbud.sbmlRole="parameter";
  kdspn.sbmlRole="parameter";
  kisbfa.sbmlRole="parameter";
  kimcm.sbmlRole="parameter";
  ksn2aa.sbmlRole="parameter";
  ksb2aa.sbmlRole="parameter";
  kdb2aa.sbmlRole="parameter";
  ksb5aa.sbmlRole="parameter";
  ksc1aa.sbmlRole="parameter";
  kdib2.sbmlRole="parameter";
  kdib5.sbmlRole="parameter";
  ks20aa.sbmlRole="parameter";
  ki20a.sbmlRole="parameter";
  kat1aa.sbmlRole="parameter";
  ksbud.sbmlRole="parameter";
  kisbfaa.sbmlRole="parameter";
  u.sbmlRole="parameter";
  kdn2.sbmlRole="parameter";
  kdb2aaa.sbmlRole="parameter";
  kdb5a.sbmlRole="parameter";
  kd1c1.sbmlRole="parameter";
  kd20.sbmlRole="parameter";
  ki20aa.sbmlRole="parameter";
  kit1a.sbmlRole="parameter";
  ksspn.sbmlRole="parameter";
  kasbf.sbmlRole="parameter";
  kamcm.sbmlRole="parameter";
  kaswi.sbmlRole="parameter";
  kiswia.sbmlRole="parameter";
  kdb5aa.sbmlRole="parameter";
  kd2c1.sbmlRole="parameter";
  kit1aa.sbmlRole="parameter";
  kiswiaa.sbmlRole="parameter";
  Cln3max.sbmlRole="parameter";
  J_spn.sbmlRole="parameter";
  Jasbf.sbmlRole="parameter";
  Jisbf.sbmlRole="parameter";
  Bck20.sbmlRole="parameter";
  Jd2c1.sbmlRole="parameter";
  Jamcm.sbmlRole="parameter";
  Jimcm.sbmlRole="parameter";
  Hct1T.sbmlRole="parameter";
  Jaswi.sbmlRole="parameter";
  Jiswi.sbmlRole="parameter";
  Jat1.sbmlRole="parameter";
  Jit1.sbmlRole="parameter";
  ec1n3.sbmlRole="parameter";
  eit1n2.sbmlRole="parameter";
  eorib2.sbmlRole="parameter";
  ec1k2.sbmlRole="parameter";
  eit1b2.sbmlRole="parameter";
  ebudb5.sbmlRole="parameter";
  ec1b2.sbmlRole="parameter";
  eit1b5.sbmlRole="parameter";
  esbfn3.sbmlRole="parameter";
  ec1b5.sbmlRole="parameter";
  esbfb5.sbmlRole="parameter";
  f.sbmlRole="parameter";
  Jn3.sbmlRole="parameter";
  Dn3.sbmlRole="parameter";
  Vi20.sbmlRole="parameter";
  s.sbmlRole="species";
  s.sbmlCompartment="compartment";
  Cln2.sbmlRole="species";
  Cln2.sbmlCompartment="compartment";
  p.sbmlRole="species";
  p.sbmlCompartment="compartment";
  Clb2T.sbmlRole="species";
  Clb2T.sbmlCompartment="compartment";
  Clb5T.sbmlRole="species";
  Clb5T.sbmlCompartment="compartment";
  Sic1T.sbmlRole="species";
  Sic1T.sbmlCompartment="compartment";
  Clb2Sic1.sbmlRole="species";
  Clb2Sic1.sbmlCompartment="compartment";
  Clb5Sic1.sbmlRole="species";
  Clb5Sic1.sbmlCompartment="compartment";
  Cdc20T.sbmlRole="species";
  Cdc20T.sbmlCompartment="compartment";
  Cdc20.sbmlRole="species";
  Cdc20.sbmlCompartment="compartment";
  Hct1.sbmlRole="species";
  Hct1.sbmlCompartment="compartment";
  mass.sbmlRole="species";
  mass.sbmlCompartment="compartment";
  ORI.sbmlRole="species";
  ORI.sbmlCompartment="compartment";
  BUD.sbmlRole="species";
  BUD.sbmlCompartment="compartment";
  SPN.sbmlRole="species";
  SPN.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";
}