// 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
// v2: x1 => x3 
// v3: x1 => x2 
// v4: x2 => p1 
// v5: x3 => p2 
// v1: s => x1 

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

  // variable definitions
  real compartment = 1 L;
  real V1 = 1;
  real V2 = 1;
  real V3 = 1;
  real V4 = 1;
  real V5 = 1;
  real K = 1;
  real Keq = 1;
  real x1(time) M;
  real x3(time) M;
  real x2(time) M;
  real p1 = 1 M;
  real p2 = 1 M;
  real s = 10 M;
  real v2(time) katal;
  real v3(time) katal;
  real v4(time) katal;
  real v5(time) katal;
  real v1(time) katal;

  // equations
  when (time=time.min) x1 = 1;
  x1:time = (-1*v2 + -1*v3 + v1)/compartment;
  when (time=time.min) x3 = 1;
  x3:time = (v2 + -1*v5)/compartment;
  when (time=time.min) x2 = 1;
  x2:time = (v3 + -1*v4)/compartment;
  v2 = V2*x1*(1-x3/(Keq*x1))/(K*(1+x1/K+x3/K));
  v3 = V3*x1*(1-x2/(Keq*x1))/(K*(1+x1/K+x2/K));
  v4 = V4*(1-p1/(Keq*x2))*x2/(K*(1+p1/K+x2/K));
  v5 = V5*(1-p2/(Keq*x3))*x3/(K*(1+p2/K+x3/K));
  v1 = s*V1*(1-x1/(Keq*s))/(K*(1+s/K+x1/K));

  // variable properties
  compartment.sbmlRole="compartment";
  V1.sbmlRole="parameter";
  V2.sbmlRole="parameter";
  V3.sbmlRole="parameter";
  V4.sbmlRole="parameter";
  V5.sbmlRole="parameter";
  K.sbmlRole="parameter";
  Keq.sbmlRole="parameter";
  x1.sbmlRole="species";
  x1.sbmlCompartment="compartment";
  x3.sbmlRole="species";
  x3.sbmlCompartment="compartment";
  x2.sbmlRole="species";
  x2.sbmlCompartment="compartment";
  p1.sbmlRole="species";
  p1.sbmlCompartment="compartment";
  p2.sbmlRole="species";
  p2.sbmlCompartment="compartment";
  s.sbmlRole="species";
  s.sbmlCompartment="compartment";
  v2.sbmlRole="rate";
  v3.sbmlRole="rate";
  v4.sbmlRole="rate";
  v5.sbmlRole="rate";
  v1.sbmlRole="rate";
}