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

  // Variable initializations:
  time_ = ;
end

model obeyesekere_1999__parameters(a_D, a_E, a_X, k, q_D, q_E, q_X, f, g, p_S, p_D, p_E, p_X, d_D, d_E, d_X, a_f, R_T, theta, GF)

  // Variable initializations:
  a_D = 0.4;
  a_E = 0.16;
  a_X = 0.08;
  k = 0.05;
  q_D = 0.6;
  q_E = 0.6;
  q_X = 0.8;
  f = 0.2;
  g = 0.528;
  p_S = 0.6;
  p_D = 0.48;
  p_E = 0.096;
  p_X = 0.48;
  d_D = 0.4;
  d_E = 0.2;
  d_X = 1.04;
  a_f = 0.9;
  R_T = 2.5;
  theta = 1.5;
  GF = 6.3;
end

model obeyesekere_1999__D(D, a_D, d_D, k, GF, E, time_)

  // Rate Rules:
  D' = (a_D * k * GF) / (1 + k * GF) - d_D * E * D;

  // Variable initializations:
  D = 0.1;
  a_D = ;
  d_D = ;
  k = ;
  GF = ;
  E = ;
  time_ = ;
end

model obeyesekere_1999__E(E, a_E, a_f, d_E, theta, R_S, X, time_)

  // Rate Rules:
  E' = a_E * (1 + a_f * (theta - R_S)) * 1 - d_E * X * E;

  // Variable initializations:
  E = 0.6;
  a_E = ;
  a_f = ;
  d_E = ;
  theta = ;
  R_S = ;
  X = ;
  time_ = ;
end

model obeyesekere_1999__R(R, p_X, p_S, X, R_T, q_X, R_S, theta, time_)

  // Rate Rules:
  R' = (p_X * X * ((R_T - R_S) - R)) / (q_X + ((R_T - R_S) - R) + X) - p_S * (theta - R_S) * R;

  // Variable initializations:
  R = 0.5;
  p_X = ;
  p_S = ;
  X = ;
  R_T = ;
  q_X = ;
  R_S = ;
  theta = ;
  time_ = ;
end

model obeyesekere_1999__R_S(R_S, p_S, p_D, p_E, theta, R, D, q_D, q_E, E, time_)

  // Rate Rules:
  R_S' = (p_S * (theta - R_S) * R - (p_D * R_S * D) / (q_D + R_S + D)) - (p_E * R_S * E) / (q_E + R_S + E);

  // Variable initializations:
  R_S = 1;
  p_S = ;
  p_D = ;
  p_E = ;
  theta = ;
  R = ;
  D = ;
  q_D = ;
  q_E = ;
  E = ;
  time_ = ;
end

model obeyesekere_1999__X(X, a_X, f, g, d_X, E, theta, R_S, time_)

  // Rate Rules:
  X' = (a_X * E + f * (theta - R_S) + g * power(X, 2) * E) - d_X * X;

  // Variable initializations:
  X = 0.7;
  a_X = ;
  f = ;
  g = ;
  d_X = ;
  E = ;
  theta = ;
  R_S = ;
  time_ = ;
end

model obeyesekere_1999__unpho_RB(unpho_RB, R, R_S)

  // Assignment Rules:
  unpho_RB := R + R_S;

  // Variable initializations:
  R = ;
  R_S = ;
end

model obeyesekere_1999__free_E2F(free_E2F, theta, R_S)

  // Assignment Rules:
  free_E2F := theta - R_S;

  // Variable initializations:
  theta = ;
  R_S = ;
end

model *obeyesekere_1999____main()

  // Sub-modules, and any changes to those submodules:
  environment: obeyesekere_1999__environment(time_);
  parameters: obeyesekere_1999__parameters(a_D, a_E, a_X, k, q_D, q_E, q_X, f, g, p_S, p_D, p_E, p_X, d_D, d_E, d_X, a_f, R_T, theta, GF);
  D: obeyesekere_1999__D(D0, a_D, d_D, k, GF, E0, time_);
  E: obeyesekere_1999__E(E0, a_E, a_f, d_E, theta, R_S0, X0, time_);
  R: obeyesekere_1999__R(R0, p_X, p_S, X0, R_T, q_X, R_S0, theta, time_);
  R_S: obeyesekere_1999__R_S(R_S0, p_S, p_D, p_E, theta, R0, D0, q_D, q_E, E0, time_);
  X: obeyesekere_1999__X(X0, a_X, f, g, d_X, E0, theta, R_S0, time_);
  unpho_RB: obeyesekere_1999__unpho_RB(unpho_RB0, R0, R_S0);
  free_E2F: obeyesekere_1999__free_E2F(free_E2F0, theta, R_S0);
end