// This model generated automatically from SBML

// unit definitions
import nsrunit;
unit conversion off;
unit item=scalar;
unit substance = mole;
unit volume = litre;
unit area = metre^2;
unit length = metre;
unit time = second;

// SBML property definitions
property sbmlRole=string;
property sbmlName=string;
property sbmlCompartment=string;

// SBML reactions
// React0: B L => BL 
// React1: BL L => BLL 
// React2: BLL => ALL 
// React3: A L => AL 
// React4: AL L => ALL 
// React5: B => A 
// React6: BL => AL 
// React7: I L => IL 
// React8: IL L => ILL 
// React9: A => I 
// React10: AL => IL 
// React11: ALL => ILL 
// React12: D L => DL 
// React13: DL L => DLL 
// React14: I => D 
// React15: IL => DL 
// React16: ILL => DLL 

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

  // variable definitions
  real comp1 = 1E-16 volume;
  real kf_0 = 3E8;
  real kr_0 = 8E3;
  real kf_1 = 1.5E8;
  real kr_1 = 16000;
  real kf_2 = 3E4;
  real kr_2 = 700;
  real kf_3 = 3E8;
  real kr_3 = 8.64;
  real kf_4 = 1.5E8;
  real kr_4 = 17.28;
  real kf_5 = .54;
  real kr_5 = 10800;
  real kf_6 = 130;
  real kr_6 = 2740;
  real kf_7 = 3E8;
  real kr_7 = 4;
  real kf_8 = 1.5E8;
  real kr_8 = 8;
  real kf_9 = 19.7;
  real kr_9 = 3.74;
  real kf_10 = 19.85;
  real kr_10 = 1.74;
  real kf_11 = 20;
  real kr_11 = .81;
  real kf_12 = 3E8;
  real kr_12 = 4;
  real kf_13 = 1.5E8;
  real kr_13 = 8;
  real kf_14 = .05;
  real kr_14 = .0012;
  real kf_15 = .05;
  real kr_15 = .0012;
  real kf_16 = .05;
  real kr_16 = .0012;
  private real BLL.amt(time) substance;
  real BLL(time) substance/volume;
  real BLL.init substance;
  private real IL.amt(time) substance;
  real IL(time) substance/volume;
  real IL.init substance;
  private real AL.amt(time) substance;
  real AL(time) substance/volume;
  real AL.init substance;
  private real A.amt(time) substance;
  real A(time) substance/volume;
  real A.init substance;
  private real BL.amt(time) substance;
  real BL(time) substance/volume;
  real BL.init substance;
  private real B.amt(time) substance;
  real B(time) substance/volume;
  real B.init substance;
  private real DLL.amt(time) substance;
  real DLL(time) substance/volume;
  real DLL.init substance;
  private real D.amt(time) substance;
  real D(time) substance/volume;
  real D.init substance;
  private real ILL.amt(time) substance;
  real ILL(time) substance/volume;
  real ILL.init substance;
  private real DL.amt(time) substance;
  real DL(time) substance/volume;
  real DL.init substance;
  private real I.amt(time) substance;
  real I(time) substance/volume;
  real I.init substance;
  private real ALL.amt(time) substance;
  real ALL(time) substance/volume;
  real ALL.init substance;
  private real L.amt(time) substance;
  real L(time) substance/volume;
  real L.init substance;
  real React0.rate(time) substance/time;
  real React1.rate(time) substance/time;
  real React2.rate(time) substance/time;
  real React3.rate(time) substance/time;
  real React4.rate(time) substance/time;
  real React5.rate(time) substance/time;
  real React6.rate(time) substance/time;
  real React7.rate(time) substance/time;
  real React8.rate(time) substance/time;
  real React9.rate(time) substance/time;
  real React10.rate(time) substance/time;
  real React11.rate(time) substance/time;
  real React12.rate(time) substance/time;
  real React13.rate(time) substance/time;
  real React14.rate(time) substance/time;
  real React15.rate(time) substance/time;
  real React16.rate(time) substance/time;

  // equations
  when (time=time.min) BLL.amt = BLL.init;
  BLL.amt:time = React1.rate + -1*React2.rate;
  BLL = BLL.amt/comp1;
  BLL.init = 0;
  when (time=time.min) IL.amt = IL.init;
  IL.amt:time = React7.rate + -1*React8.rate + React10.rate + -1*React15.rate;
  IL = IL.amt/comp1;
  IL.init = 0;
  when (time=time.min) AL.amt = AL.init;
  AL.amt:time = React3.rate + -1*React4.rate + React6.rate + -1*React10.rate;
  AL = AL.amt/comp1;
  AL.init = 0;
  when (time=time.min) A.amt = A.init;
  A.amt:time = -1*React3.rate + React5.rate + -1*React9.rate;
  A = A.amt/comp1;
  A.init = 0;
  when (time=time.min) BL.amt = BL.init;
  BL.amt:time = React0.rate + -1*React1.rate + -1*React6.rate;
  BL = BL.amt/comp1;
  BL.init = 0;
  when (time=time.min) B.amt = B.init;
  B.amt:time = -1*React0.rate + -1*React5.rate;
  B = B.amt/comp1;
  B.init = 1.6605779E-21;
  when (time=time.min) DLL.amt = DLL.init;
  DLL.amt:time = React13.rate + React16.rate;
  DLL = DLL.amt/comp1;
  DLL.init = 0;
  when (time=time.min) D.amt = D.init;
  D.amt:time = -1*React12.rate + React14.rate;
  D = D.amt/comp1;
  D.init = 0;
  when (time=time.min) ILL.amt = ILL.init;
  ILL.amt:time = React8.rate + React11.rate + -1*React16.rate;
  ILL = ILL.amt/comp1;
  ILL.init = 0;
  when (time=time.min) DL.amt = DL.init;
  DL.amt:time = React12.rate + -1*React13.rate + React15.rate;
  DL = DL.amt/comp1;
  DL.init = 0;
  when (time=time.min) I.amt = I.init;
  I.amt:time = -1*React7.rate + React9.rate + -1*React14.rate;
  I = I.amt/comp1;
  I.init = 0;
  when (time=time.min) ALL.amt = ALL.init;
  ALL.amt:time = React2.rate + React4.rate + -1*React11.rate;
  ALL = ALL.amt/comp1;
  ALL.init = 0;
  when (time=time.min) L.amt = L.init;
  L.amt:time = -1*React0.rate + -1*React1.rate + -1*React3.rate + -1*React4.rate + -1*React7.rate + -1*React8.rate + -1*React12.rate + -1*React13.rate;
  L = L.amt/comp1;
  L.init = 1.6605779E-20;
  React0.rate = comp1*(kf_0*B*L-kr_0*BL);
  React1.rate = comp1*(kf_1*BL*L-kr_1*BLL);
  React2.rate = comp1*(kf_2*BLL-kr_2*ALL);
  React3.rate = comp1*(kf_3*A*L-kr_3*AL);
  React4.rate = comp1*(kf_4*AL*L-kr_4*ALL);
  React5.rate = comp1*(kf_5*B-kr_5*A);
  React6.rate = comp1*(kf_6*BL-kr_6*AL);
  React7.rate = comp1*(kf_7*I*L-kr_7*IL);
  React8.rate = comp1*(kf_8*IL*L-kr_8*ILL);
  React9.rate = comp1*(kf_9*A-kr_9*I);
  React10.rate = comp1*(kf_10*AL-kr_10*IL);
  React11.rate = comp1*(kf_11*ALL-kr_11*ILL);
  React12.rate = comp1*(kf_12*D*L-kr_12*DL);
  React13.rate = comp1*(kf_13*DL*L-kr_13*DLL);
  React14.rate = comp1*(kf_14*I-kr_14*D);
  React15.rate = comp1*(kf_15*IL-kr_15*DL);
  React16.rate = comp1*(kf_16*ILL-kr_16*DLL);

  // variable properties
  comp1.sbmlRole="compartment";
  comp1.sbmlName="compartment1";
  kf_0.sbmlRole="parameter";
  kr_0.sbmlRole="parameter";
  kf_1.sbmlRole="parameter";
  kr_1.sbmlRole="parameter";
  kf_2.sbmlRole="parameter";
  kr_2.sbmlRole="parameter";
  kf_3.sbmlRole="parameter";
  kr_3.sbmlRole="parameter";
  kf_4.sbmlRole="parameter";
  kr_4.sbmlRole="parameter";
  kf_5.sbmlRole="parameter";
  kr_5.sbmlRole="parameter";
  kf_6.sbmlRole="parameter";
  kr_6.sbmlRole="parameter";
  kf_7.sbmlRole="parameter";
  kr_7.sbmlRole="parameter";
  kf_8.sbmlRole="parameter";
  kr_8.sbmlRole="parameter";
  kf_9.sbmlRole="parameter";
  kr_9.sbmlRole="parameter";
  kf_10.sbmlRole="parameter";
  kr_10.sbmlRole="parameter";
  kf_11.sbmlRole="parameter";
  kr_11.sbmlRole="parameter";
  kf_12.sbmlRole="parameter";
  kr_12.sbmlRole="parameter";
  kf_13.sbmlRole="parameter";
  kr_13.sbmlRole="parameter";
  kf_14.sbmlRole="parameter";
  kr_14.sbmlRole="parameter";
  kf_15.sbmlRole="parameter";
  kr_15.sbmlRole="parameter";
  kf_16.sbmlRole="parameter";
  kr_16.sbmlRole="parameter";
  BLL.amt.sbmlRole="speciesAmount";
  BLL.sbmlRole="speciesConcentration";
  BLL.sbmlName="BasalACh2";
  BLL.sbmlCompartment="comp1";
  BLL.init.sbmlRole="speciesInitialAmount";
  IL.amt.sbmlRole="speciesAmount";
  IL.sbmlRole="speciesConcentration";
  IL.sbmlName="IntermediateACh";
  IL.sbmlCompartment="comp1";
  IL.init.sbmlRole="speciesInitialAmount";
  AL.amt.sbmlRole="speciesAmount";
  AL.sbmlRole="speciesConcentration";
  AL.sbmlName="ActiveACh";
  AL.sbmlCompartment="comp1";
  AL.init.sbmlRole="speciesInitialAmount";
  A.amt.sbmlRole="speciesAmount";
  A.sbmlRole="speciesConcentration";
  A.sbmlName="Active";
  A.sbmlCompartment="comp1";
  A.init.sbmlRole="speciesInitialAmount";
  BL.amt.sbmlRole="speciesAmount";
  BL.sbmlRole="speciesConcentration";
  BL.sbmlName="BasalACh";
  BL.sbmlCompartment="comp1";
  BL.init.sbmlRole="speciesInitialAmount";
  B.amt.sbmlRole="speciesAmount";
  B.sbmlRole="speciesConcentration";
  B.sbmlName="Basal";
  B.sbmlCompartment="comp1";
  B.init.sbmlRole="speciesInitialAmount";
  DLL.amt.sbmlRole="speciesAmount";
  DLL.sbmlRole="speciesConcentration";
  DLL.sbmlName="DesensitisedACh2";
  DLL.sbmlCompartment="comp1";
  DLL.init.sbmlRole="speciesInitialAmount";
  D.amt.sbmlRole="speciesAmount";
  D.sbmlRole="speciesConcentration";
  D.sbmlName="Desensitised";
  D.sbmlCompartment="comp1";
  D.init.sbmlRole="speciesInitialAmount";
  ILL.amt.sbmlRole="speciesAmount";
  ILL.sbmlRole="speciesConcentration";
  ILL.sbmlName="IntermediateACh2";
  ILL.sbmlCompartment="comp1";
  ILL.init.sbmlRole="speciesInitialAmount";
  DL.amt.sbmlRole="speciesAmount";
  DL.sbmlRole="speciesConcentration";
  DL.sbmlName="DesensitisedACh";
  DL.sbmlCompartment="comp1";
  DL.init.sbmlRole="speciesInitialAmount";
  I.amt.sbmlRole="speciesAmount";
  I.sbmlRole="speciesConcentration";
  I.sbmlName="Intermediate";
  I.sbmlCompartment="comp1";
  I.init.sbmlRole="speciesInitialAmount";
  ALL.amt.sbmlRole="speciesAmount";
  ALL.sbmlRole="speciesConcentration";
  ALL.sbmlName="ActiveACh2";
  ALL.sbmlCompartment="comp1";
  ALL.init.sbmlRole="speciesInitialAmount";
  L.amt.sbmlRole="speciesAmount";
  L.sbmlRole="speciesConcentration";
  L.sbmlName="acetylcholine";
  L.sbmlCompartment="comp1";
  L.init.sbmlRole="speciesInitialAmount";
  React0.rate.sbmlRole="rate";
  React1.rate.sbmlRole="rate";
  React2.rate.sbmlRole="rate";
  React3.rate.sbmlRole="rate";
  React4.rate.sbmlRole="rate";
  React5.rate.sbmlRole="rate";
  React6.rate.sbmlRole="rate";
  React7.rate.sbmlRole="rate";
  React8.rate.sbmlRole="rate";
  React9.rate.sbmlRole="rate";
  React10.rate.sbmlRole="rate";
  React11.rate.sbmlRole="rate";
  React12.rate.sbmlRole="rate";
  React13.rate.sbmlRole="rate";
  React14.rate.sbmlRole="rate";
  React15.rate.sbmlRole="rate";
  React16.rate.sbmlRole="rate";
}