// 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: S1 X0 => S2 
// vR2: S2 => S3 
// vR3: S3 => S1 X4 

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

  // variable definitions
  real compartment = 1 L;
  real k1f = 2;
  real k1r = 1;
  real k2f = 1;
  real k2r = 1;
  real k3f = 5;
  real k3r = 1;
  real S1(time) M;
  real X0 = 10 M;
  real S2(time) M;
  real S3(time) M;
  real X4 = 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 + vR3)/compartment;
  when (time=time.min) S2 = 1;
  S2:time = (vR1 + -1*vR2)/compartment;
  when (time=time.min) S3 = 1;
  S3:time = (vR2 + -1*vR3)/compartment;
  vR1 = k1f*X0*S1-k1r*S2;
  vR2 = k2f*S2-k2r*S3;
  vR3 = k3f*S3-k3r*X4*S1;

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