// 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;
unit item=scalar;
unit substance = mole;
unit volume = litre;
unit area = metre^2;
unit length = metre;
unit time = second;

// 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 volume;
  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;
  private real s.amt substance;
  real s substance/volume;
  real s.init substance/volume;
  private real Cln2.amt(time) substance;
  real Cln2(time) substance/volume;
  real Cln2.init substance/volume;
  private real p.amt substance;
  real p substance/volume;
  real p.init substance/volume;
  private real Clb2T.amt(time) substance;
  real Clb2T(time) substance/volume;
  real Clb2T.init substance/volume;
  private real Clb5T.amt(time) substance;
  real Clb5T(time) substance/volume;
  real Clb5T.init substance/volume;
  private real Sic1T.amt(time) substance;
  real Sic1T(time) substance/volume;
  real Sic1T.init substance/volume;
  private real Clb2Sic1.amt(time) substance;
  real Clb2Sic1(time) substance/volume;
  real Clb2Sic1.init substance/volume;
  private real Clb5Sic1.amt(time) substance;
  real Clb5Sic1(time) substance/volume;
  real Clb5Sic1.init substance/volume;
  private real Cdc20T.amt(time) substance;
  real Cdc20T(time) substance/volume;
  real Cdc20T.init substance/volume;
  private real Cdc20.amt(time) substance;
  real Cdc20(time) substance/volume;
  real Cdc20.init substance/volume;
  private real Hct1.amt(time) substance;
  real Hct1(time) substance/volume;
  real Hct1.init substance/volume;
  private real mass.amt(time) substance;
  real mass(time) substance/volume;
  real mass.init substance/volume;
  private real ORI.amt(time) substance;
  real ORI(time) substance/volume;
  real ORI.init substance/volume;
  private real BUD.amt(time) substance;
  real BUD(time) substance/volume;
  real BUD.init substance/volume;
  private real SPN.amt(time) substance;
  real SPN(time) substance/volume;
  real SPN.init substance/volume;
  real v1.rate(time) substance/time;
  real v2.rate(time) substance/time;
  real v3.rate(time) substance/time;
  real v4.rate(time) substance/time;
  real v5.rate(time) substance/time;
  real v6.rate(time) substance/time;
  real v7.rate(time) substance/time;
  real v8.rate(time) substance/time;
  real v9.rate(time) substance/time;
  real v10.rate(time) substance/time;
  real v11.rate(time) substance/time;
  real v12.rate(time) substance/time;
  real v13.rate(time) substance/time;
  real v14.rate(time) substance/time;
  real v15.rate(time) substance/time;
  extern real v16.rate substance/time;
  real v17.rate(time) substance/time;
  real v18.rate(time) substance/time;
  real v19.rate(time) substance/time;
  real v20.rate(time) substance/time;
  real v21.rate(time) substance/time;
  real v22.rate(time) substance/time;
  real v23.rate(time) substance/time;
  real v24.rate(time) substance/time;
  real v25.rate(time) substance/time;
  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;
  s.amt = s*compartment;
  s = s.init;
  s.init = 1;
  when (time=time.min) Cln2.amt = Cln2.init*compartment;
  Cln2.amt:time = v1.rate + -1*v2.rate;
  Cln2 = Cln2.amt/compartment;
  Cln2.init = .01322672;
  p.amt = p*compartment;
  p = p.init;
  p.init = 1;
  when (time=time.min) Clb2T.amt = Clb2T.init*compartment;
  Clb2T.amt:time = v3.rate + -1*v4.rate;
  Clb2T = Clb2T.amt/compartment;
  Clb2T.init = 1.3987145;
  when (time=time.min) Clb5T.amt = Clb5T.init*compartment;
  Clb5T.amt:time = v5.rate + -1*v6.rate;
  Clb5T = Clb5T.amt/compartment;
  Clb5T.init = .07965319;
  when (time=time.min) Sic1T.amt = Sic1T.init*compartment;
  Sic1T.amt:time = v7.rate + -1*v8.rate;
  Sic1T = Sic1T.amt/compartment;
  Sic1T.init = .01595117;
  when (time=time.min) Clb2Sic1.amt = Clb2Sic1.init*compartment;
  Clb2Sic1.amt:time = v9.rate + -1*v10.rate;
  Clb2Sic1 = Clb2Sic1.amt/compartment;
  Clb2Sic1.init = .01485261;
  when (time=time.min) Clb5Sic1.amt = Clb5Sic1.init*compartment;
  Clb5Sic1.amt:time = v11.rate + -1*v12.rate;
  Clb5Sic1 = Clb5Sic1.amt/compartment;
  Clb5Sic1.init = 8.0991225E-4;
  when (time=time.min) Cdc20T.amt = Cdc20T.init*compartment;
  Cdc20T.amt:time = v13.rate + -1*v14.rate;
  Cdc20T = Cdc20T.amt/compartment;
  Cdc20T.init = .83375003;
  when (time=time.min) Cdc20.amt = Cdc20.init*compartment;
  Cdc20.amt:time = v15.rate + -1*v16.rate;
  Cdc20 = Cdc20.amt/compartment;
  Cdc20.init = .07507414;
  when (time=time.min) Hct1.amt = Hct1.init*compartment;
  Hct1.amt:time = v17.rate + -1*v18.rate;
  Hct1 = Hct1.amt/compartment;
  Hct1.init = .01220496;
  when (time=time.min) mass.amt = mass.init*compartment;
  mass.amt:time = v19.rate;
  mass = mass.amt/compartment;
  mass.init = 1.5222296;
  when (time=time.min) ORI.amt = ORI.init*compartment;
  ORI.amt:time = v20.rate + -1*v21.rate;
  ORI = ORI.amt/compartment;
  ORI.init = 15.6743543;
  when (time=time.min) BUD.amt = BUD.init*compartment;
  BUD.amt:time = v22.rate + -1*v23.rate;
  BUD = BUD.amt/compartment;
  BUD.init = 1.0339811;
  when (time=time.min) SPN.amt = SPN.init*compartment;
  SPN.amt:time = v24.rate + -1*v25.rate;
  SPN = SPN.amt/compartment;
  SPN.init = .98317565;
  v1.rate = (ksn2a+ksn2aa*SBF)*mass;
  v2.rate = kdn2*Cln2;
  v3.rate = (ksb2a+ksb2aa*Mcm1)*mass;
  v4.rate = Vdb2*Clb2T;
  v5.rate = (ksb5a+ksb5aa*MBF)*mass;
  v6.rate = Vdb5*Clb5T;
  v7.rate = ksc1a+ksc1aa*Swi5;
  v8.rate = (kd1c1+Vd2c1/(Jd2c1+Sic1T))*Sic1T;
  v9.rate = kasb2*Clb2*Sic1;
  v10.rate = (kdib2+Vdb2+kd1c1+Vd2c1/(Jd2c1+Sic1T))*Clb2Sic1;
  v11.rate = kasb5*Clb5*Sic1;
  v12.rate = (kdib5+Vdb5+kd1c1+Vd2c1/(Jd2c1+Sic1T))*Clb5Sic1;
  v13.rate = ks20a+ks20aa*Clb2;
  v14.rate = kd20*Cdc20T;
  v15.rate = ka20*(Cdc20T-Cdc20);
  v17.rate = (kat1a+kat1aa*Cdc20)*(Hct1T-Hct1)/(Jat1+Hct1T-Hct1);
  v18.rate = Vit1*Hct1/(Jit1+Hct1);
  v19.rate = u*mass;
  v20.rate = ksori*(Clb5+eorib2*Clb2);
  v21.rate = kdori*ORI;
  v22.rate = ksbud*(Cln2+Cln3max*Dn3*mass/(Jn3+Dn3*mass)+ebudb5*Clb5);
  v23.rate = kdbud*BUD;
  v24.rate = ksspn*Clb2/(J_spn+Clb2);
  v25.rate = 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));

  // 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.amt.sbmlRole="speciesAmount";
  s.sbmlRole="speciesConcentration";
  s.sbmlCompartment="compartment";
  s.init.sbmlRole="speciesInitialConcentration";
  Cln2.amt.sbmlRole="speciesAmount";
  Cln2.sbmlRole="speciesConcentration";
  Cln2.sbmlCompartment="compartment";
  Cln2.init.sbmlRole="speciesInitialConcentration";
  p.amt.sbmlRole="speciesAmount";
  p.sbmlRole="speciesConcentration";
  p.sbmlCompartment="compartment";
  p.init.sbmlRole="speciesInitialConcentration";
  Clb2T.amt.sbmlRole="speciesAmount";
  Clb2T.sbmlRole="speciesConcentration";
  Clb2T.sbmlCompartment="compartment";
  Clb2T.init.sbmlRole="speciesInitialConcentration";
  Clb5T.amt.sbmlRole="speciesAmount";
  Clb5T.sbmlRole="speciesConcentration";
  Clb5T.sbmlCompartment="compartment";
  Clb5T.init.sbmlRole="speciesInitialConcentration";
  Sic1T.amt.sbmlRole="speciesAmount";
  Sic1T.sbmlRole="speciesConcentration";
  Sic1T.sbmlCompartment="compartment";
  Sic1T.init.sbmlRole="speciesInitialConcentration";
  Clb2Sic1.amt.sbmlRole="speciesAmount";
  Clb2Sic1.sbmlRole="speciesConcentration";
  Clb2Sic1.sbmlCompartment="compartment";
  Clb2Sic1.init.sbmlRole="speciesInitialConcentration";
  Clb5Sic1.amt.sbmlRole="speciesAmount";
  Clb5Sic1.sbmlRole="speciesConcentration";
  Clb5Sic1.sbmlCompartment="compartment";
  Clb5Sic1.init.sbmlRole="speciesInitialConcentration";
  Cdc20T.amt.sbmlRole="speciesAmount";
  Cdc20T.sbmlRole="speciesConcentration";
  Cdc20T.sbmlCompartment="compartment";
  Cdc20T.init.sbmlRole="speciesInitialConcentration";
  Cdc20.amt.sbmlRole="speciesAmount";
  Cdc20.sbmlRole="speciesConcentration";
  Cdc20.sbmlCompartment="compartment";
  Cdc20.init.sbmlRole="speciesInitialConcentration";
  Hct1.amt.sbmlRole="speciesAmount";
  Hct1.sbmlRole="speciesConcentration";
  Hct1.sbmlCompartment="compartment";
  Hct1.init.sbmlRole="speciesInitialConcentration";
  mass.amt.sbmlRole="speciesAmount";
  mass.sbmlRole="speciesConcentration";
  mass.sbmlCompartment="compartment";
  mass.init.sbmlRole="speciesInitialConcentration";
  ORI.amt.sbmlRole="speciesAmount";
  ORI.sbmlRole="speciesConcentration";
  ORI.sbmlCompartment="compartment";
  ORI.init.sbmlRole="speciesInitialConcentration";
  BUD.amt.sbmlRole="speciesAmount";
  BUD.sbmlRole="speciesConcentration";
  BUD.sbmlCompartment="compartment";
  BUD.init.sbmlRole="speciesInitialConcentration";
  SPN.amt.sbmlRole="speciesAmount";
  SPN.sbmlRole="speciesConcentration";
  SPN.sbmlCompartment="compartment";
  SPN.init.sbmlRole="speciesInitialConcentration";
  v1.rate.sbmlRole="rate";
  v2.rate.sbmlRole="rate";
  v3.rate.sbmlRole="rate";
  v4.rate.sbmlRole="rate";
  v5.rate.sbmlRole="rate";
  v6.rate.sbmlRole="rate";
  v7.rate.sbmlRole="rate";
  v8.rate.sbmlRole="rate";
  v9.rate.sbmlRole="rate";
  v10.rate.sbmlRole="rate";
  v11.rate.sbmlRole="rate";
  v12.rate.sbmlRole="rate";
  v13.rate.sbmlRole="rate";
  v14.rate.sbmlRole="rate";
  v15.rate.sbmlRole="rate";
  v16.rate.sbmlRole="rate";
  v17.rate.sbmlRole="rate";
  v18.rate.sbmlRole="rate";
  v19.rate.sbmlRole="rate";
  v20.rate.sbmlRole="rate";
  v21.rate.sbmlRole="rate";
  v22.rate.sbmlRole="rate";
  v23.rate.sbmlRole="rate";
  v24.rate.sbmlRole="rate";
  v25.rate.sbmlRole="rate";
  G.call0.sbmlRole="functionCall";
  G.call1.sbmlRole="functionCall";
  G.call2.sbmlRole="functionCall";
  G.call3.sbmlRole="functionCall";
}