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

  // Variable initializations:
  time_ = ;
end

model alexander_2010__A(A, v_tilday, G, R, time_)

  // Rate Rules:
  A' = (f * v_tilday * G) / (1 + b4 + sigma4 * R) - muA * A;

  // Variable initializations:
  A = 1.0;
  v_tilday = 0.0025;
  f = 1e-4;
  muA = 0.25;
  b4 = 1.00;
  sigma4 = 1.2e-5;
  G = ;
  R = ;
  time_ = ;
end

model alexander_2010__R(R, A, E, time_)

  // Rate Rules:
  R' = (pi4 * E + beta) * A - muR * R;

  // Variable initializations:
  R = 0.0;
  pi4 = 0.016;
  beta = 200.0;
  muR = 0.25;
  A = ;
  E = ;
  time_ = ;
end

model alexander_2010__E(E, A, time_)

  // Rate Rules:
  E' = lambdaE * A - muE * E;

  // Variable initializations:
  E = 0.0;
  lambdaE = 1000.0;
  muE = 0.25;
  A = ;
  time_ = ;
end

model alexander_2010__G(G, v_tilday, E, time_)

  // Rate Rules:
  G' = gamma * E - (v_tilday * G + muG * G);

  // Variable initializations:
  G = 1e8;
  gamma = 2000.0;
  muG = 5.0;
  v_tilday = ;
  E = ;
  time_ = ;
end

model *alexander_2010____main()

  // Sub-modules, and any changes to those submodules:
  environment: alexander_2010__environment(time_);
  A: alexander_2010__A(A0, v_tilday, G0, R0, time_);
  R: alexander_2010__R(R0, A0, E0, time_);
  E: alexander_2010__E(E0, A0, time_);
  G: alexander_2010__G(G0, v_tilday, E0, time_);
end