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

  // Variable initializations:
  time_ = ;
end

model heinze_keener_midgley_model_1998_version01__GnRH(GnRH, time_)

  // Assignment Rules:
  GnRH := piecewise( 0.5 , ( geq(time_, 0.0)) && (time_ < 0.0666667 ), 0.0 , ( geq(time_, 0.0666667)) && (time_ < 0.4 ), 0.5 , ( geq(time_, 0.4)) && (time_ < 0.4666667 ), 0.0 , ( geq(time_, 0.4666667)) && (time_ < 2.4666667 ), 0.5 , ( geq(time_, 2.4666667)) && (time_ < 2.533333 ), 0.0 , ( geq(time_, 2.533333)) && (time_ < 2.6166667 ), 0.5 , ( geq(time_, 2.6166667)) && (time_ < 2.6833333 ), 0.0 , ( geq(time_, 2.6833333)) && (time_ < 4.6833333 ), 0.5 , ( geq(time_, 4.6833333)) && (time_ < 4.75 ), 0.0 , ( geq(time_, 4.75)) && (time_ < 4.916667 ), 0.5 , ( geq(time_, 4.916667)) && (time_ < 4.983333 ), 0.0 , ( geq(time_, 4.983333)) && (time_ < 6.983333 ), 0.5 , ( geq(time_, 6.983333)) && (time_ < 7.066667 ), 0.0 , ( geq(time_, 7.066667)) && (time_ < 7.733333 ), 0.5 , ( geq(time_, 7.733333)) && (time_ < 7.8 ), 0.0 , ( geq(time_, 7.8)) && (time_ < 9.8 ) );

  // Variable initializations:
  time_ = ;
end

model heinze_keener_midgley_model_1998_version01__F(F, kfb, kdf, D, GnRH, time_)

  // Rate Rules:
  F' = kdf * D - kfb * F * GnRH;

  // Variable initializations:
  F = 1.0;
  kfb = ;
  kdf = ;
  D = ;
  GnRH = ;
  time_ = ;
end

model heinze_keener_midgley_model_1998_version01__D(D, kbd, kdf, B, time_)

  // Rate Rules:
  D' = kbd * B - kdf * D;

  // Variable initializations:
  D = 0.0;
  kbd = ;
  kdf = ;
  B = ;
  time_ = ;
end

model heinze_keener_midgley_model_1998_version01__B(B, kfb, kbd, F, GnRH, time_)

  // Rate Rules:
  B' = kfb * F * GnRH - kbd * B;

  // Variable initializations:
  B = 0.0;
  kfb = ;
  kbd = ;
  F = ;
  GnRH = ;
  time_ = ;
end

model heinze_keener_midgley_model_1998_version01__R(R, s, a1, a2, B, time_)

  // Rate Rules:
  R' = s - (a1 + a2 * B) * R;

  // Variable initializations:
  R = 2115.0;
  s = ;
  a1 = ;
  a2 = ;
  B = ;
  time_ = ;
end

model heinze_keener_midgley_model_1998_version01__C(C, a1, a2, B, R, time_)

  // Rate Rules:
  C' = (a1 + a2 * B) * R;

  // Variable initializations:
  C = 0.0;
  a1 = ;
  a2 = ;
  B = ;
  R = ;
  time_ = ;
end

model heinze_keener_midgley_model_1998_version01__model_parameters(kfb, kdf, kbd, a1, a2, s)

  // Variable initializations:
  kfb = 19.35;
  kdf = 2.52;
  kbd = 9.91;
  a1 = 0.0328;
  a2 = 0.0830;
  s = 6.80;
end

model *heinze_keener_midgley_model_1998_version01____main()

  // Sub-modules, and any changes to those submodules:
  environment: heinze_keener_midgley_model_1998_version01__environment(time_);
  GnRH: heinze_keener_midgley_model_1998_version01__GnRH(GnRH0, time_);
  F: heinze_keener_midgley_model_1998_version01__F(F0, kfb, kdf, D0, GnRH0, time_);
  D: heinze_keener_midgley_model_1998_version01__D(D0, kbd, kdf, B0, time_);
  B: heinze_keener_midgley_model_1998_version01__B(B0, kfb, kbd, F0, GnRH0, time_);
  R: heinze_keener_midgley_model_1998_version01__R(R0, s, a1, a2, B0, time_);
  C: heinze_keener_midgley_model_1998_version01__C(C0, a1, a2, B0, R0, time_);
  model_parameters: heinze_keener_midgley_model_1998_version01__model_parameters(kfb, kdf, kbd, a1, a2, s);
end