// 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
// Reaction1: Br BrO3 <=> HBrO2 HOBr 
// Reaction2: Br HBrO2 <=> 2HOBr 
// Reaction3: BrO3 HBrO2 <=> Ce 2HBrO2 
// Reaction4: 2HBrO2 <=> BrO3 HOBr 
// Reaction5: Ce <=> (f)Br 

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

  // variable definitions
  real BZ = 1 L;
  real f = 1;
  real Br(time) M;
  real BrO3 = .06 M;
  real Ce(time) M;
  real HBrO2(time) M;
  real HOBr = 0 M;
  real Reaction1(time) katal;
  real k1 = 1.34;
  real Reaction2(time) katal;
  real k2 = 1.6E9;
  real Reaction3(time) katal;
  real k3 = 8E3;
  real Reaction4(time) katal;
  real k4 = 4E7;
  real Reaction5(time) katal;
  real k5 = 1;

  // equations
  when (time=time.min) Br = 1E-7;
  Br:time = (-1*Reaction1 + -1*Reaction2 + (f)*Reaction5)/BZ;
  when (time=time.min) Ce = .05;
  Ce:time = (Reaction3 + -1*Reaction5)/BZ;
  when (time=time.min) HBrO2 = 5E-11;
  HBrO2:time = (Reaction1 + -1*Reaction2 + -1*Reaction3 + 2*Reaction3 + -2*Reaction4)/BZ;
  Reaction1 = Br*BrO3*k1*BZ;
  Reaction2 = Br*HBrO2*k2*BZ;
  Reaction3 = BrO3*HBrO2*k3*BZ;
  Reaction4 = HBrO2^2*k4*BZ;
  Reaction5 = Ce*k5*BZ;

  // variable properties
  BZ.sbmlRole="compartment";
  f.sbmlRole="parameter";
  Br.sbmlRole="species";
  Br.sbmlCompartment="BZ";
  BrO3.sbmlRole="species";
  BrO3.sbmlCompartment="BZ";
  Ce.sbmlRole="species";
  Ce.sbmlCompartment="BZ";
  HBrO2.sbmlRole="species";
  HBrO2.sbmlCompartment="BZ";
  HOBr.sbmlRole="species";
  HOBr.sbmlCompartment="BZ";
  Reaction1.sbmlRole="rate";
  k1.sbmlRole="parameter";
  k1.sbmlName="kM1";
  Reaction2.sbmlRole="rate";
  k2.sbmlRole="parameter";
  k2.sbmlName="kM2";
  Reaction3.sbmlRole="rate";
  k3.sbmlRole="parameter";
  k3.sbmlName="kM3";
  Reaction4.sbmlRole="rate";
  k4.sbmlRole="parameter";
  k4.sbmlName="kM4";
  Reaction5.sbmlRole="rate";
  k5.sbmlRole="parameter";
  k5.sbmlName="kM5";
}