// 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
// R01361: space_R_space_minus_3_minus_Hydroxybutanoate => Acetoacetate 
// R00238: Acetyl_minus_CoA => Acetoacetyl_minus_CoA 
// R01360: space_S_space_minus_3_minus_Hydroxy_minus_3_minus_methylglutaryl_minus_CoA <=> Acetyl_minus_CoA Acetoacetate 

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 E4_space_1_space_3_space_4 = 0 M;
  real E2_space_3_space_1_space_9 = 0 M;
  real E1_space_1_space_1_space_30 = 0 M;
  real Acetoacetyl_minus_CoA(time) M;
  real space_S_space_minus_3_minus_Hydroxy_minus_3_minus_methylglutaryl_minus_CoA(time) M;
  real Acetoacetate(time) M;
  real space_R_space_minus_3_minus_Hydroxybutanoate(time) M;
  real Acetyl_minus_CoA(time) M;
  real R01361(time) katal;
  extern real R01361.kf;
  extern real R01361.kb;
  real R00238(time) katal;
  extern real R00238.kf;
  extern real R00238.kb;
  real R01360(time) katal;
  extern real R01360.kf;
  extern real R01360.kb;

  // equations
  when (time=time.min) Acetoacetyl_minus_CoA = 0;
  Acetoacetyl_minus_CoA:time = (R00238)/uVol;
  when (time=time.min) space_S_space_minus_3_minus_Hydroxy_minus_3_minus_methylglutaryl_minus_CoA = 0;
  space_S_space_minus_3_minus_Hydroxy_minus_3_minus_methylglutaryl_minus_CoA:time = (-1*R01360)/uVol;
  when (time=time.min) Acetoacetate = 0;
  Acetoacetate:time = (R01361 + R01360)/uVol;
  when (time=time.min) space_R_space_minus_3_minus_Hydroxybutanoate = 0;
  space_R_space_minus_3_minus_Hydroxybutanoate:time = (-1*R01361)/uVol;
  when (time=time.min) Acetyl_minus_CoA = 0;
  Acetyl_minus_CoA:time = (-1*R00238 + R01360)/uVol;
  R01361 = R01361.kf*space_R_space_minus_3_minus_Hydroxybutanoate - R01361.kb*Acetoacetate;
  R00238 = R00238.kf*Acetyl_minus_CoA - R00238.kb*Acetoacetyl_minus_CoA;
  R01360 = R01360.kf*space_S_space_minus_3_minus_Hydroxy_minus_3_minus_methylglutaryl_minus_CoA - R01360.kb*Acetyl_minus_CoA*Acetoacetate;

  // variable properties
  default.sbmlRole="compartment";
  uVol.sbmlRole="compartment";
  E4_space_1_space_3_space_4.sbmlRole="species";
  E4_space_1_space_3_space_4.sbmlCompartment="uVol";
  E2_space_3_space_1_space_9.sbmlRole="species";
  E2_space_3_space_1_space_9.sbmlCompartment="uVol";
  E1_space_1_space_1_space_30.sbmlRole="species";
  E1_space_1_space_1_space_30.sbmlCompartment="uVol";
  Acetoacetyl_minus_CoA.sbmlRole="species";
  Acetoacetyl_minus_CoA.sbmlCompartment="uVol";
  space_S_space_minus_3_minus_Hydroxy_minus_3_minus_methylglutaryl_minus_CoA.sbmlRole="species";
  space_S_space_minus_3_minus_Hydroxy_minus_3_minus_methylglutaryl_minus_CoA.sbmlCompartment="uVol";
  Acetoacetate.sbmlRole="species";
  Acetoacetate.sbmlCompartment="uVol";
  space_R_space_minus_3_minus_Hydroxybutanoate.sbmlRole="species";
  space_R_space_minus_3_minus_Hydroxybutanoate.sbmlCompartment="uVol";
  Acetyl_minus_CoA.sbmlRole="species";
  Acetyl_minus_CoA.sbmlCompartment="uVol";
  R01361.sbmlRole="rate";
  R01361.kf.sbmlRole="rate";
  R01361.kb.sbmlRole="rate";
  R00238.sbmlRole="rate";
  R00238.kf.sbmlRole="rate";
  R00238.kb.sbmlRole="rate";
  R01360.sbmlRole="rate";
  R01360.kf.sbmlRole="rate";
  R01360.kb.sbmlRole="rate";
}