//Created by libAntimony v2.4 model stilianakis_1997__environment(time_) // Variable initializations: time_ = ; end model stilianakis_1997__X(X, time_, K, V) // Rate Rules: X' = lamda - (mu0 * X + K * V * X); // Variable initializations: X = 1.0E11; mu0 = 4.0E-3; lamda = 4.0E8; time_ = ; K = ; V = ; end model stilianakis_1997__Y(Y, time_, K, V, X, Z) // Rate Rules: Y' = K * V * Y - mu1 * (1.0 + a * Z) * Y; // Variable initializations: Y = 1.0; mu1 = 0.30; a = 1.0; time_ = ; K = ; V = ; X = ; Z = ; end model stilianakis_1997__CD4(X, Y) // Assignment Rules: CD4 := (X + Y) / 1.0E11; // Variable initializations: X = ; Y = ; end model stilianakis_1997__V(V, time_, Y, Z) // Rate Rules: V' = beta * Y - mu2 * (1.0 + b * Z) * V; // Variable initializations: V = 1.0; mu2 = 1.0; beta = 1.0E3; b = 1.0; time_ = ; Y = ; Z = ; end model stilianakis_1997__Z(Z, time_, V, X) // Assignment Rules: f_X := ((1.0 + C1) * power(X / X0, 2.0)) / (C1 + power(X / X0, 2.0)); g_V := V / (C2 + V); // Rate Rules: Z' = theta * g_V + rho * (f_X - Z) * Z; // Variable initializations: Z = 0.0; theta = 1.0E-6; rho = 0.50; C1 = 0.04; C2 = 1.0E3; X0 = 1.0E11; time_ = ; V = ; X = ; end model stilianakis_1997__K(K, time_, V) // Rate Rules: K' = omega * V * (Kmax - K); // Variable initializations: K = 1.35E-14; omega = 1.0E-16; Kmax = 20.0; time_ = ; V = ; end model *stilianakis_1997____main() // Sub-modules, and any changes to those submodules: environment: stilianakis_1997__environment(time_); X: stilianakis_1997__X(X0, time_, K0, V0); Y: stilianakis_1997__Y(Y0, time_, K0, V0, X0, Z0); CD4: stilianakis_1997__CD4(X0, Y0); V: stilianakis_1997__V(V0, time_, Y0, Z0); Z: stilianakis_1997__Z(Z0, time_, V0, X0); K: stilianakis_1997__K(K0, time_, V0); end