// This model generated automatically from SBML

// unit definitions
import nsrunit;
unit conversion off;

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

// SBML reactions
// v1: s => CycB 
// v2: CycB => p 
// v3: s => Cdh1 
// v4: Cdh1 => p 
// v5: s => Cdc20T 
// v6: Cdc20T => p 
// v7: s => Cdc20A 
// v8: Cdc20A => p 
// v9: Cdc20A => p 
// v10: s => IEP 
// v11: IEP => p 
// v12: s => m 

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

  // variable definitions
  real compartment = 1 L;
  real k1 = .04;
  real k21 = .04;
  real k22 = 1;
  real k31 = 1;
  real k32 = 10;
  real k4 = 35;
  real k51 = .005;
  real k52 = .2;
  real k6 = .1;
  real k7 = 1;
  real k8 = .5;
  real k9 = .1;
  real k10 = .02;
  real J3 = .04;
  real J4 = .04;
  real J5 = .3;
  real J7 = .001;
  real J8 = .001;
  real MAD = 1;
  real mu = .01;
  real n = 4;
  real s = 1 M;
  real CycB(time) M;
  real p = 1 M;
  real Cdh1(time) M;
  real Cdc20T(time) M;
  real Cdc20A(time) M;
  real IEP(time) M;
  real m(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;

  // equations
  when (time=time.min) CycB = .63737193;
  CycB:time = (v1 + -1*v2)/compartment;
  when (time=time.min) Cdh1 = .03924449;
  Cdh1:time = (v3 + -1*v4)/compartment;
  when (time=time.min) Cdc20T = 1.6993851;
  Cdc20T:time = (v5 + -1*v6)/compartment;
  when (time=time.min) Cdc20A = .93546909;
  Cdc20A:time = (v7 + -1*v8 + -1*v9)/compartment;
  when (time=time.min) IEP = .66387991;
  IEP:time = (v10 + -1*v11)/compartment;
  when (time=time.min) m = .89928583;
  m:time = (v12)/compartment;
  v1 = k1;
  v2 = (k21+k22*Cdh1)*CycB;
  v3 = (k31+k32*Cdc20A)*(1-Cdh1)/(J3+1-Cdh1);
  v4 = k4*m*CycB*Cdh1/(J4+Cdh1);
  v5 = k51+k52*(CycB*m/J5)^n/(1+(CycB*m/J5)^n);
  v6 = k6*Cdc20T;
  v7 = k7*IEP*(Cdc20T-Cdc20A)/(J7+Cdc20T-Cdc20A);
  v8 = k8*MAD*Cdc20A/(J8+Cdc20A);
  v9 = k6*Cdc20A;
  v10 = k9*m*CycB*(1-IEP);
  v11 = k10*IEP;
  v12 = mu*m;

  // variable properties
  compartment.sbmlRole="compartment";
  k1.sbmlRole="parameter";
  k21.sbmlRole="parameter";
  k22.sbmlRole="parameter";
  k31.sbmlRole="parameter";
  k32.sbmlRole="parameter";
  k4.sbmlRole="parameter";
  k51.sbmlRole="parameter";
  k52.sbmlRole="parameter";
  k6.sbmlRole="parameter";
  k7.sbmlRole="parameter";
  k8.sbmlRole="parameter";
  k9.sbmlRole="parameter";
  k10.sbmlRole="parameter";
  J3.sbmlRole="parameter";
  J4.sbmlRole="parameter";
  J5.sbmlRole="parameter";
  J7.sbmlRole="parameter";
  J8.sbmlRole="parameter";
  MAD.sbmlRole="parameter";
  mu.sbmlRole="parameter";
  n.sbmlRole="parameter";
  s.sbmlRole="species";
  s.sbmlCompartment="compartment";
  CycB.sbmlRole="species";
  CycB.sbmlCompartment="compartment";
  p.sbmlRole="species";
  p.sbmlCompartment="compartment";
  Cdh1.sbmlRole="species";
  Cdh1.sbmlCompartment="compartment";
  Cdc20T.sbmlRole="species";
  Cdc20T.sbmlCompartment="compartment";
  Cdc20A.sbmlRole="species";
  Cdc20A.sbmlCompartment="compartment";
  IEP.sbmlRole="species";
  IEP.sbmlCompartment="compartment";
  m.sbmlRole="species";
  m.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";
}