// 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
// Reaction1: M <=> C2 YP 
// Reaction2: C2 <=> CP 
// Reaction3: CP <=> C2 
// Reaction4: CP Y <=> pM 
// Reaction5: M <=> pM 
// Reaction6: EmptySet <=> Y 
// Reaction7: Y <=> EmptySet 
// Reaction8: YP <=> EmptySet 
// Reaction9: pM <=> M 

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

  // variable definitions
  real cell = 1 volume;
  real EmptySet substance;
  private real EmptySet.conc(time) substance/volume;
  real EmptySet.init substance;
  real C2(time) substance;
  private real C2.conc(time) substance/volume;
  real C2.init substance;
  real CP(time) substance;
  private real CP.conc(time) substance/volume;
  real CP.init substance;
  real M(time) substance;
  private real M.conc(time) substance/volume;
  real M.init substance;
  real pM(time) substance;
  private real pM.conc(time) substance/volume;
  real pM.init substance;
  real Y(time) substance;
  private real Y.conc(time) substance/volume;
  real Y.init substance;
  real YP(time) substance;
  private real YP.conc(time) substance/volume;
  real YP.init substance;
  real Reaction1.rate(time) substance/time;
  real k6 = 1;
  real Reaction2.rate(time) substance/time;
  real k8notP = 1E6;
  real Reaction3.rate(time) substance/time;
  real k9 = 1E3;
  real Reaction4.rate(time) substance/time;
  real k3 = 200;
  real Reaction5.rate(time) substance/time;
  real k5notP = 0;
  real Reaction6.rate(time) substance/time;
  real k1aa = .015;
  real Reaction7.rate(time) substance/time;
  real k2 = 0;
  real Reaction8.rate(time) substance/time;
  real k7 = .6;
  real Reaction9.rate(time) substance/time;
  real k4 = 180;
  real k4prime = .018;

  // equations
  EmptySet = EmptySet.init;
  EmptySet.conc = EmptySet/cell;
  EmptySet.init = 0;
  when (time=time.min) C2 = C2.init;
  C2:time = Reaction1.rate + -1*Reaction2.rate + Reaction3.rate;
  C2.conc = C2/cell;
  C2.init = 0;
  when (time=time.min) CP = CP.init;
  CP:time = Reaction2.rate + -1*Reaction3.rate + -1*Reaction4.rate;
  CP.conc = CP/cell;
  CP.init = 1;
  when (time=time.min) M = M.init;
  M:time = -1*Reaction1.rate + -1*Reaction5.rate + Reaction9.rate;
  M.conc = M/cell;
  M.init = 0;
  when (time=time.min) pM = pM.init;
  pM:time = Reaction4.rate + Reaction5.rate + -1*Reaction9.rate;
  pM.conc = pM/cell;
  pM.init = .3;
  when (time=time.min) Y = Y.init;
  Y:time = -1*Reaction4.rate + Reaction6.rate + -1*Reaction7.rate;
  Y.conc = Y/cell;
  Y.init = 0;
  when (time=time.min) YP = YP.init;
  YP:time = Reaction1.rate + -1*Reaction8.rate;
  YP.conc = YP/cell;
  YP.init = 0;
  Reaction1.rate = k6*M;
  Reaction2.rate = C2*k8notP;
  Reaction3.rate = CP*k9;
  Reaction4.rate = CP*k3*Y;
  Reaction5.rate = k5notP*M;
  Reaction6.rate = k1aa;
  Reaction7.rate = k2*Y;
  Reaction8.rate = k7*YP;
  Reaction9.rate = pM*(k4prime+k4*M^2);

  // variable properties
  cell.sbmlRole="compartment";
  EmptySet.sbmlRole="speciesAmount";
  EmptySet.sbmlCompartment="cell";
  EmptySet.conc.sbmlRole="speciesConcentration";
  EmptySet.init.sbmlRole="speciesInitialAmount";
  C2.sbmlRole="speciesAmount";
  C2.sbmlName="cdc2k";
  C2.sbmlCompartment="cell";
  C2.conc.sbmlRole="speciesConcentration";
  C2.init.sbmlRole="speciesInitialAmount";
  CP.sbmlRole="speciesAmount";
  CP.sbmlName="cdc2k-P";
  CP.sbmlCompartment="cell";
  CP.conc.sbmlRole="speciesConcentration";
  CP.init.sbmlRole="speciesInitialAmount";
  M.sbmlRole="speciesAmount";
  M.sbmlName="p-cyclin_cdc2";
  M.sbmlCompartment="cell";
  M.conc.sbmlRole="speciesConcentration";
  M.init.sbmlRole="speciesInitialAmount";
  pM.sbmlRole="speciesAmount";
  pM.sbmlName="p-cyclin_cdc2-p";
  pM.sbmlCompartment="cell";
  pM.conc.sbmlRole="speciesConcentration";
  pM.init.sbmlRole="speciesInitialAmount";
  Y.sbmlRole="speciesAmount";
  Y.sbmlName="cyclin";
  Y.sbmlCompartment="cell";
  Y.conc.sbmlRole="speciesConcentration";
  Y.init.sbmlRole="speciesInitialAmount";
  YP.sbmlRole="speciesAmount";
  YP.sbmlName="p-cyclin";
  YP.sbmlCompartment="cell";
  YP.conc.sbmlRole="speciesConcentration";
  YP.init.sbmlRole="speciesInitialAmount";
  Reaction1.rate.sbmlRole="rate";
  k6.sbmlRole="parameter";
  Reaction2.rate.sbmlRole="rate";
  k8notP.sbmlRole="parameter";
  Reaction3.rate.sbmlRole="rate";
  k9.sbmlRole="parameter";
  Reaction4.rate.sbmlRole="rate";
  k3.sbmlRole="parameter";
  Reaction5.rate.sbmlRole="rate";
  k5notP.sbmlRole="parameter";
  Reaction6.rate.sbmlRole="rate";
  k1aa.sbmlRole="parameter";
  Reaction7.rate.sbmlRole="rate";
  k2.sbmlRole="parameter";
  Reaction8.rate.sbmlRole="rate";
  k7.sbmlRole="parameter";
  Reaction9.rate.sbmlRole="rate";
  k4.sbmlRole="parameter";
  k4prime.sbmlRole="parameter";
}