// 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
// reaction1: <=> C 
// reaction2: C 
// reaction3: C 
// reaction4: <=> M 
// reaction5: M 
// reaction6: <=> X 
// reaction7: X 

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

  // variable definitions
  real cell = 1 L;
  real V1(time);
  real V3(time);
  real VM1 = 3;
  real VM3 = 1;
  real Kc = .5;
  real C(time) M;
  real M(time) M;
  real X(time) M;
  real reaction1(time) katal;
  real vi = .025;
  real reaction2(time) katal;
  real kd = .01;
  real reaction3(time) katal;
  real vd = .25;
  real Kd = .02;
  real reaction4(time) katal;
  real K1 = .005;
  real reaction5(time) katal;
  real V2 = 1.5;
  real K2 = .005;
  real reaction6(time) katal;
  real K3 = .005;
  real reaction7(time) katal;
  real K4 = .005;
  real V4 = .5;

  // equations
  V1 = C*VM1*(C+Kc)^(-1);
  V3 = M*VM3;
  when (time=time.min) C = .01;
  C:time = (reaction1 + -1*reaction2 + -1*reaction3)/cell;
  when (time=time.min) M = .01;
  M:time = (reaction4 + -1*reaction5)/cell;
  when (time=time.min) X = .01;
  X:time = (reaction6 + -1*reaction7)/cell;
  reaction1 = cell*vi;
  reaction2 = C*cell*kd;
  reaction3 = C*cell*vd*X*(C+Kd)^(-1);
  reaction4 = cell*(1+(-1)*M)*V1*(K1+(-1)*M+1)^(-1);
  reaction5 = cell*M*V2*(K2+M)^(-1);
  reaction6 = cell*V3*(1+(-1)*X)*(K3+(-1)*X+1)^(-1);
  reaction7 = cell*V4*X*(K4+X)^(-1);

  // variable properties
  cell.sbmlRole="compartment";
  V1.sbmlRole="parameter";
  V3.sbmlRole="parameter";
  VM1.sbmlRole="parameter";
  VM3.sbmlRole="parameter";
  Kc.sbmlRole="parameter";
  C.sbmlRole="species";
  C.sbmlCompartment="cell";
  M.sbmlRole="species";
  M.sbmlCompartment="cell";
  X.sbmlRole="species";
  X.sbmlCompartment="cell";
  reaction1.sbmlRole="rate";
  vi.sbmlRole="parameter";
  reaction2.sbmlRole="rate";
  kd.sbmlRole="parameter";
  reaction3.sbmlRole="rate";
  vd.sbmlRole="parameter";
  Kd.sbmlRole="parameter";
  reaction4.sbmlRole="rate";
  K1.sbmlRole="parameter";
  reaction5.sbmlRole="rate";
  V2.sbmlRole="parameter";
  K2.sbmlRole="parameter";
  reaction6.sbmlRole="rate";
  K3.sbmlRole="parameter";
  reaction7.sbmlRole="rate";
  K4.sbmlRole="parameter";
  V4.sbmlRole="parameter";
}