// 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
// R00238: Acetyl_minus_CoA => Acetoacetyl_minus_CoA 
// R01978: space_S_space_minus_3_minus_Hydroxy_minus_3_minus_methylglutaryl_minus_CoA <=> Acetyl_minus_CoA Acetoacetyl_minus_CoA 

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 E2_space_3_space_3_space_10 = 0 M;
  real E2_space_3_space_1_space_9 = 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 Acetyl_minus_CoA(time) M;
  real R00238(time) katal;
  extern real R00238.kf;
  extern real R00238.kb;
  real R01978(time) katal;
  extern real R01978.kf;
  extern real R01978.kb;

  // equations
  when (time=time.min) Acetoacetyl_minus_CoA = 0;
  Acetoacetyl_minus_CoA:time = (R00238 + R01978)/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*R01978)/uVol;
  when (time=time.min) Acetyl_minus_CoA = 0;
  Acetyl_minus_CoA:time = (-1*R00238 + R01978)/uVol;
  R00238 = R00238.kf*Acetyl_minus_CoA - R00238.kb*Acetoacetyl_minus_CoA;
  R01978 = R01978.kf*space_S_space_minus_3_minus_Hydroxy_minus_3_minus_methylglutaryl_minus_CoA - R01978.kb*Acetyl_minus_CoA*Acetoacetyl_minus_CoA;

  // variable properties
  default.sbmlRole="compartment";
  uVol.sbmlRole="compartment";
  E2_space_3_space_3_space_10.sbmlRole="species";
  E2_space_3_space_3_space_10.sbmlCompartment="uVol";
  E2_space_3_space_1_space_9.sbmlRole="species";
  E2_space_3_space_1_space_9.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";
  Acetyl_minus_CoA.sbmlRole="species";
  Acetyl_minus_CoA.sbmlCompartment="uVol";
  R00238.sbmlRole="rate";
  R00238.kf.sbmlRole="rate";
  R00238.kb.sbmlRole="rate";
  R01978.sbmlRole="rate";
  R01978.kf.sbmlRole="rate";
  R01978.kb.sbmlRole="rate";
}