// 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: ERG => p 
// v2: DRG => p 
// v3: CD => P27 
// v4: CYCD => p 
// v5: P27 CYCE => CE 
// v6: P27 CYCA => CA 
// v7: P27 CYCD => CD 
// v8: CD => P27 CYCD 
// v9: CYCA => p 
// v10: CA => P27 
// v11: CE => P27 CYCE 
// v12: CA => P27 CYCA 
// v13: CE => P27 
// v14: CYCE => p 
// v15: P27 => p 
// v16: CE => CYCE 
// v17: CD => CYCD 
// v18: CA => CYCA 
// v19: CYCB => p 
// v20: s => CDh1 
// v21: CDh1 => p 
// v22: PPX => p 
// v23: s => IEP 
// v24: IEP => p 
// v25: CDc20T => p 
// v26: s => CDc20 
// v27: CDc20 => p 
// v28: CDc20 => p 
// v29: s => E2F 
// v30: E2F => p 
// v31: s => GM 
// v32: GM => p 
// v33: s => MASS 
// v34: s => ERG 
// v35: s => CDc20T 
// v36: s => CYCA 
// v37: s => PPX 
// v38: s => CYCE 
// v39: s => CYCD 
// v40: s => P27 
// v41: s => DRG 
// v42: s => CYCB 

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

  // variable definitions
  real compartment = 1 L;
  real PP1A(time);
  real V2(time);
  real V4(time);
  real V6(time);
  real CYCET(time);
  real CYCDT(time);
  real CYCAT(time);
  real P27T(time);
  real V8(time);
  real L(time);
  real RBH(time);
  real E2RBC(time);
  real E2FA(time);
  real k15 = .025;
  real k16 = .025;
  real J15 = .1;
  real k17a = .035;
  real k17 = 1;
  real J17 = .3;
  real k18 = 1;
  real K9 = .25;
  real K10 = .5;
  real k24 = 100;
  real k24r = 1;
  real K7a = 0;
  real K7 = .06;
  real K8a = .01;
  real K8 = .2;
  real K25 = 100;
  real K25R = 1;
  real J8 = .1;
  real YE = 1;
  real YB = .05;
  real K29 = .005;
  real K30 = 2;
  real K1a = .01;
  real K1 = .06;
  real J1 = .1;
  real K2a = .005;
  real K2 = 2;
  real K2aa = .1;
  real K5 = 2;
  real K6a = 1;
  real K6 = 10;
  real HE = .5;
  real HB = 1;
  real HA = .5;
  real RBT = 10;
  real LD = 3.3;
  real LE = 5;
  real LB = 5;
  real LA = 3;
  real K20 = 1;
  real K19a = 0;
  real K19 = 2;
  real K21 = 1;
  real PP1T = 1;
  real FE = 25;
  real FB = 2;
  real K3a = .75;
  real K3 = 14;
  real J3 = .01;
  real J4 = .01;
  real K4 = 4;
  real GE = 0;
  real GB = 1;
  real GA = .3;
  real K33 = .005;
  real K34 = .005;
  real K31 = .07;
  real K32 = .18;
  real J31 = .01;
  real J32 = .01;
  real K11a = 0;
  real K11 = .15;
  real K12 = .15;
  real K13 = .5;
  real K14 = .25;
  real J13 = .005;
  real J14 = .005;
  real E2FT = 5;
  real K22 = .1;
  real K23a = 5E-4;
  real K23 = .1;
  real K26 = 1E3;
  real K26R = 20;
  real K27 = .02;
  real K28 = .02;
  real MU = .0061;
  real eps = 1;
  real ERG(time) M;
  real p = 1 M;
  real DRG(time) M;
  real CD(time) M;
  real P27(time) M;
  real CYCD(time) M;
  real CYCE(time) M;
  real CE(time) M;
  real CYCA(time) M;
  real CA(time) M;
  real CYCB(time) M;
  real s = 1 M;
  real CDh1(time) M;
  real PPX(time) M;
  real IEP(time) M;
  real CDc20T(time) M;
  real CDc20(time) M;
  real E2F(time) M;
  real GM(time) M;
  real MASS(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;
  real v16(time) 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;
  real v26(time) katal;
  real v27(time) katal;
  real v28(time) katal;
  real v29(time) katal;
  real v30(time) katal;
  extern real v31 katal;
  real v32(time) katal;
  real v33(time) katal;
  real v34(time) katal;
  real v35(time) katal;
  real v36(time) katal;
  real v37(time) katal;
  real v38(time) katal;
  real v39(time) katal;
  real v40(time) katal;
  real v41(time) katal;
  real v42(time) katal;

  // equations
  PP1A = PP1T/(K21*(FE*(CYCE+CYCA)+FB*CYCB)+1);
  V2 = K2a*(1-CDh1)+K2*CDh1+K2aa*CDc20;
  V4 = K4*(GE*CYCE+GA*CYCA+GB*CYCB);
  V6 = K6a+K6*(HE*CYCE+HA*CYCA+HB*CYCB);
  CYCET = CYCE+CE;
  CYCDT = CYCD+CD;
  CYCAT = CYCA+CA;
  P27T = P27+CD+CE+CA;
  V8 = K8a+K8*(YE*(CYCE+CYCA)+YB*CYCB)/(J8+CYCET);
  L = (K26R+K20*(LD*CYCDT+LE*CYCE+LA*CYCA+LB*CYCB))/K26;
  RBH = RBT/(K20*(LD*CYCDT+LE*CYCE+LA*CYCA+LB*CYCB)/(K19a*(PP1T-PP1A)+K19*PP1A)+1);
  E2RBC = 2*E2FT*RBH/(E2FT+RBH+L+((E2FT+RBH+L)^2-4*E2FT*RBH)^(1/2));
  E2FA = (E2FT-E2RBC)*E2F/E2FT;
  when (time=time.min) ERG = .0121808;
  ERG:time = (-1*v1 + v34)/compartment;
  when (time=time.min) DRG = .9005328;
  DRG:time = (-1*v2 + v41)/compartment;
  when (time=time.min) CD = .0488772;
  CD:time = (-1*v3 + v7 + -1*v8 + -1*v17)/compartment;
  when (time=time.min) P27 = .011877;
  P27:time = (v3 + -1*v5 + -1*v6 + -1*v7 + v8 + v10 + v11 + v12 + v13 + -1*v15 + v40)/compartment;
  when (time=time.min) CYCD = .401389;
  CYCD:time = (-1*v4 + -1*v7 + v8 + v17 + v39)/compartment;
  when (time=time.min) CYCE = 1.09959;
  CYCE:time = (-1*v5 + v11 + -1*v14 + v16 + v38)/compartment;
  when (time=time.min) CE = .1380924;
  CE:time = (v5 + -1*v11 + -1*v13 + -1*v16)/compartment;
  when (time=time.min) CYCA = .3411;
  CYCA:time = (-1*v6 + -1*v9 + v12 + v18 + v36)/compartment;
  when (time=time.min) CA = .0434486;
  CA:time = (v6 + -1*v10 + -1*v12 + -1*v18)/compartment;
  when (time=time.min) CYCB = .00514;
  CYCB:time = (-1*v19 + v42)/compartment;
  when (time=time.min) CDh1 = .98702;
  CDh1:time = (v20 + -1*v21)/compartment;
  when (time=time.min) PPX = 1;
  PPX:time = (-1*v22 + v37)/compartment;
  when (time=time.min) IEP = 1.98E-5;
  IEP:time = (v23 + -1*v24)/compartment;
  when (time=time.min) CDc20T = .006466;
  CDc20T:time = (-1*v25 + v35)/compartment;
  when (time=time.min) CDc20 = 1E-7;
  CDc20:time = (v26 + -1*v27 + -1*v28)/compartment;
  when (time=time.min) E2F = 4.25991;
  E2F:time = (v29 + -1*v30)/compartment;
  when (time=time.min) GM = 1.032033;
  GM:time = (v31 + -1*v32)/compartment;
  when (time=time.min) MASS = 1.28486;
  MASS:time = (v33)/compartment;
  v1 = k16*ERG;
  v2 = k18*DRG;
  v3 = K10*CD;
  v4 = K10*CYCD;
  v5 = K25*P27*CYCE;
  v6 = K25*P27*CYCA;
  v7 = k24*CYCD*P27;
  v8 = k24r*CD;
  v9 = K30*CDc20*CYCA;
  v10 = K30*CDc20*CA;
  v11 = K25R*CE;
  v12 = K25R*CA;
  v13 = V8*CE;
  v14 = V8*CYCE;
  v15 = V6*P27;
  v16 = V6*CE;
  v17 = V6*CD;
  v18 = V6*CA;
  v19 = V2*CYCB;
  v20 = (K3a+K3*CDc20)*(1-CDh1)/(J3+1-CDh1);
  v21 = V4*CDh1/(J4+CDh1);
  v22 = K34*PPX;
  v23 = K31*CYCB*(1-IEP)/(J31+1-IEP);
  v24 = K32*PPX*IEP/(J32+IEP);
  v25 = K12*CDc20T;
  v26 = K13*IEP*(CDc20T-CDc20)/(J13+CDc20T-CDc20);
  v27 = K14*CDc20/(J14+CDc20);
  v28 = K12*CDc20;
  v29 = K22*(E2FT-E2F);
  v30 = (K23a+K23*(CYCA+CYCB))*E2F;
  v32 = K28*GM;
  v33 = eps*MU*GM;
  v34 = eps*k15/(1+(DRG/J15)^2);
  v35 = eps*(K11a+K11*CYCB);
  v36 = eps*K29*E2FA*MASS;
  v37 = eps*K33;
  v38 = eps*(K7a+K7*E2FA);
  v39 = eps*K9*DRG;
  v40 = eps*K5;
  v41 = eps*(k17a*ERG+k17*(DRG/J17)^2/(1+(DRG/J17)^2));
  v42 = eps*(K1a+K1*(CYCB/J1)^2/(1+(CYCB/J1)^2));

  // variable properties
  compartment.sbmlRole="compartment";
  PP1A.sbmlRole="parameter";
  V2.sbmlRole="parameter";
  V4.sbmlRole="parameter";
  V6.sbmlRole="parameter";
  CYCET.sbmlRole="parameter";
  CYCDT.sbmlRole="parameter";
  CYCAT.sbmlRole="parameter";
  P27T.sbmlRole="parameter";
  V8.sbmlRole="parameter";
  L.sbmlRole="parameter";
  RBH.sbmlRole="parameter";
  E2RBC.sbmlRole="parameter";
  E2FA.sbmlRole="parameter";
  k15.sbmlRole="parameter";
  k16.sbmlRole="parameter";
  J15.sbmlRole="parameter";
  k17a.sbmlRole="parameter";
  k17.sbmlRole="parameter";
  J17.sbmlRole="parameter";
  k18.sbmlRole="parameter";
  K9.sbmlRole="parameter";
  K10.sbmlRole="parameter";
  k24.sbmlRole="parameter";
  k24r.sbmlRole="parameter";
  K7a.sbmlRole="parameter";
  K7.sbmlRole="parameter";
  K8a.sbmlRole="parameter";
  K8.sbmlRole="parameter";
  K25.sbmlRole="parameter";
  K25R.sbmlRole="parameter";
  J8.sbmlRole="parameter";
  YE.sbmlRole="parameter";
  YB.sbmlRole="parameter";
  K29.sbmlRole="parameter";
  K30.sbmlRole="parameter";
  K1a.sbmlRole="parameter";
  K1.sbmlRole="parameter";
  J1.sbmlRole="parameter";
  K2a.sbmlRole="parameter";
  K2.sbmlRole="parameter";
  K2aa.sbmlRole="parameter";
  K5.sbmlRole="parameter";
  K6a.sbmlRole="parameter";
  K6.sbmlRole="parameter";
  HE.sbmlRole="parameter";
  HB.sbmlRole="parameter";
  HA.sbmlRole="parameter";
  RBT.sbmlRole="parameter";
  LD.sbmlRole="parameter";
  LE.sbmlRole="parameter";
  LB.sbmlRole="parameter";
  LA.sbmlRole="parameter";
  K20.sbmlRole="parameter";
  K19a.sbmlRole="parameter";
  K19.sbmlRole="parameter";
  K21.sbmlRole="parameter";
  PP1T.sbmlRole="parameter";
  FE.sbmlRole="parameter";
  FB.sbmlRole="parameter";
  K3a.sbmlRole="parameter";
  K3.sbmlRole="parameter";
  J3.sbmlRole="parameter";
  J4.sbmlRole="parameter";
  K4.sbmlRole="parameter";
  GE.sbmlRole="parameter";
  GB.sbmlRole="parameter";
  GA.sbmlRole="parameter";
  K33.sbmlRole="parameter";
  K34.sbmlRole="parameter";
  K31.sbmlRole="parameter";
  K32.sbmlRole="parameter";
  J31.sbmlRole="parameter";
  J32.sbmlRole="parameter";
  K11a.sbmlRole="parameter";
  K11.sbmlRole="parameter";
  K12.sbmlRole="parameter";
  K13.sbmlRole="parameter";
  K14.sbmlRole="parameter";
  J13.sbmlRole="parameter";
  J14.sbmlRole="parameter";
  E2FT.sbmlRole="parameter";
  K22.sbmlRole="parameter";
  K23a.sbmlRole="parameter";
  K23.sbmlRole="parameter";
  K26.sbmlRole="parameter";
  K26R.sbmlRole="parameter";
  K27.sbmlRole="parameter";
  K28.sbmlRole="parameter";
  MU.sbmlRole="parameter";
  eps.sbmlRole="parameter";
  ERG.sbmlRole="species";
  ERG.sbmlCompartment="compartment";
  p.sbmlRole="species";
  p.sbmlCompartment="compartment";
  DRG.sbmlRole="species";
  DRG.sbmlCompartment="compartment";
  CD.sbmlRole="species";
  CD.sbmlCompartment="compartment";
  P27.sbmlRole="species";
  P27.sbmlCompartment="compartment";
  CYCD.sbmlRole="species";
  CYCD.sbmlCompartment="compartment";
  CYCE.sbmlRole="species";
  CYCE.sbmlCompartment="compartment";
  CE.sbmlRole="species";
  CE.sbmlCompartment="compartment";
  CYCA.sbmlRole="species";
  CYCA.sbmlCompartment="compartment";
  CA.sbmlRole="species";
  CA.sbmlCompartment="compartment";
  CYCB.sbmlRole="species";
  CYCB.sbmlCompartment="compartment";
  s.sbmlRole="species";
  s.sbmlCompartment="compartment";
  CDh1.sbmlRole="species";
  CDh1.sbmlCompartment="compartment";
  PPX.sbmlRole="species";
  PPX.sbmlCompartment="compartment";
  IEP.sbmlRole="species";
  IEP.sbmlCompartment="compartment";
  CDc20T.sbmlRole="species";
  CDc20T.sbmlCompartment="compartment";
  CDc20.sbmlRole="species";
  CDc20.sbmlCompartment="compartment";
  E2F.sbmlRole="species";
  E2F.sbmlCompartment="compartment";
  GM.sbmlRole="species";
  GM.sbmlCompartment="compartment";
  MASS.sbmlRole="species";
  MASS.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";
  v26.sbmlRole="rate";
  v27.sbmlRole="rate";
  v28.sbmlRole="rate";
  v29.sbmlRole="rate";
  v30.sbmlRole="rate";
  v31.sbmlRole="rate";
  v32.sbmlRole="rate";
  v33.sbmlRole="rate";
  v34.sbmlRole="rate";
  v35.sbmlRole="rate";
  v36.sbmlRole="rate";
  v37.sbmlRole="rate";
  v38.sbmlRole="rate";
  v39.sbmlRole="rate";
  v40.sbmlRole="rate";
  v41.sbmlRole="rate";
  v42.sbmlRole="rate";
}