// 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
// vR1: S2 X0 => S3 S1 
// vR2: S3 X6 => X7 S2 
// vR3: S1 => X4 
// vR4: S1 => X5 

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

  // variable definitions
  real compartment = 1 L;
  real k1f = 1;
  real k1r = 1;
  real k2f = 1;
  real k2r = 1;
  real k3f = 5;
  real k3r = 1;
  real k4f = 1;
  real k4r = 1;
  real S2(time) M;
  real X0 = 10 M;
  real S3(time) M;
  real S1(time) M;
  real X6 = 5 M;
  real X7 = 1 M;
  real X4 = 2 M;
  real X5 = 1 M;
  real vR1(time) katal;
  real vR2(time) katal;
  real vR3(time) katal;
  real vR4(time) katal;

  // equations
  when (time=time.min) S2 = 1;
  S2:time = (-1*vR1 + vR2)/compartment;
  when (time=time.min) S3 = 1;
  S3:time = (vR1 + -1*vR2)/compartment;
  when (time=time.min) S1 = 1;
  S1:time = (vR1 + -1*vR3 + -1*vR4)/compartment;
  vR1 = k1f*X0*S2-k1r*S1*S3;
  vR2 = k2f*X6*S3-k2r*S2*X7;
  vR3 = k3f*S1-k3r*X4;
  vR4 = k4f*S1-k4r*X5;

  // variable properties
  compartment.sbmlRole="compartment";
  k1f.sbmlRole="parameter";
  k1r.sbmlRole="parameter";
  k2f.sbmlRole="parameter";
  k2r.sbmlRole="parameter";
  k3f.sbmlRole="parameter";
  k3r.sbmlRole="parameter";
  k4f.sbmlRole="parameter";
  k4r.sbmlRole="parameter";
  S2.sbmlRole="species";
  S2.sbmlCompartment="compartment";
  X0.sbmlRole="species";
  X0.sbmlCompartment="compartment";
  S3.sbmlRole="species";
  S3.sbmlCompartment="compartment";
  S1.sbmlRole="species";
  S1.sbmlCompartment="compartment";
  X6.sbmlRole="species";
  X6.sbmlCompartment="compartment";
  X7.sbmlRole="species";
  X7.sbmlCompartment="compartment";
  X4.sbmlRole="species";
  X4.sbmlCompartment="compartment";
  X5.sbmlRole="species";
  X5.sbmlCompartment="compartment";
  vR1.sbmlRole="rate";
  vR2.sbmlRole="rate";
  vR3.sbmlRole="rate";
  vR4.sbmlRole="rate";
}