// 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: s => a 
// vR2: a => b 
// vR3: b => p 

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

  // variable definitions
  real compartment = 1 L;
  real V1 = 200;
  real eq = 1;
  real Keq1 = 400;
  real Ks = .01;
  real Ka = 1E3;
  real Kp = .1;
  real Cp = .026;
  real n = 4;
  real L0 = 1E4;
  real s05 = 1;
  real a05 = 1E4;
  real p05 = 1;
  real h = 4;
  real alfa = .01;
  real V2 = 100;
  real Keq2 = 10;
  real K2a = 1;
  real K2b = 1;
  real V3 = 100;
  real Keq3 = 10;
  real K3b = 1;
  real K3p = 1;
  real s = 1 M;
  real a(time) M;
  real b(time) M;
  real p = 1 M;
  extern real vR1 katal;
  real vR2(time) katal;
  real vR3(time) katal;

  // equations
  when (time=time.min) a = 10;
  a:time = (vR1 + -1*vR2)/compartment;
  when (time=time.min) b = 10;
  b:time = (vR2 + -1*vR3)/compartment;
  vR2 = V2*(a-b/Keq2)/(a+K2a*(1+b/K2b));
  vR3 = V3*(b-p/Keq3)/(b+K3b*(1+p/K3p));

  // variable properties
  compartment.sbmlRole="compartment";
  V1.sbmlRole="parameter";
  eq.sbmlRole="parameter";
  Keq1.sbmlRole="parameter";
  Ks.sbmlRole="parameter";
  Ka.sbmlRole="parameter";
  Kp.sbmlRole="parameter";
  Cp.sbmlRole="parameter";
  n.sbmlRole="parameter";
  L0.sbmlRole="parameter";
  s05.sbmlRole="parameter";
  a05.sbmlRole="parameter";
  p05.sbmlRole="parameter";
  h.sbmlRole="parameter";
  alfa.sbmlRole="parameter";
  V2.sbmlRole="parameter";
  Keq2.sbmlRole="parameter";
  K2a.sbmlRole="parameter";
  K2b.sbmlRole="parameter";
  V3.sbmlRole="parameter";
  Keq3.sbmlRole="parameter";
  K3b.sbmlRole="parameter";
  K3p.sbmlRole="parameter";
  s.sbmlRole="species";
  s.sbmlCompartment="compartment";
  a.sbmlRole="species";
  a.sbmlCompartment="compartment";
  b.sbmlRole="species";
  b.sbmlCompartment="compartment";
  p.sbmlRole="species";
  p.sbmlCompartment="compartment";
  vR1.sbmlRole="rate";
  vR2.sbmlRole="rate";
  vR3.sbmlRole="rate";
}