// 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
// v1: alpha => beta 
// v2: beta => alpha 
// v3: beta => gamma 
// v4: gamma => beta 

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

  // variable definitions
  real compartment = 1 L;
  real r31 = 6;
  real r24 = 6;
  real Chi14 = 1.1;
  real Ks1 = .01;
  real Ks2 = .01;
  real Ks3 = .01;
  real Ks4 = .01;
  real Vm1 = 1;
  real alpha(time) M;
  real beta(time) M;
  real gamma(time) M;
  real v1(time) katal;
  real v2(time) katal;
  real v3(time) katal;
  real v4(time) katal;

  // equations
  when (time=time.min) alpha = .462;
  alpha:time = (-1*v1 + v2)/compartment;
  when (time=time.min) beta = .2;
  beta:time = (v1 + -1*v2 + -1*v3 + v4)/compartment;
  when (time=time.min) gamma = .338;
  gamma:time = (v3 + -1*v4)/compartment;
  v1 = Vm1*(alpha/Ks1)/(1+alpha/Ks1+beta/Ks3);
  v2 = r24*(Vm1/Chi14)*(beta/Ks2)/(1+gamma/Ks4+beta/Ks2);
  v3 = r31*Vm1*(beta/Ks3)/(1+alpha/Ks1+beta/Ks3);
  v4 = Vm1/Chi14*(gamma/Ks4)/(1+gamma/Ks4+beta/Ks2);

  // variable properties
  compartment.sbmlRole="compartment";
  r31.sbmlRole="parameter";
  r24.sbmlRole="parameter";
  Chi14.sbmlRole="parameter";
  Ks1.sbmlRole="parameter";
  Ks2.sbmlRole="parameter";
  Ks3.sbmlRole="parameter";
  Ks4.sbmlRole="parameter";
  Vm1.sbmlRole="parameter";
  alpha.sbmlRole="species";
  alpha.sbmlCompartment="compartment";
  beta.sbmlRole="species";
  beta.sbmlCompartment="compartment";
  gamma.sbmlRole="species";
  gamma.sbmlCompartment="compartment";
  v1.sbmlRole="rate";
  v2.sbmlRole="rate";
  v3.sbmlRole="rate";
  v4.sbmlRole="rate";
}