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

  // Variable initializations:
  time_ = ;
end

model perelson_neumann_markowitz_leonard_ho_1996__T(T)

  // Variable initializations:
  T = 170000;
end

model perelson_neumann_markowitz_leonard_ho_1996__T_star(T_star, T, VI, delta, time_)

  // Rate Rules:
  T_star' = k * VI * T - delta * T_star;

  // Variable initializations:
  T_star = 0;
  k = 2.4e-5;
  T = ;
  VI = ;
  delta = ;
  time_ = ;
end

model perelson_neumann_markowitz_leonard_ho_1996__VI(VI, c, time_)

  // Rate Rules:
  VI' = -(c * VI);

  // Variable initializations:
  VI = 216000.0;
  c = ;
  time_ = ;
end

model perelson_neumann_markowitz_leonard_ho_1996__VNI(VNI, T_star, delta, c, time_)

  // Rate Rules:
  VNI' = N * delta * T_star - c * VNI;

  // Variable initializations:
  VNI = 0;
  N = 774;
  T_star = ;
  delta = ;
  c = ;
  time_ = ;
end

model perelson_neumann_markowitz_leonard_ho_1996__V(V, VI, VNI)

  // Assignment Rules:
  V := VI + VNI;

  // Variable initializations:
  VI = ;
  VNI = ;
end

model perelson_neumann_markowitz_leonard_ho_1996__kinetic_parameters(c, delta)

  // Variable initializations:
  c = 3.07;
  delta = 0.49;
end

model *perelson_neumann_markowitz_leonard_ho_1996____main()

  // Sub-modules, and any changes to those submodules:
  environment: perelson_neumann_markowitz_leonard_ho_1996__environment(time_);
  T: perelson_neumann_markowitz_leonard_ho_1996__T(T0);
  T_star: perelson_neumann_markowitz_leonard_ho_1996__T_star(T_star0, T0, VI0, delta, time_);
  VI: perelson_neumann_markowitz_leonard_ho_1996__VI(VI0, c, time_);
  VNI: perelson_neumann_markowitz_leonard_ho_1996__VNI(VNI0, T_star0, delta, c, time_);
  V: perelson_neumann_markowitz_leonard_ho_1996__V(V0, VI0, VNI0);
  kinetic_parameters: perelson_neumann_markowitz_leonard_ho_1996__kinetic_parameters(c, delta);
end