// This model generated automatically from SBML

// unit definitions
import nsrunit;
unit conversion off;
unit item=scalar;
unit substance = mole;
unit volume = litre;
unit area = metre^2;
unit length = metre;
unit time = second;

// 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 volume;
  real k1f = 2;
  real k1r = 1;
  real k2f = 1;
  real k2r = 1;
  real k3f = 5;
  real k3r = 1;
  private real S1.amt(time) substance;
  real S1(time) substance/volume;
  real S1.init substance/volume;
  private real X0.amt substance;
  real X0 substance/volume;
  real X0.init substance/volume;
  private real S2.amt(time) substance;
  real S2(time) substance/volume;
  real S2.init substance/volume;
  private real S3.amt(time) substance;
  real S3(time) substance/volume;
  real S3.init substance/volume;
  private real X4.amt substance;
  real X4 substance/volume;
  real X4.init substance/volume;
  real vR1.rate(time) substance/time;
  real vR2.rate(time) substance/time;
  real vR3.rate(time) substance/time;

  // equations
  when (time=time.min) S1.amt = S1.init*compartment;
  S1.amt:time = -1*vR1.rate + vR3.rate;
  S1 = S1.amt/compartment;
  S1.init = 1;
  X0.amt = X0*compartment;
  X0 = X0.init;
  X0.init = 10;
  when (time=time.min) S2.amt = S2.init*compartment;
  S2.amt:time = vR1.rate + -1*vR2.rate;
  S2 = S2.amt/compartment;
  S2.init = 1;
  when (time=time.min) S3.amt = S3.init*compartment;
  S3.amt:time = vR2.rate + -1*vR3.rate;
  S3 = S3.amt/compartment;
  S3.init = 1;
  X4.amt = X4*compartment;
  X4 = X4.init;
  X4.init = 1;
  vR1.rate = k1f*X0*S1-k1r*S2;
  vR2.rate = k2f*S2-k2r*S3;
  vR3.rate = 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.amt.sbmlRole="speciesAmount";
  S1.sbmlRole="speciesConcentration";
  S1.sbmlCompartment="compartment";
  S1.init.sbmlRole="speciesInitialConcentration";
  X0.amt.sbmlRole="speciesAmount";
  X0.sbmlRole="speciesConcentration";
  X0.sbmlCompartment="compartment";
  X0.init.sbmlRole="speciesInitialConcentration";
  S2.amt.sbmlRole="speciesAmount";
  S2.sbmlRole="speciesConcentration";
  S2.sbmlCompartment="compartment";
  S2.init.sbmlRole="speciesInitialConcentration";
  S3.amt.sbmlRole="speciesAmount";
  S3.sbmlRole="speciesConcentration";
  S3.sbmlCompartment="compartment";
  S3.init.sbmlRole="speciesInitialConcentration";
  X4.amt.sbmlRole="speciesAmount";
  X4.sbmlRole="speciesConcentration";
  X4.sbmlCompartment="compartment";
  X4.init.sbmlRole="speciesInitialConcentration";
  vR1.rate.sbmlRole="rate";
  vR2.rate.sbmlRole="rate";
  vR3.rate.sbmlRole="rate";
}