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

  // Variable initializations:
  time_ = ;
end

model schenzle_1994__X(X, time_, K, V)

  // Assignment Rules:
  s := mu0 * 1.0E11;

  // Rate Rules:
  X' = s - (mu0 * X + K * V * X);

  // Variable initializations:
  X = 1.0E11;
  mu0 = 0.004;
  time_ = ;
  K = ;
  V = ;
end

model schenzle_1994__Y(Y, time_, K, V, X, Z)

  // Rate Rules:
  Y' = K * V * Y - (mu1 + delta1 * Z) * Y;

  // Variable initializations:
  Y = 100.0;
  mu1 = 0.30;
  delta1 = 0.30;
  time_ = ;
  K = ;
  V = ;
  X = ;
  Z = ;
end

model schenzle_1994__CD4(X, Y)

  // Assignment Rules:
  CD4 := (X + Y) / 1.0E11;

  // Variable initializations:
  X = ;
  Y = ;
end

model schenzle_1994__V(V, time_, Y, Z)

  // Rate Rules:
  V' = phi * Y - (mu2 + delta2 * Z) * V;

  // Variable initializations:
  V = 1000.0;
  mu2 = 1.0;
  phi = 1000.0;
  delta2 = 1.0;
  time_ = ;
  Y = ;
  Z = ;
end

model schenzle_1994__Z(Z, time_, V, X)

  // Assignment Rules:
  f_X := power(X / X0, 2.0) / (b + power(X / X0, 2.0));

  // Rate Rules:
  Z' = (alpha * (V / (a + V)) + beta * f_X * Z) - gamma * power(Z, 2);

  // Variable initializations:
  Z = 0.0;
  alpha = 1.0E-6;
  beta = 0.5;
  gamma = 0.5;
  a = 1000.0;
  b = 0.2;
  X0 = 1.0E11;
  time_ = ;
  V = ;
  X = ;
end

model schenzle_1994__K(K, time_, V)

  // Assignment Rules:
  Kmax := 1.35E-14 * 20.0;

  // Rate Rules:
  K' = omega * V * (Kmax - K);

  // Variable initializations:
  K = 1.35E-14;
  omega = 0.16E-15;
  time_ = ;
  V = ;
end

model *schenzle_1994____main()

  // Sub-modules, and any changes to those submodules:
  environment: schenzle_1994__environment(time_);
  X: schenzle_1994__X(X0, time_, K0, V0);
  Y: schenzle_1994__Y(Y0, time_, K0, V0, X0, Z0);
  CD4: schenzle_1994__CD4(X0, Y0);
  V: schenzle_1994__V(V0, time_, Y0, Z0);
  Z: schenzle_1994__Z(Z0, time_, V0, X0);
  K: schenzle_1994__K(K0, time_, V0);
end