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

  // Variable initializations:
  time_ = ;
end

model st_hilaire_model_2007_version01__H1(H1, A, betaIP, time_)

  // Rate Rules:
  H1' = -(betaIP * H1) + A;

  // Variable initializations:
  H1 = 0.5;
  A = ;
  betaIP = ;
  time_ = ;
end

model st_hilaire_model_2007_version01__H2(H2, H1, betaIP, betaIS, betaCP, time_)

  // Rate Rules:
  H2' = betaIP * H1 - (betaCP * H2 + betaIS * H2);

  // Variable initializations:
  H2 = 0.5;
  H1 = ;
  betaIP = ;
  betaIS = ;
  betaCP = ;
  time_ = ;
end

model st_hilaire_model_2007_version01__H3(H3, H2, betaIS, betaCS, time_)

  // Rate Rules:
  H3' = betaIS * H2 - betaCS * H3;

  // Variable initializations:
  H3 = 0.5;
  H2 = ;
  betaIS = ;
  betaCS = ;
  time_ = ;
end

model st_hilaire_model_2007_version01__A(A, A_, alpha, lamda, B_, time_)

  // Rate Rules:
  A' = piecewise( Amax * ((1.0 - exp(- lamda * (time_ - t_on))) / (1.0 - exp(- lamda * (t_off - t_on)))) , ( time_ < t_off) || (geq(time_, t_on) ), Amax * exp(- alpha * (time_ - t_off)) , geq(time_, t_off) );
  A_' = A * (1.0 - m * B_);

  // Variable initializations:
  A = ;
  A_ = ;
  Amax = 4.834;
  m = 7.0;
  t_on = 21.93;
  t_off = 29.867;
  alpha = ;
  lamda = ;
  B_ = ;
  time_ = ;
end

model st_hilaire_model_2007_version01__Process_L(B_, time_)

  // Assignment Rules:
  B_ := G * (1.0 - n) * alph;
  alph := alph_0 * power(I / I0, p) * (I / (I + 100.0));

  // Rate Rules:
  n' = 60.0 * (alph * (1.0 - n) - beta * n);
  I' = piecewise( 0.0 , ( geq(time_, 0.0)) && (time_ < 4.75 ), 9500.0 , ( geq(time_, 4.75)) && (time_ < 11.2 ), 0.0 );

  // Variable initializations:
  n = 7.0;
  p = 0.6;
  I = ;
  I0 = 9500.0;
  G = 37.0;
  alph_0 = 0.1;
  beta = 0.007;
  time_ = ;
end

model st_hilaire_model_2007_version01__Process_P(B_, time_)

  // Assignment Rules:
  B := B_ * (1.0 - 0.4 * x) * (1.0 - 0.4 * xc);

  // Rate Rules:
  x' = (pi / 12.0) * (xc + mu * (((1.0 / 3.0) * x + (4.0 / 3.0) * power(x, 3.0)) - (256.0 / 105.0) * power(x, 7.0)) + B);
  xc' = (pi / 12.0) * ((q * B * xc - x * power(24.0 / (0.99729 * tau_x), 2.0)) + k * B);

  // Variable initializations:
  x = ;
  xc = ;
  k = 0.55;
  mu = 0.13;
  q = 0.333333;
  tau_x = 24.2;
  B_ = ;
  time_ = ;
end

model st_hilaire_model_2007_version01__model_parameters(betaIP, betaCP, betaIS, betaCS, alpha, lamda)

  // Variable initializations:
  betaIP = 0.047;
  betaCP = 0.025;
  betaIS = 0.065;
  betaCS = 0.028;
  alpha = 0.067;
  lamda = 0.588;
end

model *st_hilaire_model_2007_version01____main()

  // Sub-modules, and any changes to those submodules:
  environment: st_hilaire_model_2007_version01__environment(time_);
  H1: st_hilaire_model_2007_version01__H1(H10, A0, betaIP, time_);
  H2: st_hilaire_model_2007_version01__H2(H20, H10, betaIP, betaIS, betaCP, time_);
  H3: st_hilaire_model_2007_version01__H3(H30, H20, betaIS, betaCS, time_);
  A: st_hilaire_model_2007_version01__A(A0, A_, alpha, lamda, B_, time_);
  Process_L: st_hilaire_model_2007_version01__Process_L(B_, time_);
  Process_P: st_hilaire_model_2007_version01__Process_P(B_, time_);
  model_parameters: st_hilaire_model_2007_version01__model_parameters(betaIP, betaCP, betaIS, betaCS, alpha, lamda);
end