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

  // Variable initializations:
  time_ = ;
end

model wodarz_2003__S(S, d, g, beta_1, beta_2, R_1, R_2, R_12, P_1, P_2, time_, I_1, I_2, I_12, I_21)

  // Assignment Rules:
  H := S + I_1 + R_1 + I_2 + R_2 + I_12 + I_21 + R_12;
  P := P_1 + P_2;

  // Rate Rules:
  S' = ((((r * H) / (epsilon * H + 1) - d * S) - beta_1 * S * P_1) - beta_2 * S * P_2) + g * (R_1 + R_2 + R_12);

  // Variable initializations:
  S = 100;
  r = 0.5;
  epsilon = 0.1;
  d = ;
  g = ;
  beta_1 = ;
  beta_2 = ;
  R_1 = ;
  R_2 = ;
  R_12 = ;
  P_1 = ;
  P_2 = ;
  time_ = ;
  I_1 = ;
  I_2 = ;
  I_12 = ;
  I_21 = ;
end

model wodarz_2003__I1(I_1, a_1, alpha_1, beta_1, S, P_1, time_)

  // Rate Rules:
  I_1' = (beta_1 * S * P_1 - a_1 * I_1) - alpha_1 * I_1;

  // Variable initializations:
  I_1 = 0;
  a_1 = ;
  alpha_1 = ;
  beta_1 = ;
  S = ;
  P_1 = ;
  time_ = ;
end

model wodarz_2003__I2(I_2, a_2, alpha_2, beta_2, S, P_2, time_)

  // Rate Rules:
  I_2' = (beta_2 * S * P_2 - a_2 * I_2) - alpha_2 * I_2;

  // Variable initializations:
  I_2 = 0;
  a_2 = ;
  alpha_2 = ;
  beta_2 = ;
  S = ;
  P_2 = ;
  time_ = ;
end

model wodarz_2003__I12(I_12, a_2, alpha_2, beta_2, R_1, P_2, time_)

  // Rate Rules:
  I_12' = (beta_2 * R_1 * P_2 - a_2 * I_12) - alpha_2 * I_12;

  // Variable initializations:
  I_12 = 0;
  a_2 = ;
  alpha_2 = ;
  beta_2 = ;
  R_1 = ;
  P_2 = ;
  time_ = ;
end

model wodarz_2003__I21(I_21, a_1, alpha_1, beta_1, R_2, P_1, time_)

  // Rate Rules:
  I_21' = (beta_1 * R_2 * P_1 - a_1 * I_21) - alpha_1 * I_21;

  // Variable initializations:
  I_21 = 0;
  a_1 = ;
  alpha_1 = ;
  beta_1 = ;
  R_2 = ;
  P_1 = ;
  time_ = ;
end

model wodarz_2003__R1(R_1, alpha_1, beta_2, d, g, I_1, P_2, time_)

  // Rate Rules:
  R_1' = ((alpha_1 * I_1 - d * R_1) - g * R_1) - beta_2 * R_1 * P_2;

  // Variable initializations:
  R_1 = 0;
  alpha_1 = ;
  beta_2 = ;
  d = ;
  g = ;
  I_1 = ;
  P_2 = ;
  time_ = ;
end

model wodarz_2003__R2(R_2, alpha_2, beta_1, d, g, I_2, P_1, time_)

  // Rate Rules:
  R_2' = ((alpha_2 * I_2 - d * R_2) - g * R_2) - beta_1 * R_2 * P_1;

  // Variable initializations:
  R_2 = 0;
  alpha_2 = ;
  beta_1 = ;
  d = ;
  g = ;
  I_2 = ;
  P_1 = ;
  time_ = ;
end

model wodarz_2003__R12(R_12, alpha_2, alpha_1, d, g, I_21, I_12, time_)

  // Rate Rules:
  R_12' = ((alpha_2 * I_12 + alpha_1 * I_21) - d * R_12) - g * R_12;

  // Variable initializations:
  R_12 = 0;
  alpha_2 = ;
  alpha_1 = ;
  d = ;
  g = ;
  I_21 = ;
  I_12 = ;
  time_ = ;
end

model wodarz_2003__P1(P_1, u, I_1, I_21, time_)

  // Assignment Rules:
  log_P1 := log(P_1);

  // Rate Rules:
  P_1' = k_1 * (I_1 + I_21) - u * P_1;

  // Variable initializations:
  P_1 = 1;
  k_1 = 1;
  u = ;
  I_1 = ;
  I_21 = ;
  time_ = ;
end

model wodarz_2003__P2(P_2, u, I_2, I_12, time_)

  // Assignment Rules:
  log_P2 := log(P_2);

  // Rate Rules:
  P_2' = k_2 * (I_2 + I_12) - u * P_2;

  // Variable initializations:
  P_2 = 1;
  k_2 = 1;
  u = ;
  I_2 = ;
  I_12 = ;
  time_ = ;
end

model wodarz_2003__memory_duration(g)

  // Assignment Rules:
  G := 1 / g;

  // Variable initializations:
  g = 0.01;
end

model wodarz_2003__kinetic_parameters(beta_1, alpha_1, beta_2, alpha_2, a_1, a_2, d, u)

  // Variable initializations:
  beta_1 = 1;
  alpha_1 = 0.1;
  beta_2 = 1;
  alpha_2 = 0.1;
  a_1 = 0.03;
  a_2 = 1;
  d = 0.01;
  u = 0.5;
end

model *wodarz_2003____main()

  // Sub-modules, and any changes to those submodules:
  environment: wodarz_2003__environment(time_);
  S: wodarz_2003__S(S0, d, g, beta_1, beta_2, R_1, R_2, R_12, P_1, P_2, time_, I_1, I_2, I_12, I_21);
  I1: wodarz_2003__I1(I_1, a_1, alpha_1, beta_1, S0, P_1, time_);
  I2: wodarz_2003__I2(I_2, a_2, alpha_2, beta_2, S0, P_2, time_);
  I12: wodarz_2003__I12(I_12, a_2, alpha_2, beta_2, R_1, P_2, time_);
  I21: wodarz_2003__I21(I_21, a_1, alpha_1, beta_1, R_2, P_1, time_);
  R1: wodarz_2003__R1(R_1, alpha_1, beta_2, d, g, I_1, P_2, time_);
  R2: wodarz_2003__R2(R_2, alpha_2, beta_1, d, g, I_2, P_1, time_);
  R12: wodarz_2003__R12(R_12, alpha_2, alpha_1, d, g, I_21, I_12, time_);
  P1: wodarz_2003__P1(P_1, u, I_1, I_21, time_);
  P2: wodarz_2003__P2(P_2, u, I_2, I_12, time_);
  memory_duration: wodarz_2003__memory_duration(g);
  kinetic_parameters: wodarz_2003__kinetic_parameters(beta_1, alpha_1, beta_2, alpha_2, a_1, a_2, d, u);
end