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

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