//Created by libAntimony v2.4
model sneyd_2002__environment(time_)

  // Variable initializations:
  time_ = ;
end

model sneyd_2002__R(R, phi_1, phi_2, phi_2b, k_1b, l_2b, I_1, O, p, time_)

  // Rate Rules:
  R' = (phi_2b * O + (k_1b + l_2b) * I_1) - (phi_2 * p * R + phi_1 * R);

  // Variable initializations:
  R = 1;
  phi_1 = ;
  phi_2 = ;
  phi_2b = ;
  k_1b = ;
  l_2b = ;
  I_1 = ;
  O = ;
  p = ;
  time_ = ;
end

model sneyd_2002__O(O, phi_3, phi_4, phi_4b, phi_2, phi_2b, k_3b, R, A, S, p, time_)

  // Rate Rules:
  O' = (phi_2 * p * R + phi_4b * A + k_3b * S) - (phi_2b + phi_4 + 1 * phi_3) * O;

  // Variable initializations:
  O = 0;
  phi_3 = ;
  phi_4 = ;
  phi_4b = ;
  phi_2 = ;
  phi_2b = ;
  k_3b = ;
  R = ;
  A = ;
  S = ;
  p = ;
  time_ = ;
end

model sneyd_2002__S(S, I_1, I_2, O, R, A)

  // Assignment Rules:
  S := 1 - (R + O + A + I_1 + I_2);

  // Variable initializations:
  I_1 = ;
  I_2 = ;
  O = ;
  R = ;
  A = ;
end

model sneyd_2002__I_1(I_1, phi_1, k_1b, l_2b, R, time_)

  // Rate Rules:
  I_1' = phi_1 * R - (k_1b + l_2b) * I_1;

  // Variable initializations:
  I_1 = 0;
  phi_1 = ;
  k_1b = ;
  l_2b = ;
  R = ;
  time_ = ;
end

model sneyd_2002__I_2(I_2, phi_5, k_1b, l_2b, A, time_)

  // Rate Rules:
  I_2' = phi_5 * A - (k_1b + l_2b) * I_2;

  // Variable initializations:
  I_2 = 0;
  phi_5 = ;
  k_1b = ;
  l_2b = ;
  A = ;
  time_ = ;
end

model sneyd_2002__A(A, phi_5, phi_4, phi_4b, k_1b, l_2b, I_2, O, time_)

  // Rate Rules:
  A' = (phi_4 * O + (k_1b + l_2b) * I_2) - (phi_4b * A + phi_5 * A);

  // Variable initializations:
  A = 0;
  phi_5 = ;
  phi_4 = ;
  phi_4b = ;
  k_1b = ;
  l_2b = ;
  I_2 = ;
  O = ;
  time_ = ;
end

model sneyd_2002__open_probability(O, A)

  // Assignment Rules:
  open_probability := power(0.1 * O + 0.9 * A, 4);

  // Variable initializations:
  O = ;
  A = ;
end

model sneyd_2002__reaction_rate_constants(phi_1, phi_2, phi_2b, phi_3, phi_4, phi_4b, phi_5, k_1b, k_3b, l_2b, p)

  // Assignment Rules:
  phi_1 := ((k_1a * L_1 + l_2a) * c) / (L_1 + c * (1 + L_1 / L_3));
  phi_2 := (k_2a * L_3 + l_4a * c) / (L_3 + c * (1 + L_3 / L_1));
  phi_2b := (k_2b + l_4b * c) / (1 + c / L_5);
  phi_3 := (k_3a * L_5) / (c + L_5);
  phi_4 := ((k_4a * L_5 + l_6a) * c) / (c + L_5);
  phi_4b := (L_1 * (k_4b + l_6b)) / (c + L_1);
  phi_5 := ((k_1a * L_1 + l_2a) * c) / (c + L_1);

  // Variable initializations:
  k_1a = 0.64;
  k_1b = 0.04;
  k_2a = 37.4;
  k_2b = 1.4;
  k_3a = 0.11;
  k_3b = 29.8;
  k_4a = 4;
  k_4b = 0.54;
  l_2a = 1.7;
  l_2b = 0.8;
  l_4a = 1.7;
  l_4b = 2.5;
  l_6a = 4707;
  l_6b = 11.4;
  L_1 = 0.12;
  L_3 = 0.025;
  L_5 = 54.7;
  p = 10;
  c = 1;
end

model *sneyd_2002____main()

  // Sub-modules, and any changes to those submodules:
  environment: sneyd_2002__environment(time_);
  R: sneyd_2002__R(R0, phi_1, phi_2, phi_2b, k_1b, l_2b, I_10, O0, p, time_);
  O: sneyd_2002__O(O0, phi_3, phi_4, phi_4b, phi_2, phi_2b, k_3b, R0, A0, S0, p, time_);
  S: sneyd_2002__S(S0, I_10, I_20, O0, R0, A0);
  I_1: sneyd_2002__I_1(I_10, phi_1, k_1b, l_2b, R0, time_);
  I_2: sneyd_2002__I_2(I_20, phi_5, k_1b, l_2b, A0, time_);
  A: sneyd_2002__A(A0, phi_5, phi_4, phi_4b, k_1b, l_2b, I_20, O0, time_);
  open_probability: sneyd_2002__open_probability(O0, A0);
  reaction_rate_constants: sneyd_2002__reaction_rate_constants(phi_1, phi_2, phi_2b, phi_3, phi_4, phi_4b, phi_5, k_1b, k_3b, l_2b, p);
end