// 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
// vgk: glc => F6P 
// vpfk: F6P => FBP 
// vfba: FBP => g3P 

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

  // variable definitions
  real compartment = 1 L;
  real Vgk = 10;
  real hgk = 1.7;
  real keqpgi = .3;
  real Vpfk = 100;
  real Vfba = 25;
  real ksgk = .008;
  real kspfk = .004;
  real ksfba = 5E-6;
  real kxpfk = 1E-5;
  real alpha = 5;
  real hx = 2.5;
  real hpfk = 2.5;
  real hact = 1;
  real sigcorr = 1.5;
  real fac = 1.8E5;
  real glc = .01 M;
  real F6P(time) M;
  real FBP(time) M;
  real g3P = 0 M;
  real vgk(time) katal;
  real vpfk(time) katal;
  real vfba(time) katal;

  // equations
  when (time=time.min) F6P = .00403814;
  F6P:time = (vgk + -1*vpfk)/compartment;
  when (time=time.min) FBP = 9.1172373E-7;
  FBP:time = (vpfk + -1*vfba)/compartment;
  vgk = Vgk/fac*(glc/ksgk)^hgk/(1+(glc/ksgk)^hgk);
  vpfk = Vpfk/fac*(keqpgi/(1+keqpgi)*F6P/kspfk)^(hpfk-(hpfk-hact)*(FBP/ksfba/(1+FBP/ksfba)))/((keqpgi/(1+keqpgi)*F6P/kspfk)^(hpfk-(hpfk-hact)*(FBP/ksfba/(1+FBP/ksfba)))+(1+(FBP/kxpfk)^hx)/(1+alpha^(hpfk-(hpfk-hact)*(FBP/ksfba/(1+FBP/ksfba)))*(FBP/kxpfk)^hx));
  vfba = Vfba/fac*(FBP/ksfba)/(FBP/ksfba+1);

  // variable properties
  compartment.sbmlRole="compartment";
  Vgk.sbmlRole="parameter";
  hgk.sbmlRole="parameter";
  keqpgi.sbmlRole="parameter";
  Vpfk.sbmlRole="parameter";
  Vfba.sbmlRole="parameter";
  ksgk.sbmlRole="parameter";
  kspfk.sbmlRole="parameter";
  ksfba.sbmlRole="parameter";
  kxpfk.sbmlRole="parameter";
  alpha.sbmlRole="parameter";
  hx.sbmlRole="parameter";
  hpfk.sbmlRole="parameter";
  hact.sbmlRole="parameter";
  sigcorr.sbmlRole="parameter";
  fac.sbmlRole="parameter";
  glc.sbmlRole="species";
  glc.sbmlCompartment="compartment";
  F6P.sbmlRole="species";
  F6P.sbmlCompartment="compartment";
  FBP.sbmlRole="species";
  FBP.sbmlCompartment="compartment";
  g3P.sbmlRole="species";
  g3P.sbmlCompartment="compartment";
  vgk.sbmlRole="rate";
  vpfk.sbmlRole="rate";
  vfba.sbmlRole="rate";
}