//Created by libAntimony v2.4 model kroll_model_2000__environment(time_) // Variable initializations: time_ = ; end model kroll_model_2000__Y(Y, C1, k1, P, K, C, C2, k2, ky, time_) // Rate Rules: Y' = k1 * C1 * (P / (K + P)) * C - (k2 * C2 * (1.0 - P / (K + P)) * Y + ky * Y); // Variable initializations: Y = 10.0; C1 = ; k1 = ; P = ; K = ; C = ; C2 = ; k2 = ; ky = ; time_ = ; end model kroll_model_2000__X(X, C2, k2, k3, P, Y, K, time_) // Rate Rules: X' = k2 * C2 * (1.0 - P / (K + P)) * Y - k3 * X; // Variable initializations: X = 500.0; C2 = ; k2 = ; k3 = ; P = ; Y = ; K = ; time_ = ; end model kroll_model_2000__P(P, k4, time_) // Rate Rules: P' = piecewise( 10.0 - k4 * P , ( geq(time_, 0.0)) && (time_ < 6.0 ), -(k4 * P) , ( geq(time_, 6.0)) && (time_ < 12.0 ), 10.0 - k4 * P , ( geq(time_, 12.0)) && (time_ < 18.0 ), -(k4 * P) , ( geq(time_, 18.0)) && (time_ < 24.0 ), 10.0 - k4 * P , ( geq(time_, 24.0)) && (time_ < 30.0 ), -(k4 * P) , ( geq(time_, 30.0)) && (time_ < 36.0 ), 10.0 - k4 * P , ( geq(time_, 36.0)) && (time_ < 42.0 ), -(k4 * P) , ( geq(time_, 42.0)) && (time_ < 48.0 ), 10.0 - k4 * P , ( geq(time_, 48.0)) && (time_ < 54.0 ), -(k4 * P) , ( geq(time_, 54.0)) && (time_ < 60.0 ), 10.0 - k4 * P , ( geq(time_, 60.0)) && (time_ < 66.0 ), -(k4 * P) , ( geq(time_, 66.0)) && (time_ < 72.0 ), 10.0 - k4 * P , ( geq(time_, 72.0)) && (time_ < 78.0 ), -(k4 * P) , ( geq(time_, 78.0)) && (time_ < 84.0 ), 10.0 - k4 * P , ( geq(time_, 84.0)) && (time_ < 90.0 ), -(k4 * P) , ( geq(time_, 90.0)) && (time_ < 96.0 ) ); // Variable initializations: P = 10.0; k4 = ; time_ = ; end model kroll_model_2000__Z(Z, C3, k5, K2, IL6, k6, time_) // Rate Rules: Z' = k5 * C3 * (IL6 / (K2 + IL6)) - k6 * Z; // Variable initializations: Z = 200.0; C3 = ; k5 = ; K2 = ; IL6 = ; k6 = ; time_ = ; end model kroll_model_2000__IL6(IL6, X, time_) // Rate Rules: IL6' = 0.1 * X - 10.0 * IL6; // Variable initializations: IL6 = 1.9; X = ; time_ = ; end model kroll_model_2000__model_parameters(C, k1, k2, k3, k4, k5, k6, ky, K, K2, C1, C2, C3) // Variable initializations: C = 50.0; k1 = 1.0; k2 = 1.3; k3 = 0.05; k4 = 0.9; k5 = 5.0; k6 = 0.02; ky = 0.01; K = 5.0; K2 = 2.0; C1 = 1.0; C2 = 1.0; C3 = 1.0; end model *kroll_model_2000____main() // Sub-modules, and any changes to those submodules: environment: kroll_model_2000__environment(time_); Y: kroll_model_2000__Y(Y0, C1, k1, P0, K, C, C2, k2, ky, time_); X: kroll_model_2000__X(X0, C2, k2, k3, P0, Y0, K, time_); P: kroll_model_2000__P(P0, k4, time_); Z: kroll_model_2000__Z(Z0, C3, k5, K2, IL60, k6, time_); IL6: kroll_model_2000__IL6(IL60, X0, time_); model_parameters: kroll_model_2000__model_parameters(C, k1, k2, k3, k4, k5, k6, ky, K, K2, C1, C2, C3); end