//Created by libAntimony v2.4 model rattanakul_model_2003__environment(time_) // Variable initializations: time_ = ; end model rattanakul_model_2003__x(x, a1, k1, b1, kp, y, time_) // Rate Rules: x' = piecewise( a1 / (k1 + y) - b1 * x , ( geq(time_, 0.0)) && (time_ < 10000.0 ), (a1 / (k1 + y) - b1 * x) + kp , ( geq(time_, 10000.0)) && (time_ < 10360.0 ), a1 / (k1 + y) - b1 * x , ( geq(time_, 10360.0)) && (time_ < 11440.0 ), (a1 / (k1 + y) - b1 * x) + kp , ( geq(time_, 11440.0)) && (time_ < 11800.0 ), a1 / (k1 + y) - b1 * x , ( geq(time_, 11800.0)) && (time_ < 12880.0 ), (a1 / (k1 + y) - b1 * x) + kp , ( geq(time_, 12880.0)) && (time_ < 13240.0 ), a1 / (k1 + y) - b1 * x , ( geq(time_, 13240.0)) && (time_ < 14320.0 ), (a1 / (k1 + y) - b1 * x) + kp , ( geq(time_, 14320.0)) && (time_ < 14680.0 ), a1 / (k1 + y) - b1 * x , ( geq(time_, 14680.0)) && (time_ < 15760.0 ), (a1 / (k1 + y) - b1 * x) + kp , ( geq(time_, 15760.0)) && (time_ < 16120.0 ) ); // Variable initializations: x = 2.0; a1 = ; k1 = ; b1 = ; kp = ; y = ; time_ = ; end model rattanakul_model_2003__y(y, epsilon, a2, a3, k2, b2, x, z, time_) // Rate Rules: y' = epsilon * (((a2 + a3 * x) * y * z) / (k2 + power(x, 2.0)) - b2 * y); // Variable initializations: y = 1.0; epsilon = ; a2 = ; a3 = ; k2 = ; b2 = ; x = ; z = ; time_ = ; end model rattanakul_model_2003__z(z, epsilon, delta, a4, a5, k3, b3, x, time_) // Rate Rules: z' = epsilon * delta * (a4 * x - ((a5 * x * z) / (k3 + x) + b3 * z)); // Variable initializations: z = 0.15; epsilon = ; delta = ; a4 = ; a5 = ; k3 = ; b3 = ; x = ; time_ = ; end model rattanakul_model_2003__model_parameters(epsilon, delta, a1, a2, a3, a4, a5, b1, b2, b3, k1, k2, k3, kp) // Variable initializations: epsilon = 0.1; delta = 0.9; a1 = 0.05; a2 = 0.009; a3 = 0.675; a4 = 0.01; a5 = 0.005; b1 = 0.1; b2 = 0.3; b3 = 0.01; k1 = 0.1; k2 = 0.5; k3 = 0.025; kp = 0.5; end model *rattanakul_model_2003____main() // Sub-modules, and any changes to those submodules: environment: rattanakul_model_2003__environment(time_); x: rattanakul_model_2003__x(x0, a1, k1, b1, kp, y0, time_); y: rattanakul_model_2003__y(y0, epsilon, a2, a3, k2, b2, x0, z0, time_); z: rattanakul_model_2003__z(z0, epsilon, delta, a4, a5, k3, b3, x0, time_); model_parameters: rattanakul_model_2003__model_parameters(epsilon, delta, a1, a2, a3, a4, a5, b1, b2, b3, k1, k2, k3, kp); end