// 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: c3 => c3a 
// v2: c8 => c8a 
// v3: iap c3a => c3aiap 
// v4: iap => breakdown 
// v5: c8a => breakdown 
// v6: c3a => breakdown 
// v7: c3aiap => breakdown 
// v8: iap => breakdown 
// v9: c8 => breakdown 
// v10: c3 => breakdown 
// v11: c8a bar => c8abar 
// v12: bar => breakdown 
// v13: c8abar => breakdown 

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

  // variable definitions
  real compartment = 1 L;
  real external(time);
  real k1 = 5.8E-5;
  real k2 = 1E-5;
  real k3f = 5E-4;
  real k3r = .21;
  real k4 = 3E-4;
  real k5 = .0058;
  real k6 = .0058;
  real k7 = .0173;
  real k8f = .0116;
  real k8r = 464;
  real k9f = .0039;
  real k9r = 507;
  real k10f = .0039;
  real k10r = 81.9;
  real k11f = 5E-4;
  real k11r = .21;
  real k12f = .001;
  real k12r = 40;
  real k13 = .0116;
  real c3(time) M;
  real c3a(time) M;
  real c8(time) M;
  real c8a(time) M;
  real iap(time) M;
  real c3aiap(time) M;
  real breakdown = 1 M;
  real bar(time) M;
  real c8abar(time) 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 v13(time) katal;

  // equations
  external = breakdown;
  when (time=time.min) c3 = 21000;
  c3:time = (-1*v1 + -1*v10)/compartment;
  when (time=time.min) c3a = 0;
  c3a:time = (v1 + -1*v3 + -1*v6)/compartment;
  when (time=time.min) c8 = 1.3E5;
  c8:time = (-1*v2 + -1*v9)/compartment;
  when (time=time.min) c8a = 5E3;
  c8a:time = (v2 + -1*v5 + -1*v11)/compartment;
  when (time=time.min) iap = 4E4;
  iap:time = (-1*v3 + -1*v4 + -1*v8)/compartment;
  when (time=time.min) c3aiap = 0;
  c3aiap:time = (v3 + -1*v7)/compartment;
  when (time=time.min) bar = 4E4;
  bar:time = (-1*v11 + -1*v12)/compartment;
  when (time=time.min) c8abar = 0;
  c8abar:time = (v11 + -1*v13)/compartment;
  v1 = k1*c8a*c3;
  v2 = k2*c3a*c8;
  v3 = k3f*c3a*iap-k3r*c3aiap;
  v4 = k4*c3a*iap;
  v5 = k5*c8a;
  v6 = k6*c3a;
  v7 = k7*c3aiap;
  v8 = k8f*iap-k8r;
  v9 = k9f*c8-k9r;
  v10 = k10f*c3-k10r;
  v11 = k11f*c8a*bar-k11r*c8abar;
  v12 = k12f*bar-k12r;
  v13 = k13*c8abar;

  // variable properties
  compartment.sbmlRole="compartment";
  external.sbmlRole="parameter";
  k1.sbmlRole="parameter";
  k2.sbmlRole="parameter";
  k3f.sbmlRole="parameter";
  k3r.sbmlRole="parameter";
  k4.sbmlRole="parameter";
  k5.sbmlRole="parameter";
  k6.sbmlRole="parameter";
  k7.sbmlRole="parameter";
  k8f.sbmlRole="parameter";
  k8r.sbmlRole="parameter";
  k9f.sbmlRole="parameter";
  k9r.sbmlRole="parameter";
  k10f.sbmlRole="parameter";
  k10r.sbmlRole="parameter";
  k11f.sbmlRole="parameter";
  k11r.sbmlRole="parameter";
  k12f.sbmlRole="parameter";
  k12r.sbmlRole="parameter";
  k13.sbmlRole="parameter";
  c3.sbmlRole="species";
  c3.sbmlCompartment="compartment";
  c3a.sbmlRole="species";
  c3a.sbmlCompartment="compartment";
  c8.sbmlRole="species";
  c8.sbmlCompartment="compartment";
  c8a.sbmlRole="species";
  c8a.sbmlCompartment="compartment";
  iap.sbmlRole="species";
  iap.sbmlCompartment="compartment";
  c3aiap.sbmlRole="species";
  c3aiap.sbmlCompartment="compartment";
  breakdown.sbmlRole="species";
  breakdown.sbmlCompartment="compartment";
  bar.sbmlRole="species";
  bar.sbmlCompartment="compartment";
  c8abar.sbmlRole="species";
  c8abar.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";
  v13.sbmlRole="rate";
}