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

  // Variable initializations:
  time_ = ;
end

model nelson_1995a__uninfected_CD4(time_, mu_T, k_1, k_1_, T_1, V, V_, T)

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

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

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

model nelson_1995a__latently_infected_CD4(time_, k_1, k_1_, k_2, mu_T, V, V_, T, T_1)

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

  // Variable initializations:
  time_ = ;
  k_1 = ;
  k_1_ = ;
  k_2 = 0.017;
  mu_T = ;
  V = ;
  V_ = ;
  T = ;
  T_1 = 0;
end

model nelson_1995a__actively_infected_CD4(time_, k_2, k_1D, mu_b, D, T_1, T_2)

  // Rate Rules:
  T_2' = (k_2 * T_1 - k_1D * D * T_2) - mu_b * T_2;

  // Variable initializations:
  time_ = ;
  k_2 = ;
  k_1D = 2.4E-6;
  mu_b = 0.24;
  D = ;
  T_1 = ;
  T_2 = 0;
end

model nelson_1995a__actively_coinfected_CD4(time_, k_1D, k_s, mu_bD, D, T_2, T_D2)

  // Rate Rules:
  T_D2' = (k_1D * D * T_2 - k_s * T_D2) - mu_bD * T_D2;

  // Variable initializations:
  time_ = ;
  k_1D = ;
  k_s = 0.24;
  mu_bD = 0.17;
  D = ;
  T_2 = ;
  T_D2 = 0;
end

model nelson_1995a__stably_coinfected_CD4(time_, k_s, T_D2, T_D1)

  // Rate Rules:
  T_D1' = k_s * T_D2 - mu_TD * T_D1;

  // Variable initializations:
  time_ = ;
  mu_TD = 0.02;
  k_s = ;
  T_D2 = ;
  T_D1 = 0;
end

model nelson_1995a__DIV(time_, mu_bD, N_t, mu_b, T_D2, T_D1, D)

  // Assignment Rules:
  N_D_t := 0.2 * N_t;
  pi_D_t := 0.3 * mu_b * N_t;

  // Rate Rules:
  D' = (N_D_t * mu_bD * T_D2 + pi_D_t * T_D1) - mu_D * D;

  // Variable initializations:
  time_ = ;
  mu_bD = ;
  mu_D = 2.4;
  N_t = ;
  mu_b = ;
  T_D2 = ;
  T_D1 = ;
  D = 0;
end

model nelson_1995a__HIV1(time_, mu_b, mu_bD, k_1, mu_V, N_t, T, T_2, T_D2, V)

  // Assignment Rules:
  N_2_t := 0.6 * N_t;

  // Rate Rules:
  V' = ((N_t * mu_b * T_2 + N_2_t * mu_bD * T_D2) - k_1 * V * T) - mu_V * V;

  // Variable initializations:
  time_ = ;
  mu_b = ;
  mu_bD = ;
  k_1 = 2.4E-5;
  mu_V = 2.4;
  N_t = ;
  T = ;
  T_2 = ;
  T_D2 = ;
  V = 1E-3;
end

model nelson_1995a__hybrid_HIV1(time_, mu_bD, mu_b, mu_V, k_1_, N_t, T, T_D2, T_D1, V_)

  // Assignment Rules:
  N_t_ := 0.2 * N_t;
  pi_t_ := 0.1 * mu_b * N_t;

  // Rate Rules:
  V_' = ((N_t_ * mu_bD * T_D2 + pi_t_ * T_D1) - mu_V * V_) - k_1_ * T * V_;

  // Variable initializations:
  time_ = ;
  mu_bD = ;
  mu_b = ;
  mu_V = ;
  k_1_ = 2.4E-6;
  N_t = ;
  T = ;
  T_D2 = ;
  T_D1 = ;
  V_ = 0;
end

model nelson_1995a__production_function(time_, N_t)

  // Assignment Rules:
  N_t := N_0 * (1 + gamma * (power(time_, 2) / (power(time_, 2) + power(t_c, 2))));

  // Variable initializations:
  time_ = ;
  N_0 = 300;
  gamma = 25;
  t_c = 7305;
end

model nelson_1995a__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 *nelson_1995a____main()

  // Sub-modules, and any changes to those submodules:
  environment: nelson_1995a__environment(time_);
  uninfected_CD4: nelson_1995a__uninfected_CD4(time_, mu_T, k_1, k_1_, T_1, V, V_, T);
  latently_infected_CD4: nelson_1995a__latently_infected_CD4(time_, k_1, k_1_, k_2, mu_T, V, V_, T, T_1);
  actively_infected_CD4: nelson_1995a__actively_infected_CD4(time_, k_2, k_1D, mu_b, D, T_1, T_2);
  actively_coinfected_CD4: nelson_1995a__actively_coinfected_CD4(time_, k_1D, k_s, mu_bD, D, T_2, T_D2);
  stably_coinfected_CD4: nelson_1995a__stably_coinfected_CD4(time_, k_s, T_D2, T_D1);
  DIV: nelson_1995a__DIV(time_, mu_bD, N_t, mu_b, T_D2, T_D1, D);
  HIV1: nelson_1995a__HIV1(time_, mu_b, mu_bD, k_1, mu_V, N_t, T, T_2, T_D2, V);
  hybrid_HIV1: nelson_1995a__hybrid_HIV1(time_, mu_bD, mu_b, mu_V, k_1_, N_t, T, T_D2, T_D1, V_);
  production_function: nelson_1995a__production_function(time_, N_t);
  cell_population: nelson_1995a__cell_population(T, T_1, T_2);
end