Generated Code
The following is python code generated by the CellML API from this CellML file. (Back to language selection)
The raw code is available.
# Size of variable arrays:
sizeAlgebraic = 0
sizeStates = 66
sizeConstants = 255
from math import *
from numpy import *
def createLegends():
legend_states = [""] * sizeStates
legend_rates = [""] * sizeStates
legend_algebraic = [""] * sizeAlgebraic
legend_voi = ""
legend_constants = [""] * sizeConstants
legend_voi = "time in component environment (hour)"
legend_states[0] = "AP1 in component AP1 (micromolar)"
legend_constants[0] = "v_sap1 in component model_parameters (flux)"
legend_constants[1] = "GF in component model_parameters (micromolar)"
legend_constants[2] = "K_agf in component model_parameters (micromolar)"
legend_constants[3] = "k_dap1 in component model_parameters (first_order_rate_constant)"
legend_constants[4] = "eps in component model_parameters (dimensionless)"
legend_states[1] = "pRB in component pRB (micromolar)"
legend_constants[5] = "v_sprb in component model_parameters (flux)"
legend_constants[6] = "k_pc1 in component model_parameters (second_order_rate_constant)"
legend_states[2] = "E2F in component E2F (micromolar)"
legend_constants[7] = "k_pc2 in component model_parameters (first_order_rate_constant)"
legend_states[3] = "pRBc1 in component pRBc1 (micromolar)"
legend_constants[8] = "V_1 in component model_parameters (first_order_rate_constant)"
legend_constants[9] = "K_1 in component model_parameters (micromolar)"
legend_states[4] = "Md in component Md (micromolar)"
legend_states[5] = "Mdp27 in component Mdp27 (micromolar)"
legend_constants[10] = "V_2 in component model_parameters (flux)"
legend_states[6] = "pRBp in component pRBp (micromolar)"
legend_constants[11] = "K_2 in component model_parameters (micromolar)"
legend_constants[12] = "k_dprb in component model_parameters (first_order_rate_constant)"
legend_constants[13] = "V_3 in component model_parameters (first_order_rate_constant)"
legend_constants[14] = "K_3 in component model_parameters (micromolar)"
legend_states[7] = "Me in component Me (micromolar)"
legend_constants[15] = "V_4 in component model_parameters (flux)"
legend_states[8] = "pRBpp in component pRBpp (micromolar)"
legend_constants[16] = "K_4 in component model_parameters (micromolar)"
legend_constants[17] = "k_pc3 in component model_parameters (second_order_rate_constant)"
legend_constants[18] = "k_pc4 in component model_parameters (first_order_rate_constant)"
legend_states[9] = "pRBc2 in component pRBc2 (micromolar)"
legend_constants[19] = "k_dpRBp in component model_parameters (first_order_rate_constant)"
legend_constants[20] = "k_dpRBpp in component model_parameters (first_order_rate_constant)"
legend_constants[21] = "v_se2f in component model_parameters (flux)"
legend_constants[22] = "V_1e2f in component model_parameters (first_order_rate_constant)"
legend_states[10] = "Ma in component Ma (micromolar)"
legend_constants[23] = "K_1e2f in component model_parameters (micromolar)"
legend_constants[24] = "V_2e2f in component model_parameters (flux)"
legend_states[11] = "E2Fp in component E2Fp (micromolar)"
legend_constants[25] = "K_2e2f in component model_parameters (micromolar)"
legend_constants[26] = "k_de2f in component model_parameters (first_order_rate_constant)"
legend_constants[27] = "k_de2fp in component model_parameters (first_order_rate_constant)"
legend_states[12] = "Cd in component Cd (micromolar)"
legend_constants[28] = "k_cd1 in component model_parameters (first_order_rate_constant)"
legend_constants[29] = "k_cd2 in component model_parameters (first_order_rate_constant)"
legend_constants[30] = "K_i7 in component model_parameters (micromolar)"
legend_constants[31] = "K_i8 in component model_parameters (micromolar)"
legend_constants[32] = "k_com1 in component model_parameters (second_order_rate_constant)"
legend_constants[33] = "Cdk4_tot in component model_parameters (micromolar)"
legend_states[13] = "Mdi in component Mdi (micromolar)"
legend_constants[34] = "k_decom1 in component model_parameters (first_order_rate_constant)"
legend_constants[35] = "V_dd in component model_parameters (flux)"
legend_constants[36] = "K_dd in component model_parameters (micromolar)"
legend_constants[37] = "k_ddd in component model_parameters (first_order_rate_constant)"
legend_constants[38] = "V_m2d in component model_parameters (flux)"
legend_constants[39] = "K_2d in component model_parameters (micromolar)"
legend_constants[40] = "V_m1d in component model_parameters (flux)"
legend_constants[41] = "K_1d in component model_parameters (micromolar)"
legend_constants[42] = "k_c1 in component model_parameters (second_order_rate_constant)"
legend_states[14] = "p27 in component p27 (micromolar)"
legend_constants[43] = "k_c2 in component model_parameters (first_order_rate_constant)"
legend_states[15] = "Ce in component Ce (micromolar)"
legend_constants[44] = "k_ce in component model_parameters (first_order_rate_constant)"
legend_constants[45] = "K_i9 in component model_parameters (micromolar)"
legend_constants[46] = "K_i10 in component model_parameters (micromolar)"
legend_constants[47] = "k_com2 in component model_parameters (second_order_rate_constant)"
legend_constants[48] = "Cdk2_tot in component model_parameters (micromolar)"
legend_states[16] = "Mei in component Mei (micromolar)"
legend_states[17] = "Mep27 in component Mep27 (micromolar)"
legend_states[18] = "Mai in component Mai (micromolar)"
legend_states[19] = "Map27 in component Map27 (micromolar)"
legend_constants[49] = "k_decom2 in component model_parameters (first_order_rate_constant)"
legend_constants[50] = "V_de in component model_parameters (flux)"
legend_states[20] = "Skp2 in component Skp2 (micromolar)"
legend_constants[51] = "K_dceskp2 in component model_parameters (micromolar)"
legend_constants[52] = "K_de in component model_parameters (micromolar)"
legend_constants[53] = "k_dde in component model_parameters (first_order_rate_constant)"
legend_constants[54] = "V_m2e in component model_parameters (first_order_rate_constant)"
legend_states[21] = "Wee1 in component Wee1 (micromolar)"
legend_constants[55] = "i_b1 in component model_parameters (micromolar)"
legend_constants[56] = "K_2e in component model_parameters (micromolar)"
legend_constants[57] = "V_m1e in component model_parameters (first_order_rate_constant)"
legend_states[22] = "Pe in component Pe (micromolar)"
legend_constants[58] = "K_1e in component model_parameters (micromolar)"
legend_constants[59] = "k_c3 in component model_parameters (second_order_rate_constant)"
legend_constants[60] = "k_c4 in component model_parameters (first_order_rate_constant)"
legend_constants[61] = "v_sskp2 in component model_parameters (flux)"
legend_constants[62] = "V_dskp2 in component model_parameters (flux)"
legend_constants[63] = "K_dskp2 in component model_parameters (micromolar)"
legend_states[23] = "Cdh1a in component Cdh1a (micromolar)"
legend_constants[64] = "K_cdh1 in component model_parameters (micromolar)"
legend_constants[65] = "k_ddskp2 in component model_parameters (first_order_rate_constant)"
legend_states[24] = "Pei in component Pei (micromolar)"
legend_constants[66] = "v_spei in component model_parameters (flux)"
legend_constants[67] = "V_6e in component model_parameters (flux)"
legend_constants[68] = "x_e1 in component model_parameters (dimensionless)"
legend_constants[69] = "x_e2 in component model_parameters (per_micromolar)"
legend_states[25] = "Chk1 in component Chk1 (micromolar)"
legend_constants[70] = "K_6e in component model_parameters (micromolar)"
legend_constants[71] = "V_m5e in component model_parameters (first_order_rate_constant)"
legend_constants[72] = "a_e in component model_parameters (micromolar)"
legend_constants[73] = "K_5e in component model_parameters (micromolar)"
legend_constants[74] = "k_dpei in component model_parameters (first_order_rate_constant)"
legend_constants[75] = "k_dpe in component model_parameters (first_order_rate_constant)"
legend_states[26] = "Ca in component Ca (micromolar)"
legend_constants[76] = "k_ca in component model_parameters (first_order_rate_constant)"
legend_constants[77] = "K_i11 in component model_parameters (micromolar)"
legend_constants[78] = "K_i12 in component model_parameters (micromolar)"
legend_constants[79] = "k_com3 in component model_parameters (second_order_rate_constant)"
legend_constants[80] = "k_decom3 in component model_parameters (first_order_rate_constant)"
legend_constants[81] = "V_da in component model_parameters (flux)"
legend_constants[82] = "K_da in component model_parameters (micromolar)"
legend_states[27] = "Cdc20a in component Cdc20a (micromolar)"
legend_constants[83] = "K_acdc20 in component model_parameters (micromolar)"
legend_constants[84] = "k_dda in component model_parameters (first_order_rate_constant)"
legend_constants[85] = "V_m2a in component model_parameters (first_order_rate_constant)"
legend_constants[86] = "i_b2 in component model_parameters (micromolar)"
legend_constants[87] = "K_2a in component model_parameters (micromolar)"
legend_constants[88] = "V_m1a in component model_parameters (first_order_rate_constant)"
legend_states[28] = "Pa in component Pa (micromolar)"
legend_constants[89] = "K_1a in component model_parameters (micromolar)"
legend_constants[90] = "k_c5 in component model_parameters (second_order_rate_constant)"
legend_constants[91] = "k_c6 in component model_parameters (first_order_rate_constant)"
legend_constants[92] = "v_s1p27 in component model_parameters (flux)"
legend_constants[93] = "v_s2p27 in component model_parameters (first_order_rate_constant)"
legend_constants[94] = "K_i13 in component model_parameters (micromolar)"
legend_constants[95] = "K_i14 in component model_parameters (micromolar)"
legend_constants[96] = "k_c7 in component model_parameters (second_order_rate_constant)"
legend_states[29] = "Mb in component Mb (micromolar)"
legend_constants[97] = "k_c8 in component model_parameters (first_order_rate_constant)"
legend_states[30] = "Mbp27 in component Mbp27 (micromolar)"
legend_constants[98] = "V_1p27 in component model_parameters (first_order_rate_constant)"
legend_constants[99] = "K_1p27 in component model_parameters (micromolar)"
legend_constants[100] = "V_2p27 in component model_parameters (flux)"
legend_constants[101] = "K_2p27 in component model_parameters (micromolar)"
legend_states[31] = "p27p in component p27p (micromolar)"
legend_constants[102] = "k_ddp27 in component model_parameters (first_order_rate_constant)"
legend_constants[103] = "V_dp27p in component model_parameters (flux)"
legend_constants[104] = "K_dp27skp2 in component model_parameters (micromolar)"
legend_constants[105] = "K_dp27p in component model_parameters (micromolar)"
legend_constants[106] = "k_ddp27p in component model_parameters (first_order_rate_constant)"
legend_states[32] = "Cdh1i in component Cdh1i (micromolar)"
legend_constants[107] = "V_2cdh1 in component model_parameters (first_order_rate_constant)"
legend_constants[108] = "K_2cdh1 in component model_parameters (micromolar)"
legend_constants[109] = "V_1cdh1 in component model_parameters (flux)"
legend_constants[110] = "K_1cdh1 in component model_parameters (micromolar)"
legend_constants[111] = "k_dcdh1i in component model_parameters (first_order_rate_constant)"
legend_constants[112] = "v_scdh1a in component model_parameters (flux)"
legend_constants[113] = "k_dcdh1a in component model_parameters (first_order_rate_constant)"
legend_states[33] = "Pai in component Pai (micromolar)"
legend_constants[114] = "v_spai in component model_parameters (flux)"
legend_constants[115] = "V_6a in component model_parameters (flux)"
legend_constants[116] = "x_a1 in component model_parameters (dimensionless)"
legend_constants[117] = "x_a2 in component model_parameters (per_micromolar)"
legend_constants[118] = "K_6a in component model_parameters (micromolar)"
legend_constants[119] = "V_m5a in component model_parameters (first_order_rate_constant)"
legend_constants[120] = "a_a in component model_parameters (micromolar)"
legend_constants[121] = "K_5a in component model_parameters (micromolar)"
legend_constants[122] = "k_dpai in component model_parameters (first_order_rate_constant)"
legend_constants[123] = "k_dpa in component model_parameters (first_order_rate_constant)"
legend_states[34] = "Cb in component Cb (micromolar)"
legend_constants[124] = "v_cb in component model_parameters (flux)"
legend_constants[125] = "k_com4 in component model_parameters (second_order_rate_constant)"
legend_constants[126] = "Cdk1_tot in component model_parameters (micromolar)"
legend_states[35] = "Mbi in component Mbi (micromolar)"
legend_constants[127] = "k_decom4 in component model_parameters (first_order_rate_constant)"
legend_constants[128] = "V_db in component model_parameters (flux)"
legend_constants[129] = "K_db in component model_parameters (micromolar)"
legend_constants[130] = "K_dbcdc20 in component model_parameters (micromolar)"
legend_constants[131] = "K_dbcdh1 in component model_parameters (micromolar)"
legend_constants[132] = "k_ddb in component model_parameters (first_order_rate_constant)"
legend_constants[133] = "V_m2b in component model_parameters (first_order_rate_constant)"
legend_constants[134] = "i_b3 in component model_parameters (micromolar)"
legend_constants[135] = "K_2b in component model_parameters (micromolar)"
legend_constants[136] = "V_m1b in component model_parameters (first_order_rate_constant)"
legend_states[36] = "Pb in component Pb (micromolar)"
legend_constants[137] = "K_1b in component model_parameters (micromolar)"
legend_states[37] = "Cdc20i in component Cdc20i (micromolar)"
legend_constants[138] = "v_scdc20i in component model_parameters (flux)"
legend_constants[139] = "V_m3b in component model_parameters (first_order_rate_constant)"
legend_constants[140] = "K_3b in component model_parameters (micromolar)"
legend_constants[141] = "V_m4b in component model_parameters (flux)"
legend_constants[142] = "K_4b in component model_parameters (micromolar)"
legend_constants[143] = "k_dcdc20i in component model_parameters (first_order_rate_constant)"
legend_constants[144] = "k_dcdc20a in component model_parameters (first_order_rate_constant)"
legend_states[38] = "Pbi in component Pbi (micromolar)"
legend_constants[145] = "v_spbi in component model_parameters (flux)"
legend_constants[146] = "V_6b in component model_parameters (flux)"
legend_constants[147] = "x_b1 in component model_parameters (dimensionless)"
legend_constants[148] = "x_b2 in component model_parameters (per_micromolar)"
legend_constants[149] = "K_6b in component model_parameters (micromolar)"
legend_constants[150] = "V_m5b in component model_parameters (first_order_rate_constant)"
legend_constants[151] = "a_b in component model_parameters (micromolar)"
legend_constants[152] = "K_5b in component model_parameters (micromolar)"
legend_constants[153] = "k_dpbi in component model_parameters (first_order_rate_constant)"
legend_constants[154] = "k_dpb in component model_parameters (first_order_rate_constant)"
legend_constants[155] = "v_swee1 in component model_parameters (flux)"
legend_constants[156] = "k_sw in component model_parameters (first_order_rate_constant)"
legend_states[39] = "Mw in component Mw (micromolar)"
legend_constants[157] = "V_m7b in component model_parameters (first_order_rate_constant)"
legend_constants[158] = "i_b in component model_parameters (micromolar)"
legend_constants[159] = "K_7b in component model_parameters (micromolar)"
legend_constants[160] = "V_m8b in component model_parameters (flux)"
legend_states[40] = "Wee1p in component Wee1p (micromolar)"
legend_constants[161] = "K_8b in component model_parameters (micromolar)"
legend_constants[162] = "k_dwee1 in component model_parameters (first_order_rate_constant)"
legend_constants[163] = "k_dwee1p in component model_parameters (first_order_rate_constant)"
legend_states[41] = "Cdc45 in component Cdc45 (micromolar)"
legend_constants[164] = "V_1cdc45 in component model_parameters (first_order_rate_constant)"
legend_constants[165] = "Cdc45_tot in component model_parameters (micromolar)"
legend_constants[166] = "K_1cdc45 in component model_parameters (micromolar)"
legend_constants[167] = "V_2cdc45 in component model_parameters (flux)"
legend_constants[168] = "K_2cdc45 in component model_parameters (micromolar)"
legend_constants[169] = "k_spol in component model_parameters (second_order_rate_constant)"
legend_constants[170] = "Pol_tot in component model_parameters (micromolar)"
legend_states[42] = "Pol in component Pol (micromolar)"
legend_constants[171] = "k_dpol in component model_parameters (first_order_rate_constant)"
legend_states[43] = "Primer in component Primer (micromolar)"
legend_constants[172] = "k_sprim in component model_parameters (first_order_rate_constant)"
legend_constants[173] = "k_dprim in component model_parameters (first_order_rate_constant)"
legend_constants[174] = "k_aatr in component model_parameters (second_order_rate_constant)"
legend_constants[175] = "ATR_tot in component model_parameters (micromolar)"
legend_states[44] = "ATR in component ATR (micromolar)"
legend_constants[176] = "k_datr in component model_parameters (first_order_rate_constant)"
legend_constants[177] = "V_1chk in component model_parameters (first_order_rate_constant)"
legend_constants[178] = "Chk1_tot in component model_parameters (micromolar)"
legend_constants[179] = "K_1chk in component model_parameters (micromolar)"
legend_constants[180] = "V_2chk in component model_parameters (flux)"
legend_constants[181] = "K_2chk in component model_parameters (micromolar)"
legend_constants[182] = "v_sw in component model_parameters (flux)"
legend_states[45] = "BN in component BN (nanomolar)"
legend_constants[183] = "n_gerard in component model_parameters (dimensionless)"
legend_constants[184] = "K_iw in component model_parameters (nanomolar)"
legend_constants[185] = "v_dw in component model_parameters (flux)"
legend_constants[186] = "K_dw in component model_parameters (micromolar)"
legend_states[46] = "X in component X (micromolar)"
legend_constants[187] = "V_1x in component model_parameters (first_order_rate_constant)"
legend_constants[188] = "X_tot in component model_parameters (micromolar)"
legend_constants[189] = "K_1x in component model_parameters (micromolar)"
legend_constants[190] = "V_2x in component model_parameters (flux)"
legend_constants[191] = "K_2x in component model_parameters (micromolar)"
legend_states[47] = "CbA in component CbA (micromolar)"
legend_states[48] = "MP in component MP (nanomolar)"
legend_constants[192] = "vsP in component MP (nano_flux)"
legend_constants[193] = "vmP in component MP (nano_flux)"
legend_constants[194] = "kdmp in component MP (first_order_rate_constant_nano)"
legend_constants[195] = "KAP in component MP (nanomolar)"
legend_constants[196] = "KmP in component MP (nanomolar)"
legend_constants[197] = "n in component model_parameters (dimensionless)"
legend_states[49] = "MC in component MC (nanomolar)"
legend_constants[198] = "vsC in component MC (nano_flux)"
legend_constants[199] = "vmC in component MC (nano_flux)"
legend_constants[200] = "kdmc in component MC (first_order_rate_constant_nano)"
legend_constants[201] = "KAC in component MC (nanomolar)"
legend_constants[202] = "KmC in component MC (nanomolar)"
legend_states[50] = "MB in component MB (nanomolar)"
legend_constants[203] = "vsB in component MB (nano_flux)"
legend_constants[204] = "vmB in component MB (nano_flux)"
legend_constants[205] = "kdmb in component MB (first_order_rate_constant_nano)"
legend_constants[206] = "KIB in component MB (nanomolar)"
legend_constants[207] = "KmB in component MB (nanomolar)"
legend_constants[208] = "m in component model_parameters (dimensionless)"
legend_states[51] = "RN in component RN (nanomolar)"
legend_states[52] = "MR in component MR (nanomolar)"
legend_constants[209] = "vsR in component MR (nano_flux)"
legend_constants[210] = "vmR in component MR (nano_flux)"
legend_constants[211] = "kdmr in component MR (first_order_rate_constant_nano)"
legend_constants[212] = "KAR in component MR (nanomolar)"
legend_constants[213] = "KmR in component MR (nanomolar)"
legend_constants[214] = "h in component model_parameters (dimensionless)"
legend_states[53] = "PC in component PC (nanomolar)"
legend_constants[215] = "ksP in component model_parameters (first_order_rate_constant_nano)"
legend_constants[216] = "Kp in component model_parameters (nanomolar)"
legend_constants[217] = "Kdp in component model_parameters (nanomolar)"
legend_constants[218] = "k3 in component model_parameters (second_order_rate_constant_nano)"
legend_constants[219] = "k4 in component model_parameters (first_order_rate_constant_nano)"
legend_constants[220] = "kdn in component model_parameters (first_order_rate_constant_nano)"
legend_constants[221] = "V1P in component model_parameters (nano_flux)"
legend_constants[222] = "V2P in component model_parameters (nano_flux)"
legend_states[54] = "PCP in component PCP (nanomolar)"
legend_states[55] = "PCC in component PCC (nanomolar)"
legend_states[56] = "CC in component CC (nanomolar)"
legend_constants[223] = "ksC in component model_parameters (first_order_rate_constant_nano)"
legend_constants[224] = "kdnc in component model_parameters (first_order_rate_constant_nano)"
legend_constants[225] = "V1C in component model_parameters (nano_flux)"
legend_constants[226] = "V2C in component model_parameters (nano_flux)"
legend_states[57] = "CCP in component CCP (nanomolar)"
legend_states[58] = "RC in component RC (nanomolar)"
legend_constants[227] = "ksR in component model_parameters (first_order_rate_constant_nano)"
legend_constants[228] = "Kd in component model_parameters (nanomolar)"
legend_constants[229] = "k9 in component model_parameters (first_order_rate_constant_nano)"
legend_constants[230] = "k10 in component model_parameters (first_order_rate_constant_nano)"
legend_constants[231] = "vdRC in component model_parameters (nano_flux)"
legend_constants[232] = "vdPC in component model_parameters (nano_flux)"
legend_constants[233] = "vdCC in component model_parameters (nano_flux)"
legend_constants[234] = "k1 in component model_parameters (first_order_rate_constant_nano)"
legend_constants[235] = "k2 in component model_parameters (first_order_rate_constant_nano)"
legend_constants[236] = "V1PC in component model_parameters (nano_flux)"
legend_constants[237] = "V2PC in component model_parameters (nano_flux)"
legend_states[59] = "PCCP in component PCCP (nanomolar)"
legend_states[60] = "PCN in component PCN (nanomolar)"
legend_constants[238] = "k7 in component model_parameters (second_order_rate_constant_nano)"
legend_constants[239] = "k8 in component model_parameters (first_order_rate_constant_nano)"
legend_constants[240] = "V3PC in component model_parameters (nano_flux)"
legend_constants[241] = "V4PC in component model_parameters (nano_flux)"
legend_states[61] = "PCNP in component PCNP (nanomolar)"
legend_states[62] = "IN in component IN (nanomolar)"
legend_constants[242] = "vdRN in component model_parameters (nano_flux)"
legend_constants[243] = "vdPCC in component model_parameters (nano_flux)"
legend_constants[244] = "vdPCN in component model_parameters (nano_flux)"
legend_states[63] = "BC in component BC (nanomolar)"
legend_constants[245] = "ksB in component model_parameters (first_order_rate_constant_nano)"
legend_constants[246] = "k5 in component model_parameters (first_order_rate_constant_nano)"
legend_constants[247] = "k6 in component model_parameters (first_order_rate_constant_nano)"
legend_constants[248] = "V1B in component model_parameters (nano_flux)"
legend_constants[249] = "V2B in component model_parameters (nano_flux)"
legend_states[64] = "BCP in component BCP (nanomolar)"
legend_constants[250] = "vdBC in component model_parameters (nano_flux)"
legend_constants[251] = "V3B in component model_parameters (nano_flux)"
legend_constants[252] = "V4B in component model_parameters (nano_flux)"
legend_states[65] = "BNP in component BNP (nanomolar)"
legend_constants[253] = "vdBN in component model_parameters (nano_flux)"
legend_constants[254] = "vdIN in component model_parameters (nano_flux)"
legend_rates[0] = "d/dt AP1 in component AP1 (micromolar)"
legend_rates[1] = "d/dt pRB in component pRB (micromolar)"
legend_rates[3] = "d/dt pRBc1 in component pRBc1 (micromolar)"
legend_rates[6] = "d/dt pRBp in component pRBp (micromolar)"
legend_rates[9] = "d/dt pRBc2 in component pRBc2 (micromolar)"
legend_rates[8] = "d/dt pRBpp in component pRBpp (micromolar)"
legend_rates[2] = "d/dt E2F in component E2F (micromolar)"
legend_rates[11] = "d/dt E2Fp in component E2Fp (micromolar)"
legend_rates[12] = "d/dt Cd in component Cd (micromolar)"
legend_rates[13] = "d/dt Mdi in component Mdi (micromolar)"
legend_rates[4] = "d/dt Md in component Md (micromolar)"
legend_rates[5] = "d/dt Mdp27 in component Mdp27 (micromolar)"
legend_rates[15] = "d/dt Ce in component Ce (micromolar)"
legend_rates[16] = "d/dt Mei in component Mei (micromolar)"
legend_rates[7] = "d/dt Me in component Me (micromolar)"
legend_rates[20] = "d/dt Skp2 in component Skp2 (micromolar)"
legend_rates[17] = "d/dt Mep27 in component Mep27 (micromolar)"
legend_rates[24] = "d/dt Pei in component Pei (micromolar)"
legend_rates[22] = "d/dt Pe in component Pe (micromolar)"
legend_rates[26] = "d/dt Ca in component Ca (micromolar)"
legend_rates[18] = "d/dt Mai in component Mai (micromolar)"
legend_rates[10] = "d/dt Ma in component Ma (micromolar)"
legend_rates[19] = "d/dt Map27 in component Map27 (micromolar)"
legend_rates[14] = "d/dt p27 in component p27 (micromolar)"
legend_rates[31] = "d/dt p27p in component p27p (micromolar)"
legend_rates[32] = "d/dt Cdh1i in component Cdh1i (micromolar)"
legend_rates[23] = "d/dt Cdh1a in component Cdh1a (micromolar)"
legend_rates[33] = "d/dt Pai in component Pai (micromolar)"
legend_rates[28] = "d/dt Pa in component Pa (micromolar)"
legend_rates[34] = "d/dt Cb in component Cb (micromolar)"
legend_rates[35] = "d/dt Mbi in component Mbi (micromolar)"
legend_rates[29] = "d/dt Mb in component Mb (micromolar)"
legend_rates[30] = "d/dt Mbp27 in component Mbp27 (micromolar)"
legend_rates[37] = "d/dt Cdc20i in component Cdc20i (micromolar)"
legend_rates[27] = "d/dt Cdc20a in component Cdc20a (micromolar)"
legend_rates[38] = "d/dt Pbi in component Pbi (micromolar)"
legend_rates[36] = "d/dt Pb in component Pb (micromolar)"
legend_rates[21] = "d/dt Wee1 in component Wee1 (micromolar)"
legend_rates[40] = "d/dt Wee1p in component Wee1p (micromolar)"
legend_rates[41] = "d/dt Cdc45 in component Cdc45 (micromolar)"
legend_rates[42] = "d/dt Pol in component Pol (micromolar)"
legend_rates[43] = "d/dt Primer in component Primer (micromolar)"
legend_rates[44] = "d/dt ATR in component ATR (micromolar)"
legend_rates[25] = "d/dt Chk1 in component Chk1 (micromolar)"
legend_rates[39] = "d/dt Mw in component Mw (micromolar)"
legend_rates[46] = "d/dt X in component X (micromolar)"
legend_rates[47] = "d/dt CbA in component CbA (micromolar)"
legend_rates[48] = "d/dt MP in component MP (nanomolar)"
legend_rates[49] = "d/dt MC in component MC (nanomolar)"
legend_rates[50] = "d/dt MB in component MB (nanomolar)"
legend_rates[52] = "d/dt MR in component MR (nanomolar)"
legend_rates[53] = "d/dt PC in component PC (nanomolar)"
legend_rates[56] = "d/dt CC in component CC (nanomolar)"
legend_rates[58] = "d/dt RC in component RC (nanomolar)"
legend_rates[54] = "d/dt PCP in component PCP (nanomolar)"
legend_rates[57] = "d/dt CCP in component CCP (nanomolar)"
legend_rates[55] = "d/dt PCC in component PCC (nanomolar)"
legend_rates[60] = "d/dt PCN in component PCN (nanomolar)"
legend_rates[51] = "d/dt RN in component RN (nanomolar)"
legend_rates[59] = "d/dt PCCP in component PCCP (nanomolar)"
legend_rates[61] = "d/dt PCNP in component PCNP (nanomolar)"
legend_rates[63] = "d/dt BC in component BC (nanomolar)"
legend_rates[64] = "d/dt BCP in component BCP (nanomolar)"
legend_rates[45] = "d/dt BN in component BN (nanomolar)"
legend_rates[65] = "d/dt BNP in component BNP (nanomolar)"
legend_rates[62] = "d/dt IN in component IN (nanomolar)"
return (legend_states, legend_algebraic, legend_voi, legend_constants)
def initConsts():
constants = [0.0] * sizeConstants; states = [0.0] * sizeStates;
states[0] = 0.01
constants[0] = 1
constants[1] = 1
constants[2] = 0.1
constants[3] = 0.15
constants[4] = 17
states[1] = 1
constants[5] = 0.8
constants[6] = 0.05
states[2] = 0.01
constants[7] = 0.5
states[3] = 0.1
constants[8] = 2.2
constants[9] = 0.1
states[4] = 0.01
states[5] = 0.01
constants[10] = 2
states[6] = 0.25
constants[11] = 0.1
constants[12] = 0.01
constants[13] = 1
constants[14] = 0.1
states[7] = 0.01
constants[15] = 2
states[8] = 0.1
constants[16] = 0.1
constants[17] = 0.025
constants[18] = 0.5
states[9] = 0.05
constants[19] = 0.06
constants[20] = 0.04
constants[21] = 0.17
constants[22] = 4
states[10] = 0.01
constants[23] = 5
constants[24] = 0.75
states[11] = 0.01
constants[25] = 5
constants[26] = 0.002
constants[27] = 1.1
states[12] = 0.01
constants[28] = 0.4
constants[29] = 0.005
constants[30] = 0.1
constants[31] = 2
constants[32] = 0.175
constants[33] = 1.5
states[13] = 0.01
constants[34] = 0.1
constants[35] = 5
constants[36] = 0.1
constants[37] = 0.005
constants[38] = 0.2
constants[39] = 0.1
constants[40] = 1
constants[41] = 0.1
constants[42] = 0.15
states[14] = 0.01
constants[43] = 0.05
states[15] = 0.01
constants[44] = 0.29
constants[45] = 0.1
constants[46] = 2
constants[47] = 0.2
constants[48] = 2
states[16] = 0.01
states[17] = 0.01
states[18] = 0.01
states[19] = 0.01
constants[49] = 0.1
constants[50] = 3
states[20] = 0.01
constants[51] = 2
constants[52] = 0.1
constants[53] = 0.005
constants[54] = 1.4
states[21] = 0.1
constants[55] = 0.5
constants[56] = 0.1
constants[57] = 2
states[22] = 0.01
constants[58] = 0.1
constants[59] = 0.2
constants[60] = 0.1
constants[61] = 0.15
constants[62] = 1.1
constants[63] = 0.5
states[23] = 0.01
constants[64] = 0.4
constants[65] = 0.005
states[24] = 0.01
constants[66] = 0.13
constants[67] = 0.8
constants[68] = 1
constants[69] = 1
states[25] = 0.01
constants[70] = 0.1
constants[71] = 5
constants[72] = 0.25
constants[73] = 0.1
constants[74] = 0.15
constants[75] = 0.075
states[26] = 0.01
constants[76] = 0.0375
constants[77] = 0.1
constants[78] = 2
constants[79] = 0.2
constants[80] = 0.1
constants[81] = 2.5
constants[82] = 1.1
states[27] = 0.01
constants[83] = 2
constants[84] = 0.005
constants[85] = 1.85
constants[86] = 0.5
constants[87] = 0.1
constants[88] = 2
states[28] = 0.01
constants[89] = 0.1
constants[90] = 0.15
constants[91] = 0.125
constants[92] = 0.8
constants[93] = 0.1
constants[94] = 0.1
constants[95] = 2
constants[96] = 0.12
states[29] = 0.01
constants[97] = 0.2
states[30] = 0.01
constants[98] = 100
constants[99] = 0.5
constants[100] = 0.1
constants[101] = 0.5
states[31] = 0.01
constants[102] = 0.06
constants[103] = 5
constants[104] = 0.1
constants[105] = 0.1
constants[106] = 0.01
states[32] = 0.01
constants[107] = 8
constants[108] = 0.01
constants[109] = 1.25
constants[110] = 0.01
constants[111] = 0.2
constants[112] = 0.11
constants[113] = 0.1
states[33] = 0.01
constants[114] = 0.105
constants[115] = 1
constants[116] = 1
constants[117] = 1
constants[118] = 0.1
constants[119] = 4
constants[120] = 0.2
constants[121] = 0.1
constants[122] = 0.15
constants[123] = 0.075
states[34] = 0.01
constants[124] = 0.05
constants[125] = 0.25
constants[126] = 0.5
states[35] = 0.01
constants[127] = 0.1
constants[128] = 0.06
constants[129] = 0.005
constants[130] = 0.2
constants[131] = 0.1
constants[132] = 0.005
constants[133] = 2.1
constants[134] = 0.5
constants[135] = 0.1
constants[136] = 3.9
states[36] = 0.01
constants[137] = 0.1
states[37] = 0.01
constants[138] = 0.1
constants[139] = 8
constants[140] = 0.1
constants[141] = 0.7
constants[142] = 0.1
constants[143] = 0.14
constants[144] = 0.05
states[38] = 0.01
constants[145] = 0.12
constants[146] = 1
constants[147] = 1
constants[148] = 1
constants[149] = 0.1
constants[150] = 5
constants[151] = 0.11
constants[152] = 0.1
constants[153] = 0.2
constants[154] = 0.1
constants[155] = 0.06
constants[156] = 5
states[39] = 0
constants[157] = 1.2
constants[158] = 0.75
constants[159] = 0.1
constants[160] = 1
states[40] = 0.01
constants[161] = 0.1
constants[162] = 0.1
constants[163] = 0.2
states[41] = 0.01
constants[164] = 0.8
constants[165] = 0.5
constants[166] = 0.02
constants[167] = 0.12
constants[168] = 0.02
constants[169] = 0.8
constants[170] = 0.5
states[42] = 0.01
constants[171] = 0.2
states[43] = 0.01
constants[172] = 0.05
constants[173] = 0.15
constants[174] = 0.022
constants[175] = 0.5
states[44] = 0.01
constants[176] = 0.15
constants[177] = 4
constants[178] = 0.5
constants[179] = 0.5
constants[180] = 0.1
constants[181] = 0.5
constants[182] = 0
states[45] = 0.1
constants[183] = 4
constants[184] = 0.5
constants[185] = 0.12
constants[186] = 0.5
states[46] = 0.01
constants[187] = 10
constants[188] = 1
constants[189] = 0.1
constants[190] = 2
constants[191] = 0.1
states[47] = 0.01
states[48] = 0.1
constants[192] = 2.4
constants[193] = 2.2
constants[194] = 0.02
constants[195] = 0.6
constants[196] = 0.3
constants[197] = 2
states[49] = 1.2
constants[198] = 2.2
constants[199] = 2
constants[200] = 0.02
constants[201] = 0.6
constants[202] = 0.4
states[50] = 9
constants[203] = 1.8
constants[204] = 1.3
constants[205] = 0.02
constants[206] = 2.2
constants[207] = 0.4
constants[208] = 2
states[51] = 0.1
states[52] = 1.5
constants[209] = 1.6
constants[210] = 1.6
constants[211] = 0.02
constants[212] = 0.6
constants[213] = 0.4
constants[214] = 2
states[53] = 0.1
constants[215] = 1.2
constants[216] = 1.006
constants[217] = 0.1
constants[218] = 0.8
constants[219] = 0.4
constants[220] = 0.02
constants[221] = 9.6
constants[222] = 0.6
states[54] = 0.1
states[55] = 0.1
states[56] = 0.1
constants[223] = 3.2
constants[224] = 0.02
constants[225] = 1.2
constants[226] = 0.2
states[57] = 0.1
states[58] = 0.1
constants[227] = 1.7
constants[228] = 0.3
constants[229] = 0.8
constants[230] = 0.4
constants[231] = 4.4
constants[232] = 3.4
constants[233] = 1.4
constants[234] = 0.8
constants[235] = 0.4
constants[236] = 2.4
constants[237] = 0.2
states[59] = 0.1
states[60] = 0.1
constants[238] = 1
constants[239] = 0.2
constants[240] = 2.4
constants[241] = 0.2
states[61] = 0.1
states[62] = 0.1
constants[242] = 0.8
constants[243] = 1.4
constants[244] = 1.4
states[63] = 0.1
constants[245] = 0.32
constants[246] = 0.8
constants[247] = 0.4
constants[248] = 1.4
constants[249] = 0.2
states[64] = 0.1
constants[250] = 3
constants[251] = 1.4
constants[252] = 0.4
states[65] = 0.1
constants[253] = 3
constants[254] = 1.6
return (states, constants)
def computeRates(voi, states, constants):
rates = [0.0] * sizeStates; algebraic = [0.0] * sizeAlgebraic
rates[0] = ((constants[0]*constants[1])/(constants[2]+constants[1])-constants[3]*states[0])*constants[4]
rates[1] = (((((constants[5]-constants[6]*states[1]*states[2])+constants[7]*states[3])-((constants[8]*states[1])/(constants[9]+states[1]))*(states[4]+states[5]))+(constants[10]*states[6])/(constants[11]+states[6]))-constants[12]*states[1])*constants[4]
rates[3] = (constants[6]*states[1]*states[2]-constants[7]*states[3])*constants[4]
rates[6] = ((((((((constants[8]*states[1])/(constants[9]+states[1]))*(states[4]+states[5])-(constants[10]*states[6])/(constants[11]+states[6]))-((constants[13]*states[6])/(constants[14]+states[6]))*states[7])+(constants[15]*states[8])/(constants[16]+states[8]))-constants[17]*states[6]*states[2])+constants[18]*states[9])-constants[19]*states[6])*constants[4]
rates[9] = (constants[17]*states[6]*states[2]-constants[18]*states[9])*constants[4]
rates[8] = ((((constants[13]*states[6])/(constants[14]+states[6]))*states[7]-(constants[15]*states[8])/(constants[16]+states[8]))-constants[20]*states[8])*constants[4]
rates[2] = (((((((constants[21]-constants[6]*states[1]*states[2])+constants[7]*states[3])-constants[17]*states[6]*states[2])+constants[18]*states[9])-(constants[22]*states[10]*states[2])/(constants[23]+states[2]))+(constants[24]*states[11])/(constants[25]+states[11]))-constants[26]*states[2])*constants[4]
rates[11] = (((constants[22]*states[10]*states[2])/(constants[23]+states[2])-(constants[24]*states[11])/(constants[25]+states[11]))-constants[27]*states[11])*constants[4]
rates[12] = (((((constants[28]*states[0]+(((constants[29]*states[2]*constants[30])/(constants[30]+states[1]))*constants[31])/(constants[31]+states[6]))-constants[32]*states[12]*(constants[33]-(states[13]+states[4]+states[5])))+constants[34]*states[13])-(constants[35]*states[12])/(constants[36]+states[12]))-constants[37]*states[12])*constants[4]
rates[13] = (((constants[32]*states[12]*(constants[33]-(states[13]+states[4]+states[5]))-constants[34]*states[13])+(constants[38]*states[4])/(constants[39]+states[4]))-(constants[40]*states[13])/(constants[41]+states[13]))*constants[4]
rates[4] = ((((constants[40]*states[13])/(constants[41]+states[13])-(constants[38]*states[4])/(constants[39]+states[4]))-constants[42]*states[4]*states[14])+constants[43]*states[5])*constants[4]
rates[5] = (constants[42]*states[4]*states[14]-constants[43]*states[5])*constants[4]
rates[15] = (((((((constants[44]*states[2]*constants[45])/(constants[45]+states[1]))*constants[46])/(constants[46]+states[6])-constants[47]*states[15]*(constants[48]-(states[16]+states[7]+states[17]+states[18]+states[10]+states[19])))+constants[49]*states[16])-(((constants[50]*states[20])/(constants[51]+states[20]))*states[15])/(constants[52]+states[15]))-constants[53]*states[15])*constants[4]
rates[16] = (((constants[47]*states[15]*(constants[48]-(states[16]+states[7]+states[17]+states[18]+states[10]+states[19]))-constants[49]*states[16])+(constants[54]*(states[21]+constants[55])*states[7])/(constants[56]+states[7]))-(constants[57]*states[22]*states[16])/(constants[58]+states[16]))*constants[4]
rates[7] = ((((constants[57]*states[22]*states[16])/(constants[58]+states[16])-(constants[54]*(states[21]+constants[55])*states[7])/(constants[56]+states[7]))-constants[59]*states[7]*states[14])+constants[60]*states[17])*constants[4]
rates[20] = ((constants[61]-(((constants[62]*states[20])/(constants[63]+states[20]))*states[23])/(constants[64]+states[23]))-constants[65]*states[20])*constants[4]
rates[17] = (constants[59]*states[7]*states[14]-constants[60]*states[17])*constants[4]
rates[24] = (((constants[66]+(constants[67]*(constants[68]+constants[69]*states[25])*states[22])/(constants[70]+states[22]))-(constants[71]*(states[7]+constants[72])*states[24])/(constants[73]+states[24]))-constants[74]*states[24])*constants[4]
rates[22] = ((constants[71]*(states[7]+constants[72])*(states[24]/(constants[73]+states[24]))-constants[67]*(constants[68]+constants[69]*states[25])*(states[22]/(constants[70]+states[22])))-constants[75]*states[22])*constants[4]
rates[26] = (((((((constants[76]*states[2]*constants[77])/(constants[77]+states[1]))*constants[78])/(constants[78]+states[6])-constants[79]*states[26]*(constants[48]-(states[16]+states[7]+states[17]+states[18]+states[10]+states[19])))+constants[80]*states[18])-(((constants[81]*states[26])/(constants[82]+states[26]))*states[27])/(constants[83]+states[27]))-constants[84]*states[26])*constants[4]
rates[18] = (((constants[79]*states[26]*(constants[48]-(states[16]+states[7]+states[17]+states[18]+states[10]+states[19]))-constants[80]*states[18])+(constants[85]*(states[21]+constants[86])*states[10])/(constants[87]+states[10]))-(constants[88]*states[28]*states[18])/(constants[89]+states[18]))*constants[4]
rates[10] = ((((constants[88]*states[28]*states[18])/(constants[89]+states[18])-(constants[85]*(states[21]+constants[86])*states[10])/(constants[87]+states[10]))-constants[90]*states[10]*states[14])+constants[91]*states[19])*constants[4]
rates[19] = (constants[90]*states[10]*states[14]-constants[91]*states[19])*constants[4]
rates[14] = ((((((((((((constants[92]+(((constants[93]*states[2]*constants[94])/(constants[94]+states[1]))*constants[95])/(constants[95]+states[6]))-constants[42]*states[4]*states[14])+constants[43]*states[5])-constants[59]*states[7]*states[14])+constants[60]*states[17])-constants[90]*states[10]*states[14])+constants[91]*states[19])-constants[96]*states[29]*states[14])+constants[97]*states[30])-(constants[98]*states[7]*states[14])/(constants[99]+states[14]))+(constants[100]*states[31])/(constants[101]+states[31]))-constants[102]*states[14])*constants[4]
rates[31] = ((((constants[98]*states[7]*states[14])/(constants[99]+states[14])-(constants[100]*states[31])/(constants[101]+states[31]))-(((constants[103]*states[20])/(constants[104]+states[20]))*states[31])/(constants[105]+states[31]))-constants[106]*states[31])*constants[4]
rates[32] = ((((constants[107]*states[23])/(constants[108]+states[23]))*(states[10]+states[29])-(constants[109]*states[32])/(constants[110]+states[32]))-constants[111]*states[32])*constants[4]
rates[23] = (((constants[112]+(constants[109]*states[32])/(constants[110]+states[32]))-((constants[107]*states[23])/(constants[108]+states[23]))*(states[10]+states[29]))-constants[113]*states[23])*constants[4]
rates[33] = (((constants[114]+(constants[115]*(constants[116]+constants[117]*states[25])*states[28])/(constants[118]+states[28]))-(constants[119]*(states[10]+constants[120])*states[33])/(constants[121]+states[33]))-constants[122]*states[33])*constants[4]
rates[28] = (((constants[119]*(states[10]+constants[120])*states[33])/(constants[121]+states[33])-(constants[115]*(constants[116]+constants[117]*states[25])*states[28])/(constants[118]+states[28]))-constants[123]*states[28])*constants[4]
rates[34] = ((((constants[124]-constants[125]*states[34]*(constants[126]-(states[35]+states[29]+states[30])))+constants[127]*states[35])-((constants[128]*states[34])/(constants[129]+states[34]))*(states[27]/(constants[130]+states[27])+states[23]/(constants[131]+states[23])))-constants[132]*states[34])*constants[4]
rates[35] = (((constants[125]*states[34]*(constants[126]-(states[35]+states[29]+states[30]))-constants[127]*states[35])+(constants[133]*(states[21]+constants[134])*states[29])/(constants[135]+states[29]))-(constants[136]*states[36]*states[35])/(constants[137]+states[35]))*constants[4]
rates[29] = ((((constants[136]*states[36]*states[35])/(constants[137]+states[35])-(constants[133]*(states[21]+constants[134])*states[29])/(constants[135]+states[29]))-constants[96]*states[29]*states[14])+constants[97]*states[30])*constants[4]
rates[30] = (constants[96]*states[29]*states[14]-constants[97]*states[30])*constants[4]
rates[37] = (((constants[138]-(constants[139]*states[29]*states[37])/(constants[140]+states[37]))+(constants[141]*states[27])/(constants[142]+states[27]))-constants[143]*states[37])*constants[4]
rates[27] = (((constants[139]*states[29]*states[37])/(constants[140]+states[37])-(constants[141]*states[27])/(constants[142]+states[27]))-constants[144]*states[27])*constants[4]
rates[38] = (((constants[145]+(constants[146]*(constants[147]+constants[148]*states[25])*states[36])/(constants[149]+states[36]))-(constants[150]*(states[29]+constants[151])*states[38])/(constants[152]+states[38]))-constants[153]*states[38])*constants[4]
rates[36] = (((constants[150]*(states[29]+constants[151])*states[38])/(constants[152]+states[38])-(constants[146]*(constants[147]+constants[148]*states[25])*states[36])/(constants[149]+states[36]))-constants[154]*states[36])*constants[4]
rates[21] = ((((constants[155]+constants[156]*states[39])-(constants[157]*(states[29]+constants[158])*states[21])/(constants[159]+states[21]))+(constants[160]*states[40])/(constants[161]+states[40]))-constants[162]*states[21])*constants[4]
rates[40] = (((constants[157]*(states[29]+constants[158])*states[21])/(constants[159]+states[21])-(constants[160]*states[40])/(constants[161]+states[40]))-constants[163]*states[40])*constants[4]
rates[41] = ((((constants[164]*states[7]*(constants[165]-states[41]))/(constants[166]+(constants[165]-states[41]))-(constants[167]*states[41])/(constants[168]+states[41]))-constants[169]*(constants[170]-states[42])*states[41])+constants[171]*states[42])*constants[4]
rates[42] = (constants[169]*(constants[170]-states[42])*states[41]-constants[171]*states[42])*constants[4]
rates[43] = (((constants[172]*states[42]-constants[173]*states[43])-constants[174]*(constants[175]-states[44])*states[43])+constants[176]*states[44])*constants[4]
rates[44] = (constants[174]*(constants[175]-states[44])*states[43]-constants[176]*states[44])*constants[4]
rates[25] = ((constants[177]*states[44]*(constants[178]-states[25]))/(constants[179]+(constants[178]-states[25]))-(constants[180]*states[25])/(constants[181]+states[25]))*constants[4]
rates[39] = (constants[182]*(power(states[45], constants[183])))/(power(constants[184], constants[183])+power(states[45], constants[183]))-(constants[185]*states[39])/(constants[186]+states[39])
rates[46] = (constants[187]*states[10]*((constants[188]-states[46])/(constants[189]+(constants[188]-states[46])))-constants[190]*(states[46]/(constants[191]+states[46])))*constants[4]
rates[47] = ((((constants[124]*states[46]*1.00000-constants[125]*states[47]*(constants[126]-(states[35]+states[29]+states[30])))+constants[127]*states[35])-constants[128]*(states[47]/(constants[129]+states[47]))*(states[27]/(constants[130]+states[27])+states[23]/(constants[131]+states[23])))-constants[132]*states[47])*constants[4]
rates[48] = (constants[192]*(power(states[45], constants[197])))/(power(constants[195], constants[197])+power(states[45], constants[197]))-((constants[193]*states[48])/(constants[196]+states[48])+constants[194]*states[48])
rates[49] = (constants[198]*(power(states[45], constants[197])))/(power(constants[201], constants[197])+power(states[45], constants[197]))-((constants[199]*states[49])/(constants[202]+states[49])+constants[200]*states[49])
rates[50] = (constants[203]*(power(constants[206], constants[208])))/(power(constants[206], constants[208])+power(states[51], constants[208]))-((constants[204]*states[50])/(constants[207]+states[50])+constants[205]*states[50])
rates[52] = (constants[209]*(power(states[45], constants[214])))/(power(constants[212], constants[214])+power(states[45], constants[214]))-((constants[210]*states[52])/(constants[213]+states[52])+constants[211]*states[52])
rates[53] = (constants[215]*states[48]+(constants[222]*states[54])/(constants[217]+states[54])+constants[219]*states[55])-((constants[221]*states[53])/(constants[216]+states[53])+constants[218]*states[53]*states[56]+constants[220]*states[53])
rates[56] = (constants[223]*states[49]+(constants[226]*states[57])/(constants[217]+states[57])+constants[219]*states[55])-((constants[225]*states[56])/(constants[216]+states[56])+constants[218]*states[53]*states[56]+constants[224]*states[56])
rates[58] = (constants[227]*states[52]+constants[230]*states[51])-(constants[229]*states[58]+(constants[231]*states[58])/(constants[228]+states[58])+constants[220]*states[58])
rates[54] = (constants[221]*states[53])/(constants[216]+states[53])-((constants[222]*states[54])/(constants[217]+states[54])+(constants[232]*states[54])/(constants[228]+states[54])+constants[220]*states[54])
rates[57] = (constants[225]*states[56])/(constants[216]+states[56])-((constants[226]*states[57])/(constants[217]+states[57])+(constants[233]*states[57])/(constants[228]+states[57])+constants[220]*states[57])
rates[55] = ((constants[237]*states[59])/(constants[217]+states[59])+constants[218]*states[53]*states[56]+constants[235]*states[60])-((constants[236]*states[55])/(constants[216]+states[55])+constants[219]*states[55]+constants[234]*states[55]+constants[220]*states[55])
rates[60] = ((constants[241]*states[61])/(constants[217]+states[61])+constants[234]*states[55]+constants[239]*states[62])-((constants[240]*states[60])/(constants[216]+states[60])+constants[235]*states[60]+constants[238]*states[45]*states[60]+constants[220]*states[60])
rates[51] = constants[229]*states[58]-(constants[230]*states[51]+(constants[242]*states[51])/(constants[228]+states[51])+constants[220]*states[51])
rates[59] = (constants[236]*states[55])/(constants[216]+states[55])-((constants[237]*states[59])/(constants[217]+states[59])+(constants[243]*states[59])/(constants[228]+states[59])+constants[220]*states[59])
rates[61] = (constants[240]*states[60])/(constants[216]+states[60])-((constants[241]*states[61])/(constants[217]+states[61])+(constants[244]*states[61])/(constants[228]+states[61])+constants[220]*states[61])
rates[63] = ((constants[249]*states[64])/(constants[217]+states[64])+constants[247]*states[45]+constants[245]*states[50])-((constants[248]*states[63])/(constants[216]+states[63])+constants[246]*states[63]+constants[220]*states[63])
rates[64] = (constants[248]*states[63])/(constants[216]+states[63])-((constants[249]*states[64])/(constants[217]+states[64])+(constants[250]*states[64])/(constants[228]+states[64])+constants[220]*states[64])
rates[45] = ((constants[252]*states[65])/(constants[217]+states[65])+constants[246]*states[63]+constants[239]*states[62])-((constants[251]*states[45])/(constants[216]+states[45])+constants[247]*states[45]+constants[238]*states[45]*states[60]+constants[220]*states[45])
rates[65] = (constants[251]*states[45])/(constants[216]+states[45])-((constants[252]*states[65])/(constants[217]+states[65])+(constants[253]*states[65])/(constants[228]+states[65])+constants[220]*states[65])
rates[62] = constants[238]*states[45]*states[60]-(constants[239]*states[62]+(constants[254]*states[62])/(constants[228]+states[62])+constants[220]*states[62])
return(rates)
def computeAlgebraic(constants, states, voi):
algebraic = array([[0.0] * len(voi)] * sizeAlgebraic)
states = array(states)
voi = array(voi)
return algebraic
def solve_model():
"""Solve model with ODE solver"""
from scipy.integrate import ode
# Initialise constants and state variables
(init_states, constants) = initConsts()
# Set timespan to solve over
voi = linspace(0, 10, 500)
# Construct ODE object to solve
r = ode(computeRates)
r.set_integrator('vode', method='bdf', atol=1e-06, rtol=1e-06, max_step=1)
r.set_initial_value(init_states, voi[0])
r.set_f_params(constants)
# Solve model
states = array([[0.0] * len(voi)] * sizeStates)
states[:,0] = init_states
for (i,t) in enumerate(voi[1:]):
if r.successful():
r.integrate(t)
states[:,i+1] = r.y
else:
break
# Compute algebraic variables
algebraic = computeAlgebraic(constants, states, voi)
return (voi, states, algebraic)
def plot_model(voi, states, algebraic):
"""Plot variables against variable of integration"""
import pylab
(legend_states, legend_algebraic, legend_voi, legend_constants) = createLegends()
pylab.figure(1)
pylab.plot(voi,vstack((states,algebraic)).T)
pylab.xlabel(legend_voi)
pylab.legend(legend_states + legend_algebraic, loc='best')
pylab.show()
if __name__ == "__main__":
(voi, states, algebraic) = solve_model()
plot_model(voi, states, algebraic)