// 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
// v1: NADH O2 => H2O2 NAD 
// v2: per3 H2O2 => coI 
// v3: coI ArH => coII Ar 
// v4: coII ArH => per3 Ar 
// v5: O2 NADrad => super NAD 
// v6: per3 super => coIII 
// v7: 2super => O2 H2O2 
// v8: NADrad coIII => coI NAD 
// v9: 2NADrad => NAD2 
// v10: per3 NADrad => per2 NAD 
// v11: O2 per2 => coIII 
// v12: NADHres => NADH 
// v131: O2g => O2 
// v132: O2 => O2g 
// v14: NADH Ar => NADrad ArH 

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

  // variable definitions
  real compartment = 1 L;
  real k1 = 3E-6;
  real k2 = 18;
  real k3 = .15;
  real k4 = .0052;
  real k5 = 20;
  real k6 = 17;
  real k7 = 20;
  real k8 = 40;
  real k9 = 60;
  real k10 = 1.8;
  real k11 = .1;
  real k12 = .08;
  real k13f = .072;
  real k13b = .006;
  real k14 = .7;
  real NADH(time) M;
  real O2(time) M;
  real H2O2(time) M;
  real NAD = 0 M;
  real per3(time) M;
  real coI(time) M;
  real ArH(time) M;
  real coII(time) M;
  real Ar(time) M;
  real NADrad(time) M;
  real super(time) M;
  real coIII(time) M;
  real NAD2 = 0 M;
  real per2(time) M;
  real NADHres = 1 M;
  real O2g = 1 M;
  real v1(time) katal;
  real v2(time) katal;
  real v3(time) katal;
  real v4(time) katal;
  real v5(time) katal;
  real v6(time) katal;
  real v7(time) katal;
  real v8(time) katal;
  real v9(time) katal;
  real v10(time) katal;
  real v11(time) katal;
  real v12(time) katal;
  real v131(time) katal;
  real v132(time) katal;
  real v14(time) katal;

  // equations
  when (time=time.min) NADH = 0;
  NADH:time = (-1*v1 + v12 + -1*v14)/compartment;
  when (time=time.min) O2 = 12;
  O2:time = (-1*v1 + -1*v5 + v7 + -1*v11 + v131 + -1*v132)/compartment;
  when (time=time.min) H2O2 = .7;
  H2O2:time = (v1 + -1*v2 + v7)/compartment;
  when (time=time.min) per3 = 1.2;
  per3:time = (-1*v2 + v4 + -1*v6 + -1*v10)/compartment;
  when (time=time.min) coI = 0;
  coI:time = (v2 + -1*v3 + v8)/compartment;
  when (time=time.min) ArH = 300;
  ArH:time = (-1*v3 + -1*v4 + v14)/compartment;
  when (time=time.min) coII = 0;
  coII:time = (v3 + -1*v4)/compartment;
  when (time=time.min) Ar = 0;
  Ar:time = (v3 + v4 + -1*v14)/compartment;
  when (time=time.min) NADrad = 0;
  NADrad:time = (-1*v5 + -1*v8 + -2*v9 + -1*v10 + v14)/compartment;
  when (time=time.min) super = 0;
  super:time = (v5 + -1*v6 + -2*v7)/compartment;
  when (time=time.min) coIII = 0;
  coIII:time = (v6 + -1*v8 + v11)/compartment;
  when (time=time.min) per2 = 0;
  per2:time = (v10 + -1*v11)/compartment;
  v1 = k1*NADH*O2;
  v2 = k2*H2O2*per3;
  v3 = k3*coI*ArH;
  v4 = k4*coII*ArH;
  v5 = k5*NADrad*O2;
  v6 = k6*super*per3;
  v7 = k7*super*super;
  v8 = k8*coIII*NADrad;
  v9 = k9*NADrad*NADrad;
  v10 = k10*per3*NADrad;
  v11 = k11*per2*O2;
  v12 = k12;
  v131 = k13f;
  v132 = k13b*O2;
  v14 = k14*Ar*NADH;

  // variable properties
  compartment.sbmlRole="compartment";
  k1.sbmlRole="parameter";
  k2.sbmlRole="parameter";
  k3.sbmlRole="parameter";
  k4.sbmlRole="parameter";
  k5.sbmlRole="parameter";
  k6.sbmlRole="parameter";
  k7.sbmlRole="parameter";
  k8.sbmlRole="parameter";
  k9.sbmlRole="parameter";
  k10.sbmlRole="parameter";
  k11.sbmlRole="parameter";
  k12.sbmlRole="parameter";
  k13f.sbmlRole="parameter";
  k13b.sbmlRole="parameter";
  k14.sbmlRole="parameter";
  NADH.sbmlRole="species";
  NADH.sbmlCompartment="compartment";
  O2.sbmlRole="species";
  O2.sbmlCompartment="compartment";
  H2O2.sbmlRole="species";
  H2O2.sbmlCompartment="compartment";
  NAD.sbmlRole="species";
  NAD.sbmlCompartment="compartment";
  per3.sbmlRole="species";
  per3.sbmlCompartment="compartment";
  coI.sbmlRole="species";
  coI.sbmlCompartment="compartment";
  ArH.sbmlRole="species";
  ArH.sbmlCompartment="compartment";
  coII.sbmlRole="species";
  coII.sbmlCompartment="compartment";
  Ar.sbmlRole="species";
  Ar.sbmlCompartment="compartment";
  NADrad.sbmlRole="species";
  NADrad.sbmlCompartment="compartment";
  super.sbmlRole="species";
  super.sbmlCompartment="compartment";
  coIII.sbmlRole="species";
  coIII.sbmlCompartment="compartment";
  NAD2.sbmlRole="species";
  NAD2.sbmlCompartment="compartment";
  per2.sbmlRole="species";
  per2.sbmlCompartment="compartment";
  NADHres.sbmlRole="species";
  NADHres.sbmlCompartment="compartment";
  O2g.sbmlRole="species";
  O2g.sbmlCompartment="compartment";
  v1.sbmlRole="rate";
  v2.sbmlRole="rate";
  v3.sbmlRole="rate";
  v4.sbmlRole="rate";
  v5.sbmlRole="rate";
  v6.sbmlRole="rate";
  v7.sbmlRole="rate";
  v8.sbmlRole="rate";
  v9.sbmlRole="rate";
  v10.sbmlRole="rate";
  v11.sbmlRole="rate";
  v12.sbmlRole="rate";
  v131.sbmlRole="rate";
  v132.sbmlRole="rate";
  v14.sbmlRole="rate";
}