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

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 = 1;
  real k3r = 1;
  real X0 = 10 M;
  real S1(time) M;
  real S2(time) M;
  real X5 = 1 M;
  real S3(time) M;
  real S4(time) M;
  real X6 = 10 M;
  real X7 = 1 M;
  real vR1(time) katal;
  real vR2(time) katal;
  real vR3(time) katal;

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

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