// 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: s => Enz 
// v2: Enz => s 
// v3: s => Ts 
// v4: Ts => s 
// v5: Ts => s 

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

  // variable definitions
  real compartment = 1 L;
  real To(time);
  real Tt(time);
  real ki2 = 820;
  real ki1 = 3.53;
  real nH = 1.92;
  real mu = .01;
  real gval = 25;
  real kg = .2;
  real k1 = 65;
  real k2 = 25;
  real Ot = .0033;
  real Tex = .1;
  real eval = .9;
  real fval = 380;
  real Tomax = 100;
  real s = 0 M;
  real Enz(time) M;
  real Ts(time) M;
  real v1(time) katal;
  real v2(time) katal;
  real v3(time) katal;
  real v4(time) katal;
  real v5(time) katal;

  // equations
  To = (-1)/(2*Tex)*(eval*fval+fval*Tex+Tex*Ts-(4*fval*Tex^2*Tomax+((-1)*eval*fval-fval*Tex-Tex*Ts)^2)^(1/2));
  Tt = To+Ts;
  when (time=time.min) Enz = 0;
  Enz:time = (v1 + -1*v2)/compartment;
  when (time=time.min) Ts = 0;
  Ts:time = (v3 + -1*v4 + -1*v5)/compartment;
  v1 = k1*Ot*ki1^nH/(ki1^nH+Tt^nH);
  v2 = mu*Enz;
  v3 = k2*Enz*ki2/(ki2+Tt);
  v4 = gval*Ts/(kg+Ts);
  v5 = mu*Ts;

  // variable properties
  compartment.sbmlRole="compartment";
  To.sbmlRole="parameter";
  Tt.sbmlRole="parameter";
  ki2.sbmlRole="parameter";
  ki1.sbmlRole="parameter";
  nH.sbmlRole="parameter";
  mu.sbmlRole="parameter";
  gval.sbmlRole="parameter";
  kg.sbmlRole="parameter";
  k1.sbmlRole="parameter";
  k2.sbmlRole="parameter";
  Ot.sbmlRole="parameter";
  Tex.sbmlRole="parameter";
  eval.sbmlRole="parameter";
  fval.sbmlRole="parameter";
  Tomax.sbmlRole="parameter";
  s.sbmlRole="species";
  s.sbmlCompartment="compartment";
  Enz.sbmlRole="species";
  Enz.sbmlCompartment="compartment";
  Ts.sbmlRole="species";
  Ts.sbmlCompartment="compartment";
  v1.sbmlRole="rate";
  v2.sbmlRole="rate";
  v3.sbmlRole="rate";
  v4.sbmlRole="rate";
  v5.sbmlRole="rate";
}