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

  // Variable initializations:
  time_ = ;
end

model perelson_1993b__uninfected(time_, mu_T, k_1, T_1, T_2, V, T)

  // Assignment Rules:
  s_V := (theta * s) / (theta + V);

  // Rate Rules:
  T' = ((s_V - mu_T * T) + r * T * (1 - (T + T_1 + T_2) / T_max)) - k_1 * V * T;

  // Variable initializations:
  time_ = ;
  s = 10;
  r = 0.03;
  T_max = 1500;
  mu_T = 0.02;
  theta = 1;
  k_1 = ;
  T_1 = ;
  T_2 = ;
  V = ;
  T = 1000;
end

model perelson_1993b__latently_infected(time_, k_1, k_2, mu_T, V, T, T_1)

  // Rate Rules:
  T_1' = (k_1 * V * T - mu_T * T_1) - k_2 * T_1;

  // Variable initializations:
  time_ = ;
  k_1 = 2.4E-5;
  k_2 = ;
  mu_T = ;
  V = ;
  T = ;
  T_1 = 0;
end

model perelson_1993b__actively_infected(time_, k_2, mu_b, T_1, T_2)

  // Rate Rules:
  T_2' = k_2 * T_1 - mu_b * T_2;

  // Variable initializations:
  time_ = ;
  k_2 = 3E-3;
  mu_b = 0.24;
  T_1 = ;
  T_2 = 0;
end

model perelson_1993b__free_virus_particle(time_, mu_b, k_1, T, T_2, V)

  // Rate Rules:
  V' = (N * mu_b * T_2 - k_1 * V * T) - mu_V * V;

  // Variable initializations:
  time_ = ;
  mu_b = ;
  k_1 = ;
  mu_V = 2.4;
  N = 1400;
  T = ;
  T_2 = ;
  V = 1.0E-3;
end

model perelson_1993b__T_cell_population(T, T_1, T_2)

  // Assignment Rules:
  T_tot := T + T_1 + T_2;

  // Variable initializations:
  T = ;
  T_1 = ;
  T_2 = ;
end

model *perelson_1993b____main()

  // Sub-modules, and any changes to those submodules:
  environment: perelson_1993b__environment(time_);
  uninfected: perelson_1993b__uninfected(time_, mu_T, k_1, T_1, T_2, V, T);
  latently_infected: perelson_1993b__latently_infected(time_, k_1, k_2, mu_T, V, T, T_1);
  actively_infected: perelson_1993b__actively_infected(time_, k_2, mu_b, T_1, T_2);
  free_virus_particle: perelson_1993b__free_virus_particle(time_, mu_b, k_1, T, T_2, V);
  T_cell_population: perelson_1993b__T_cell_population(T, T_1, T_2);
end