// 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
// R00286: UDP_minus_glucose <=> UDP_minus_glucuronate 
// R02957: D_minus_Glucuronolactone => D_minus_Glucarate 
// R02933: L_minus_Gulono_minus_1_space_4_minus_lactone => L_minus_Gulonate 

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

  // variable definitions
  real default = 1 L;
  real uVol = 1 L;
  real E3_space_1_space_1_space_17 = 0 M;
  real E1_space_1_space_1_space_22 = 0 M;
  real E1_space_2_space_1_space_3 = 0 M;
  real UDP_minus_glucuronate(time) M;
  real L_minus_Gulonate(time) M;
  real L_minus_Gulono_minus_1_space_4_minus_lactone(time) M;
  real UDP_minus_glucose(time) M;
  real D_minus_Glucuronolactone(time) M;
  real D_minus_Glucarate(time) M;
  real R00286(time) katal;
  extern real R00286.kf;
  extern real R00286.kb;
  real R02957(time) katal;
  extern real R02957.kf;
  extern real R02957.kb;
  real R02933(time) katal;
  extern real R02933.kf;
  extern real R02933.kb;

  // equations
  when (time=time.min) UDP_minus_glucuronate = 0;
  UDP_minus_glucuronate:time = (R00286)/uVol;
  when (time=time.min) L_minus_Gulonate = 0;
  L_minus_Gulonate:time = (R02933)/uVol;
  when (time=time.min) L_minus_Gulono_minus_1_space_4_minus_lactone = 0;
  L_minus_Gulono_minus_1_space_4_minus_lactone:time = (-1*R02933)/uVol;
  when (time=time.min) UDP_minus_glucose = 0;
  UDP_minus_glucose:time = (-1*R00286)/uVol;
  when (time=time.min) D_minus_Glucuronolactone = 0;
  D_minus_Glucuronolactone:time = (-1*R02957)/uVol;
  when (time=time.min) D_minus_Glucarate = 0;
  D_minus_Glucarate:time = (R02957)/uVol;
  R00286 = R00286.kf*UDP_minus_glucose - R00286.kb*UDP_minus_glucuronate;
  R02957 = R02957.kf*D_minus_Glucuronolactone - R02957.kb*D_minus_Glucarate;
  R02933 = R02933.kf*L_minus_Gulono_minus_1_space_4_minus_lactone - R02933.kb*L_minus_Gulonate;

  // variable properties
  default.sbmlRole="compartment";
  uVol.sbmlRole="compartment";
  E3_space_1_space_1_space_17.sbmlRole="species";
  E3_space_1_space_1_space_17.sbmlCompartment="uVol";
  E1_space_1_space_1_space_22.sbmlRole="species";
  E1_space_1_space_1_space_22.sbmlCompartment="uVol";
  E1_space_2_space_1_space_3.sbmlRole="species";
  E1_space_2_space_1_space_3.sbmlCompartment="uVol";
  UDP_minus_glucuronate.sbmlRole="species";
  UDP_minus_glucuronate.sbmlCompartment="uVol";
  L_minus_Gulonate.sbmlRole="species";
  L_minus_Gulonate.sbmlCompartment="uVol";
  L_minus_Gulono_minus_1_space_4_minus_lactone.sbmlRole="species";
  L_minus_Gulono_minus_1_space_4_minus_lactone.sbmlCompartment="uVol";
  UDP_minus_glucose.sbmlRole="species";
  UDP_minus_glucose.sbmlCompartment="uVol";
  D_minus_Glucuronolactone.sbmlRole="species";
  D_minus_Glucuronolactone.sbmlCompartment="uVol";
  D_minus_Glucarate.sbmlRole="species";
  D_minus_Glucarate.sbmlCompartment="uVol";
  R00286.sbmlRole="rate";
  R00286.kf.sbmlRole="rate";
  R00286.kb.sbmlRole="rate";
  R02957.sbmlRole="rate";
  R02957.kf.sbmlRole="rate";
  R02957.kb.sbmlRole="rate";
  R02933.sbmlRole="rate";
  R02933.kf.sbmlRole="rate";
  R02933.kb.sbmlRole="rate";
}