//Created by libAntimony v2.4
model jelic_cupic_kolaranic_2005_Fig5__environment(time_, tau, C_0)

  // Assignment Rules:
  time_ := tau / C_0;

  // Variable initializations:
  tau = ;
  C_0 = ;
end

model jelic_cupic_kolaranic_2005_Fig5__a(a, alpha, beta, KP, C_1, g, tau)

  // Assignment Rules:
  A := C_1 * a;

  // Rate Rules:
  a' = KP - ((1 + alpha + beta) * a + a * power(g, 2));

  // Variable initializations:
  a = 4.39927;
  alpha = ;
  beta = ;
  KP = ;
  C_1 = ;
  g = ;
  tau = ;
end

model jelic_cupic_kolaranic_2005_Fig5__g(g, alpha, gamma, L, C_2, a, tau)

  // Assignment Rules:
  G := C_2 * g;

  // Rate Rules:
  g' = ((1 - alpha) * a + a * power(g, 2)) - (L + gamma * g);

  // Variable initializations:
  g = 1.96477;
  alpha = ;
  gamma = ;
  L = ;
  C_2 = ;
  a = ;
  tau = ;
end

model jelic_cupic_kolaranic_2005_Fig5__reaction_constants(alpha, beta, gamma, L, KP, time_, C_0, C_1, C_2)

  // Assignment Rules:
  KP := K * P;
  P := 0.735 + 0.1563 * sin((2 * 3.141592653 * time_) / (24 * 3600)) + 0.3127 * abs(sin((3.141592653 * time_) / (24 * 3600)));
  C_0 := k2;
  C_1 := root(k2 / k4);
  C_2 := root(k2 / k4);

  // Variable initializations:
  alpha = 0.008;
  beta = 1.485;
  gamma = 11.385;
  L = 0.1;
  k2 = 6e-4;
  k3 = 0.0000048;
  k6 = 0.000891;
  k7 = 0.006831;
  k0 = 8.7831e-11;
  k4 = 2.1e12;
  km = 6.9001e-14;
  K = 30;
  time_ = ;
end

model *jelic_cupic_kolaranic_2005_Fig5____main()

  // Sub-modules, and any changes to those submodules:
  environment: jelic_cupic_kolaranic_2005_Fig5__environment(time_, tau, C_0);
  a: jelic_cupic_kolaranic_2005_Fig5__a(a0, alpha, beta, KP, C_1, g0, tau);
  g: jelic_cupic_kolaranic_2005_Fig5__g(g0, alpha, gamma, L, C_2, a0, tau);
  reaction_constants: jelic_cupic_kolaranic_2005_Fig5__reaction_constants(alpha, beta, gamma, L, KP, time_, C_0, C_1, C_2);
end