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

  // Variable initializations:
  time_ = ;
end

model bonhoeffer_2000__T(T, b, I, time_)

  // Rate Rules:
  T' = s - (dT * T + b * T * I);

  // Variable initializations:
  T = 1.0;
  s = 10.0;
  dT = 0.01;
  b = ;
  I = ;
  time_ = ;
end

model bonhoeffer_2000__I(I, b, ql, qa, T, E, time_)

  // Rate Rules:
  I' = (b * T * I - (dI * I + p * E * I + ql * I)) + qa * I;

  // Variable initializations:
  I = 1.0;
  p = 0.05;
  dI = 0.3;
  b = ;
  ql = ;
  qa = ;
  T = ;
  E = ;
  time_ = ;
end

model bonhoeffer_2000__Il(Il, ql, qa, I, time_)

  // Rate Rules:
  Il' = ql * I - (al * Il + qa * Il);

  // Variable initializations:
  Il = 1.0;
  al = 0.01;
  ql = ;
  qa = ;
  I = ;
  time_ = ;
end

model bonhoeffer_2000__E(E, I, time_)

  // Rate Rules:
  E' = (c * E * I) / (I + K) - (dE * E + (dE_ * E * I) / (K_ + I));

  // Variable initializations:
  E = 1.0;
  c = 0.3;
  dE = 0.1;
  K = 0.1;
  dE_ = 0.25;
  K_ = 05.0;
  I = ;
  time_ = ;
end

model bonhoeffer_2000__kinetic_parameters(b, qa, ql)

  // Variable initializations:
  b = 0.001;
  qa = 0.001;
  ql = 0.001;
end

model *bonhoeffer_2000____main()

  // Sub-modules, and any changes to those submodules:
  environment: bonhoeffer_2000__environment(time_);
  T: bonhoeffer_2000__T(T0, b, I0, time_);
  I: bonhoeffer_2000__I(I0, b, ql, qa, T0, E0, time_);
  Il: bonhoeffer_2000__Il(Il0, ql, qa, I0, time_);
  E: bonhoeffer_2000__E(E0, I0, time_);
  kinetic_parameters: bonhoeffer_2000__kinetic_parameters(b, qa, ql);
end