# 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)