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

  // Variable initializations:
  time_ = ;
end

model mittler_1998__T(T, k, p, c, delta)

  // Assignment Rules:
  T := (c * delta) / (k * p);

  // Variable initializations:
  k = ;
  p = ;
  c = ;
  delta = ;
end

model mittler_1998__I(I, delta, k_, T, E4, VI_0, c, p, time_)

  // Assignment Rules:
  I_0 := (c / p) * VI_0;

  // Rate Rules:
  I' = k_ * T * E4 - delta * I;

  // Variable initializations:
  I = 0.1;
  delta = ;
  k_ = ;
  T = ;
  E4 = ;
  VI_0 = ;
  c = ;
  p = ;
  time_ = ;
end

model mittler_1998__VI(VI, VI_0, h, I, c, p, time_)

  // Rate Rules:
  VI' = (1 - h) * p * I - c * VI;

  // Variable initializations:
  VI = 200000;
  VI_0 = 200000;
  h = ;
  I = ;
  c = ;
  p = ;
  time_ = ;
end

model mittler_1998__VNI(VNI, h, I, c, p, time_)

  // Rate Rules:
  VNI' = h * p * I - c * VNI;

  // Variable initializations:
  VNI = 0;
  h = ;
  I = ;
  c = ;
  p = ;
  time_ = ;
end

model mittler_1998__virus_total(VI, VNI)

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

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

model mittler_1998__E1(E1, b_, VI, time_)

  // Rate Rules:
  E1' = (VI - E1) / b_;

  // Variable initializations:
  E1 = 0;
  b_ = ;
  VI = ;
  time_ = ;
end

model mittler_1998__E2(E2, b_, E1, time_)

  // Rate Rules:
  E2' = (E1 - E2) / b_;

  // Variable initializations:
  E2 = 0;
  b_ = ;
  E1 = ;
  time_ = ;
end

model mittler_1998__E3(E3, b_, E2, time_)

  // Rate Rules:
  E3' = (E2 - E3) / b_;

  // Variable initializations:
  E3 = 0;
  b_ = ;
  E2 = ;
  time_ = ;
end

model mittler_1998__E4(E4, b_, E3, time_)

  // Rate Rules:
  E4' = (E3 - E4) / b_;

  // Variable initializations:
  E4 = 0;
  b_ = ;
  E3 = ;
  time_ = ;
end

model mittler_1998__Heavyside_function(h, time_)

  // Assignment Rules:
  h := piecewise( 0 , time_ < tau_p , 1 );

  // Variable initializations:
  tau_p = 0;
  time_ = ;
end

model mittler_1998__kinetic_parameters(c, delta, p, b_, k_, k)

  // Assignment Rules:
  b_ := b / (1 + m * b);
  k_ := k / power(1 + m * b, n);

  // Variable initializations:
  c = 3;
  delta = 0.5;
  p = 774;
  k = 2.4e-5;
  b = 0.25;
  m = 0.01;
  n = 4;
end

model *mittler_1998____main()

  // Sub-modules, and any changes to those submodules:
  environment: mittler_1998__environment(time_);
  T: mittler_1998__T(T0, k, p, c, delta);
  I: mittler_1998__I(I0, delta, k_, T0, E40, VI_0, c, p, time_);
  VI: mittler_1998__VI(VI0, VI_0, h, I0, c, p, time_);
  VNI: mittler_1998__VNI(VNI0, h, I0, c, p, time_);
  virus_total: mittler_1998__virus_total(VI0, VNI0);
  E1: mittler_1998__E1(E10, b_, VI0, time_);
  E2: mittler_1998__E2(E20, b_, E10, time_);
  E3: mittler_1998__E3(E30, b_, E20, time_);
  E4: mittler_1998__E4(E40, b_, E30, time_);
  Heavyside_function: mittler_1998__Heavyside_function(h, time_);
  kinetic_parameters: mittler_1998__kinetic_parameters(c, delta, p, b_, k_, k);
end