// This model generated automatically from SBML

// 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 => MPF MPFtotal 
// v2: MPF => p 
// v3: s => MPF 
// v4: MPF => p 
// v5: MPFtotal => p 
// v6: s => Cdc20 
// v7: Cdc20 => p 
// v8: s => Cdc20 
// v9: Cdc20 APCP => APCPCdc20 
// v10: s => Mad2Cdc20P 
// v11: Mad2Cdc20P => p 
// v12: APCPCdc20 => p 
// v13: s => APCP 
// v14: APCP => p 
// v15: s => APCP 
// v16: s => Wee 
// v17: Wee => p 
// v18: s => Cdc25P 
// v19: Cdc25P => p 

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

  // variable definitions
  real compartment = 1 volume;
  real kwee(time);
  real k25(time);
  real APC(time);
  real Cdc20P(time);
  real Mad2(time);
  real kcycsyn = .4;
  real kcycdeg = 2;
  real kacycdeg = .08;
  real kapcph = .04;
  real kapcdeph = .04;
  real k20deph = 8;
  real k20ph = 10;
  real kass = 80;
  real kdiss = 20;
  real kass2 = 10;
  real kdiss2 = .1;
  real ka25 = .2;
  real kaa25 = 10;
  real Vi25 = .2;
  real Vapa25 = 1;
  real kawee = .1;
  real kaawee = 4;
  real Vawee = .2;
  real Vaaiwee = 1;
  real Mad2active = .01;
  real APCtotal = 1;
  real Cdc20total = 1;
  real Cdc25total = 1;
  real Weetotal = 1;
  real Jwee = .01;
  real Jcycdeg = .01;
  real J25 = .01;
  private real s.amt substance;
  real s substance/volume;
  real s.init substance/volume;
  private real MPF.amt(time) substance;
  real MPF(time) substance/volume;
  real MPF.init substance/volume;
  private real MPFtotal.amt(time) substance;
  real MPFtotal(time) substance/volume;
  real MPFtotal.init substance/volume;
  private real p.amt substance;
  real p substance/volume;
  real p.init substance/volume;
  private real Cdc20.amt(time) substance;
  real Cdc20(time) substance/volume;
  real Cdc20.init substance/volume;
  private real APCP.amt(time) substance;
  real APCP(time) substance/volume;
  real APCP.init substance/volume;
  private real APCPCdc20.amt(time) substance;
  real APCPCdc20(time) substance/volume;
  real APCPCdc20.init substance/volume;
  private real Mad2Cdc20P.amt(time) substance;
  real Mad2Cdc20P(time) substance/volume;
  real Mad2Cdc20P.init substance/volume;
  private real Wee.amt(time) substance;
  real Wee(time) substance/volume;
  real Wee.init substance/volume;
  private real Cdc25P.amt(time) substance;
  real Cdc25P(time) substance/volume;
  real Cdc25P.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;
  real v16.rate(time) substance/time;
  real v17.rate(time) substance/time;
  real v18.rate(time) substance/time;
  real v19.rate(time) substance/time;

  // equations
  kwee = Weetotal*kawee+(kaawee-kawee)*Wee;
  k25 = Cdc25total*ka25+(kaa25-ka25)*Cdc25P;
  APC = APCtotal-APCP-APCPCdc20;
  Cdc20P = Cdc20total-Cdc20-APCPCdc20-Mad2Cdc20P;
  Mad2 = Mad2active-Mad2Cdc20P;
  s.amt = s*compartment;
  s = s.init;
  s.init = 1;
  when (time=time.min) MPF.amt = MPF.init*compartment;
  MPF.amt:time = v1.rate + -1*v2.rate + v3.rate + -1*v4.rate;
  MPF = MPF.amt/compartment;
  MPF.init = .015;
  when (time=time.min) MPFtotal.amt = MPFtotal.init*compartment;
  MPFtotal.amt:time = v1.rate + -1*v5.rate;
  MPFtotal = MPFtotal.amt/compartment;
  MPFtotal.init = .016;
  p.amt = p*compartment;
  p = p.init;
  p.init = 1;
  when (time=time.min) Cdc20.amt = Cdc20.init*compartment;
  Cdc20.amt:time = v6.rate + -1*v7.rate + v8.rate + -1*v9.rate;
  Cdc20 = Cdc20.amt/compartment;
  Cdc20.init = .652;
  when (time=time.min) APCP.amt = APCP.init*compartment;
  APCP.amt:time = -1*v9.rate + v13.rate + -1*v14.rate + v15.rate;
  APCP = APCP.amt/compartment;
  APCP.init = .124;
  when (time=time.min) APCPCdc20.amt = APCPCdc20.init*compartment;
  APCPCdc20.amt:time = v9.rate + -1*v12.rate;
  APCPCdc20 = APCPCdc20.amt/compartment;
  APCPCdc20.init = .32;
  when (time=time.min) Mad2Cdc20P.amt = Mad2Cdc20P.init*compartment;
  Mad2Cdc20P.amt:time = v10.rate + -1*v11.rate;
  Mad2Cdc20P = Mad2Cdc20P.amt/compartment;
  Mad2Cdc20P.init = .009;
  when (time=time.min) Wee.amt = Wee.init*compartment;
  Wee.amt:time = v16.rate + -1*v17.rate;
  Wee = Wee.amt/compartment;
  Wee.init = .718;
  when (time=time.min) Cdc25P.amt = Cdc25P.init*compartment;
  Cdc25P.amt:time = v18.rate + -1*v19.rate;
  Cdc25P = Cdc25P.amt/compartment;
  Cdc25P.init = .281;
  v1.rate = kcycsyn;
  v2.rate = kwee*MPF;
  v3.rate = k25*(MPFtotal-MPF);
  v4.rate = (kacycdeg+kcycdeg*(APCPCdc20/(Jcycdeg+MPFtotal)))*MPF;
  v5.rate = (kacycdeg+kcycdeg*(APCPCdc20/(Jcycdeg+MPFtotal)))*MPFtotal;
  v6.rate = k20deph*Cdc20P;
  v7.rate = k20ph*MPF*Cdc20;
  v8.rate = (kapcdeph+kdiss)*APCPCdc20;
  v9.rate = kass*Cdc20*APCP;
  v10.rate = kass2*Mad2*Cdc20P;
  v11.rate = kdiss2*Mad2Cdc20P;
  v12.rate = (kdiss+kapcdeph+k20ph*MPF)*APCPCdc20;
  v13.rate = kapcph*MPF*APC;
  v14.rate = kapcdeph*APCP;
  v15.rate = (kdiss+k20ph*MPF)*APCPCdc20;
  v16.rate = Vawee*((Weetotal-Wee)/(Jwee+Weetotal-Wee));
  v17.rate = Vaaiwee*MPF*Wee;
  v18.rate = Vapa25*MPF*(Cdc25total-Cdc25P);
  v19.rate = Vi25*Cdc25P/(J25+Cdc25P);

  // variable properties
  compartment.sbmlRole="compartment";
  kwee.sbmlRole="parameter";
  k25.sbmlRole="parameter";
  APC.sbmlRole="parameter";
  Cdc20P.sbmlRole="parameter";
  Mad2.sbmlRole="parameter";
  kcycsyn.sbmlRole="parameter";
  kcycdeg.sbmlRole="parameter";
  kacycdeg.sbmlRole="parameter";
  kapcph.sbmlRole="parameter";
  kapcdeph.sbmlRole="parameter";
  k20deph.sbmlRole="parameter";
  k20ph.sbmlRole="parameter";
  kass.sbmlRole="parameter";
  kdiss.sbmlRole="parameter";
  kass2.sbmlRole="parameter";
  kdiss2.sbmlRole="parameter";
  ka25.sbmlRole="parameter";
  kaa25.sbmlRole="parameter";
  Vi25.sbmlRole="parameter";
  Vapa25.sbmlRole="parameter";
  kawee.sbmlRole="parameter";
  kaawee.sbmlRole="parameter";
  Vawee.sbmlRole="parameter";
  Vaaiwee.sbmlRole="parameter";
  Mad2active.sbmlRole="parameter";
  APCtotal.sbmlRole="parameter";
  Cdc20total.sbmlRole="parameter";
  Cdc25total.sbmlRole="parameter";
  Weetotal.sbmlRole="parameter";
  Jwee.sbmlRole="parameter";
  Jcycdeg.sbmlRole="parameter";
  J25.sbmlRole="parameter";
  s.amt.sbmlRole="speciesAmount";
  s.sbmlRole="speciesConcentration";
  s.sbmlCompartment="compartment";
  s.init.sbmlRole="speciesInitialConcentration";
  MPF.amt.sbmlRole="speciesAmount";
  MPF.sbmlRole="speciesConcentration";
  MPF.sbmlCompartment="compartment";
  MPF.init.sbmlRole="speciesInitialConcentration";
  MPFtotal.amt.sbmlRole="speciesAmount";
  MPFtotal.sbmlRole="speciesConcentration";
  MPFtotal.sbmlCompartment="compartment";
  MPFtotal.init.sbmlRole="speciesInitialConcentration";
  p.amt.sbmlRole="speciesAmount";
  p.sbmlRole="speciesConcentration";
  p.sbmlCompartment="compartment";
  p.init.sbmlRole="speciesInitialConcentration";
  Cdc20.amt.sbmlRole="speciesAmount";
  Cdc20.sbmlRole="speciesConcentration";
  Cdc20.sbmlCompartment="compartment";
  Cdc20.init.sbmlRole="speciesInitialConcentration";
  APCP.amt.sbmlRole="speciesAmount";
  APCP.sbmlRole="speciesConcentration";
  APCP.sbmlCompartment="compartment";
  APCP.init.sbmlRole="speciesInitialConcentration";
  APCPCdc20.amt.sbmlRole="speciesAmount";
  APCPCdc20.sbmlRole="speciesConcentration";
  APCPCdc20.sbmlCompartment="compartment";
  APCPCdc20.init.sbmlRole="speciesInitialConcentration";
  Mad2Cdc20P.amt.sbmlRole="speciesAmount";
  Mad2Cdc20P.sbmlRole="speciesConcentration";
  Mad2Cdc20P.sbmlCompartment="compartment";
  Mad2Cdc20P.init.sbmlRole="speciesInitialConcentration";
  Wee.amt.sbmlRole="speciesAmount";
  Wee.sbmlRole="speciesConcentration";
  Wee.sbmlCompartment="compartment";
  Wee.init.sbmlRole="speciesInitialConcentration";
  Cdc25P.amt.sbmlRole="speciesAmount";
  Cdc25P.sbmlRole="speciesConcentration";
  Cdc25P.sbmlCompartment="compartment";
  Cdc25P.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";
}