# Size of variable arrays: sizeAlgebraic = 288 sizeStates = 78 sizeConstants = 730 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 (second)" legend_algebraic[5] = "I_CaT in component I_CaT (picoampere)" legend_constants[0] = "g_CaT in component I_CaT (picosiemens)" legend_constants[1] = "ECaT in component I_CaT (millivolt)" legend_constants[2] = "alpha_OT in component I_CaT (first_order_rate_constant)" legend_constants[3] = "beta_OT in component I_CaT (first_order_rate_constant)" legend_states[0] = "OCaT in component I_CaT (dimensionless)" legend_states[1] = "Vm in component membrane (millivolt)" legend_states[2] = "dT in component I_CaT_d_gate (dimensionless)" legend_states[3] = "fT in component I_CaT_f_gate (dimensionless)" legend_algebraic[0] = "infinity_dT in component I_CaT_d_gate (dimensionless)" legend_constants[648] = "tau_dT1 in component I_CaT_d_gate (second)" legend_constants[4] = "A_dT1 in component I_CaT_d_gate (second)" legend_constants[5] = "k_dT in component I_CaT_d_gate (per_millivolt)" legend_constants[6] = "V_dT in component I_CaT_d_gate (millivolt)" legend_algebraic[1] = "infinity_fT in component I_CaT_f_gate (dimensionless)" legend_algebraic[17] = "tau_fT1 in component I_CaT_f_gate (second)" legend_constants[7] = "A_fT1 in component I_CaT_f_gate (second)" legend_constants[8] = "A_fT2 in component I_CaT_f_gate (second)" legend_constants[9] = "A_fT3 in component I_CaT_f_gate (per_millivolt)" legend_constants[10] = "A_fT4 in component I_CaT_f_gate (millivolt)" legend_constants[11] = "k_fT in component I_CaT_f_gate (per_millivolt)" legend_constants[12] = "V_fT in component I_CaT_f_gate (millivolt)" legend_algebraic[21] = "I_CaExt in component I_CaExt (picoampere)" legend_constants[13] = "gCaExt in component I_CaExt (femtoampere)" legend_constants[14] = "KCaExt in component I_CaExt (micromolar)" legend_states[4] = "CCy in component CCy (micromolar)" legend_algebraic[22] = "I_Kv1_1 in component I_Kv1_1 (picoampere)" legend_constants[15] = "g_Kv1_1 in component I_Kv1_1 (picosiemens)" legend_constants[649] = "EK in component Nerst_potentials (millivolt)" legend_states[5] = "dv1_1 in component I_Kv1_1_d_gate (dimensionless)" legend_states[6] = "fv1_1 in component I_Kv1_1_f_gate (dimensionless)" legend_algebraic[2] = "alpha_dv1_1 in component I_Kv1_1_d_gate (first_order_rate_constant)" legend_constants[16] = "A_dv1_11 in component I_Kv1_1_d_gate (first_order_rate_constant)" legend_constants[17] = "A_dv1_12 in component I_Kv1_1_d_gate (dimensionless)" legend_constants[18] = "A_dv1_13 in component I_Kv1_1_d_gate (per_millivolt)" legend_constants[19] = "A_dv1_14 in component I_Kv1_1_d_gate (millivolt)" legend_algebraic[18] = "beta_dv1_1 in component I_Kv1_1_d_gate (first_order_rate_constant)" legend_algebraic[3] = "alpha_fv1_1 in component I_Kv1_1_f_gate (first_order_rate_constant)" legend_constants[20] = "A_fv1_11 in component I_Kv1_1_f_gate (first_order_rate_constant)" legend_constants[21] = "A_fv1_12 in component I_Kv1_1_f_gate (dimensionless)" legend_constants[22] = "A_fv1_13 in component I_Kv1_1_f_gate (per_millivolt)" legend_constants[23] = "A_fv1_14 in component I_Kv1_1_f_gate (millivolt)" legend_algebraic[19] = "beta_fv1_1 in component I_Kv1_1_f_gate (first_order_rate_constant)" legend_algebraic[23] = "I_KERG in component I_KERG (picoampere)" legend_constants[24] = "g_KERG in component I_KERG (picosiemens)" legend_states[7] = "dERG in component I_KERG_d_gate (dimensionless)" legend_algebraic[4] = "alpha_dERG in component I_KERG_d_gate (first_order_rate_constant)" legend_constants[25] = "A_dERG1 in component I_KERG_d_gate (first_order_rate_constant)" legend_constants[26] = "A_dERG2 in component I_KERG_d_gate (dimensionless)" legend_constants[27] = "A_dERG3 in component I_KERG_d_gate (per_millivolt)" legend_constants[28] = "A_dERG4 in component I_KERG_d_gate (millivolt)" legend_algebraic[20] = "beta_dERG in component I_KERG_d_gate (first_order_rate_constant)" legend_algebraic[24] = "I_KB in component I_KB (picoampere)" legend_constants[29] = "g_KB in component I_KB (picosiemens)" legend_constants[30] = "EKB in component I_KB (millivolt)" legend_algebraic[25] = "I_L in component I_L (picoampere)" legend_constants[31] = "g_L in component I_L (picosiemens)" legend_constants[32] = "EL in component I_L (millivolt)" legend_algebraic[26] = "I_ion_Cy in component membrane (picoampere)" legend_constants[33] = "Cm in component membrane (picofarad)" legend_algebraic[64] = "I_ion_PU in component PU_membrane (picoampere)" legend_algebraic[88] = "I_ion_PU in component PU_membrane (picoampere)" legend_algebraic[112] = "I_ion_PU in component PU_membrane (picoampere)" legend_algebraic[136] = "I_ion_PU in component PU_membrane (picoampere)" legend_algebraic[160] = "I_ion_PU in component PU_membrane (picoampere)" legend_algebraic[184] = "I_ion_PU in component PU_membrane (picoampere)" legend_algebraic[208] = "I_ion_PU in component PU_membrane (picoampere)" legend_algebraic[232] = "I_ion_PU in component PU_membrane (picoampere)" legend_algebraic[256] = "I_ion_PU in component PU_membrane (picoampere)" legend_algebraic[279] = "I_ion_PU in component PU_membrane (picoampere)" legend_algebraic[6] = "ECa_Cy in component Nerst_potentials (millivolt)" legend_constants[34] = "T in component model_parameters (kelvin)" legend_constants[35] = "R in component model_parameters (attojoule_per_zeptomole_kelvin)" legend_constants[36] = "F in component model_parameters (femtocoulomb_per_zeptomole)" legend_constants[37] = "CO in component model_parameters (micromolar)" legend_constants[38] = "KO in component model_parameters (micromolar)" legend_constants[39] = "Ki in component model_parameters (micromolar)" legend_constants[40] = "lambda_S2_Cy in component CCy (dimensionless)" legend_constants[41] = "delta_SCy in component CCy (micromolar_per_coulomb)" legend_algebraic[37] = "JCy in component JCy (flux)" legend_algebraic[27] = "JS2Cy1 in component JS2Cy (flux)" legend_algebraic[28] = "JS2Cy2 in component JS2Cy (flux)" legend_algebraic[29] = "JS2Cy3 in component JS2Cy (flux)" legend_algebraic[30] = "JS2Cy4 in component JS2Cy (flux)" legend_algebraic[31] = "JS2Cy5 in component JS2Cy (flux)" legend_algebraic[32] = "JS2Cy6 in component JS2Cy (flux)" legend_algebraic[33] = "JS2Cy7 in component JS2Cy (flux)" legend_algebraic[34] = "JS2Cy8 in component JS2Cy (flux)" legend_algebraic[35] = "JS2Cy9 in component JS2Cy (flux)" legend_algebraic[36] = "JS2Cy10 in component JS2Cy (flux)" legend_constants[42] = "mu_S2Cy1 in component JS2Cy (first_order_rate_constant)" legend_constants[43] = "mu_S2Cy2 in component JS2Cy (first_order_rate_constant)" legend_constants[44] = "mu_S2Cy3 in component JS2Cy (first_order_rate_constant)" legend_constants[45] = "mu_S2Cy4 in component JS2Cy (first_order_rate_constant)" legend_constants[46] = "mu_S2Cy5 in component JS2Cy (first_order_rate_constant)" legend_constants[47] = "mu_S2Cy6 in component JS2Cy (first_order_rate_constant)" legend_constants[48] = "mu_S2Cy7 in component JS2Cy (first_order_rate_constant)" legend_constants[49] = "mu_S2Cy8 in component JS2Cy (first_order_rate_constant)" legend_constants[50] = "mu_S2Cy9 in component JS2Cy (first_order_rate_constant)" legend_constants[51] = "mu_S2Cy10 in component JS2Cy (first_order_rate_constant)" legend_states[8] = "CS2 in component CS2 (micromolar)" legend_states[9] = "CS2 in component CS2 (micromolar)" legend_states[10] = "CS2 in component CS2 (micromolar)" legend_states[11] = "CS2 in component CS2 (micromolar)" legend_states[12] = "CS2 in component CS2 (micromolar)" legend_states[13] = "CS2 in component CS2 (micromolar)" legend_states[14] = "CS2 in component CS2 (micromolar)" legend_states[15] = "CS2 in component CS2 (micromolar)" legend_states[16] = "CS2 in component CS2 (micromolar)" legend_states[17] = "CS2 in component CS2 (micromolar)" legend_algebraic[38] = "JS1S2_1 in component JS1S2 (flux)" legend_algebraic[39] = "JS1S2_2 in component JS1S2 (flux)" legend_algebraic[40] = "JS1S2_3 in component JS1S2 (flux)" legend_algebraic[41] = "JS1S2_4 in component JS1S2 (flux)" legend_algebraic[42] = "JS1S2_5 in component JS1S2 (flux)" legend_algebraic[43] = "JS1S2_6 in component JS1S2 (flux)" legend_algebraic[44] = "JS1S2_7 in component JS1S2 (flux)" legend_algebraic[45] = "JS1S2_8 in component JS1S2 (flux)" legend_algebraic[46] = "JS1S2_9 in component JS1S2 (flux)" legend_algebraic[47] = "JS1S2_10 in component JS1S2 (flux)" legend_constants[651] = "mu_S1S2_1 in component JS1S2 (first_order_rate_constant)" legend_constants[652] = "mu_S1S2_2 in component JS1S2 (first_order_rate_constant)" legend_constants[653] = "mu_S1S2_3 in component JS1S2 (first_order_rate_constant)" legend_constants[654] = "mu_S1S2_4 in component JS1S2 (first_order_rate_constant)" legend_constants[655] = "mu_S1S2_5 in component JS1S2 (first_order_rate_constant)" legend_constants[656] = "mu_S1S2_6 in component JS1S2 (first_order_rate_constant)" legend_constants[657] = "mu_S1S2_7 in component JS1S2 (first_order_rate_constant)" legend_constants[658] = "mu_S1S2_8 in component JS1S2 (first_order_rate_constant)" legend_constants[659] = "mu_S1S2_9 in component JS1S2 (first_order_rate_constant)" legend_constants[660] = "mu_S1S2_10 in component JS1S2 (first_order_rate_constant)" legend_constants[52] = "mu_A in component JS1S2 (first_order_rate_constant)" legend_constants[53] = "mu_B in component JS1S2 (first_order_rate_constant)" legend_constants[54] = "n_PU in component model_parameters (dimensionless)" legend_states[18] = "CS1 in component CS1 (micromolar)" legend_states[19] = "CS1 in component CS1 (micromolar)" legend_states[20] = "CS1 in component CS1 (micromolar)" legend_states[21] = "CS1 in component CS1 (micromolar)" legend_states[22] = "CS1 in component CS1 (micromolar)" legend_states[23] = "CS1 in component CS1 (micromolar)" legend_states[24] = "CS1 in component CS1 (micromolar)" legend_states[25] = "CS1 in component CS1 (micromolar)" legend_states[26] = "CS1 in component CS1 (micromolar)" legend_states[27] = "CS1 in component CS1 (micromolar)" legend_constants[55] = "mu_Cy in component JCy (first_order_rate_constant)" legend_constants[56] = "C_infinity in component JCy (micromolar)" legend_constants[57] = "P in component model_parameters (micromolar)" legend_algebraic[58] = "I_iCa in component PU_membrane (picoampere)" legend_algebraic[62] = "I_iNa in component PU_membrane (picoampere)" legend_constants[675] = "alpha_scale in component model_parameters (dimensionless)" legend_algebraic[50] = "I_Ca in component I_Ca (picoampere)" legend_algebraic[52] = "I_Na in component I_Na (picoampere)" legend_algebraic[60] = "I_NaP in component I_NaP (picoampere)" legend_algebraic[54] = "I_NSCC_Ca in component I_NSCC_Ca (picoampere)" legend_algebraic[57] = "I_PM in component I_PM (picoampere)" legend_algebraic[56] = "I_NSCC_Na in component I_NSCC_Na (picoampere)" legend_constants[58] = "gCa_ in component I_Ca (picosiemens)" legend_constants[59] = "kCa in component I_Ca (per_millivolt)" legend_algebraic[48] = "gCa in component I_Ca (picosiemens)" legend_algebraic[49] = "ECa_PU in component I_Ca (millivolt)" legend_constants[60] = "T in component model_parameters (kelvin)" legend_constants[61] = "R in component model_parameters (attojoule_per_zeptomole_kelvin)" legend_constants[62] = "F in component model_parameters (femtocoulomb_per_zeptomole)" legend_constants[63] = "CO in component model_parameters (micromolar)" legend_constants[64] = "gNa in component I_Na (picosiemens)" legend_algebraic[51] = "ENa_PU in component I_Na (millivolt)" legend_constants[65] = "NO in component model_parameters (micromolar)" legend_states[28] = "NS1 in component NS1 (micromolar)" legend_constants[66] = "gNSCC_Ca_ in component I_NSCC_Ca (picosiemens)" legend_algebraic[53] = "gNSCC_Ca in component I_NSCC_Ca (picosiemens)" legend_constants[67] = "hNSCC in component model_parameters (dimensionless)" legend_constants[68] = "ENSCC in component model_parameters (millivolt)" legend_constants[69] = "KNSCC in component model_parameters (micromolar)" legend_constants[70] = "gNSCC_Na_ in component I_NSCC_Na (picosiemens)" legend_algebraic[55] = "gNSCC_Na in component I_NSCC_Na (picosiemens)" legend_constants[71] = "gPM in component I_PM (femtoampere)" legend_constants[72] = "KPM in component I_PM (micromolar)" legend_constants[73] = "gNaP in component I_NaP (picosiemens)" legend_constants[74] = "KNaP in component I_NaP (micromolar)" legend_constants[75] = "hNaP in component I_NaP (dimensionless)" legend_constants[76] = "ENaP in component I_NaP (millivolt)" legend_algebraic[59] = "JSERCA in component JSERCA (flux)" legend_constants[77] = "VSERCA in component JSERCA (first_order_rate_constant)" legend_constants[78] = "A2 in component JSERCA (dimensionless)" legend_constants[79] = "A4 in component JSERCA (per_micromolar)" legend_constants[80] = "A5 in component JSERCA (per_micromolar)" legend_constants[81] = "A6 in component JSERCA (per_micromolar2)" legend_states[29] = "CER in component CER (micromolar)" legend_algebraic[63] = "JMCU in component JMCU (flux)" legend_constants[82] = "VMCU in component JMCU (flux)" legend_constants[83] = "KMCU in component JMCU (micromolar)" legend_algebraic[61] = "epsilon_INH in component JMCU (dimensionless)" legend_constants[84] = "KINH in component JMCU (micromolar)" legend_constants[85] = "hINH in component JMCU (dimensionless)" legend_states[30] = "CMT in component CMT (micromolar)" legend_algebraic[71] = "JIPR in component JIPR (flux)" legend_constants[86] = "kIPR in component JIPR (first_order_rate_constant)" legend_constants[87] = "k_1 in component JIPR (flux)" legend_constants[88] = "k1 in component JIPR (first_order_rate_constant)" legend_constants[89] = "k2 in component JIPR (first_order_rate_constant)" legend_constants[90] = "r2 in component JIPR (first_order_rate_constant)" legend_constants[91] = "r_2 in component JIPR (flux)" legend_constants[92] = "r4 in component JIPR (first_order_rate_constant)" legend_constants[93] = "R1 in component JIPR (micromolar)" legend_constants[94] = "R3 in component JIPR (micromolar)" legend_algebraic[65] = "phi1 in component JIPR (first_order_rate_constant)" legend_algebraic[68] = "phi_1 in component JIPR (flux)" legend_algebraic[72] = "phi2 in component JIPR (first_order_rate_constant)" legend_algebraic[74] = "phi3 in component JIPR (first_order_rate_constant)" legend_states[31] = "H in component JIPR (dimensionless)" legend_constants[95] = "g_beta in component JIPR (first_order_rate_constant)" legend_constants[96] = "h_beta in component JIPR (dimensionless)" legend_constants[97] = "g_alpha in component JIPR (first_order_rate_constant)" legend_constants[98] = "K_beta in component JIPR (micromolar)" legend_constants[661] = "alpha_zeta in component JIPR (first_order_rate_constant)" legend_algebraic[7] = "beta_zeta in component JIPR (first_order_rate_constant)" legend_states[32] = "zeta in component JIPR (dimensionless)" legend_constants[99] = "K_phi3_act in component JIPR (micromolar)" legend_constants[100] = "K_phi3_inh in component JIPR (micromolar)" legend_constants[101] = "h_phi3_act in component JIPR (dimensionless)" legend_constants[102] = "h_phi3_inh in component JIPR (dimensionless)" legend_constants[103] = "g_phi3 in component JIPR (first_order_rate_constant)" legend_constants[104] = "P in component model_parameters (micromolar)" legend_algebraic[66] = "JNCX in component JNCX (flux)" legend_constants[105] = "VNCX in component JNCX (flux)" legend_constants[106] = "KNCX in component JNCX (micromolar)" legend_constants[650] = "lambda_MT_S1 in component CS1 (dimensionless)" legend_constants[671] = "lambda_ER_S1 in component CS1 (dimensionless)" legend_constants[107] = "delta_SPU in component model_parameters (micromolar_per_coulomb)" legend_constants[108] = "gamma_S1 in component model_parameters (dimensionless)" legend_constants[109] = "gamma_MT in component model_parameters (dimensionless)" legend_constants[110] = "gamma_ER in component model_parameters (dimensionless)" legend_constants[672] = "lambda_MT_S2 in component CS2 (dimensionless)" legend_constants[673] = "lambda_ER_S2 in component CS2 (dimensionless)" legend_constants[674] = "lambda_S1_S2 in component CS2 (dimensionless)" legend_constants[111] = "gamma_S2 in component model_parameters (dimensionless)" legend_algebraic[69] = "fm in component CMT (dimensionless)" legend_constants[112] = "Km in component CMT (micromolar)" legend_constants[113] = "Bm in component CMT (micromolar)" legend_constants[114] = "KO in component model_parameters (micromolar)" legend_constants[115] = "Ki in component model_parameters (micromolar)" legend_constants[116] = "n_PU_base in component model_parameters (dimensionless)" legend_algebraic[82] = "I_iCa in component PU_membrane (picoampere)" legend_algebraic[86] = "I_iNa in component PU_membrane (picoampere)" legend_constants[681] = "alpha_scale in component model_parameters (dimensionless)" legend_algebraic[73] = "I_Ca in component I_Ca (picoampere)" legend_algebraic[76] = "I_Na in component I_Na (picoampere)" legend_algebraic[84] = "I_NaP in component I_NaP (picoampere)" legend_algebraic[78] = "I_NSCC_Ca in component I_NSCC_Ca (picoampere)" legend_algebraic[81] = "I_PM in component I_PM (picoampere)" legend_algebraic[80] = "I_NSCC_Na in component I_NSCC_Na (picoampere)" legend_constants[117] = "gCa_ in component I_Ca (picosiemens)" legend_constants[118] = "kCa in component I_Ca (per_millivolt)" legend_algebraic[67] = "gCa in component I_Ca (picosiemens)" legend_algebraic[70] = "ECa_PU in component I_Ca (millivolt)" legend_constants[119] = "T in component model_parameters (kelvin)" legend_constants[120] = "R in component model_parameters (attojoule_per_zeptomole_kelvin)" legend_constants[121] = "F in component model_parameters (femtocoulomb_per_zeptomole)" legend_constants[122] = "CO in component model_parameters (micromolar)" legend_constants[123] = "gNa in component I_Na (picosiemens)" legend_algebraic[75] = "ENa_PU in component I_Na (millivolt)" legend_constants[124] = "NO in component model_parameters (micromolar)" legend_states[33] = "NS1 in component NS1 (micromolar)" legend_constants[125] = "gNSCC_Ca_ in component I_NSCC_Ca (picosiemens)" legend_algebraic[77] = "gNSCC_Ca in component I_NSCC_Ca (picosiemens)" legend_constants[126] = "hNSCC in component model_parameters (dimensionless)" legend_constants[127] = "ENSCC in component model_parameters (millivolt)" legend_constants[128] = "KNSCC in component model_parameters (micromolar)" legend_constants[129] = "gNSCC_Na_ in component I_NSCC_Na (picosiemens)" legend_algebraic[79] = "gNSCC_Na in component I_NSCC_Na (picosiemens)" legend_constants[130] = "gPM in component I_PM (femtoampere)" legend_constants[131] = "KPM in component I_PM (micromolar)" legend_constants[132] = "gNaP in component I_NaP (picosiemens)" legend_constants[133] = "KNaP in component I_NaP (micromolar)" legend_constants[134] = "hNaP in component I_NaP (dimensionless)" legend_constants[135] = "ENaP in component I_NaP (millivolt)" legend_algebraic[83] = "JSERCA in component JSERCA (flux)" legend_constants[136] = "VSERCA in component JSERCA (first_order_rate_constant)" legend_constants[137] = "A2 in component JSERCA (dimensionless)" legend_constants[138] = "A4 in component JSERCA (per_micromolar)" legend_constants[139] = "A5 in component JSERCA (per_micromolar)" legend_constants[140] = "A6 in component JSERCA (per_micromolar2)" legend_states[34] = "CER in component CER (micromolar)" legend_algebraic[87] = "JMCU in component JMCU (flux)" legend_constants[141] = "VMCU in component JMCU (flux)" legend_constants[142] = "KMCU in component JMCU (micromolar)" legend_algebraic[85] = "epsilon_INH in component JMCU (dimensionless)" legend_constants[143] = "KINH in component JMCU (micromolar)" legend_constants[144] = "hINH in component JMCU (dimensionless)" legend_states[35] = "CMT in component CMT (micromolar)" legend_algebraic[95] = "JIPR in component JIPR (flux)" legend_constants[145] = "kIPR in component JIPR (first_order_rate_constant)" legend_constants[146] = "k_1 in component JIPR (flux)" legend_constants[147] = "k1 in component JIPR (first_order_rate_constant)" legend_constants[148] = "k2 in component JIPR (first_order_rate_constant)" legend_constants[149] = "r2 in component JIPR (first_order_rate_constant)" legend_constants[150] = "r_2 in component JIPR (flux)" legend_constants[151] = "r4 in component JIPR (first_order_rate_constant)" legend_constants[152] = "R1 in component JIPR (micromolar)" legend_constants[153] = "R3 in component JIPR (micromolar)" legend_algebraic[89] = "phi1 in component JIPR (first_order_rate_constant)" legend_algebraic[92] = "phi_1 in component JIPR (flux)" legend_algebraic[96] = "phi2 in component JIPR (first_order_rate_constant)" legend_algebraic[98] = "phi3 in component JIPR (first_order_rate_constant)" legend_states[36] = "H in component JIPR (dimensionless)" legend_constants[154] = "g_beta in component JIPR (first_order_rate_constant)" legend_constants[155] = "h_beta in component JIPR (dimensionless)" legend_constants[156] = "g_alpha in component JIPR (first_order_rate_constant)" legend_constants[157] = "K_beta in component JIPR (micromolar)" legend_constants[662] = "alpha_zeta in component JIPR (first_order_rate_constant)" legend_algebraic[8] = "beta_zeta in component JIPR (first_order_rate_constant)" legend_states[37] = "zeta in component JIPR (dimensionless)" legend_constants[158] = "K_phi3_act in component JIPR (micromolar)" legend_constants[159] = "K_phi3_inh in component JIPR (micromolar)" legend_constants[160] = "h_phi3_act in component JIPR (dimensionless)" legend_constants[161] = "h_phi3_inh in component JIPR (dimensionless)" legend_constants[162] = "g_phi3 in component JIPR (first_order_rate_constant)" legend_constants[163] = "P in component model_parameters (micromolar)" legend_algebraic[90] = "JNCX in component JNCX (flux)" legend_constants[164] = "VNCX in component JNCX (flux)" legend_constants[165] = "KNCX in component JNCX (micromolar)" legend_constants[676] = "lambda_MT_S1 in component CS1 (dimensionless)" legend_constants[677] = "lambda_ER_S1 in component CS1 (dimensionless)" legend_constants[166] = "delta_SPU in component model_parameters (micromolar_per_coulomb)" legend_constants[167] = "gamma_S1 in component model_parameters (dimensionless)" legend_constants[168] = "gamma_MT in component model_parameters (dimensionless)" legend_constants[169] = "gamma_ER in component model_parameters (dimensionless)" legend_constants[678] = "lambda_MT_S2 in component CS2 (dimensionless)" legend_constants[679] = "lambda_ER_S2 in component CS2 (dimensionless)" legend_constants[680] = "lambda_S1_S2 in component CS2 (dimensionless)" legend_constants[170] = "gamma_S2 in component model_parameters (dimensionless)" legend_algebraic[93] = "fm in component CMT (dimensionless)" legend_constants[171] = "Km in component CMT (micromolar)" legend_constants[172] = "Bm in component CMT (micromolar)" legend_constants[173] = "KO in component model_parameters (micromolar)" legend_constants[174] = "Ki in component model_parameters (micromolar)" legend_constants[175] = "n_PU_base in component model_parameters (dimensionless)" legend_algebraic[106] = "I_iCa in component PU_membrane (picoampere)" legend_algebraic[110] = "I_iNa in component PU_membrane (picoampere)" legend_constants[687] = "alpha_scale in component model_parameters (dimensionless)" legend_algebraic[97] = "I_Ca in component I_Ca (picoampere)" legend_algebraic[100] = "I_Na in component I_Na (picoampere)" legend_algebraic[108] = "I_NaP in component I_NaP (picoampere)" legend_algebraic[102] = "I_NSCC_Ca in component I_NSCC_Ca (picoampere)" legend_algebraic[105] = "I_PM in component I_PM (picoampere)" legend_algebraic[104] = "I_NSCC_Na in component I_NSCC_Na (picoampere)" legend_constants[176] = "gCa_ in component I_Ca (picosiemens)" legend_constants[177] = "kCa in component I_Ca (per_millivolt)" legend_algebraic[91] = "gCa in component I_Ca (picosiemens)" legend_algebraic[94] = "ECa_PU in component I_Ca (millivolt)" legend_constants[178] = "T in component model_parameters (kelvin)" legend_constants[179] = "R in component model_parameters (attojoule_per_zeptomole_kelvin)" legend_constants[180] = "F in component model_parameters (femtocoulomb_per_zeptomole)" legend_constants[181] = "CO in component model_parameters (micromolar)" legend_constants[182] = "gNa in component I_Na (picosiemens)" legend_algebraic[99] = "ENa_PU in component I_Na (millivolt)" legend_constants[183] = "NO in component model_parameters (micromolar)" legend_states[38] = "NS1 in component NS1 (micromolar)" legend_constants[184] = "gNSCC_Ca_ in component I_NSCC_Ca (picosiemens)" legend_algebraic[101] = "gNSCC_Ca in component I_NSCC_Ca (picosiemens)" legend_constants[185] = "hNSCC in component model_parameters (dimensionless)" legend_constants[186] = "ENSCC in component model_parameters (millivolt)" legend_constants[187] = "KNSCC in component model_parameters (micromolar)" legend_constants[188] = "gNSCC_Na_ in component I_NSCC_Na (picosiemens)" legend_algebraic[103] = "gNSCC_Na in component I_NSCC_Na (picosiemens)" legend_constants[189] = "gPM in component I_PM (femtoampere)" legend_constants[190] = "KPM in component I_PM (micromolar)" legend_constants[191] = "gNaP in component I_NaP (picosiemens)" legend_constants[192] = "KNaP in component I_NaP (micromolar)" legend_constants[193] = "hNaP in component I_NaP (dimensionless)" legend_constants[194] = "ENaP in component I_NaP (millivolt)" legend_algebraic[107] = "JSERCA in component JSERCA (flux)" legend_constants[195] = "VSERCA in component JSERCA (first_order_rate_constant)" legend_constants[196] = "A2 in component JSERCA (dimensionless)" legend_constants[197] = "A4 in component JSERCA (per_micromolar)" legend_constants[198] = "A5 in component JSERCA (per_micromolar)" legend_constants[199] = "A6 in component JSERCA (per_micromolar2)" legend_states[39] = "CER in component CER (micromolar)" legend_algebraic[111] = "JMCU in component JMCU (flux)" legend_constants[200] = "VMCU in component JMCU (flux)" legend_constants[201] = "KMCU in component JMCU (micromolar)" legend_algebraic[109] = "epsilon_INH in component JMCU (dimensionless)" legend_constants[202] = "KINH in component JMCU (micromolar)" legend_constants[203] = "hINH in component JMCU (dimensionless)" legend_states[40] = "CMT in component CMT (micromolar)" legend_algebraic[119] = "JIPR in component JIPR (flux)" legend_constants[204] = "kIPR in component JIPR (first_order_rate_constant)" legend_constants[205] = "k_1 in component JIPR (flux)" legend_constants[206] = "k1 in component JIPR (first_order_rate_constant)" legend_constants[207] = "k2 in component JIPR (first_order_rate_constant)" legend_constants[208] = "r2 in component JIPR (first_order_rate_constant)" legend_constants[209] = "r_2 in component JIPR (flux)" legend_constants[210] = "r4 in component JIPR (first_order_rate_constant)" legend_constants[211] = "R1 in component JIPR (micromolar)" legend_constants[212] = "R3 in component JIPR (micromolar)" legend_algebraic[113] = "phi1 in component JIPR (first_order_rate_constant)" legend_algebraic[116] = "phi_1 in component JIPR (flux)" legend_algebraic[120] = "phi2 in component JIPR (first_order_rate_constant)" legend_algebraic[122] = "phi3 in component JIPR (first_order_rate_constant)" legend_states[41] = "H in component JIPR (dimensionless)" legend_constants[213] = "g_beta in component JIPR (first_order_rate_constant)" legend_constants[214] = "h_beta in component JIPR (dimensionless)" legend_constants[215] = "g_alpha in component JIPR (first_order_rate_constant)" legend_constants[216] = "K_beta in component JIPR (micromolar)" legend_constants[663] = "alpha_zeta in component JIPR (first_order_rate_constant)" legend_algebraic[9] = "beta_zeta in component JIPR (first_order_rate_constant)" legend_states[42] = "zeta in component JIPR (dimensionless)" legend_constants[217] = "K_phi3_act in component JIPR (micromolar)" legend_constants[218] = "K_phi3_inh in component JIPR (micromolar)" legend_constants[219] = "h_phi3_act in component JIPR (dimensionless)" legend_constants[220] = "h_phi3_inh in component JIPR (dimensionless)" legend_constants[221] = "g_phi3 in component JIPR (first_order_rate_constant)" legend_constants[222] = "P in component model_parameters (micromolar)" legend_algebraic[114] = "JNCX in component JNCX (flux)" legend_constants[223] = "VNCX in component JNCX (flux)" legend_constants[224] = "KNCX in component JNCX (micromolar)" legend_constants[682] = "lambda_MT_S1 in component CS1 (dimensionless)" legend_constants[683] = "lambda_ER_S1 in component CS1 (dimensionless)" legend_constants[225] = "delta_SPU in component model_parameters (micromolar_per_coulomb)" legend_constants[226] = "gamma_S1 in component model_parameters (dimensionless)" legend_constants[227] = "gamma_MT in component model_parameters (dimensionless)" legend_constants[228] = "gamma_ER in component model_parameters (dimensionless)" legend_constants[684] = "lambda_MT_S2 in component CS2 (dimensionless)" legend_constants[685] = "lambda_ER_S2 in component CS2 (dimensionless)" legend_constants[686] = "lambda_S1_S2 in component CS2 (dimensionless)" legend_constants[229] = "gamma_S2 in component model_parameters (dimensionless)" legend_algebraic[117] = "fm in component CMT (dimensionless)" legend_constants[230] = "Km in component CMT (micromolar)" legend_constants[231] = "Bm in component CMT (micromolar)" legend_constants[232] = "KO in component model_parameters (micromolar)" legend_constants[233] = "Ki in component model_parameters (micromolar)" legend_constants[234] = "n_PU_base in component model_parameters (dimensionless)" legend_algebraic[130] = "I_iCa in component PU_membrane (picoampere)" legend_algebraic[134] = "I_iNa in component PU_membrane (picoampere)" legend_constants[693] = "alpha_scale in component model_parameters (dimensionless)" legend_algebraic[121] = "I_Ca in component I_Ca (picoampere)" legend_algebraic[124] = "I_Na in component I_Na (picoampere)" legend_algebraic[132] = "I_NaP in component I_NaP (picoampere)" legend_algebraic[126] = "I_NSCC_Ca in component I_NSCC_Ca (picoampere)" legend_algebraic[129] = "I_PM in component I_PM (picoampere)" legend_algebraic[128] = "I_NSCC_Na in component I_NSCC_Na (picoampere)" legend_constants[235] = "gCa_ in component I_Ca (picosiemens)" legend_constants[236] = "kCa in component I_Ca (per_millivolt)" legend_algebraic[115] = "gCa in component I_Ca (picosiemens)" legend_algebraic[118] = "ECa_PU in component I_Ca (millivolt)" legend_constants[237] = "T in component model_parameters (kelvin)" legend_constants[238] = "R in component model_parameters (attojoule_per_zeptomole_kelvin)" legend_constants[239] = "F in component model_parameters (femtocoulomb_per_zeptomole)" legend_constants[240] = "CO in component model_parameters (micromolar)" legend_constants[241] = "gNa in component I_Na (picosiemens)" legend_algebraic[123] = "ENa_PU in component I_Na (millivolt)" legend_constants[242] = "NO in component model_parameters (micromolar)" legend_states[43] = "NS1 in component NS1 (micromolar)" legend_constants[243] = "gNSCC_Ca_ in component I_NSCC_Ca (picosiemens)" legend_algebraic[125] = "gNSCC_Ca in component I_NSCC_Ca (picosiemens)" legend_constants[244] = "hNSCC in component model_parameters (dimensionless)" legend_constants[245] = "ENSCC in component model_parameters (millivolt)" legend_constants[246] = "KNSCC in component model_parameters (micromolar)" legend_constants[247] = "gNSCC_Na_ in component I_NSCC_Na (picosiemens)" legend_algebraic[127] = "gNSCC_Na in component I_NSCC_Na (picosiemens)" legend_constants[248] = "gPM in component I_PM (femtoampere)" legend_constants[249] = "KPM in component I_PM (micromolar)" legend_constants[250] = "gNaP in component I_NaP (picosiemens)" legend_constants[251] = "KNaP in component I_NaP (micromolar)" legend_constants[252] = "hNaP in component I_NaP (dimensionless)" legend_constants[253] = "ENaP in component I_NaP (millivolt)" legend_algebraic[131] = "JSERCA in component JSERCA (flux)" legend_constants[254] = "VSERCA in component JSERCA (first_order_rate_constant)" legend_constants[255] = "A2 in component JSERCA (dimensionless)" legend_constants[256] = "A4 in component JSERCA (per_micromolar)" legend_constants[257] = "A5 in component JSERCA (per_micromolar)" legend_constants[258] = "A6 in component JSERCA (per_micromolar2)" legend_states[44] = "CER in component CER (micromolar)" legend_algebraic[135] = "JMCU in component JMCU (flux)" legend_constants[259] = "VMCU in component JMCU (flux)" legend_constants[260] = "KMCU in component JMCU (micromolar)" legend_algebraic[133] = "epsilon_INH in component JMCU (dimensionless)" legend_constants[261] = "KINH in component JMCU (micromolar)" legend_constants[262] = "hINH in component JMCU (dimensionless)" legend_states[45] = "CMT in component CMT (micromolar)" legend_algebraic[143] = "JIPR in component JIPR (flux)" legend_constants[263] = "kIPR in component JIPR (first_order_rate_constant)" legend_constants[264] = "k_1 in component JIPR (flux)" legend_constants[265] = "k1 in component JIPR (first_order_rate_constant)" legend_constants[266] = "k2 in component JIPR (first_order_rate_constant)" legend_constants[267] = "r2 in component JIPR (first_order_rate_constant)" legend_constants[268] = "r_2 in component JIPR (flux)" legend_constants[269] = "r4 in component JIPR (first_order_rate_constant)" legend_constants[270] = "R1 in component JIPR (micromolar)" legend_constants[271] = "R3 in component JIPR (micromolar)" legend_algebraic[137] = "phi1 in component JIPR (first_order_rate_constant)" legend_algebraic[140] = "phi_1 in component JIPR (flux)" legend_algebraic[144] = "phi2 in component JIPR (first_order_rate_constant)" legend_algebraic[146] = "phi3 in component JIPR (first_order_rate_constant)" legend_states[46] = "H in component JIPR (dimensionless)" legend_constants[272] = "g_beta in component JIPR (first_order_rate_constant)" legend_constants[273] = "h_beta in component JIPR (dimensionless)" legend_constants[274] = "g_alpha in component JIPR (first_order_rate_constant)" legend_constants[275] = "K_beta in component JIPR (micromolar)" legend_constants[664] = "alpha_zeta in component JIPR (first_order_rate_constant)" legend_algebraic[10] = "beta_zeta in component JIPR (first_order_rate_constant)" legend_states[47] = "zeta in component JIPR (dimensionless)" legend_constants[276] = "K_phi3_act in component JIPR (micromolar)" legend_constants[277] = "K_phi3_inh in component JIPR (micromolar)" legend_constants[278] = "h_phi3_act in component JIPR (dimensionless)" legend_constants[279] = "h_phi3_inh in component JIPR (dimensionless)" legend_constants[280] = "g_phi3 in component JIPR (first_order_rate_constant)" legend_constants[281] = "P in component model_parameters (micromolar)" legend_algebraic[138] = "JNCX in component JNCX (flux)" legend_constants[282] = "VNCX in component JNCX (flux)" legend_constants[283] = "KNCX in component JNCX (micromolar)" legend_constants[688] = "lambda_MT_S1 in component CS1 (dimensionless)" legend_constants[689] = "lambda_ER_S1 in component CS1 (dimensionless)" legend_constants[284] = "delta_SPU in component model_parameters (micromolar_per_coulomb)" legend_constants[285] = "gamma_S1 in component model_parameters (dimensionless)" legend_constants[286] = "gamma_MT in component model_parameters (dimensionless)" legend_constants[287] = "gamma_ER in component model_parameters (dimensionless)" legend_constants[690] = "lambda_MT_S2 in component CS2 (dimensionless)" legend_constants[691] = "lambda_ER_S2 in component CS2 (dimensionless)" legend_constants[692] = "lambda_S1_S2 in component CS2 (dimensionless)" legend_constants[288] = "gamma_S2 in component model_parameters (dimensionless)" legend_algebraic[141] = "fm in component CMT (dimensionless)" legend_constants[289] = "Km in component CMT (micromolar)" legend_constants[290] = "Bm in component CMT (micromolar)" legend_constants[291] = "KO in component model_parameters (micromolar)" legend_constants[292] = "Ki in component model_parameters (micromolar)" legend_constants[293] = "n_PU_base in component model_parameters (dimensionless)" legend_algebraic[154] = "I_iCa in component PU_membrane (picoampere)" legend_algebraic[158] = "I_iNa in component PU_membrane (picoampere)" legend_constants[699] = "alpha_scale in component model_parameters (dimensionless)" legend_algebraic[145] = "I_Ca in component I_Ca (picoampere)" legend_algebraic[148] = "I_Na in component I_Na (picoampere)" legend_algebraic[156] = "I_NaP in component I_NaP (picoampere)" legend_algebraic[150] = "I_NSCC_Ca in component I_NSCC_Ca (picoampere)" legend_algebraic[153] = "I_PM in component I_PM (picoampere)" legend_algebraic[152] = "I_NSCC_Na in component I_NSCC_Na (picoampere)" legend_constants[294] = "gCa_ in component I_Ca (picosiemens)" legend_constants[295] = "kCa in component I_Ca (per_millivolt)" legend_algebraic[139] = "gCa in component I_Ca (picosiemens)" legend_algebraic[142] = "ECa_PU in component I_Ca (millivolt)" legend_constants[296] = "T in component model_parameters (kelvin)" legend_constants[297] = "R in component model_parameters (attojoule_per_zeptomole_kelvin)" legend_constants[298] = "F in component model_parameters (femtocoulomb_per_zeptomole)" legend_constants[299] = "CO in component model_parameters (micromolar)" legend_constants[300] = "gNa in component I_Na (picosiemens)" legend_algebraic[147] = "ENa_PU in component I_Na (millivolt)" legend_constants[301] = "NO in component model_parameters (micromolar)" legend_states[48] = "NS1 in component NS1 (micromolar)" legend_constants[302] = "gNSCC_Ca_ in component I_NSCC_Ca (picosiemens)" legend_algebraic[149] = "gNSCC_Ca in component I_NSCC_Ca (picosiemens)" legend_constants[303] = "hNSCC in component model_parameters (dimensionless)" legend_constants[304] = "ENSCC in component model_parameters (millivolt)" legend_constants[305] = "KNSCC in component model_parameters (micromolar)" legend_constants[306] = "gNSCC_Na_ in component I_NSCC_Na (picosiemens)" legend_algebraic[151] = "gNSCC_Na in component I_NSCC_Na (picosiemens)" legend_constants[307] = "gPM in component I_PM (femtoampere)" legend_constants[308] = "KPM in component I_PM (micromolar)" legend_constants[309] = "gNaP in component I_NaP (picosiemens)" legend_constants[310] = "KNaP in component I_NaP (micromolar)" legend_constants[311] = "hNaP in component I_NaP (dimensionless)" legend_constants[312] = "ENaP in component I_NaP (millivolt)" legend_algebraic[155] = "JSERCA in component JSERCA (flux)" legend_constants[313] = "VSERCA in component JSERCA (first_order_rate_constant)" legend_constants[314] = "A2 in component JSERCA (dimensionless)" legend_constants[315] = "A4 in component JSERCA (per_micromolar)" legend_constants[316] = "A5 in component JSERCA (per_micromolar)" legend_constants[317] = "A6 in component JSERCA (per_micromolar2)" legend_states[49] = "CER in component CER (micromolar)" legend_algebraic[159] = "JMCU in component JMCU (flux)" legend_constants[318] = "VMCU in component JMCU (flux)" legend_constants[319] = "KMCU in component JMCU (micromolar)" legend_algebraic[157] = "epsilon_INH in component JMCU (dimensionless)" legend_constants[320] = "KINH in component JMCU (micromolar)" legend_constants[321] = "hINH in component JMCU (dimensionless)" legend_states[50] = "CMT in component CMT (micromolar)" legend_algebraic[167] = "JIPR in component JIPR (flux)" legend_constants[322] = "kIPR in component JIPR (first_order_rate_constant)" legend_constants[323] = "k_1 in component JIPR (flux)" legend_constants[324] = "k1 in component JIPR (first_order_rate_constant)" legend_constants[325] = "k2 in component JIPR (first_order_rate_constant)" legend_constants[326] = "r2 in component JIPR (first_order_rate_constant)" legend_constants[327] = "r_2 in component JIPR (flux)" legend_constants[328] = "r4 in component JIPR (first_order_rate_constant)" legend_constants[329] = "R1 in component JIPR (micromolar)" legend_constants[330] = "R3 in component JIPR (micromolar)" legend_algebraic[161] = "phi1 in component JIPR (first_order_rate_constant)" legend_algebraic[164] = "phi_1 in component JIPR (flux)" legend_algebraic[168] = "phi2 in component JIPR (first_order_rate_constant)" legend_algebraic[170] = "phi3 in component JIPR (first_order_rate_constant)" legend_states[51] = "H in component JIPR (dimensionless)" legend_constants[331] = "g_beta in component JIPR (first_order_rate_constant)" legend_constants[332] = "h_beta in component JIPR (dimensionless)" legend_constants[333] = "g_alpha in component JIPR (first_order_rate_constant)" legend_constants[334] = "K_beta in component JIPR (micromolar)" legend_constants[665] = "alpha_zeta in component JIPR (first_order_rate_constant)" legend_algebraic[11] = "beta_zeta in component JIPR (first_order_rate_constant)" legend_states[52] = "zeta in component JIPR (dimensionless)" legend_constants[335] = "K_phi3_act in component JIPR (micromolar)" legend_constants[336] = "K_phi3_inh in component JIPR (micromolar)" legend_constants[337] = "h_phi3_act in component JIPR (dimensionless)" legend_constants[338] = "h_phi3_inh in component JIPR (dimensionless)" legend_constants[339] = "g_phi3 in component JIPR (first_order_rate_constant)" legend_constants[340] = "P in component model_parameters (micromolar)" legend_algebraic[162] = "JNCX in component JNCX (flux)" legend_constants[341] = "VNCX in component JNCX (flux)" legend_constants[342] = "KNCX in component JNCX (micromolar)" legend_constants[694] = "lambda_MT_S1 in component CS1 (dimensionless)" legend_constants[695] = "lambda_ER_S1 in component CS1 (dimensionless)" legend_constants[343] = "delta_SPU in component model_parameters (micromolar_per_coulomb)" legend_constants[344] = "gamma_S1 in component model_parameters (dimensionless)" legend_constants[345] = "gamma_MT in component model_parameters (dimensionless)" legend_constants[346] = "gamma_ER in component model_parameters (dimensionless)" legend_constants[696] = "lambda_MT_S2 in component CS2 (dimensionless)" legend_constants[697] = "lambda_ER_S2 in component CS2 (dimensionless)" legend_constants[698] = "lambda_S1_S2 in component CS2 (dimensionless)" legend_constants[347] = "gamma_S2 in component model_parameters (dimensionless)" legend_algebraic[165] = "fm in component CMT (dimensionless)" legend_constants[348] = "Km in component CMT (micromolar)" legend_constants[349] = "Bm in component CMT (micromolar)" legend_constants[350] = "KO in component model_parameters (micromolar)" legend_constants[351] = "Ki in component model_parameters (micromolar)" legend_constants[352] = "n_PU_base in component model_parameters (dimensionless)" legend_algebraic[178] = "I_iCa in component PU_membrane (picoampere)" legend_algebraic[182] = "I_iNa in component PU_membrane (picoampere)" legend_constants[705] = "alpha_scale in component model_parameters (dimensionless)" legend_algebraic[169] = "I_Ca in component I_Ca (picoampere)" legend_algebraic[172] = "I_Na in component I_Na (picoampere)" legend_algebraic[180] = "I_NaP in component I_NaP (picoampere)" legend_algebraic[174] = "I_NSCC_Ca in component I_NSCC_Ca (picoampere)" legend_algebraic[177] = "I_PM in component I_PM (picoampere)" legend_algebraic[176] = "I_NSCC_Na in component I_NSCC_Na (picoampere)" legend_constants[353] = "gCa_ in component I_Ca (picosiemens)" legend_constants[354] = "kCa in component I_Ca (per_millivolt)" legend_algebraic[163] = "gCa in component I_Ca (picosiemens)" legend_algebraic[166] = "ECa_PU in component I_Ca (millivolt)" legend_constants[355] = "T in component model_parameters (kelvin)" legend_constants[356] = "R in component model_parameters (attojoule_per_zeptomole_kelvin)" legend_constants[357] = "F in component model_parameters (femtocoulomb_per_zeptomole)" legend_constants[358] = "CO in component model_parameters (micromolar)" legend_constants[359] = "gNa in component I_Na (picosiemens)" legend_algebraic[171] = "ENa_PU in component I_Na (millivolt)" legend_constants[360] = "NO in component model_parameters (micromolar)" legend_states[53] = "NS1 in component NS1 (micromolar)" legend_constants[361] = "gNSCC_Ca_ in component I_NSCC_Ca (picosiemens)" legend_algebraic[173] = "gNSCC_Ca in component I_NSCC_Ca (picosiemens)" legend_constants[362] = "hNSCC in component model_parameters (dimensionless)" legend_constants[363] = "ENSCC in component model_parameters (millivolt)" legend_constants[364] = "KNSCC in component model_parameters (micromolar)" legend_constants[365] = "gNSCC_Na_ in component I_NSCC_Na (picosiemens)" legend_algebraic[175] = "gNSCC_Na in component I_NSCC_Na (picosiemens)" legend_constants[366] = "gPM in component I_PM (femtoampere)" legend_constants[367] = "KPM in component I_PM (micromolar)" legend_constants[368] = "gNaP in component I_NaP (picosiemens)" legend_constants[369] = "KNaP in component I_NaP (micromolar)" legend_constants[370] = "hNaP in component I_NaP (dimensionless)" legend_constants[371] = "ENaP in component I_NaP (millivolt)" legend_algebraic[179] = "JSERCA in component JSERCA (flux)" legend_constants[372] = "VSERCA in component JSERCA (first_order_rate_constant)" legend_constants[373] = "A2 in component JSERCA (dimensionless)" legend_constants[374] = "A4 in component JSERCA (per_micromolar)" legend_constants[375] = "A5 in component JSERCA (per_micromolar)" legend_constants[376] = "A6 in component JSERCA (per_micromolar2)" legend_states[54] = "CER in component CER (micromolar)" legend_algebraic[183] = "JMCU in component JMCU (flux)" legend_constants[377] = "VMCU in component JMCU (flux)" legend_constants[378] = "KMCU in component JMCU (micromolar)" legend_algebraic[181] = "epsilon_INH in component JMCU (dimensionless)" legend_constants[379] = "KINH in component JMCU (micromolar)" legend_constants[380] = "hINH in component JMCU (dimensionless)" legend_states[55] = "CMT in component CMT (micromolar)" legend_algebraic[191] = "JIPR in component JIPR (flux)" legend_constants[381] = "kIPR in component JIPR (first_order_rate_constant)" legend_constants[382] = "k_1 in component JIPR (flux)" legend_constants[383] = "k1 in component JIPR (first_order_rate_constant)" legend_constants[384] = "k2 in component JIPR (first_order_rate_constant)" legend_constants[385] = "r2 in component JIPR (first_order_rate_constant)" legend_constants[386] = "r_2 in component JIPR (flux)" legend_constants[387] = "r4 in component JIPR (first_order_rate_constant)" legend_constants[388] = "R1 in component JIPR (micromolar)" legend_constants[389] = "R3 in component JIPR (micromolar)" legend_algebraic[185] = "phi1 in component JIPR (first_order_rate_constant)" legend_algebraic[188] = "phi_1 in component JIPR (flux)" legend_algebraic[192] = "phi2 in component JIPR (first_order_rate_constant)" legend_algebraic[194] = "phi3 in component JIPR (first_order_rate_constant)" legend_states[56] = "H in component JIPR (dimensionless)" legend_constants[390] = "g_beta in component JIPR (first_order_rate_constant)" legend_constants[391] = "h_beta in component JIPR (dimensionless)" legend_constants[392] = "g_alpha in component JIPR (first_order_rate_constant)" legend_constants[393] = "K_beta in component JIPR (micromolar)" legend_constants[666] = "alpha_zeta in component JIPR (first_order_rate_constant)" legend_algebraic[12] = "beta_zeta in component JIPR (first_order_rate_constant)" legend_states[57] = "zeta in component JIPR (dimensionless)" legend_constants[394] = "K_phi3_act in component JIPR (micromolar)" legend_constants[395] = "K_phi3_inh in component JIPR (micromolar)" legend_constants[396] = "h_phi3_act in component JIPR (dimensionless)" legend_constants[397] = "h_phi3_inh in component JIPR (dimensionless)" legend_constants[398] = "g_phi3 in component JIPR (first_order_rate_constant)" legend_constants[399] = "P in component model_parameters (micromolar)" legend_algebraic[186] = "JNCX in component JNCX (flux)" legend_constants[400] = "VNCX in component JNCX (flux)" legend_constants[401] = "KNCX in component JNCX (micromolar)" legend_constants[700] = "lambda_MT_S1 in component CS1 (dimensionless)" legend_constants[701] = "lambda_ER_S1 in component CS1 (dimensionless)" legend_constants[402] = "delta_SPU in component model_parameters (micromolar_per_coulomb)" legend_constants[403] = "gamma_S1 in component model_parameters (dimensionless)" legend_constants[404] = "gamma_MT in component model_parameters (dimensionless)" legend_constants[405] = "gamma_ER in component model_parameters (dimensionless)" legend_constants[702] = "lambda_MT_S2 in component CS2 (dimensionless)" legend_constants[703] = "lambda_ER_S2 in component CS2 (dimensionless)" legend_constants[704] = "lambda_S1_S2 in component CS2 (dimensionless)" legend_constants[406] = "gamma_S2 in component model_parameters (dimensionless)" legend_algebraic[189] = "fm in component CMT (dimensionless)" legend_constants[407] = "Km in component CMT (micromolar)" legend_constants[408] = "Bm in component CMT (micromolar)" legend_constants[409] = "KO in component model_parameters (micromolar)" legend_constants[410] = "Ki in component model_parameters (micromolar)" legend_constants[411] = "n_PU_base in component model_parameters (dimensionless)" legend_algebraic[202] = "I_iCa in component PU_membrane (picoampere)" legend_algebraic[206] = "I_iNa in component PU_membrane (picoampere)" legend_constants[711] = "alpha_scale in component model_parameters (dimensionless)" legend_algebraic[193] = "I_Ca in component I_Ca (picoampere)" legend_algebraic[196] = "I_Na in component I_Na (picoampere)" legend_algebraic[204] = "I_NaP in component I_NaP (picoampere)" legend_algebraic[198] = "I_NSCC_Ca in component I_NSCC_Ca (picoampere)" legend_algebraic[201] = "I_PM in component I_PM (picoampere)" legend_algebraic[200] = "I_NSCC_Na in component I_NSCC_Na (picoampere)" legend_constants[412] = "gCa_ in component I_Ca (picosiemens)" legend_constants[413] = "kCa in component I_Ca (per_millivolt)" legend_algebraic[187] = "gCa in component I_Ca (picosiemens)" legend_algebraic[190] = "ECa_PU in component I_Ca (millivolt)" legend_constants[414] = "T in component model_parameters (kelvin)" legend_constants[415] = "R in component model_parameters (attojoule_per_zeptomole_kelvin)" legend_constants[416] = "F in component model_parameters (femtocoulomb_per_zeptomole)" legend_constants[417] = "CO in component model_parameters (micromolar)" legend_constants[418] = "gNa in component I_Na (picosiemens)" legend_algebraic[195] = "ENa_PU in component I_Na (millivolt)" legend_constants[419] = "NO in component model_parameters (micromolar)" legend_states[58] = "NS1 in component NS1 (micromolar)" legend_constants[420] = "gNSCC_Ca_ in component I_NSCC_Ca (picosiemens)" legend_algebraic[197] = "gNSCC_Ca in component I_NSCC_Ca (picosiemens)" legend_constants[421] = "hNSCC in component model_parameters (dimensionless)" legend_constants[422] = "ENSCC in component model_parameters (millivolt)" legend_constants[423] = "KNSCC in component model_parameters (micromolar)" legend_constants[424] = "gNSCC_Na_ in component I_NSCC_Na (picosiemens)" legend_algebraic[199] = "gNSCC_Na in component I_NSCC_Na (picosiemens)" legend_constants[425] = "gPM in component I_PM (femtoampere)" legend_constants[426] = "KPM in component I_PM (micromolar)" legend_constants[427] = "gNaP in component I_NaP (picosiemens)" legend_constants[428] = "KNaP in component I_NaP (micromolar)" legend_constants[429] = "hNaP in component I_NaP (dimensionless)" legend_constants[430] = "ENaP in component I_NaP (millivolt)" legend_algebraic[203] = "JSERCA in component JSERCA (flux)" legend_constants[431] = "VSERCA in component JSERCA (first_order_rate_constant)" legend_constants[432] = "A2 in component JSERCA (dimensionless)" legend_constants[433] = "A4 in component JSERCA (per_micromolar)" legend_constants[434] = "A5 in component JSERCA (per_micromolar)" legend_constants[435] = "A6 in component JSERCA (per_micromolar2)" legend_states[59] = "CER in component CER (micromolar)" legend_algebraic[207] = "JMCU in component JMCU (flux)" legend_constants[436] = "VMCU in component JMCU (flux)" legend_constants[437] = "KMCU in component JMCU (micromolar)" legend_algebraic[205] = "epsilon_INH in component JMCU (dimensionless)" legend_constants[438] = "KINH in component JMCU (micromolar)" legend_constants[439] = "hINH in component JMCU (dimensionless)" legend_states[60] = "CMT in component CMT (micromolar)" legend_algebraic[215] = "JIPR in component JIPR (flux)" legend_constants[440] = "kIPR in component JIPR (first_order_rate_constant)" legend_constants[441] = "k_1 in component JIPR (flux)" legend_constants[442] = "k1 in component JIPR (first_order_rate_constant)" legend_constants[443] = "k2 in component JIPR (first_order_rate_constant)" legend_constants[444] = "r2 in component JIPR (first_order_rate_constant)" legend_constants[445] = "r_2 in component JIPR (flux)" legend_constants[446] = "r4 in component JIPR (first_order_rate_constant)" legend_constants[447] = "R1 in component JIPR (micromolar)" legend_constants[448] = "R3 in component JIPR (micromolar)" legend_algebraic[209] = "phi1 in component JIPR (first_order_rate_constant)" legend_algebraic[212] = "phi_1 in component JIPR (flux)" legend_algebraic[216] = "phi2 in component JIPR (first_order_rate_constant)" legend_algebraic[218] = "phi3 in component JIPR (first_order_rate_constant)" legend_states[61] = "H in component JIPR (dimensionless)" legend_constants[449] = "g_beta in component JIPR (first_order_rate_constant)" legend_constants[450] = "h_beta in component JIPR (dimensionless)" legend_constants[451] = "g_alpha in component JIPR (first_order_rate_constant)" legend_constants[452] = "K_beta in component JIPR (micromolar)" legend_constants[667] = "alpha_zeta in component JIPR (first_order_rate_constant)" legend_algebraic[13] = "beta_zeta in component JIPR (first_order_rate_constant)" legend_states[62] = "zeta in component JIPR (dimensionless)" legend_constants[453] = "K_phi3_act in component JIPR (micromolar)" legend_constants[454] = "K_phi3_inh in component JIPR (micromolar)" legend_constants[455] = "h_phi3_act in component JIPR (dimensionless)" legend_constants[456] = "h_phi3_inh in component JIPR (dimensionless)" legend_constants[457] = "g_phi3 in component JIPR (first_order_rate_constant)" legend_constants[458] = "P in component model_parameters (micromolar)" legend_algebraic[210] = "JNCX in component JNCX (flux)" legend_constants[459] = "VNCX in component JNCX (flux)" legend_constants[460] = "KNCX in component JNCX (micromolar)" legend_constants[706] = "lambda_MT_S1 in component CS1 (dimensionless)" legend_constants[707] = "lambda_ER_S1 in component CS1 (dimensionless)" legend_constants[461] = "delta_SPU in component model_parameters (micromolar_per_coulomb)" legend_constants[462] = "gamma_S1 in component model_parameters (dimensionless)" legend_constants[463] = "gamma_MT in component model_parameters (dimensionless)" legend_constants[464] = "gamma_ER in component model_parameters (dimensionless)" legend_constants[708] = "lambda_MT_S2 in component CS2 (dimensionless)" legend_constants[709] = "lambda_ER_S2 in component CS2 (dimensionless)" legend_constants[710] = "lambda_S1_S2 in component CS2 (dimensionless)" legend_constants[465] = "gamma_S2 in component model_parameters (dimensionless)" legend_algebraic[213] = "fm in component CMT (dimensionless)" legend_constants[466] = "Km in component CMT (micromolar)" legend_constants[467] = "Bm in component CMT (micromolar)" legend_constants[468] = "KO in component model_parameters (micromolar)" legend_constants[469] = "Ki in component model_parameters (micromolar)" legend_constants[470] = "n_PU_base in component model_parameters (dimensionless)" legend_algebraic[226] = "I_iCa in component PU_membrane (picoampere)" legend_algebraic[230] = "I_iNa in component PU_membrane (picoampere)" legend_constants[717] = "alpha_scale in component model_parameters (dimensionless)" legend_algebraic[217] = "I_Ca in component I_Ca (picoampere)" legend_algebraic[220] = "I_Na in component I_Na (picoampere)" legend_algebraic[228] = "I_NaP in component I_NaP (picoampere)" legend_algebraic[222] = "I_NSCC_Ca in component I_NSCC_Ca (picoampere)" legend_algebraic[225] = "I_PM in component I_PM (picoampere)" legend_algebraic[224] = "I_NSCC_Na in component I_NSCC_Na (picoampere)" legend_constants[471] = "gCa_ in component I_Ca (picosiemens)" legend_constants[472] = "kCa in component I_Ca (per_millivolt)" legend_algebraic[211] = "gCa in component I_Ca (picosiemens)" legend_algebraic[214] = "ECa_PU in component I_Ca (millivolt)" legend_constants[473] = "T in component model_parameters (kelvin)" legend_constants[474] = "R in component model_parameters (attojoule_per_zeptomole_kelvin)" legend_constants[475] = "F in component model_parameters (femtocoulomb_per_zeptomole)" legend_constants[476] = "CO in component model_parameters (micromolar)" legend_constants[477] = "gNa in component I_Na (picosiemens)" legend_algebraic[219] = "ENa_PU in component I_Na (millivolt)" legend_constants[478] = "NO in component model_parameters (micromolar)" legend_states[63] = "NS1 in component NS1 (micromolar)" legend_constants[479] = "gNSCC_Ca_ in component I_NSCC_Ca (picosiemens)" legend_algebraic[221] = "gNSCC_Ca in component I_NSCC_Ca (picosiemens)" legend_constants[480] = "hNSCC in component model_parameters (dimensionless)" legend_constants[481] = "ENSCC in component model_parameters (millivolt)" legend_constants[482] = "KNSCC in component model_parameters (micromolar)" legend_constants[483] = "gNSCC_Na_ in component I_NSCC_Na (picosiemens)" legend_algebraic[223] = "gNSCC_Na in component I_NSCC_Na (picosiemens)" legend_constants[484] = "gPM in component I_PM (femtoampere)" legend_constants[485] = "KPM in component I_PM (micromolar)" legend_constants[486] = "gNaP in component I_NaP (picosiemens)" legend_constants[487] = "KNaP in component I_NaP (micromolar)" legend_constants[488] = "hNaP in component I_NaP (dimensionless)" legend_constants[489] = "ENaP in component I_NaP (millivolt)" legend_algebraic[227] = "JSERCA in component JSERCA (flux)" legend_constants[490] = "VSERCA in component JSERCA (first_order_rate_constant)" legend_constants[491] = "A2 in component JSERCA (dimensionless)" legend_constants[492] = "A4 in component JSERCA (per_micromolar)" legend_constants[493] = "A5 in component JSERCA (per_micromolar)" legend_constants[494] = "A6 in component JSERCA (per_micromolar2)" legend_states[64] = "CER in component CER (micromolar)" legend_algebraic[231] = "JMCU in component JMCU (flux)" legend_constants[495] = "VMCU in component JMCU (flux)" legend_constants[496] = "KMCU in component JMCU (micromolar)" legend_algebraic[229] = "epsilon_INH in component JMCU (dimensionless)" legend_constants[497] = "KINH in component JMCU (micromolar)" legend_constants[498] = "hINH in component JMCU (dimensionless)" legend_states[65] = "CMT in component CMT (micromolar)" legend_algebraic[239] = "JIPR in component JIPR (flux)" legend_constants[499] = "kIPR in component JIPR (first_order_rate_constant)" legend_constants[500] = "k_1 in component JIPR (flux)" legend_constants[501] = "k1 in component JIPR (first_order_rate_constant)" legend_constants[502] = "k2 in component JIPR (first_order_rate_constant)" legend_constants[503] = "r2 in component JIPR (first_order_rate_constant)" legend_constants[504] = "r_2 in component JIPR (flux)" legend_constants[505] = "r4 in component JIPR (first_order_rate_constant)" legend_constants[506] = "R1 in component JIPR (micromolar)" legend_constants[507] = "R3 in component JIPR (micromolar)" legend_algebraic[233] = "phi1 in component JIPR (first_order_rate_constant)" legend_algebraic[236] = "phi_1 in component JIPR (flux)" legend_algebraic[240] = "phi2 in component JIPR (first_order_rate_constant)" legend_algebraic[242] = "phi3 in component JIPR (first_order_rate_constant)" legend_states[66] = "H in component JIPR (dimensionless)" legend_constants[508] = "g_beta in component JIPR (first_order_rate_constant)" legend_constants[509] = "h_beta in component JIPR (dimensionless)" legend_constants[510] = "g_alpha in component JIPR (first_order_rate_constant)" legend_constants[511] = "K_beta in component JIPR (micromolar)" legend_constants[668] = "alpha_zeta in component JIPR (first_order_rate_constant)" legend_algebraic[14] = "beta_zeta in component JIPR (first_order_rate_constant)" legend_states[67] = "zeta in component JIPR (dimensionless)" legend_constants[512] = "K_phi3_act in component JIPR (micromolar)" legend_constants[513] = "K_phi3_inh in component JIPR (micromolar)" legend_constants[514] = "h_phi3_act in component JIPR (dimensionless)" legend_constants[515] = "h_phi3_inh in component JIPR (dimensionless)" legend_constants[516] = "g_phi3 in component JIPR (first_order_rate_constant)" legend_constants[517] = "P in component model_parameters (micromolar)" legend_algebraic[234] = "JNCX in component JNCX (flux)" legend_constants[518] = "VNCX in component JNCX (flux)" legend_constants[519] = "KNCX in component JNCX (micromolar)" legend_constants[712] = "lambda_MT_S1 in component CS1 (dimensionless)" legend_constants[713] = "lambda_ER_S1 in component CS1 (dimensionless)" legend_constants[520] = "delta_SPU in component model_parameters (micromolar_per_coulomb)" legend_constants[521] = "gamma_S1 in component model_parameters (dimensionless)" legend_constants[522] = "gamma_MT in component model_parameters (dimensionless)" legend_constants[523] = "gamma_ER in component model_parameters (dimensionless)" legend_constants[714] = "lambda_MT_S2 in component CS2 (dimensionless)" legend_constants[715] = "lambda_ER_S2 in component CS2 (dimensionless)" legend_constants[716] = "lambda_S1_S2 in component CS2 (dimensionless)" legend_constants[524] = "gamma_S2 in component model_parameters (dimensionless)" legend_algebraic[237] = "fm in component CMT (dimensionless)" legend_constants[525] = "Km in component CMT (micromolar)" legend_constants[526] = "Bm in component CMT (micromolar)" legend_constants[527] = "KO in component model_parameters (micromolar)" legend_constants[528] = "Ki in component model_parameters (micromolar)" legend_constants[529] = "n_PU_base in component model_parameters (dimensionless)" legend_algebraic[250] = "I_iCa in component PU_membrane (picoampere)" legend_algebraic[254] = "I_iNa in component PU_membrane (picoampere)" legend_constants[723] = "alpha_scale in component model_parameters (dimensionless)" legend_algebraic[241] = "I_Ca in component I_Ca (picoampere)" legend_algebraic[244] = "I_Na in component I_Na (picoampere)" legend_algebraic[252] = "I_NaP in component I_NaP (picoampere)" legend_algebraic[246] = "I_NSCC_Ca in component I_NSCC_Ca (picoampere)" legend_algebraic[249] = "I_PM in component I_PM (picoampere)" legend_algebraic[248] = "I_NSCC_Na in component I_NSCC_Na (picoampere)" legend_constants[530] = "gCa_ in component I_Ca (picosiemens)" legend_constants[531] = "kCa in component I_Ca (per_millivolt)" legend_algebraic[235] = "gCa in component I_Ca (picosiemens)" legend_algebraic[238] = "ECa_PU in component I_Ca (millivolt)" legend_constants[532] = "T in component model_parameters (kelvin)" legend_constants[533] = "R in component model_parameters (attojoule_per_zeptomole_kelvin)" legend_constants[534] = "F in component model_parameters (femtocoulomb_per_zeptomole)" legend_constants[535] = "CO in component model_parameters (micromolar)" legend_constants[536] = "gNa in component I_Na (picosiemens)" legend_algebraic[243] = "ENa_PU in component I_Na (millivolt)" legend_constants[537] = "NO in component model_parameters (micromolar)" legend_states[68] = "NS1 in component NS1 (micromolar)" legend_constants[538] = "gNSCC_Ca_ in component I_NSCC_Ca (picosiemens)" legend_algebraic[245] = "gNSCC_Ca in component I_NSCC_Ca (picosiemens)" legend_constants[539] = "hNSCC in component model_parameters (dimensionless)" legend_constants[540] = "ENSCC in component model_parameters (millivolt)" legend_constants[541] = "KNSCC in component model_parameters (micromolar)" legend_constants[542] = "gNSCC_Na_ in component I_NSCC_Na (picosiemens)" legend_algebraic[247] = "gNSCC_Na in component I_NSCC_Na (picosiemens)" legend_constants[543] = "gPM in component I_PM (femtoampere)" legend_constants[544] = "KPM in component I_PM (micromolar)" legend_constants[545] = "gNaP in component I_NaP (picosiemens)" legend_constants[546] = "KNaP in component I_NaP (micromolar)" legend_constants[547] = "hNaP in component I_NaP (dimensionless)" legend_constants[548] = "ENaP in component I_NaP (millivolt)" legend_algebraic[251] = "JSERCA in component JSERCA (flux)" legend_constants[549] = "VSERCA in component JSERCA (first_order_rate_constant)" legend_constants[550] = "A2 in component JSERCA (dimensionless)" legend_constants[551] = "A4 in component JSERCA (per_micromolar)" legend_constants[552] = "A5 in component JSERCA (per_micromolar)" legend_constants[553] = "A6 in component JSERCA (per_micromolar2)" legend_states[69] = "CER in component CER (micromolar)" legend_algebraic[255] = "JMCU in component JMCU (flux)" legend_constants[554] = "VMCU in component JMCU (flux)" legend_constants[555] = "KMCU in component JMCU (micromolar)" legend_algebraic[253] = "epsilon_INH in component JMCU (dimensionless)" legend_constants[556] = "KINH in component JMCU (micromolar)" legend_constants[557] = "hINH in component JMCU (dimensionless)" legend_states[70] = "CMT in component CMT (micromolar)" legend_algebraic[263] = "JIPR in component JIPR (flux)" legend_constants[558] = "kIPR in component JIPR (first_order_rate_constant)" legend_constants[559] = "k_1 in component JIPR (flux)" legend_constants[560] = "k1 in component JIPR (first_order_rate_constant)" legend_constants[561] = "k2 in component JIPR (first_order_rate_constant)" legend_constants[562] = "r2 in component JIPR (first_order_rate_constant)" legend_constants[563] = "r_2 in component JIPR (flux)" legend_constants[564] = "r4 in component JIPR (first_order_rate_constant)" legend_constants[565] = "R1 in component JIPR (micromolar)" legend_constants[566] = "R3 in component JIPR (micromolar)" legend_algebraic[257] = "phi1 in component JIPR (first_order_rate_constant)" legend_algebraic[260] = "phi_1 in component JIPR (flux)" legend_algebraic[264] = "phi2 in component JIPR (first_order_rate_constant)" legend_algebraic[266] = "phi3 in component JIPR (first_order_rate_constant)" legend_states[71] = "H in component JIPR (dimensionless)" legend_constants[567] = "g_beta in component JIPR (first_order_rate_constant)" legend_constants[568] = "h_beta in component JIPR (dimensionless)" legend_constants[569] = "g_alpha in component JIPR (first_order_rate_constant)" legend_constants[570] = "K_beta in component JIPR (micromolar)" legend_constants[669] = "alpha_zeta in component JIPR (first_order_rate_constant)" legend_algebraic[15] = "beta_zeta in component JIPR (first_order_rate_constant)" legend_states[72] = "zeta in component JIPR (dimensionless)" legend_constants[571] = "K_phi3_act in component JIPR (micromolar)" legend_constants[572] = "K_phi3_inh in component JIPR (micromolar)" legend_constants[573] = "h_phi3_act in component JIPR (dimensionless)" legend_constants[574] = "h_phi3_inh in component JIPR (dimensionless)" legend_constants[575] = "g_phi3 in component JIPR (first_order_rate_constant)" legend_constants[576] = "P in component model_parameters (micromolar)" legend_algebraic[258] = "JNCX in component JNCX (flux)" legend_constants[577] = "VNCX in component JNCX (flux)" legend_constants[578] = "KNCX in component JNCX (micromolar)" legend_constants[718] = "lambda_MT_S1 in component CS1 (dimensionless)" legend_constants[719] = "lambda_ER_S1 in component CS1 (dimensionless)" legend_constants[579] = "delta_SPU in component model_parameters (micromolar_per_coulomb)" legend_constants[580] = "gamma_S1 in component model_parameters (dimensionless)" legend_constants[581] = "gamma_MT in component model_parameters (dimensionless)" legend_constants[582] = "gamma_ER in component model_parameters (dimensionless)" legend_constants[720] = "lambda_MT_S2 in component CS2 (dimensionless)" legend_constants[721] = "lambda_ER_S2 in component CS2 (dimensionless)" legend_constants[722] = "lambda_S1_S2 in component CS2 (dimensionless)" legend_constants[583] = "gamma_S2 in component model_parameters (dimensionless)" legend_algebraic[261] = "fm in component CMT (dimensionless)" legend_constants[584] = "Km in component CMT (micromolar)" legend_constants[585] = "Bm in component CMT (micromolar)" legend_constants[586] = "KO in component model_parameters (micromolar)" legend_constants[587] = "Ki in component model_parameters (micromolar)" legend_constants[588] = "n_PU_base in component model_parameters (dimensionless)" legend_algebraic[274] = "I_iCa in component PU_membrane (picoampere)" legend_algebraic[277] = "I_iNa in component PU_membrane (picoampere)" legend_constants[729] = "alpha_scale in component model_parameters (dimensionless)" legend_algebraic[265] = "I_Ca in component I_Ca (picoampere)" legend_algebraic[268] = "I_Na in component I_Na (picoampere)" legend_algebraic[275] = "I_NaP in component I_NaP (picoampere)" legend_algebraic[270] = "I_NSCC_Ca in component I_NSCC_Ca (picoampere)" legend_algebraic[273] = "I_PM in component I_PM (picoampere)" legend_algebraic[272] = "I_NSCC_Na in component I_NSCC_Na (picoampere)" legend_constants[589] = "gCa_ in component I_Ca (picosiemens)" legend_constants[590] = "kCa in component I_Ca (per_millivolt)" legend_algebraic[259] = "gCa in component I_Ca (picosiemens)" legend_algebraic[262] = "ECa_PU in component I_Ca (millivolt)" legend_constants[591] = "T in component model_parameters (kelvin)" legend_constants[592] = "R in component model_parameters (attojoule_per_zeptomole_kelvin)" legend_constants[593] = "F in component model_parameters (femtocoulomb_per_zeptomole)" legend_constants[594] = "CO in component model_parameters (micromolar)" legend_constants[595] = "gNa in component I_Na (picosiemens)" legend_algebraic[267] = "ENa_PU in component I_Na (millivolt)" legend_constants[596] = "NO in component model_parameters (micromolar)" legend_states[73] = "NS1 in component NS1 (micromolar)" legend_constants[597] = "gNSCC_Ca_ in component I_NSCC_Ca (picosiemens)" legend_algebraic[269] = "gNSCC_Ca in component I_NSCC_Ca (picosiemens)" legend_constants[598] = "hNSCC in component model_parameters (dimensionless)" legend_constants[599] = "ENSCC in component model_parameters (millivolt)" legend_constants[600] = "KNSCC in component model_parameters (micromolar)" legend_constants[601] = "gNSCC_Na_ in component I_NSCC_Na (picosiemens)" legend_algebraic[271] = "gNSCC_Na in component I_NSCC_Na (picosiemens)" legend_constants[602] = "gPM in component I_PM (femtoampere)" legend_constants[603] = "KPM in component I_PM (micromolar)" legend_constants[604] = "gNaP in component I_NaP (picosiemens)" legend_constants[605] = "KNaP in component I_NaP (micromolar)" legend_constants[606] = "hNaP in component I_NaP (dimensionless)" legend_constants[607] = "ENaP in component I_NaP (millivolt)" legend_algebraic[276] = "JSERCA in component JSERCA (flux)" legend_constants[608] = "VSERCA in component JSERCA (first_order_rate_constant)" legend_constants[609] = "A2 in component JSERCA (dimensionless)" legend_constants[610] = "A4 in component JSERCA (per_micromolar)" legend_constants[611] = "A5 in component JSERCA (per_micromolar)" legend_constants[612] = "A6 in component JSERCA (per_micromolar2)" legend_states[74] = "CER in component CER (micromolar)" legend_algebraic[280] = "JMCU in component JMCU (flux)" legend_constants[613] = "VMCU in component JMCU (flux)" legend_constants[614] = "KMCU in component JMCU (micromolar)" legend_algebraic[278] = "epsilon_INH in component JMCU (dimensionless)" legend_constants[615] = "KINH in component JMCU (micromolar)" legend_constants[616] = "hINH in component JMCU (dimensionless)" legend_states[75] = "CMT in component CMT (micromolar)" legend_algebraic[285] = "JIPR in component JIPR (flux)" legend_constants[617] = "kIPR in component JIPR (first_order_rate_constant)" legend_constants[618] = "k_1 in component JIPR (flux)" legend_constants[619] = "k1 in component JIPR (first_order_rate_constant)" legend_constants[620] = "k2 in component JIPR (first_order_rate_constant)" legend_constants[621] = "r2 in component JIPR (first_order_rate_constant)" legend_constants[622] = "r_2 in component JIPR (flux)" legend_constants[623] = "r4 in component JIPR (first_order_rate_constant)" legend_constants[624] = "R1 in component JIPR (micromolar)" legend_constants[625] = "R3 in component JIPR (micromolar)" legend_algebraic[281] = "phi1 in component JIPR (first_order_rate_constant)" legend_algebraic[283] = "phi_1 in component JIPR (flux)" legend_algebraic[286] = "phi2 in component JIPR (first_order_rate_constant)" legend_algebraic[287] = "phi3 in component JIPR (first_order_rate_constant)" legend_states[76] = "H in component JIPR (dimensionless)" legend_constants[626] = "g_beta in component JIPR (first_order_rate_constant)" legend_constants[627] = "h_beta in component JIPR (dimensionless)" legend_constants[628] = "g_alpha in component JIPR (first_order_rate_constant)" legend_constants[629] = "K_beta in component JIPR (micromolar)" legend_constants[670] = "alpha_zeta in component JIPR (first_order_rate_constant)" legend_algebraic[16] = "beta_zeta in component JIPR (first_order_rate_constant)" legend_states[77] = "zeta in component JIPR (dimensionless)" legend_constants[630] = "K_phi3_act in component JIPR (micromolar)" legend_constants[631] = "K_phi3_inh in component JIPR (micromolar)" legend_constants[632] = "h_phi3_act in component JIPR (dimensionless)" legend_constants[633] = "h_phi3_inh in component JIPR (dimensionless)" legend_constants[634] = "g_phi3 in component JIPR (first_order_rate_constant)" legend_constants[635] = "P in component model_parameters (micromolar)" legend_algebraic[282] = "JNCX in component JNCX (flux)" legend_constants[636] = "VNCX in component JNCX (flux)" legend_constants[637] = "KNCX in component JNCX (micromolar)" legend_constants[724] = "lambda_MT_S1 in component CS1 (dimensionless)" legend_constants[725] = "lambda_ER_S1 in component CS1 (dimensionless)" legend_constants[638] = "delta_SPU in component model_parameters (micromolar_per_coulomb)" legend_constants[639] = "gamma_S1 in component model_parameters (dimensionless)" legend_constants[640] = "gamma_MT in component model_parameters (dimensionless)" legend_constants[641] = "gamma_ER in component model_parameters (dimensionless)" legend_constants[726] = "lambda_MT_S2 in component CS2 (dimensionless)" legend_constants[727] = "lambda_ER_S2 in component CS2 (dimensionless)" legend_constants[728] = "lambda_S1_S2 in component CS2 (dimensionless)" legend_constants[642] = "gamma_S2 in component model_parameters (dimensionless)" legend_algebraic[284] = "fm in component CMT (dimensionless)" legend_constants[643] = "Km in component CMT (micromolar)" legend_constants[644] = "Bm in component CMT (micromolar)" legend_constants[645] = "KO in component model_parameters (micromolar)" legend_constants[646] = "Ki in component model_parameters (micromolar)" legend_constants[647] = "n_PU_base in component model_parameters (dimensionless)" legend_rates[0] = "d/dt OCaT in component I_CaT (dimensionless)" legend_rates[2] = "d/dt dT in component I_CaT_d_gate (dimensionless)" legend_rates[3] = "d/dt fT in component I_CaT_f_gate (dimensionless)" legend_rates[5] = "d/dt dv1_1 in component I_Kv1_1_d_gate (dimensionless)" legend_rates[6] = "d/dt fv1_1 in component I_Kv1_1_f_gate (dimensionless)" legend_rates[7] = "d/dt dERG in component I_KERG_d_gate (dimensionless)" legend_rates[1] = "d/dt Vm in component membrane (millivolt)" legend_rates[4] = "d/dt CCy in component CCy (micromolar)" legend_rates[31] = "d/dt H in component JIPR (dimensionless)" legend_rates[32] = "d/dt zeta in component JIPR (dimensionless)" legend_rates[18] = "d/dt CS1 in component CS1 (micromolar)" legend_rates[8] = "d/dt CS2 in component CS2 (micromolar)" legend_rates[29] = "d/dt CER in component CER (micromolar)" legend_rates[30] = "d/dt CMT in component CMT (micromolar)" legend_rates[28] = "d/dt NS1 in component NS1 (micromolar)" legend_rates[36] = "d/dt H in component JIPR (dimensionless)" legend_rates[37] = "d/dt zeta in component JIPR (dimensionless)" legend_rates[19] = "d/dt CS1 in component CS1 (micromolar)" legend_rates[9] = "d/dt CS2 in component CS2 (micromolar)" legend_rates[34] = "d/dt CER in component CER (micromolar)" legend_rates[35] = "d/dt CMT in component CMT (micromolar)" legend_rates[33] = "d/dt NS1 in component NS1 (micromolar)" legend_rates[41] = "d/dt H in component JIPR (dimensionless)" legend_rates[42] = "d/dt zeta in component JIPR (dimensionless)" legend_rates[20] = "d/dt CS1 in component CS1 (micromolar)" legend_rates[10] = "d/dt CS2 in component CS2 (micromolar)" legend_rates[39] = "d/dt CER in component CER (micromolar)" legend_rates[40] = "d/dt CMT in component CMT (micromolar)" legend_rates[38] = "d/dt NS1 in component NS1 (micromolar)" legend_rates[46] = "d/dt H in component JIPR (dimensionless)" legend_rates[47] = "d/dt zeta in component JIPR (dimensionless)" legend_rates[21] = "d/dt CS1 in component CS1 (micromolar)" legend_rates[11] = "d/dt CS2 in component CS2 (micromolar)" legend_rates[44] = "d/dt CER in component CER (micromolar)" legend_rates[45] = "d/dt CMT in component CMT (micromolar)" legend_rates[43] = "d/dt NS1 in component NS1 (micromolar)" legend_rates[51] = "d/dt H in component JIPR (dimensionless)" legend_rates[52] = "d/dt zeta in component JIPR (dimensionless)" legend_rates[22] = "d/dt CS1 in component CS1 (micromolar)" legend_rates[12] = "d/dt CS2 in component CS2 (micromolar)" legend_rates[49] = "d/dt CER in component CER (micromolar)" legend_rates[50] = "d/dt CMT in component CMT (micromolar)" legend_rates[48] = "d/dt NS1 in component NS1 (micromolar)" legend_rates[56] = "d/dt H in component JIPR (dimensionless)" legend_rates[57] = "d/dt zeta in component JIPR (dimensionless)" legend_rates[23] = "d/dt CS1 in component CS1 (micromolar)" legend_rates[13] = "d/dt CS2 in component CS2 (micromolar)" legend_rates[54] = "d/dt CER in component CER (micromolar)" legend_rates[55] = "d/dt CMT in component CMT (micromolar)" legend_rates[53] = "d/dt NS1 in component NS1 (micromolar)" legend_rates[61] = "d/dt H in component JIPR (dimensionless)" legend_rates[62] = "d/dt zeta in component JIPR (dimensionless)" legend_rates[24] = "d/dt CS1 in component CS1 (micromolar)" legend_rates[14] = "d/dt CS2 in component CS2 (micromolar)" legend_rates[59] = "d/dt CER in component CER (micromolar)" legend_rates[60] = "d/dt CMT in component CMT (micromolar)" legend_rates[58] = "d/dt NS1 in component NS1 (micromolar)" legend_rates[66] = "d/dt H in component JIPR (dimensionless)" legend_rates[67] = "d/dt zeta in component JIPR (dimensionless)" legend_rates[25] = "d/dt CS1 in component CS1 (micromolar)" legend_rates[15] = "d/dt CS2 in component CS2 (micromolar)" legend_rates[64] = "d/dt CER in component CER (micromolar)" legend_rates[65] = "d/dt CMT in component CMT (micromolar)" legend_rates[63] = "d/dt NS1 in component NS1 (micromolar)" legend_rates[71] = "d/dt H in component JIPR (dimensionless)" legend_rates[72] = "d/dt zeta in component JIPR (dimensionless)" legend_rates[26] = "d/dt CS1 in component CS1 (micromolar)" legend_rates[16] = "d/dt CS2 in component CS2 (micromolar)" legend_rates[69] = "d/dt CER in component CER (micromolar)" legend_rates[70] = "d/dt CMT in component CMT (micromolar)" legend_rates[68] = "d/dt NS1 in component NS1 (micromolar)" legend_rates[76] = "d/dt H in component JIPR (dimensionless)" legend_rates[77] = "d/dt zeta in component JIPR (dimensionless)" legend_rates[27] = "d/dt CS1 in component CS1 (micromolar)" legend_rates[17] = "d/dt CS2 in component CS2 (micromolar)" legend_rates[74] = "d/dt CER in component CER (micromolar)" legend_rates[75] = "d/dt CMT in component CMT (micromolar)" legend_rates[73] = "d/dt NS1 in component NS1 (micromolar)" return (legend_states, legend_algebraic, legend_voi, legend_constants) def initConsts(): constants = [0.0] * sizeConstants; states = [0.0] * sizeStates; constants[0] = 800.0 constants[1] = 17.0 constants[2] = 240 constants[3] = 72 states[0] = 0 states[1] = -70.0 states[2] = 0.010 states[3] = 0.001 constants[4] = 0.0025 constants[5] = -0.60 constants[6] = -53 constants[7] = 0.019 constants[8] = 6.75 constants[9] = 2 constants[10] = -40 constants[11] = 1 constants[12] = -65 constants[13] = 100.0 constants[14] = 1.0 states[4] = 0.12 constants[15] = 10.0 states[5] = 0.000 states[6] = 1.0 constants[16] = 1000.0 constants[17] = 0.80 constants[18] = -0.13 constants[19] = 25 constants[20] = 333.0 constants[21] = 0.10 constants[22] = 0.23 constants[23] = 44.8 constants[24] = 6 states[7] = 0.000 constants[25] = 1000.0 constants[26] = 0.70 constants[27] = -0.56 constants[28] = 30 constants[29] = 13.5 constants[30] = -70.0 constants[31] = 0.80 constants[32] = 0.0 constants[33] = 20.0 constants[34] = 310.16 constants[35] = 8.314E-3 constants[36] = 0.09649 constants[37] = 1.8E3 constants[38] = 5400.0 constants[39] = 145E3 constants[40] = 1.0E-3 constants[41] = 0.002 constants[42] = 0.000 constants[43] = 0.3361 constants[44] = 0.3658 constants[45] = 0.4357 constants[46] = 0.5351 constants[47] = 0.6366 constants[48] = 0.7386 constants[49] = 0.8361 constants[50] = 0.9320 constants[51] = 1.0185 states[8] = 0.023 states[9] = 0.023 states[10] = 0.023 states[11] = 0.023 states[12] = 0.023 states[13] = 0.023 states[14] = 0.023 states[15] = 0.023 states[16] = 0.023 states[17] = 0.023 constants[52] = 0.30 constants[53] = 0.24 constants[54] = 10.0 states[18] = 0.120 states[19] = 0.120 states[20] = 0.120 states[21] = 0.120 states[22] = 0.120 states[23] = 0.120 states[24] = 0.120 states[25] = 0.120 states[26] = 0.120 states[27] = 0.120 constants[55] = 1.30 constants[56] = 0.12 constants[57] = 1.0 constants[58] = 0.074 constants[59] = 0.013 constants[60] = 310.16 constants[61] = 8.314E-3 constants[62] = 0.09649 constants[63] = 1.8E3 constants[64] = 13.5 constants[65] = 140.0E3 states[28] = 1.01E4 constants[66] = 0.10 constants[67] = 3.0 constants[68] = 0.0 constants[69] = 0.12 constants[70] = 160.0 constants[71] = 675.0 constants[72] = 1.0 constants[73] = 187.5 constants[74] = 1.0E4 constants[75] = 4.0 constants[76] = 10.0 constants[77] = 1.0E5 constants[78] = 6E-4 constants[79] = 3.57 constants[80] = 2.7E-5 constants[81] = 2.31E-5 states[29] = 200.0 constants[82] = 800.0 constants[83] = 10.0 constants[84] = 10.0 constants[85] = 4.0 states[30] = 0.200 constants[86] = 2000.0 constants[87] = 6.4 constants[88] = 0.0 constants[89] = 4.0 constants[90] = 250.0 constants[91] = 0.0 constants[92] = 750.0 constants[93] = 36.0 constants[94] = 300.0 states[31] = 0.200 constants[95] = 1500.0 constants[96] = 4.0 constants[97] = 0.85 constants[98] = 0.35 states[32] = 0.300 constants[99] = 0.1 constants[100] = 0.5 constants[101] = 3.0 constants[102] = 3.0 constants[103] = 4.5 constants[104] = 1.0 constants[105] = 3.0 constants[106] = 0.3 constants[107] = 18.5 constants[108] = 100.0 constants[109] = 200.0 constants[110] = 20.0 constants[111] = 1.0 constants[112] = 0.01 constants[113] = 100.0 constants[114] = 5400.0 constants[115] = 145E3 constants[116] = 50.0 constants[117] = 0.074 constants[118] = 0.013 constants[119] = 310.16 constants[120] = 8.314E-3 constants[121] = 0.09649 constants[122] = 1.8E3 constants[123] = 13.5 constants[124] = 140.0E3 states[33] = 1.01E4 constants[125] = 0.10 constants[126] = 3.0 constants[127] = 0.0 constants[128] = 0.12 constants[129] = 160.0 constants[130] = 675.0 constants[131] = 1.0 constants[132] = 187.5 constants[133] = 1.0E4 constants[134] = 4.0 constants[135] = 10.0 constants[136] = 1.0E5 constants[137] = 6E-4 constants[138] = 3.57 constants[139] = 2.7E-5 constants[140] = 2.31E-5 states[34] = 200.0 constants[141] = 800.0 constants[142] = 10.0 constants[143] = 10.0 constants[144] = 4.0 states[35] = 0.200 constants[145] = 2000.0 constants[146] = 6.4 constants[147] = 0.0 constants[148] = 4.0 constants[149] = 250.0 constants[150] = 0.0 constants[151] = 750.0 constants[152] = 36.0 constants[153] = 300.0 states[36] = 0.200 constants[154] = 1500.0 constants[155] = 4.0 constants[156] = 0.85 constants[157] = 0.35 states[37] = 0.300 constants[158] = 0.1 constants[159] = 0.5 constants[160] = 3.0 constants[161] = 3.0 constants[162] = 4.5 constants[163] = 1.0 constants[164] = 3.0 constants[165] = 0.3 constants[166] = 18.5 constants[167] = 100.0 constants[168] = 200.0 constants[169] = 20.0 constants[170] = 1.0 constants[171] = 0.01 constants[172] = 100.0 constants[173] = 5400.0 constants[174] = 145E3 constants[175] = 50.0 constants[176] = 0.074 constants[177] = 0.013 constants[178] = 310.16 constants[179] = 8.314E-3 constants[180] = 0.09649 constants[181] = 1.8E3 constants[182] = 13.5 constants[183] = 140.0E3 states[38] = 1.01E4 constants[184] = 0.10 constants[185] = 3.0 constants[186] = 0.0 constants[187] = 0.12 constants[188] = 160.0 constants[189] = 675.0 constants[190] = 1.0 constants[191] = 187.5 constants[192] = 1.0E4 constants[193] = 4.0 constants[194] = 10.0 constants[195] = 1.0E5 constants[196] = 6E-4 constants[197] = 3.57 constants[198] = 2.7E-5 constants[199] = 2.31E-5 states[39] = 200.0 constants[200] = 800.0 constants[201] = 10.0 constants[202] = 10.0 constants[203] = 4.0 states[40] = 0.200 constants[204] = 2000.0 constants[205] = 6.4 constants[206] = 0.0 constants[207] = 4.0 constants[208] = 250.0 constants[209] = 0.0 constants[210] = 750.0 constants[211] = 36.0 constants[212] = 300.0 states[41] = 0.200 constants[213] = 1500.0 constants[214] = 4.0 constants[215] = 0.85 constants[216] = 0.35 states[42] = 0.300 constants[217] = 0.1 constants[218] = 0.5 constants[219] = 3.0 constants[220] = 3.0 constants[221] = 4.5 constants[222] = 1.0 constants[223] = 3.0 constants[224] = 0.3 constants[225] = 18.5 constants[226] = 100.0 constants[227] = 200.0 constants[228] = 20.0 constants[229] = 1.0 constants[230] = 0.01 constants[231] = 100.0 constants[232] = 5400.0 constants[233] = 145E3 constants[234] = 50.0 constants[235] = 0.074 constants[236] = 0.013 constants[237] = 310.16 constants[238] = 8.314E-3 constants[239] = 0.09649 constants[240] = 1.8E3 constants[241] = 13.5 constants[242] = 140.0E3 states[43] = 1.01E4 constants[243] = 0.10 constants[244] = 3.0 constants[245] = 0.0 constants[246] = 0.12 constants[247] = 160.0 constants[248] = 675.0 constants[249] = 1.0 constants[250] = 187.5 constants[251] = 1.0E4 constants[252] = 4.0 constants[253] = 10.0 constants[254] = 1.0E5 constants[255] = 6E-4 constants[256] = 3.57 constants[257] = 2.7E-5 constants[258] = 2.31E-5 states[44] = 200.0 constants[259] = 800.0 constants[260] = 10.0 constants[261] = 10.0 constants[262] = 4.0 states[45] = 0.200 constants[263] = 2000.0 constants[264] = 6.4 constants[265] = 0.0 constants[266] = 4.0 constants[267] = 250.0 constants[268] = 0.0 constants[269] = 750.0 constants[270] = 36.0 constants[271] = 300.0 states[46] = 0.200 constants[272] = 1500.0 constants[273] = 4.0 constants[274] = 0.85 constants[275] = 0.35 states[47] = 0.300 constants[276] = 0.1 constants[277] = 0.5 constants[278] = 3.0 constants[279] = 3.0 constants[280] = 4.5 constants[281] = 1.0 constants[282] = 3.0 constants[283] = 0.3 constants[284] = 18.5 constants[285] = 100.0 constants[286] = 200.0 constants[287] = 20.0 constants[288] = 1.0 constants[289] = 0.01 constants[290] = 100.0 constants[291] = 5400.0 constants[292] = 145E3 constants[293] = 50.0 constants[294] = 0.074 constants[295] = 0.013 constants[296] = 310.16 constants[297] = 8.314E-3 constants[298] = 0.09649 constants[299] = 1.8E3 constants[300] = 13.5 constants[301] = 140.0E3 states[48] = 1.01E4 constants[302] = 0.10 constants[303] = 3.0 constants[304] = 0.0 constants[305] = 0.12 constants[306] = 160.0 constants[307] = 675.0 constants[308] = 1.0 constants[309] = 187.5 constants[310] = 1.0E4 constants[311] = 4.0 constants[312] = 10.0 constants[313] = 1.0E5 constants[314] = 6E-4 constants[315] = 3.57 constants[316] = 2.7E-5 constants[317] = 2.31E-5 states[49] = 200.0 constants[318] = 800.0 constants[319] = 10.0 constants[320] = 10.0 constants[321] = 4.0 states[50] = 0.200 constants[322] = 2000.0 constants[323] = 6.4 constants[324] = 0.0 constants[325] = 4.0 constants[326] = 250.0 constants[327] = 0.0 constants[328] = 750.0 constants[329] = 36.0 constants[330] = 300.0 states[51] = 0.200 constants[331] = 1500.0 constants[332] = 4.0 constants[333] = 0.85 constants[334] = 0.35 states[52] = 0.300 constants[335] = 0.1 constants[336] = 0.5 constants[337] = 3.0 constants[338] = 3.0 constants[339] = 4.5 constants[340] = 1.0 constants[341] = 3.0 constants[342] = 0.3 constants[343] = 18.5 constants[344] = 100.0 constants[345] = 200.0 constants[346] = 20.0 constants[347] = 1.0 constants[348] = 0.01 constants[349] = 100.0 constants[350] = 5400.0 constants[351] = 145E3 constants[352] = 50.0 constants[353] = 0.074 constants[354] = 0.013 constants[355] = 310.16 constants[356] = 8.314E-3 constants[357] = 0.09649 constants[358] = 1.8E3 constants[359] = 13.5 constants[360] = 140.0E3 states[53] = 1.01E4 constants[361] = 0.10 constants[362] = 3.0 constants[363] = 0.0 constants[364] = 0.12 constants[365] = 160.0 constants[366] = 675.0 constants[367] = 1.0 constants[368] = 187.5 constants[369] = 1.0E4 constants[370] = 4.0 constants[371] = 10.0 constants[372] = 1.0E5 constants[373] = 6E-4 constants[374] = 3.57 constants[375] = 2.7E-5 constants[376] = 2.31E-5 states[54] = 200.0 constants[377] = 800.0 constants[378] = 10.0 constants[379] = 10.0 constants[380] = 4.0 states[55] = 0.200 constants[381] = 2000.0 constants[382] = 6.4 constants[383] = 0.0 constants[384] = 4.0 constants[385] = 250.0 constants[386] = 0.0 constants[387] = 750.0 constants[388] = 36.0 constants[389] = 300.0 states[56] = 0.200 constants[390] = 1500.0 constants[391] = 4.0 constants[392] = 0.85 constants[393] = 0.35 states[57] = 0.300 constants[394] = 0.1 constants[395] = 0.5 constants[396] = 3.0 constants[397] = 3.0 constants[398] = 4.5 constants[399] = 1.0 constants[400] = 3.0 constants[401] = 0.3 constants[402] = 18.5 constants[403] = 100.0 constants[404] = 200.0 constants[405] = 20.0 constants[406] = 1.0 constants[407] = 0.01 constants[408] = 100.0 constants[409] = 5400.0 constants[410] = 145E3 constants[411] = 50.0 constants[412] = 0.074 constants[413] = 0.013 constants[414] = 310.16 constants[415] = 8.314E-3 constants[416] = 0.09649 constants[417] = 1.8E3 constants[418] = 13.5 constants[419] = 140.0E3 states[58] = 1.01E4 constants[420] = 0.10 constants[421] = 3.0 constants[422] = 0.0 constants[423] = 0.12 constants[424] = 160.0 constants[425] = 675.0 constants[426] = 1.0 constants[427] = 187.5 constants[428] = 1.0E4 constants[429] = 4.0 constants[430] = 10.0 constants[431] = 1.0E5 constants[432] = 6E-4 constants[433] = 3.57 constants[434] = 2.7E-5 constants[435] = 2.31E-5 states[59] = 200.0 constants[436] = 800.0 constants[437] = 10.0 constants[438] = 10.0 constants[439] = 4.0 states[60] = 0.200 constants[440] = 2000.0 constants[441] = 6.4 constants[442] = 0.0 constants[443] = 4.0 constants[444] = 250.0 constants[445] = 0.0 constants[446] = 750.0 constants[447] = 36.0 constants[448] = 300.0 states[61] = 0.200 constants[449] = 1500.0 constants[450] = 4.0 constants[451] = 0.85 constants[452] = 0.35 states[62] = 0.300 constants[453] = 0.1 constants[454] = 0.5 constants[455] = 3.0 constants[456] = 3.0 constants[457] = 4.5 constants[458] = 1.0 constants[459] = 3.0 constants[460] = 0.3 constants[461] = 18.5 constants[462] = 100.0 constants[463] = 200.0 constants[464] = 20.0 constants[465] = 1.0 constants[466] = 0.01 constants[467] = 100.0 constants[468] = 5400.0 constants[469] = 145E3 constants[470] = 50.0 constants[471] = 0.074 constants[472] = 0.013 constants[473] = 310.16 constants[474] = 8.314E-3 constants[475] = 0.09649 constants[476] = 1.8E3 constants[477] = 13.5 constants[478] = 140.0E3 states[63] = 1.01E4 constants[479] = 0.10 constants[480] = 3.0 constants[481] = 0.0 constants[482] = 0.12 constants[483] = 160.0 constants[484] = 675.0 constants[485] = 1.0 constants[486] = 187.5 constants[487] = 1.0E4 constants[488] = 4.0 constants[489] = 10.0 constants[490] = 1.0E5 constants[491] = 6E-4 constants[492] = 3.57 constants[493] = 2.7E-5 constants[494] = 2.31E-5 states[64] = 200.0 constants[495] = 800.0 constants[496] = 10.0 constants[497] = 10.0 constants[498] = 4.0 states[65] = 0.200 constants[499] = 2000.0 constants[500] = 6.4 constants[501] = 0.0 constants[502] = 4.0 constants[503] = 250.0 constants[504] = 0.0 constants[505] = 750.0 constants[506] = 36.0 constants[507] = 300.0 states[66] = 0.200 constants[508] = 1500.0 constants[509] = 4.0 constants[510] = 0.85 constants[511] = 0.35 states[67] = 0.300 constants[512] = 0.1 constants[513] = 0.5 constants[514] = 3.0 constants[515] = 3.0 constants[516] = 4.5 constants[517] = 1.0 constants[518] = 3.0 constants[519] = 0.3 constants[520] = 18.5 constants[521] = 100.0 constants[522] = 200.0 constants[523] = 20.0 constants[524] = 1.0 constants[525] = 0.01 constants[526] = 100.0 constants[527] = 5400.0 constants[528] = 145E3 constants[529] = 50.0 constants[530] = 0.074 constants[531] = 0.013 constants[532] = 310.16 constants[533] = 8.314E-3 constants[534] = 0.09649 constants[535] = 1.8E3 constants[536] = 13.5 constants[537] = 140.0E3 states[68] = 1.01E4 constants[538] = 0.10 constants[539] = 3.0 constants[540] = 0.0 constants[541] = 0.12 constants[542] = 160.0 constants[543] = 675.0 constants[544] = 1.0 constants[545] = 187.5 constants[546] = 1.0E4 constants[547] = 4.0 constants[548] = 10.0 constants[549] = 1.0E5 constants[550] = 6E-4 constants[551] = 3.57 constants[552] = 2.7E-5 constants[553] = 2.31E-5 states[69] = 200.0 constants[554] = 800.0 constants[555] = 10.0 constants[556] = 10.0 constants[557] = 4.0 states[70] = 0.200 constants[558] = 2000.0 constants[559] = 6.4 constants[560] = 0.0 constants[561] = 4.0 constants[562] = 250.0 constants[563] = 0.0 constants[564] = 750.0 constants[565] = 36.0 constants[566] = 300.0 states[71] = 0.200 constants[567] = 1500.0 constants[568] = 4.0 constants[569] = 0.85 constants[570] = 0.35 states[72] = 0.300 constants[571] = 0.1 constants[572] = 0.5 constants[573] = 3.0 constants[574] = 3.0 constants[575] = 4.5 constants[576] = 1.0 constants[577] = 3.0 constants[578] = 0.3 constants[579] = 18.5 constants[580] = 100.0 constants[581] = 200.0 constants[582] = 20.0 constants[583] = 1.0 constants[584] = 0.01 constants[585] = 100.0 constants[586] = 5400.0 constants[587] = 145E3 constants[588] = 50.0 constants[589] = 0.074 constants[590] = 0.013 constants[591] = 310.16 constants[592] = 8.314E-3 constants[593] = 0.09649 constants[594] = 1.8E3 constants[595] = 13.5 constants[596] = 140.0E3 states[73] = 1.01E4 constants[597] = 0.10 constants[598] = 3.0 constants[599] = 0.0 constants[600] = 0.12 constants[601] = 160.0 constants[602] = 675.0 constants[603] = 1.0 constants[604] = 187.5 constants[605] = 1.0E4 constants[606] = 4.0 constants[607] = 10.0 constants[608] = 1.0E5 constants[609] = 6E-4 constants[610] = 3.57 constants[611] = 2.7E-5 constants[612] = 2.31E-5 states[74] = 200.0 constants[613] = 800.0 constants[614] = 10.0 constants[615] = 10.0 constants[616] = 4.0 states[75] = 0.200 constants[617] = 2000.0 constants[618] = 6.4 constants[619] = 0.0 constants[620] = 4.0 constants[621] = 250.0 constants[622] = 0.0 constants[623] = 750.0 constants[624] = 36.0 constants[625] = 300.0 states[76] = 0.200 constants[626] = 1500.0 constants[627] = 4.0 constants[628] = 0.85 constants[629] = 0.35 states[77] = 0.300 constants[630] = 0.1 constants[631] = 0.5 constants[632] = 3.0 constants[633] = 3.0 constants[634] = 4.5 constants[635] = 1.0 constants[636] = 3.0 constants[637] = 0.3 constants[638] = 18.5 constants[639] = 100.0 constants[640] = 200.0 constants[641] = 20.0 constants[642] = 1.0 constants[643] = 0.01 constants[644] = 100.0 constants[645] = 5400.0 constants[646] = 145E3 constants[647] = 50.0 constants[648] = constants[4] constants[649] = ((constants[35]*constants[34])/constants[36])*log(constants[38]/constants[39]) constants[650] = constants[109]/constants[108] constants[651] = constants[52]+(constants[53]-constants[52])*((1.00000-1.00000)/(constants[54]-1.00000)) constants[652] = constants[52]+(constants[53]-constants[52])*((2.00000-1.00000)/(constants[54]-1.00000)) constants[653] = constants[52]+(constants[53]-constants[52])*((3.00000-1.00000)/(constants[54]-1.00000)) constants[654] = constants[52]+(constants[53]-constants[52])*((4.00000-1.00000)/(constants[54]-1.00000)) constants[655] = constants[52]+(constants[53]-constants[52])*((5.00000-1.00000)/(constants[54]-1.00000)) constants[656] = constants[52]+(constants[53]-constants[52])*((6.00000-1.00000)/(constants[54]-1.00000)) constants[657] = constants[52]+(constants[53]-constants[52])*((7.00000-1.00000)/(constants[54]-1.00000)) constants[658] = constants[52]+(constants[53]-constants[52])*((8.00000-1.00000)/(constants[54]-1.00000)) constants[659] = constants[52]+(constants[53]-constants[52])*((9.00000-1.00000)/(constants[54]-1.00000)) constants[660] = constants[52]+(constants[53]-constants[52])*((10.0000-1.00000)/(constants[54]-1.00000)) constants[661] = constants[97] constants[662] = constants[156] constants[663] = constants[215] constants[664] = constants[274] constants[665] = constants[333] constants[666] = constants[392] constants[667] = constants[451] constants[668] = constants[510] constants[669] = constants[569] constants[670] = constants[628] constants[671] = constants[110]/constants[108] constants[672] = constants[109]/constants[111] constants[673] = constants[110]/constants[111] constants[674] = constants[108]/constants[111] constants[675] = constants[116]/constants[54] constants[676] = constants[168]/constants[167] constants[677] = constants[169]/constants[167] constants[678] = constants[168]/constants[170] constants[679] = constants[169]/constants[170] constants[680] = constants[167]/constants[170] constants[681] = constants[175]/constants[54] constants[682] = constants[227]/constants[226] constants[683] = constants[228]/constants[226] constants[684] = constants[227]/constants[229] constants[685] = constants[228]/constants[229] constants[686] = constants[226]/constants[229] constants[687] = constants[234]/constants[54] constants[688] = constants[286]/constants[285] constants[689] = constants[287]/constants[285] constants[690] = constants[286]/constants[288] constants[691] = constants[287]/constants[288] constants[692] = constants[285]/constants[288] constants[693] = constants[293]/constants[54] constants[694] = constants[345]/constants[344] constants[695] = constants[346]/constants[344] constants[696] = constants[345]/constants[347] constants[697] = constants[346]/constants[347] constants[698] = constants[344]/constants[347] constants[699] = constants[352]/constants[54] constants[700] = constants[404]/constants[403] constants[701] = constants[405]/constants[403] constants[702] = constants[404]/constants[406] constants[703] = constants[405]/constants[406] constants[704] = constants[403]/constants[406] constants[705] = constants[411]/constants[54] constants[706] = constants[463]/constants[462] constants[707] = constants[464]/constants[462] constants[708] = constants[463]/constants[465] constants[709] = constants[464]/constants[465] constants[710] = constants[462]/constants[465] constants[711] = constants[470]/constants[54] constants[712] = constants[522]/constants[521] constants[713] = constants[523]/constants[521] constants[714] = constants[522]/constants[524] constants[715] = constants[523]/constants[524] constants[716] = constants[521]/constants[524] constants[717] = constants[529]/constants[54] constants[718] = constants[581]/constants[580] constants[719] = constants[582]/constants[580] constants[720] = constants[581]/constants[583] constants[721] = constants[582]/constants[583] constants[722] = constants[580]/constants[583] constants[723] = constants[588]/constants[54] constants[724] = constants[640]/constants[639] constants[725] = constants[641]/constants[639] constants[726] = constants[640]/constants[642] constants[727] = constants[641]/constants[642] constants[728] = constants[639]/constants[642] constants[729] = constants[647]/constants[54] return (states, constants) def computeRates(voi, states, constants): rates = [0.0] * sizeStates; algebraic = [0.0] * sizeAlgebraic rates[0] = constants[2]*states[2]*states[3]-constants[3]*states[0] algebraic[0] = 1.00000/(1.00000+exp(constants[5]*(states[1]-constants[6]))) rates[2] = (algebraic[0]-states[2])/constants[648] algebraic[7] = constants[95]/(1.00000+power(constants[98]/states[8], constants[96])) rates[32] = constants[661]*(1.00000-states[32])-algebraic[7]*states[32] algebraic[8] = constants[154]/(1.00000+power(constants[157]/states[9], constants[155])) rates[37] = constants[662]*(1.00000-states[37])-algebraic[8]*states[37] algebraic[9] = constants[213]/(1.00000+power(constants[216]/states[10], constants[214])) rates[42] = constants[663]*(1.00000-states[42])-algebraic[9]*states[42] algebraic[10] = constants[272]/(1.00000+power(constants[275]/states[11], constants[273])) rates[47] = constants[664]*(1.00000-states[47])-algebraic[10]*states[47] algebraic[11] = constants[331]/(1.00000+power(constants[334]/states[12], constants[332])) rates[52] = constants[665]*(1.00000-states[52])-algebraic[11]*states[52] algebraic[12] = constants[390]/(1.00000+power(constants[393]/states[13], constants[391])) rates[57] = constants[666]*(1.00000-states[57])-algebraic[12]*states[57] algebraic[13] = constants[449]/(1.00000+power(constants[452]/states[14], constants[450])) rates[62] = constants[667]*(1.00000-states[62])-algebraic[13]*states[62] algebraic[14] = constants[508]/(1.00000+power(constants[511]/states[15], constants[509])) rates[67] = constants[668]*(1.00000-states[67])-algebraic[14]*states[67] algebraic[15] = constants[567]/(1.00000+power(constants[570]/states[16], constants[568])) rates[72] = constants[669]*(1.00000-states[72])-algebraic[15]*states[72] algebraic[16] = constants[626]/(1.00000+power(constants[629]/states[17], constants[627])) rates[77] = constants[670]*(1.00000-states[77])-algebraic[16]*states[77] algebraic[1] = 1.00000/(1.00000+exp(constants[11]*(states[1]-constants[12]))) algebraic[17] = constants[7]+(constants[8]-constants[7])/(1.00000+exp(constants[9]*(states[1]-constants[10]))) rates[3] = (algebraic[1]-states[3])/algebraic[17] algebraic[2] = constants[16]*(constants[17]/(1.00000+exp(constants[18]*(states[1]-constants[19])))+(1.00000-constants[17])) algebraic[18] = constants[16]*(constants[17]-constants[17]/(1.00000+exp(constants[18]*(states[1]-constants[19])))) rates[5] = algebraic[2]*(1.00000-states[5])-algebraic[18]*states[5] algebraic[3] = constants[20]*(constants[21]/(1.00000+exp(constants[22]*(states[1]-constants[23])))+(1.00000-constants[21])) algebraic[19] = constants[20]*(constants[21]-constants[21]/(1.00000+exp(constants[22]*(states[1]-constants[23])))) rates[6] = algebraic[3]*(1.00000-states[6])-algebraic[19]*states[6] algebraic[4] = constants[25]*(constants[26]/(1.00000+exp(constants[27]*(states[1]-constants[28])))+(1.00000-constants[26])) algebraic[20] = constants[25]*(constants[26]-constants[26]/(1.00000+exp(constants[27]*(states[1]-constants[28])))) rates[7] = algebraic[4]*(1.00000-states[7])-algebraic[20]*states[7] algebraic[5] = constants[0]*states[0]*(states[1]-constants[1]) algebraic[21] = constants[13]*(states[4]/(constants[14]+states[4])) algebraic[37] = constants[55]*(constants[56]-states[4]) algebraic[27] = constants[42]*(states[4]-states[8]) algebraic[28] = constants[43]*(states[4]-states[9]) algebraic[29] = constants[44]*(states[4]-states[10]) algebraic[30] = constants[45]*(states[4]-states[11]) algebraic[31] = constants[46]*(states[4]-states[12]) algebraic[32] = constants[47]*(states[4]-states[13]) algebraic[33] = constants[48]*(states[4]-states[14]) algebraic[34] = constants[49]*(states[4]-states[15]) algebraic[35] = constants[50]*(states[4]-states[16]) algebraic[36] = constants[51]*(states[4]-states[17]) rates[4] = algebraic[37]-(constants[40]*algebraic[27]+constants[40]*algebraic[28]+constants[40]*algebraic[29]+constants[40]*algebraic[30]+constants[40]*algebraic[31]+constants[40]*algebraic[32]+constants[40]*algebraic[33]+constants[40]*algebraic[34]+constants[40]*algebraic[35]+constants[40]*algebraic[36]+(constants[41]/2.00000)*(algebraic[5]+algebraic[21])) algebraic[51] = ((constants[61]*constants[60])/constants[62])*log(constants[65]/states[28]) algebraic[52] = constants[64]*(states[1]-algebraic[51]) algebraic[60] = constants[73]*((power(states[28], constants[75]))/(power(constants[74], constants[75])+power(states[28], constants[75])))*(constants[76]-states[1]) algebraic[55] = constants[70]*((power(constants[69], constants[67]))/(power(constants[69], constants[67])+power(states[18], constants[67]))) algebraic[56] = algebraic[55]*(states[1]-constants[68]) algebraic[62] = constants[675]*(algebraic[56]+algebraic[52]+algebraic[60]) rates[28] = -(constants[107]/(1.00000*constants[675]))*algebraic[62] algebraic[38] = constants[651]*(states[8]-states[18]) algebraic[48] = constants[58]*exp(constants[59]*states[1]) algebraic[49] = ((constants[61]*constants[60])/(2.00000*constants[62]))*log(constants[63]/states[18]) algebraic[50] = algebraic[48]*(states[1]-algebraic[49]) algebraic[53] = constants[66]*((power(constants[69], constants[67]))/(power(constants[69], constants[67])+power(states[18], constants[67]))) algebraic[54] = algebraic[53]*(states[1]-constants[68]) algebraic[57] = constants[71]*((power(states[18], 2.00000))/(power(constants[72], 2.00000)+power(states[18], 2.00000))) algebraic[58] = constants[675]*(algebraic[50]+algebraic[54]+algebraic[57]) algebraic[59] = (constants[77]*(states[18]-constants[78]*states[29]))/(1.00000+constants[79]*states[18]+constants[80]*states[29]+constants[81]*states[18]*states[29]) algebraic[66] = constants[105]*(states[30]/(states[30]+constants[106])) rates[18] = (algebraic[38]+constants[650]*algebraic[66])-((constants[107]/(constants[675]*2.00000))*algebraic[58]+constants[671]*algebraic[59]) algebraic[61] = (power(constants[84], constants[85]))/(power(constants[84], constants[85])+power(states[30], constants[85])) algebraic[63] = constants[82]*((power(states[8], 2.00000))/(power(constants[83], 2.00000)+power(states[8], 2.00000)))*algebraic[61] algebraic[69] = 1.00000/(1.00000+(constants[112]*constants[113])/(power(constants[112]+states[30], 2.00000))) rates[30] = algebraic[69]*(algebraic[63]-algebraic[66]) algebraic[65] = (constants[88]*constants[93]+constants[90]*states[8])/(constants[93]+states[8]) algebraic[68] = ((constants[87]+constants[91])*constants[94])/(constants[94]+states[8]) algebraic[71] = constants[86]*(power((constants[104]*algebraic[65]*states[31])/(constants[104]*algebraic[65]+algebraic[68]), 4.00000))*(states[29]-states[8]) rates[8] = (algebraic[27]+constants[673]*algebraic[71])-(constants[674]*algebraic[38]+constants[672]*algebraic[63]) rates[29] = algebraic[59]-algebraic[71] algebraic[72] = (constants[89]*constants[94]+constants[92]*states[8])/(constants[94]+states[8]) algebraic[74] = (constants[103]*states[32])/((1.00000+power(constants[99]/states[8], constants[101]))*(1.00000+power(states[8]/constants[100], constants[102]))) rates[31] = algebraic[74]*(1.00000-states[31])-((constants[104]*algebraic[65]*algebraic[72])/(constants[104]*algebraic[65]+algebraic[68]))*states[31] algebraic[75] = ((constants[120]*constants[119])/constants[121])*log(constants[124]/states[33]) algebraic[76] = constants[123]*(states[1]-algebraic[75]) algebraic[84] = constants[132]*((power(states[33], constants[134]))/(power(constants[133], constants[134])+power(states[33], constants[134])))*(constants[135]-states[1]) algebraic[79] = constants[129]*((power(constants[128], constants[126]))/(power(constants[128], constants[126])+power(states[19], constants[126]))) algebraic[80] = algebraic[79]*(states[1]-constants[127]) algebraic[86] = constants[681]*(algebraic[80]+algebraic[76]+algebraic[84]) rates[33] = -(constants[166]/(1.00000*constants[681]))*algebraic[86] algebraic[39] = constants[652]*(states[9]-states[19]) algebraic[67] = constants[117]*exp(constants[118]*states[1]) algebraic[70] = ((constants[120]*constants[119])/(2.00000*constants[121]))*log(constants[122]/states[19]) algebraic[73] = algebraic[67]*(states[1]-algebraic[70]) algebraic[77] = constants[125]*((power(constants[128], constants[126]))/(power(constants[128], constants[126])+power(states[19], constants[126]))) algebraic[78] = algebraic[77]*(states[1]-constants[127]) algebraic[81] = constants[130]*((power(states[19], 2.00000))/(power(constants[131], 2.00000)+power(states[19], 2.00000))) algebraic[82] = constants[681]*(algebraic[73]+algebraic[78]+algebraic[81]) algebraic[83] = (constants[136]*(states[19]-constants[137]*states[34]))/(1.00000+constants[138]*states[19]+constants[139]*states[34]+constants[140]*states[19]*states[34]) algebraic[90] = constants[164]*(states[35]/(states[35]+constants[165])) rates[19] = (algebraic[39]+constants[676]*algebraic[90])-((constants[166]/(constants[681]*2.00000))*algebraic[82]+constants[677]*algebraic[83]) algebraic[85] = (power(constants[143], constants[144]))/(power(constants[143], constants[144])+power(states[35], constants[144])) algebraic[87] = constants[141]*((power(states[9], 2.00000))/(power(constants[142], 2.00000)+power(states[9], 2.00000)))*algebraic[85] algebraic[93] = 1.00000/(1.00000+(constants[171]*constants[172])/(power(constants[171]+states[35], 2.00000))) rates[35] = algebraic[93]*(algebraic[87]-algebraic[90]) algebraic[89] = (constants[147]*constants[152]+constants[149]*states[9])/(constants[152]+states[9]) algebraic[92] = ((constants[146]+constants[150])*constants[153])/(constants[153]+states[9]) algebraic[95] = constants[145]*(power((constants[163]*algebraic[89]*states[36])/(constants[163]*algebraic[89]+algebraic[92]), 4.00000))*(states[34]-states[9]) rates[9] = (algebraic[28]+constants[679]*algebraic[95])-(constants[680]*algebraic[39]+constants[678]*algebraic[87]) rates[34] = algebraic[83]-algebraic[95] algebraic[96] = (constants[148]*constants[153]+constants[151]*states[9])/(constants[153]+states[9]) algebraic[98] = (constants[162]*states[37])/((1.00000+power(constants[158]/states[9], constants[160]))*(1.00000+power(states[9]/constants[159], constants[161]))) rates[36] = algebraic[98]*(1.00000-states[36])-((constants[163]*algebraic[89]*algebraic[96])/(constants[163]*algebraic[89]+algebraic[92]))*states[36] algebraic[99] = ((constants[179]*constants[178])/constants[180])*log(constants[183]/states[38]) algebraic[100] = constants[182]*(states[1]-algebraic[99]) algebraic[108] = constants[191]*((power(states[38], constants[193]))/(power(constants[192], constants[193])+power(states[38], constants[193])))*(constants[194]-states[1]) algebraic[103] = constants[188]*((power(constants[187], constants[185]))/(power(constants[187], constants[185])+power(states[20], constants[185]))) algebraic[104] = algebraic[103]*(states[1]-constants[186]) algebraic[110] = constants[687]*(algebraic[104]+algebraic[100]+algebraic[108]) rates[38] = -(constants[225]/(1.00000*constants[687]))*algebraic[110] algebraic[40] = constants[653]*(states[10]-states[20]) algebraic[91] = constants[176]*exp(constants[177]*states[1]) algebraic[94] = ((constants[179]*constants[178])/(2.00000*constants[180]))*log(constants[181]/states[20]) algebraic[97] = algebraic[91]*(states[1]-algebraic[94]) algebraic[101] = constants[184]*((power(constants[187], constants[185]))/(power(constants[187], constants[185])+power(states[20], constants[185]))) algebraic[102] = algebraic[101]*(states[1]-constants[186]) algebraic[105] = constants[189]*((power(states[20], 2.00000))/(power(constants[190], 2.00000)+power(states[20], 2.00000))) algebraic[106] = constants[687]*(algebraic[97]+algebraic[102]+algebraic[105]) algebraic[107] = (constants[195]*(states[20]-constants[196]*states[39]))/(1.00000+constants[197]*states[20]+constants[198]*states[39]+constants[199]*states[20]*states[39]) algebraic[114] = constants[223]*(states[40]/(states[40]+constants[224])) rates[20] = (algebraic[40]+constants[682]*algebraic[114])-((constants[225]/(constants[687]*2.00000))*algebraic[106]+constants[683]*algebraic[107]) algebraic[109] = (power(constants[202], constants[203]))/(power(constants[202], constants[203])+power(states[40], constants[203])) algebraic[111] = constants[200]*((power(states[10], 2.00000))/(power(constants[201], 2.00000)+power(states[10], 2.00000)))*algebraic[109] algebraic[117] = 1.00000/(1.00000+(constants[230]*constants[231])/(power(constants[230]+states[40], 2.00000))) rates[40] = algebraic[117]*(algebraic[111]-algebraic[114]) algebraic[113] = (constants[206]*constants[211]+constants[208]*states[10])/(constants[211]+states[10]) algebraic[116] = ((constants[205]+constants[209])*constants[212])/(constants[212]+states[10]) algebraic[119] = constants[204]*(power((constants[222]*algebraic[113]*states[41])/(constants[222]*algebraic[113]+algebraic[116]), 4.00000))*(states[39]-states[10]) rates[10] = (algebraic[29]+constants[685]*algebraic[119])-(constants[686]*algebraic[40]+constants[684]*algebraic[111]) rates[39] = algebraic[107]-algebraic[119] algebraic[120] = (constants[207]*constants[212]+constants[210]*states[10])/(constants[212]+states[10]) algebraic[122] = (constants[221]*states[42])/((1.00000+power(constants[217]/states[10], constants[219]))*(1.00000+power(states[10]/constants[218], constants[220]))) rates[41] = algebraic[122]*(1.00000-states[41])-((constants[222]*algebraic[113]*algebraic[120])/(constants[222]*algebraic[113]+algebraic[116]))*states[41] algebraic[123] = ((constants[238]*constants[237])/constants[239])*log(constants[242]/states[43]) algebraic[124] = constants[241]*(states[1]-algebraic[123]) algebraic[132] = constants[250]*((power(states[43], constants[252]))/(power(constants[251], constants[252])+power(states[43], constants[252])))*(constants[253]-states[1]) algebraic[127] = constants[247]*((power(constants[246], constants[244]))/(power(constants[246], constants[244])+power(states[21], constants[244]))) algebraic[128] = algebraic[127]*(states[1]-constants[245]) algebraic[134] = constants[693]*(algebraic[128]+algebraic[124]+algebraic[132]) rates[43] = -(constants[284]/(1.00000*constants[693]))*algebraic[134] algebraic[41] = constants[654]*(states[11]-states[21]) algebraic[115] = constants[235]*exp(constants[236]*states[1]) algebraic[118] = ((constants[238]*constants[237])/(2.00000*constants[239]))*log(constants[240]/states[21]) algebraic[121] = algebraic[115]*(states[1]-algebraic[118]) algebraic[125] = constants[243]*((power(constants[246], constants[244]))/(power(constants[246], constants[244])+power(states[21], constants[244]))) algebraic[126] = algebraic[125]*(states[1]-constants[245]) algebraic[129] = constants[248]*((power(states[21], 2.00000))/(power(constants[249], 2.00000)+power(states[21], 2.00000))) algebraic[130] = constants[693]*(algebraic[121]+algebraic[126]+algebraic[129]) algebraic[131] = (constants[254]*(states[21]-constants[255]*states[44]))/(1.00000+constants[256]*states[21]+constants[257]*states[44]+constants[258]*states[21]*states[44]) algebraic[138] = constants[282]*(states[45]/(states[45]+constants[283])) rates[21] = (algebraic[41]+constants[688]*algebraic[138])-((constants[284]/(constants[693]*2.00000))*algebraic[130]+constants[689]*algebraic[131]) algebraic[133] = (power(constants[261], constants[262]))/(power(constants[261], constants[262])+power(states[45], constants[262])) algebraic[135] = constants[259]*((power(states[11], 2.00000))/(power(constants[260], 2.00000)+power(states[11], 2.00000)))*algebraic[133] algebraic[141] = 1.00000/(1.00000+(constants[289]*constants[290])/(power(constants[289]+states[45], 2.00000))) rates[45] = algebraic[141]*(algebraic[135]-algebraic[138]) algebraic[137] = (constants[265]*constants[270]+constants[267]*states[11])/(constants[270]+states[11]) algebraic[140] = ((constants[264]+constants[268])*constants[271])/(constants[271]+states[11]) algebraic[143] = constants[263]*(power((constants[281]*algebraic[137]*states[46])/(constants[281]*algebraic[137]+algebraic[140]), 4.00000))*(states[44]-states[11]) rates[11] = (algebraic[30]+constants[691]*algebraic[143])-(constants[692]*algebraic[41]+constants[690]*algebraic[135]) rates[44] = algebraic[131]-algebraic[143] algebraic[144] = (constants[266]*constants[271]+constants[269]*states[11])/(constants[271]+states[11]) algebraic[146] = (constants[280]*states[47])/((1.00000+power(constants[276]/states[11], constants[278]))*(1.00000+power(states[11]/constants[277], constants[279]))) rates[46] = algebraic[146]*(1.00000-states[46])-((constants[281]*algebraic[137]*algebraic[144])/(constants[281]*algebraic[137]+algebraic[140]))*states[46] algebraic[147] = ((constants[297]*constants[296])/constants[298])*log(constants[301]/states[48]) algebraic[148] = constants[300]*(states[1]-algebraic[147]) algebraic[156] = constants[309]*((power(states[48], constants[311]))/(power(constants[310], constants[311])+power(states[48], constants[311])))*(constants[312]-states[1]) algebraic[151] = constants[306]*((power(constants[305], constants[303]))/(power(constants[305], constants[303])+power(states[22], constants[303]))) algebraic[152] = algebraic[151]*(states[1]-constants[304]) algebraic[158] = constants[699]*(algebraic[152]+algebraic[148]+algebraic[156]) rates[48] = -(constants[343]/(1.00000*constants[699]))*algebraic[158] algebraic[42] = constants[655]*(states[12]-states[22]) algebraic[139] = constants[294]*exp(constants[295]*states[1]) algebraic[142] = ((constants[297]*constants[296])/(2.00000*constants[298]))*log(constants[299]/states[22]) algebraic[145] = algebraic[139]*(states[1]-algebraic[142]) algebraic[149] = constants[302]*((power(constants[305], constants[303]))/(power(constants[305], constants[303])+power(states[22], constants[303]))) algebraic[150] = algebraic[149]*(states[1]-constants[304]) algebraic[153] = constants[307]*((power(states[22], 2.00000))/(power(constants[308], 2.00000)+power(states[22], 2.00000))) algebraic[154] = constants[699]*(algebraic[145]+algebraic[150]+algebraic[153]) algebraic[155] = (constants[313]*(states[22]-constants[314]*states[49]))/(1.00000+constants[315]*states[22]+constants[316]*states[49]+constants[317]*states[22]*states[49]) algebraic[162] = constants[341]*(states[50]/(states[50]+constants[342])) rates[22] = (algebraic[42]+constants[694]*algebraic[162])-((constants[343]/(constants[699]*2.00000))*algebraic[154]+constants[695]*algebraic[155]) algebraic[157] = (power(constants[320], constants[321]))/(power(constants[320], constants[321])+power(states[50], constants[321])) algebraic[159] = constants[318]*((power(states[12], 2.00000))/(power(constants[319], 2.00000)+power(states[12], 2.00000)))*algebraic[157] algebraic[165] = 1.00000/(1.00000+(constants[348]*constants[349])/(power(constants[348]+states[50], 2.00000))) rates[50] = algebraic[165]*(algebraic[159]-algebraic[162]) algebraic[161] = (constants[324]*constants[329]+constants[326]*states[12])/(constants[329]+states[12]) algebraic[164] = ((constants[323]+constants[327])*constants[330])/(constants[330]+states[12]) algebraic[167] = constants[322]*(power((constants[340]*algebraic[161]*states[51])/(constants[340]*algebraic[161]+algebraic[164]), 4.00000))*(states[49]-states[12]) rates[12] = (algebraic[31]+constants[697]*algebraic[167])-(constants[698]*algebraic[42]+constants[696]*algebraic[159]) rates[49] = algebraic[155]-algebraic[167] algebraic[168] = (constants[325]*constants[330]+constants[328]*states[12])/(constants[330]+states[12]) algebraic[170] = (constants[339]*states[52])/((1.00000+power(constants[335]/states[12], constants[337]))*(1.00000+power(states[12]/constants[336], constants[338]))) rates[51] = algebraic[170]*(1.00000-states[51])-((constants[340]*algebraic[161]*algebraic[168])/(constants[340]*algebraic[161]+algebraic[164]))*states[51] algebraic[171] = ((constants[356]*constants[355])/constants[357])*log(constants[360]/states[53]) algebraic[172] = constants[359]*(states[1]-algebraic[171]) algebraic[180] = constants[368]*((power(states[53], constants[370]))/(power(constants[369], constants[370])+power(states[53], constants[370])))*(constants[371]-states[1]) algebraic[175] = constants[365]*((power(constants[364], constants[362]))/(power(constants[364], constants[362])+power(states[23], constants[362]))) algebraic[176] = algebraic[175]*(states[1]-constants[363]) algebraic[182] = constants[705]*(algebraic[176]+algebraic[172]+algebraic[180]) rates[53] = -(constants[402]/(1.00000*constants[705]))*algebraic[182] algebraic[43] = constants[656]*(states[13]-states[23]) algebraic[163] = constants[353]*exp(constants[354]*states[1]) algebraic[166] = ((constants[356]*constants[355])/(2.00000*constants[357]))*log(constants[358]/states[23]) algebraic[169] = algebraic[163]*(states[1]-algebraic[166]) algebraic[173] = constants[361]*((power(constants[364], constants[362]))/(power(constants[364], constants[362])+power(states[23], constants[362]))) algebraic[174] = algebraic[173]*(states[1]-constants[363]) algebraic[177] = constants[366]*((power(states[23], 2.00000))/(power(constants[367], 2.00000)+power(states[23], 2.00000))) algebraic[178] = constants[705]*(algebraic[169]+algebraic[174]+algebraic[177]) algebraic[179] = (constants[372]*(states[23]-constants[373]*states[54]))/(1.00000+constants[374]*states[23]+constants[375]*states[54]+constants[376]*states[23]*states[54]) algebraic[186] = constants[400]*(states[55]/(states[55]+constants[401])) rates[23] = (algebraic[43]+constants[700]*algebraic[186])-((constants[402]/(constants[705]*2.00000))*algebraic[178]+constants[701]*algebraic[179]) algebraic[181] = (power(constants[379], constants[380]))/(power(constants[379], constants[380])+power(states[55], constants[380])) algebraic[183] = constants[377]*((power(states[13], 2.00000))/(power(constants[378], 2.00000)+power(states[13], 2.00000)))*algebraic[181] algebraic[189] = 1.00000/(1.00000+(constants[407]*constants[408])/(power(constants[407]+states[55], 2.00000))) rates[55] = algebraic[189]*(algebraic[183]-algebraic[186]) algebraic[185] = (constants[383]*constants[388]+constants[385]*states[13])/(constants[388]+states[13]) algebraic[188] = ((constants[382]+constants[386])*constants[389])/(constants[389]+states[13]) algebraic[191] = constants[381]*(power((constants[399]*algebraic[185]*states[56])/(constants[399]*algebraic[185]+algebraic[188]), 4.00000))*(states[54]-states[13]) rates[13] = (algebraic[32]+constants[703]*algebraic[191])-(constants[704]*algebraic[43]+constants[702]*algebraic[183]) rates[54] = algebraic[179]-algebraic[191] algebraic[192] = (constants[384]*constants[389]+constants[387]*states[13])/(constants[389]+states[13]) algebraic[194] = (constants[398]*states[57])/((1.00000+power(constants[394]/states[13], constants[396]))*(1.00000+power(states[13]/constants[395], constants[397]))) rates[56] = algebraic[194]*(1.00000-states[56])-((constants[399]*algebraic[185]*algebraic[192])/(constants[399]*algebraic[185]+algebraic[188]))*states[56] algebraic[195] = ((constants[415]*constants[414])/constants[416])*log(constants[419]/states[58]) algebraic[196] = constants[418]*(states[1]-algebraic[195]) algebraic[204] = constants[427]*((power(states[58], constants[429]))/(power(constants[428], constants[429])+power(states[58], constants[429])))*(constants[430]-states[1]) algebraic[199] = constants[424]*((power(constants[423], constants[421]))/(power(constants[423], constants[421])+power(states[24], constants[421]))) algebraic[200] = algebraic[199]*(states[1]-constants[422]) algebraic[206] = constants[711]*(algebraic[200]+algebraic[196]+algebraic[204]) rates[58] = -(constants[461]/(1.00000*constants[711]))*algebraic[206] algebraic[44] = constants[657]*(states[14]-states[24]) algebraic[187] = constants[412]*exp(constants[413]*states[1]) algebraic[190] = ((constants[415]*constants[414])/(2.00000*constants[416]))*log(constants[417]/states[24]) algebraic[193] = algebraic[187]*(states[1]-algebraic[190]) algebraic[197] = constants[420]*((power(constants[423], constants[421]))/(power(constants[423], constants[421])+power(states[24], constants[421]))) algebraic[198] = algebraic[197]*(states[1]-constants[422]) algebraic[201] = constants[425]*((power(states[24], 2.00000))/(power(constants[426], 2.00000)+power(states[24], 2.00000))) algebraic[202] = constants[711]*(algebraic[193]+algebraic[198]+algebraic[201]) algebraic[203] = (constants[431]*(states[24]-constants[432]*states[59]))/(1.00000+constants[433]*states[24]+constants[434]*states[59]+constants[435]*states[24]*states[59]) algebraic[210] = constants[459]*(states[60]/(states[60]+constants[460])) rates[24] = (algebraic[44]+constants[706]*algebraic[210])-((constants[461]/(constants[711]*2.00000))*algebraic[202]+constants[707]*algebraic[203]) algebraic[205] = (power(constants[438], constants[439]))/(power(constants[438], constants[439])+power(states[60], constants[439])) algebraic[207] = constants[436]*((power(states[14], 2.00000))/(power(constants[437], 2.00000)+power(states[14], 2.00000)))*algebraic[205] algebraic[213] = 1.00000/(1.00000+(constants[466]*constants[467])/(power(constants[466]+states[60], 2.00000))) rates[60] = algebraic[213]*(algebraic[207]-algebraic[210]) algebraic[209] = (constants[442]*constants[447]+constants[444]*states[14])/(constants[447]+states[14]) algebraic[212] = ((constants[441]+constants[445])*constants[448])/(constants[448]+states[14]) algebraic[215] = constants[440]*(power((constants[458]*algebraic[209]*states[61])/(constants[458]*algebraic[209]+algebraic[212]), 4.00000))*(states[59]-states[14]) rates[14] = (algebraic[33]+constants[709]*algebraic[215])-(constants[710]*algebraic[44]+constants[708]*algebraic[207]) rates[59] = algebraic[203]-algebraic[215] algebraic[216] = (constants[443]*constants[448]+constants[446]*states[14])/(constants[448]+states[14]) algebraic[218] = (constants[457]*states[62])/((1.00000+power(constants[453]/states[14], constants[455]))*(1.00000+power(states[14]/constants[454], constants[456]))) rates[61] = algebraic[218]*(1.00000-states[61])-((constants[458]*algebraic[209]*algebraic[216])/(constants[458]*algebraic[209]+algebraic[212]))*states[61] algebraic[219] = ((constants[474]*constants[473])/constants[475])*log(constants[478]/states[63]) algebraic[220] = constants[477]*(states[1]-algebraic[219]) algebraic[228] = constants[486]*((power(states[63], constants[488]))/(power(constants[487], constants[488])+power(states[63], constants[488])))*(constants[489]-states[1]) algebraic[223] = constants[483]*((power(constants[482], constants[480]))/(power(constants[482], constants[480])+power(states[25], constants[480]))) algebraic[224] = algebraic[223]*(states[1]-constants[481]) algebraic[230] = constants[717]*(algebraic[224]+algebraic[220]+algebraic[228]) rates[63] = -(constants[520]/(1.00000*constants[717]))*algebraic[230] algebraic[45] = constants[658]*(states[15]-states[25]) algebraic[211] = constants[471]*exp(constants[472]*states[1]) algebraic[214] = ((constants[474]*constants[473])/(2.00000*constants[475]))*log(constants[476]/states[25]) algebraic[217] = algebraic[211]*(states[1]-algebraic[214]) algebraic[221] = constants[479]*((power(constants[482], constants[480]))/(power(constants[482], constants[480])+power(states[25], constants[480]))) algebraic[222] = algebraic[221]*(states[1]-constants[481]) algebraic[225] = constants[484]*((power(states[25], 2.00000))/(power(constants[485], 2.00000)+power(states[25], 2.00000))) algebraic[226] = constants[717]*(algebraic[217]+algebraic[222]+algebraic[225]) algebraic[227] = (constants[490]*(states[25]-constants[491]*states[64]))/(1.00000+constants[492]*states[25]+constants[493]*states[64]+constants[494]*states[25]*states[64]) algebraic[234] = constants[518]*(states[65]/(states[65]+constants[519])) rates[25] = (algebraic[45]+constants[712]*algebraic[234])-((constants[520]/(constants[717]*2.00000))*algebraic[226]+constants[713]*algebraic[227]) algebraic[229] = (power(constants[497], constants[498]))/(power(constants[497], constants[498])+power(states[65], constants[498])) algebraic[231] = constants[495]*((power(states[15], 2.00000))/(power(constants[496], 2.00000)+power(states[15], 2.00000)))*algebraic[229] algebraic[237] = 1.00000/(1.00000+(constants[525]*constants[526])/(power(constants[525]+states[65], 2.00000))) rates[65] = algebraic[237]*(algebraic[231]-algebraic[234]) algebraic[233] = (constants[501]*constants[506]+constants[503]*states[15])/(constants[506]+states[15]) algebraic[236] = ((constants[500]+constants[504])*constants[507])/(constants[507]+states[15]) algebraic[239] = constants[499]*(power((constants[517]*algebraic[233]*states[66])/(constants[517]*algebraic[233]+algebraic[236]), 4.00000))*(states[64]-states[15]) rates[15] = (algebraic[34]+constants[715]*algebraic[239])-(constants[716]*algebraic[45]+constants[714]*algebraic[231]) rates[64] = algebraic[227]-algebraic[239] algebraic[240] = (constants[502]*constants[507]+constants[505]*states[15])/(constants[507]+states[15]) algebraic[242] = (constants[516]*states[67])/((1.00000+power(constants[512]/states[15], constants[514]))*(1.00000+power(states[15]/constants[513], constants[515]))) rates[66] = algebraic[242]*(1.00000-states[66])-((constants[517]*algebraic[233]*algebraic[240])/(constants[517]*algebraic[233]+algebraic[236]))*states[66] algebraic[243] = ((constants[533]*constants[532])/constants[534])*log(constants[537]/states[68]) algebraic[244] = constants[536]*(states[1]-algebraic[243]) algebraic[252] = constants[545]*((power(states[68], constants[547]))/(power(constants[546], constants[547])+power(states[68], constants[547])))*(constants[548]-states[1]) algebraic[247] = constants[542]*((power(constants[541], constants[539]))/(power(constants[541], constants[539])+power(states[26], constants[539]))) algebraic[248] = algebraic[247]*(states[1]-constants[540]) algebraic[254] = constants[723]*(algebraic[248]+algebraic[244]+algebraic[252]) rates[68] = -(constants[579]/(1.00000*constants[723]))*algebraic[254] algebraic[46] = constants[659]*(states[16]-states[26]) algebraic[235] = constants[530]*exp(constants[531]*states[1]) algebraic[238] = ((constants[533]*constants[532])/(2.00000*constants[534]))*log(constants[535]/states[26]) algebraic[241] = algebraic[235]*(states[1]-algebraic[238]) algebraic[245] = constants[538]*((power(constants[541], constants[539]))/(power(constants[541], constants[539])+power(states[26], constants[539]))) algebraic[246] = algebraic[245]*(states[1]-constants[540]) algebraic[249] = constants[543]*((power(states[26], 2.00000))/(power(constants[544], 2.00000)+power(states[26], 2.00000))) algebraic[250] = constants[723]*(algebraic[241]+algebraic[246]+algebraic[249]) algebraic[251] = (constants[549]*(states[26]-constants[550]*states[69]))/(1.00000+constants[551]*states[26]+constants[552]*states[69]+constants[553]*states[26]*states[69]) algebraic[258] = constants[577]*(states[70]/(states[70]+constants[578])) rates[26] = (algebraic[46]+constants[718]*algebraic[258])-((constants[579]/(constants[723]*2.00000))*algebraic[250]+constants[719]*algebraic[251]) algebraic[253] = (power(constants[556], constants[557]))/(power(constants[556], constants[557])+power(states[70], constants[557])) algebraic[255] = constants[554]*((power(states[16], 2.00000))/(power(constants[555], 2.00000)+power(states[16], 2.00000)))*algebraic[253] algebraic[261] = 1.00000/(1.00000+(constants[584]*constants[585])/(power(constants[584]+states[70], 2.00000))) rates[70] = algebraic[261]*(algebraic[255]-algebraic[258]) algebraic[257] = (constants[560]*constants[565]+constants[562]*states[16])/(constants[565]+states[16]) algebraic[260] = ((constants[559]+constants[563])*constants[566])/(constants[566]+states[16]) algebraic[263] = constants[558]*(power((constants[576]*algebraic[257]*states[71])/(constants[576]*algebraic[257]+algebraic[260]), 4.00000))*(states[69]-states[16]) rates[16] = (algebraic[35]+constants[721]*algebraic[263])-(constants[722]*algebraic[46]+constants[720]*algebraic[255]) rates[69] = algebraic[251]-algebraic[263] algebraic[264] = (constants[561]*constants[566]+constants[564]*states[16])/(constants[566]+states[16]) algebraic[266] = (constants[575]*states[72])/((1.00000+power(constants[571]/states[16], constants[573]))*(1.00000+power(states[16]/constants[572], constants[574]))) rates[71] = algebraic[266]*(1.00000-states[71])-((constants[576]*algebraic[257]*algebraic[264])/(constants[576]*algebraic[257]+algebraic[260]))*states[71] algebraic[267] = ((constants[592]*constants[591])/constants[593])*log(constants[596]/states[73]) algebraic[268] = constants[595]*(states[1]-algebraic[267]) algebraic[275] = constants[604]*((power(states[73], constants[606]))/(power(constants[605], constants[606])+power(states[73], constants[606])))*(constants[607]-states[1]) algebraic[271] = constants[601]*((power(constants[600], constants[598]))/(power(constants[600], constants[598])+power(states[27], constants[598]))) algebraic[272] = algebraic[271]*(states[1]-constants[599]) algebraic[277] = constants[729]*(algebraic[272]+algebraic[268]+algebraic[275]) rates[73] = -(constants[638]/(1.00000*constants[729]))*algebraic[277] algebraic[22] = constants[15]*states[5]*states[6]*(states[1]-constants[649]) algebraic[23] = constants[24]*states[7]*(states[1]-constants[649]) algebraic[24] = constants[29]*(states[1]-constants[30]) algebraic[25] = constants[31]*(states[1]-constants[32]) algebraic[26] = algebraic[5]+algebraic[21]+algebraic[22]+algebraic[23]+algebraic[24]+algebraic[25] algebraic[64] = algebraic[58]+algebraic[62] algebraic[88] = algebraic[82]+algebraic[86] algebraic[112] = algebraic[106]+algebraic[110] algebraic[136] = algebraic[130]+algebraic[134] algebraic[160] = algebraic[154]+algebraic[158] algebraic[184] = algebraic[178]+algebraic[182] algebraic[208] = algebraic[202]+algebraic[206] algebraic[232] = algebraic[226]+algebraic[230] algebraic[256] = algebraic[250]+algebraic[254] algebraic[259] = constants[589]*exp(constants[590]*states[1]) algebraic[262] = ((constants[592]*constants[591])/(2.00000*constants[593]))*log(constants[594]/states[27]) algebraic[265] = algebraic[259]*(states[1]-algebraic[262]) algebraic[269] = constants[597]*((power(constants[600], constants[598]))/(power(constants[600], constants[598])+power(states[27], constants[598]))) algebraic[270] = algebraic[269]*(states[1]-constants[599]) algebraic[273] = constants[602]*((power(states[27], 2.00000))/(power(constants[603], 2.00000)+power(states[27], 2.00000))) algebraic[274] = constants[729]*(algebraic[265]+algebraic[270]+algebraic[273]) algebraic[279] = algebraic[274]+algebraic[277] rates[1] = -(1.00000/constants[33])*(algebraic[26]+algebraic[64]+algebraic[88]+algebraic[112]+algebraic[136]+algebraic[160]+algebraic[184]+algebraic[208]+algebraic[232]+algebraic[256]+algebraic[279]) algebraic[47] = constants[660]*(states[17]-states[27]) algebraic[276] = (constants[608]*(states[27]-constants[609]*states[74]))/(1.00000+constants[610]*states[27]+constants[611]*states[74]+constants[612]*states[27]*states[74]) algebraic[282] = constants[636]*(states[75]/(states[75]+constants[637])) rates[27] = (algebraic[47]+constants[724]*algebraic[282])-((constants[638]/(constants[729]*2.00000))*algebraic[274]+constants[725]*algebraic[276]) algebraic[278] = (power(constants[615], constants[616]))/(power(constants[615], constants[616])+power(states[75], constants[616])) algebraic[280] = constants[613]*((power(states[17], 2.00000))/(power(constants[614], 2.00000)+power(states[17], 2.00000)))*algebraic[278] algebraic[284] = 1.00000/(1.00000+(constants[643]*constants[644])/(power(constants[643]+states[75], 2.00000))) rates[75] = algebraic[284]*(algebraic[280]-algebraic[282]) algebraic[281] = (constants[619]*constants[624]+constants[621]*states[17])/(constants[624]+states[17]) algebraic[283] = ((constants[618]+constants[622])*constants[625])/(constants[625]+states[17]) algebraic[285] = constants[617]*(power((constants[635]*algebraic[281]*states[76])/(constants[635]*algebraic[281]+algebraic[283]), 4.00000))*(states[74]-states[17]) rates[17] = (algebraic[36]+constants[727]*algebraic[285])-(constants[728]*algebraic[47]+constants[726]*algebraic[280]) rates[74] = algebraic[276]-algebraic[285] algebraic[286] = (constants[620]*constants[625]+constants[623]*states[17])/(constants[625]+states[17]) algebraic[287] = (constants[634]*states[77])/((1.00000+power(constants[630]/states[17], constants[632]))*(1.00000+power(states[17]/constants[631], constants[633]))) rates[76] = algebraic[287]*(1.00000-states[76])-((constants[635]*algebraic[281]*algebraic[286])/(constants[635]*algebraic[281]+algebraic[283]))*states[76] return(rates) def computeAlgebraic(constants, states, voi): algebraic = array([[0.0] * len(voi)] * sizeAlgebraic) states = array(states) voi = array(voi) algebraic[0] = 1.00000/(1.00000+exp(constants[5]*(states[1]-constants[6]))) algebraic[7] = constants[95]/(1.00000+power(constants[98]/states[8], constants[96])) algebraic[8] = constants[154]/(1.00000+power(constants[157]/states[9], constants[155])) algebraic[9] = constants[213]/(1.00000+power(constants[216]/states[10], constants[214])) algebraic[10] = constants[272]/(1.00000+power(constants[275]/states[11], constants[273])) algebraic[11] = constants[331]/(1.00000+power(constants[334]/states[12], constants[332])) algebraic[12] = constants[390]/(1.00000+power(constants[393]/states[13], constants[391])) algebraic[13] = constants[449]/(1.00000+power(constants[452]/states[14], constants[450])) algebraic[14] = constants[508]/(1.00000+power(constants[511]/states[15], constants[509])) algebraic[15] = constants[567]/(1.00000+power(constants[570]/states[16], constants[568])) algebraic[16] = constants[626]/(1.00000+power(constants[629]/states[17], constants[627])) algebraic[1] = 1.00000/(1.00000+exp(constants[11]*(states[1]-constants[12]))) algebraic[17] = constants[7]+(constants[8]-constants[7])/(1.00000+exp(constants[9]*(states[1]-constants[10]))) algebraic[2] = constants[16]*(constants[17]/(1.00000+exp(constants[18]*(states[1]-constants[19])))+(1.00000-constants[17])) algebraic[18] = constants[16]*(constants[17]-constants[17]/(1.00000+exp(constants[18]*(states[1]-constants[19])))) algebraic[3] = constants[20]*(constants[21]/(1.00000+exp(constants[22]*(states[1]-constants[23])))+(1.00000-constants[21])) algebraic[19] = constants[20]*(constants[21]-constants[21]/(1.00000+exp(constants[22]*(states[1]-constants[23])))) algebraic[4] = constants[25]*(constants[26]/(1.00000+exp(constants[27]*(states[1]-constants[28])))+(1.00000-constants[26])) algebraic[20] = constants[25]*(constants[26]-constants[26]/(1.00000+exp(constants[27]*(states[1]-constants[28])))) algebraic[5] = constants[0]*states[0]*(states[1]-constants[1]) algebraic[21] = constants[13]*(states[4]/(constants[14]+states[4])) algebraic[37] = constants[55]*(constants[56]-states[4]) algebraic[27] = constants[42]*(states[4]-states[8]) algebraic[28] = constants[43]*(states[4]-states[9]) algebraic[29] = constants[44]*(states[4]-states[10]) algebraic[30] = constants[45]*(states[4]-states[11]) algebraic[31] = constants[46]*(states[4]-states[12]) algebraic[32] = constants[47]*(states[4]-states[13]) algebraic[33] = constants[48]*(states[4]-states[14]) algebraic[34] = constants[49]*(states[4]-states[15]) algebraic[35] = constants[50]*(states[4]-states[16]) algebraic[36] = constants[51]*(states[4]-states[17]) algebraic[51] = ((constants[61]*constants[60])/constants[62])*log(constants[65]/states[28]) algebraic[52] = constants[64]*(states[1]-algebraic[51]) algebraic[60] = constants[73]*((power(states[28], constants[75]))/(power(constants[74], constants[75])+power(states[28], constants[75])))*(constants[76]-states[1]) algebraic[55] = constants[70]*((power(constants[69], constants[67]))/(power(constants[69], constants[67])+power(states[18], constants[67]))) algebraic[56] = algebraic[55]*(states[1]-constants[68]) algebraic[62] = constants[675]*(algebraic[56]+algebraic[52]+algebraic[60]) algebraic[38] = constants[651]*(states[8]-states[18]) algebraic[48] = constants[58]*exp(constants[59]*states[1]) algebraic[49] = ((constants[61]*constants[60])/(2.00000*constants[62]))*log(constants[63]/states[18]) algebraic[50] = algebraic[48]*(states[1]-algebraic[49]) algebraic[53] = constants[66]*((power(constants[69], constants[67]))/(power(constants[69], constants[67])+power(states[18], constants[67]))) algebraic[54] = algebraic[53]*(states[1]-constants[68]) algebraic[57] = constants[71]*((power(states[18], 2.00000))/(power(constants[72], 2.00000)+power(states[18], 2.00000))) algebraic[58] = constants[675]*(algebraic[50]+algebraic[54]+algebraic[57]) algebraic[59] = (constants[77]*(states[18]-constants[78]*states[29]))/(1.00000+constants[79]*states[18]+constants[80]*states[29]+constants[81]*states[18]*states[29]) algebraic[66] = constants[105]*(states[30]/(states[30]+constants[106])) algebraic[61] = (power(constants[84], constants[85]))/(power(constants[84], constants[85])+power(states[30], constants[85])) algebraic[63] = constants[82]*((power(states[8], 2.00000))/(power(constants[83], 2.00000)+power(states[8], 2.00000)))*algebraic[61] algebraic[69] = 1.00000/(1.00000+(constants[112]*constants[113])/(power(constants[112]+states[30], 2.00000))) algebraic[65] = (constants[88]*constants[93]+constants[90]*states[8])/(constants[93]+states[8]) algebraic[68] = ((constants[87]+constants[91])*constants[94])/(constants[94]+states[8]) algebraic[71] = constants[86]*(power((constants[104]*algebraic[65]*states[31])/(constants[104]*algebraic[65]+algebraic[68]), 4.00000))*(states[29]-states[8]) algebraic[72] = (constants[89]*constants[94]+constants[92]*states[8])/(constants[94]+states[8]) algebraic[74] = (constants[103]*states[32])/((1.00000+power(constants[99]/states[8], constants[101]))*(1.00000+power(states[8]/constants[100], constants[102]))) algebraic[75] = ((constants[120]*constants[119])/constants[121])*log(constants[124]/states[33]) algebraic[76] = constants[123]*(states[1]-algebraic[75]) algebraic[84] = constants[132]*((power(states[33], constants[134]))/(power(constants[133], constants[134])+power(states[33], constants[134])))*(constants[135]-states[1]) algebraic[79] = constants[129]*((power(constants[128], constants[126]))/(power(constants[128], constants[126])+power(states[19], constants[126]))) algebraic[80] = algebraic[79]*(states[1]-constants[127]) algebraic[86] = constants[681]*(algebraic[80]+algebraic[76]+algebraic[84]) algebraic[39] = constants[652]*(states[9]-states[19]) algebraic[67] = constants[117]*exp(constants[118]*states[1]) algebraic[70] = ((constants[120]*constants[119])/(2.00000*constants[121]))*log(constants[122]/states[19]) algebraic[73] = algebraic[67]*(states[1]-algebraic[70]) algebraic[77] = constants[125]*((power(constants[128], constants[126]))/(power(constants[128], constants[126])+power(states[19], constants[126]))) algebraic[78] = algebraic[77]*(states[1]-constants[127]) algebraic[81] = constants[130]*((power(states[19], 2.00000))/(power(constants[131], 2.00000)+power(states[19], 2.00000))) algebraic[82] = constants[681]*(algebraic[73]+algebraic[78]+algebraic[81]) algebraic[83] = (constants[136]*(states[19]-constants[137]*states[34]))/(1.00000+constants[138]*states[19]+constants[139]*states[34]+constants[140]*states[19]*states[34]) algebraic[90] = constants[164]*(states[35]/(states[35]+constants[165])) algebraic[85] = (power(constants[143], constants[144]))/(power(constants[143], constants[144])+power(states[35], constants[144])) algebraic[87] = constants[141]*((power(states[9], 2.00000))/(power(constants[142], 2.00000)+power(states[9], 2.00000)))*algebraic[85] algebraic[93] = 1.00000/(1.00000+(constants[171]*constants[172])/(power(constants[171]+states[35], 2.00000))) algebraic[89] = (constants[147]*constants[152]+constants[149]*states[9])/(constants[152]+states[9]) algebraic[92] = ((constants[146]+constants[150])*constants[153])/(constants[153]+states[9]) algebraic[95] = constants[145]*(power((constants[163]*algebraic[89]*states[36])/(constants[163]*algebraic[89]+algebraic[92]), 4.00000))*(states[34]-states[9]) algebraic[96] = (constants[148]*constants[153]+constants[151]*states[9])/(constants[153]+states[9]) algebraic[98] = (constants[162]*states[37])/((1.00000+power(constants[158]/states[9], constants[160]))*(1.00000+power(states[9]/constants[159], constants[161]))) algebraic[99] = ((constants[179]*constants[178])/constants[180])*log(constants[183]/states[38]) algebraic[100] = constants[182]*(states[1]-algebraic[99]) algebraic[108] = constants[191]*((power(states[38], constants[193]))/(power(constants[192], constants[193])+power(states[38], constants[193])))*(constants[194]-states[1]) algebraic[103] = constants[188]*((power(constants[187], constants[185]))/(power(constants[187], constants[185])+power(states[20], constants[185]))) algebraic[104] = algebraic[103]*(states[1]-constants[186]) algebraic[110] = constants[687]*(algebraic[104]+algebraic[100]+algebraic[108]) algebraic[40] = constants[653]*(states[10]-states[20]) algebraic[91] = constants[176]*exp(constants[177]*states[1]) algebraic[94] = ((constants[179]*constants[178])/(2.00000*constants[180]))*log(constants[181]/states[20]) algebraic[97] = algebraic[91]*(states[1]-algebraic[94]) algebraic[101] = constants[184]*((power(constants[187], constants[185]))/(power(constants[187], constants[185])+power(states[20], constants[185]))) algebraic[102] = algebraic[101]*(states[1]-constants[186]) algebraic[105] = constants[189]*((power(states[20], 2.00000))/(power(constants[190], 2.00000)+power(states[20], 2.00000))) algebraic[106] = constants[687]*(algebraic[97]+algebraic[102]+algebraic[105]) algebraic[107] = (constants[195]*(states[20]-constants[196]*states[39]))/(1.00000+constants[197]*states[20]+constants[198]*states[39]+constants[199]*states[20]*states[39]) algebraic[114] = constants[223]*(states[40]/(states[40]+constants[224])) algebraic[109] = (power(constants[202], constants[203]))/(power(constants[202], constants[203])+power(states[40], constants[203])) algebraic[111] = constants[200]*((power(states[10], 2.00000))/(power(constants[201], 2.00000)+power(states[10], 2.00000)))*algebraic[109] algebraic[117] = 1.00000/(1.00000+(constants[230]*constants[231])/(power(constants[230]+states[40], 2.00000))) algebraic[113] = (constants[206]*constants[211]+constants[208]*states[10])/(constants[211]+states[10]) algebraic[116] = ((constants[205]+constants[209])*constants[212])/(constants[212]+states[10]) algebraic[119] = constants[204]*(power((constants[222]*algebraic[113]*states[41])/(constants[222]*algebraic[113]+algebraic[116]), 4.00000))*(states[39]-states[10]) algebraic[120] = (constants[207]*constants[212]+constants[210]*states[10])/(constants[212]+states[10]) algebraic[122] = (constants[221]*states[42])/((1.00000+power(constants[217]/states[10], constants[219]))*(1.00000+power(states[10]/constants[218], constants[220]))) algebraic[123] = ((constants[238]*constants[237])/constants[239])*log(constants[242]/states[43]) algebraic[124] = constants[241]*(states[1]-algebraic[123]) algebraic[132] = constants[250]*((power(states[43], constants[252]))/(power(constants[251], constants[252])+power(states[43], constants[252])))*(constants[253]-states[1]) algebraic[127] = constants[247]*((power(constants[246], constants[244]))/(power(constants[246], constants[244])+power(states[21], constants[244]))) algebraic[128] = algebraic[127]*(states[1]-constants[245]) algebraic[134] = constants[693]*(algebraic[128]+algebraic[124]+algebraic[132]) algebraic[41] = constants[654]*(states[11]-states[21]) algebraic[115] = constants[235]*exp(constants[236]*states[1]) algebraic[118] = ((constants[238]*constants[237])/(2.00000*constants[239]))*log(constants[240]/states[21]) algebraic[121] = algebraic[115]*(states[1]-algebraic[118]) algebraic[125] = constants[243]*((power(constants[246], constants[244]))/(power(constants[246], constants[244])+power(states[21], constants[244]))) algebraic[126] = algebraic[125]*(states[1]-constants[245]) algebraic[129] = constants[248]*((power(states[21], 2.00000))/(power(constants[249], 2.00000)+power(states[21], 2.00000))) algebraic[130] = constants[693]*(algebraic[121]+algebraic[126]+algebraic[129]) algebraic[131] = (constants[254]*(states[21]-constants[255]*states[44]))/(1.00000+constants[256]*states[21]+constants[257]*states[44]+constants[258]*states[21]*states[44]) algebraic[138] = constants[282]*(states[45]/(states[45]+constants[283])) algebraic[133] = (power(constants[261], constants[262]))/(power(constants[261], constants[262])+power(states[45], constants[262])) algebraic[135] = constants[259]*((power(states[11], 2.00000))/(power(constants[260], 2.00000)+power(states[11], 2.00000)))*algebraic[133] algebraic[141] = 1.00000/(1.00000+(constants[289]*constants[290])/(power(constants[289]+states[45], 2.00000))) algebraic[137] = (constants[265]*constants[270]+constants[267]*states[11])/(constants[270]+states[11]) algebraic[140] = ((constants[264]+constants[268])*constants[271])/(constants[271]+states[11]) algebraic[143] = constants[263]*(power((constants[281]*algebraic[137]*states[46])/(constants[281]*algebraic[137]+algebraic[140]), 4.00000))*(states[44]-states[11]) algebraic[144] = (constants[266]*constants[271]+constants[269]*states[11])/(constants[271]+states[11]) algebraic[146] = (constants[280]*states[47])/((1.00000+power(constants[276]/states[11], constants[278]))*(1.00000+power(states[11]/constants[277], constants[279]))) algebraic[147] = ((constants[297]*constants[296])/constants[298])*log(constants[301]/states[48]) algebraic[148] = constants[300]*(states[1]-algebraic[147]) algebraic[156] = constants[309]*((power(states[48], constants[311]))/(power(constants[310], constants[311])+power(states[48], constants[311])))*(constants[312]-states[1]) algebraic[151] = constants[306]*((power(constants[305], constants[303]))/(power(constants[305], constants[303])+power(states[22], constants[303]))) algebraic[152] = algebraic[151]*(states[1]-constants[304]) algebraic[158] = constants[699]*(algebraic[152]+algebraic[148]+algebraic[156]) algebraic[42] = constants[655]*(states[12]-states[22]) algebraic[139] = constants[294]*exp(constants[295]*states[1]) algebraic[142] = ((constants[297]*constants[296])/(2.00000*constants[298]))*log(constants[299]/states[22]) algebraic[145] = algebraic[139]*(states[1]-algebraic[142]) algebraic[149] = constants[302]*((power(constants[305], constants[303]))/(power(constants[305], constants[303])+power(states[22], constants[303]))) algebraic[150] = algebraic[149]*(states[1]-constants[304]) algebraic[153] = constants[307]*((power(states[22], 2.00000))/(power(constants[308], 2.00000)+power(states[22], 2.00000))) algebraic[154] = constants[699]*(algebraic[145]+algebraic[150]+algebraic[153]) algebraic[155] = (constants[313]*(states[22]-constants[314]*states[49]))/(1.00000+constants[315]*states[22]+constants[316]*states[49]+constants[317]*states[22]*states[49]) algebraic[162] = constants[341]*(states[50]/(states[50]+constants[342])) algebraic[157] = (power(constants[320], constants[321]))/(power(constants[320], constants[321])+power(states[50], constants[321])) algebraic[159] = constants[318]*((power(states[12], 2.00000))/(power(constants[319], 2.00000)+power(states[12], 2.00000)))*algebraic[157] algebraic[165] = 1.00000/(1.00000+(constants[348]*constants[349])/(power(constants[348]+states[50], 2.00000))) algebraic[161] = (constants[324]*constants[329]+constants[326]*states[12])/(constants[329]+states[12]) algebraic[164] = ((constants[323]+constants[327])*constants[330])/(constants[330]+states[12]) algebraic[167] = constants[322]*(power((constants[340]*algebraic[161]*states[51])/(constants[340]*algebraic[161]+algebraic[164]), 4.00000))*(states[49]-states[12]) algebraic[168] = (constants[325]*constants[330]+constants[328]*states[12])/(constants[330]+states[12]) algebraic[170] = (constants[339]*states[52])/((1.00000+power(constants[335]/states[12], constants[337]))*(1.00000+power(states[12]/constants[336], constants[338]))) algebraic[171] = ((constants[356]*constants[355])/constants[357])*log(constants[360]/states[53]) algebraic[172] = constants[359]*(states[1]-algebraic[171]) algebraic[180] = constants[368]*((power(states[53], constants[370]))/(power(constants[369], constants[370])+power(states[53], constants[370])))*(constants[371]-states[1]) algebraic[175] = constants[365]*((power(constants[364], constants[362]))/(power(constants[364], constants[362])+power(states[23], constants[362]))) algebraic[176] = algebraic[175]*(states[1]-constants[363]) algebraic[182] = constants[705]*(algebraic[176]+algebraic[172]+algebraic[180]) algebraic[43] = constants[656]*(states[13]-states[23]) algebraic[163] = constants[353]*exp(constants[354]*states[1]) algebraic[166] = ((constants[356]*constants[355])/(2.00000*constants[357]))*log(constants[358]/states[23]) algebraic[169] = algebraic[163]*(states[1]-algebraic[166]) algebraic[173] = constants[361]*((power(constants[364], constants[362]))/(power(constants[364], constants[362])+power(states[23], constants[362]))) algebraic[174] = algebraic[173]*(states[1]-constants[363]) algebraic[177] = constants[366]*((power(states[23], 2.00000))/(power(constants[367], 2.00000)+power(states[23], 2.00000))) algebraic[178] = constants[705]*(algebraic[169]+algebraic[174]+algebraic[177]) algebraic[179] = (constants[372]*(states[23]-constants[373]*states[54]))/(1.00000+constants[374]*states[23]+constants[375]*states[54]+constants[376]*states[23]*states[54]) algebraic[186] = constants[400]*(states[55]/(states[55]+constants[401])) algebraic[181] = (power(constants[379], constants[380]))/(power(constants[379], constants[380])+power(states[55], constants[380])) algebraic[183] = constants[377]*((power(states[13], 2.00000))/(power(constants[378], 2.00000)+power(states[13], 2.00000)))*algebraic[181] algebraic[189] = 1.00000/(1.00000+(constants[407]*constants[408])/(power(constants[407]+states[55], 2.00000))) algebraic[185] = (constants[383]*constants[388]+constants[385]*states[13])/(constants[388]+states[13]) algebraic[188] = ((constants[382]+constants[386])*constants[389])/(constants[389]+states[13]) algebraic[191] = constants[381]*(power((constants[399]*algebraic[185]*states[56])/(constants[399]*algebraic[185]+algebraic[188]), 4.00000))*(states[54]-states[13]) algebraic[192] = (constants[384]*constants[389]+constants[387]*states[13])/(constants[389]+states[13]) algebraic[194] = (constants[398]*states[57])/((1.00000+power(constants[394]/states[13], constants[396]))*(1.00000+power(states[13]/constants[395], constants[397]))) algebraic[195] = ((constants[415]*constants[414])/constants[416])*log(constants[419]/states[58]) algebraic[196] = constants[418]*(states[1]-algebraic[195]) algebraic[204] = constants[427]*((power(states[58], constants[429]))/(power(constants[428], constants[429])+power(states[58], constants[429])))*(constants[430]-states[1]) algebraic[199] = constants[424]*((power(constants[423], constants[421]))/(power(constants[423], constants[421])+power(states[24], constants[421]))) algebraic[200] = algebraic[199]*(states[1]-constants[422]) algebraic[206] = constants[711]*(algebraic[200]+algebraic[196]+algebraic[204]) algebraic[44] = constants[657]*(states[14]-states[24]) algebraic[187] = constants[412]*exp(constants[413]*states[1]) algebraic[190] = ((constants[415]*constants[414])/(2.00000*constants[416]))*log(constants[417]/states[24]) algebraic[193] = algebraic[187]*(states[1]-algebraic[190]) algebraic[197] = constants[420]*((power(constants[423], constants[421]))/(power(constants[423], constants[421])+power(states[24], constants[421]))) algebraic[198] = algebraic[197]*(states[1]-constants[422]) algebraic[201] = constants[425]*((power(states[24], 2.00000))/(power(constants[426], 2.00000)+power(states[24], 2.00000))) algebraic[202] = constants[711]*(algebraic[193]+algebraic[198]+algebraic[201]) algebraic[203] = (constants[431]*(states[24]-constants[432]*states[59]))/(1.00000+constants[433]*states[24]+constants[434]*states[59]+constants[435]*states[24]*states[59]) algebraic[210] = constants[459]*(states[60]/(states[60]+constants[460])) algebraic[205] = (power(constants[438], constants[439]))/(power(constants[438], constants[439])+power(states[60], constants[439])) algebraic[207] = constants[436]*((power(states[14], 2.00000))/(power(constants[437], 2.00000)+power(states[14], 2.00000)))*algebraic[205] algebraic[213] = 1.00000/(1.00000+(constants[466]*constants[467])/(power(constants[466]+states[60], 2.00000))) algebraic[209] = (constants[442]*constants[447]+constants[444]*states[14])/(constants[447]+states[14]) algebraic[212] = ((constants[441]+constants[445])*constants[448])/(constants[448]+states[14]) algebraic[215] = constants[440]*(power((constants[458]*algebraic[209]*states[61])/(constants[458]*algebraic[209]+algebraic[212]), 4.00000))*(states[59]-states[14]) algebraic[216] = (constants[443]*constants[448]+constants[446]*states[14])/(constants[448]+states[14]) algebraic[218] = (constants[457]*states[62])/((1.00000+power(constants[453]/states[14], constants[455]))*(1.00000+power(states[14]/constants[454], constants[456]))) algebraic[219] = ((constants[474]*constants[473])/constants[475])*log(constants[478]/states[63]) algebraic[220] = constants[477]*(states[1]-algebraic[219]) algebraic[228] = constants[486]*((power(states[63], constants[488]))/(power(constants[487], constants[488])+power(states[63], constants[488])))*(constants[489]-states[1]) algebraic[223] = constants[483]*((power(constants[482], constants[480]))/(power(constants[482], constants[480])+power(states[25], constants[480]))) algebraic[224] = algebraic[223]*(states[1]-constants[481]) algebraic[230] = constants[717]*(algebraic[224]+algebraic[220]+algebraic[228]) algebraic[45] = constants[658]*(states[15]-states[25]) algebraic[211] = constants[471]*exp(constants[472]*states[1]) algebraic[214] = ((constants[474]*constants[473])/(2.00000*constants[475]))*log(constants[476]/states[25]) algebraic[217] = algebraic[211]*(states[1]-algebraic[214]) algebraic[221] = constants[479]*((power(constants[482], constants[480]))/(power(constants[482], constants[480])+power(states[25], constants[480]))) algebraic[222] = algebraic[221]*(states[1]-constants[481]) algebraic[225] = constants[484]*((power(states[25], 2.00000))/(power(constants[485], 2.00000)+power(states[25], 2.00000))) algebraic[226] = constants[717]*(algebraic[217]+algebraic[222]+algebraic[225]) algebraic[227] = (constants[490]*(states[25]-constants[491]*states[64]))/(1.00000+constants[492]*states[25]+constants[493]*states[64]+constants[494]*states[25]*states[64]) algebraic[234] = constants[518]*(states[65]/(states[65]+constants[519])) algebraic[229] = (power(constants[497], constants[498]))/(power(constants[497], constants[498])+power(states[65], constants[498])) algebraic[231] = constants[495]*((power(states[15], 2.00000))/(power(constants[496], 2.00000)+power(states[15], 2.00000)))*algebraic[229] algebraic[237] = 1.00000/(1.00000+(constants[525]*constants[526])/(power(constants[525]+states[65], 2.00000))) algebraic[233] = (constants[501]*constants[506]+constants[503]*states[15])/(constants[506]+states[15]) algebraic[236] = ((constants[500]+constants[504])*constants[507])/(constants[507]+states[15]) algebraic[239] = constants[499]*(power((constants[517]*algebraic[233]*states[66])/(constants[517]*algebraic[233]+algebraic[236]), 4.00000))*(states[64]-states[15]) algebraic[240] = (constants[502]*constants[507]+constants[505]*states[15])/(constants[507]+states[15]) algebraic[242] = (constants[516]*states[67])/((1.00000+power(constants[512]/states[15], constants[514]))*(1.00000+power(states[15]/constants[513], constants[515]))) algebraic[243] = ((constants[533]*constants[532])/constants[534])*log(constants[537]/states[68]) algebraic[244] = constants[536]*(states[1]-algebraic[243]) algebraic[252] = constants[545]*((power(states[68], constants[547]))/(power(constants[546], constants[547])+power(states[68], constants[547])))*(constants[548]-states[1]) algebraic[247] = constants[542]*((power(constants[541], constants[539]))/(power(constants[541], constants[539])+power(states[26], constants[539]))) algebraic[248] = algebraic[247]*(states[1]-constants[540]) algebraic[254] = constants[723]*(algebraic[248]+algebraic[244]+algebraic[252]) algebraic[46] = constants[659]*(states[16]-states[26]) algebraic[235] = constants[530]*exp(constants[531]*states[1]) algebraic[238] = ((constants[533]*constants[532])/(2.00000*constants[534]))*log(constants[535]/states[26]) algebraic[241] = algebraic[235]*(states[1]-algebraic[238]) algebraic[245] = constants[538]*((power(constants[541], constants[539]))/(power(constants[541], constants[539])+power(states[26], constants[539]))) algebraic[246] = algebraic[245]*(states[1]-constants[540]) algebraic[249] = constants[543]*((power(states[26], 2.00000))/(power(constants[544], 2.00000)+power(states[26], 2.00000))) algebraic[250] = constants[723]*(algebraic[241]+algebraic[246]+algebraic[249]) algebraic[251] = (constants[549]*(states[26]-constants[550]*states[69]))/(1.00000+constants[551]*states[26]+constants[552]*states[69]+constants[553]*states[26]*states[69]) algebraic[258] = constants[577]*(states[70]/(states[70]+constants[578])) algebraic[253] = (power(constants[556], constants[557]))/(power(constants[556], constants[557])+power(states[70], constants[557])) algebraic[255] = constants[554]*((power(states[16], 2.00000))/(power(constants[555], 2.00000)+power(states[16], 2.00000)))*algebraic[253] algebraic[261] = 1.00000/(1.00000+(constants[584]*constants[585])/(power(constants[584]+states[70], 2.00000))) algebraic[257] = (constants[560]*constants[565]+constants[562]*states[16])/(constants[565]+states[16]) algebraic[260] = ((constants[559]+constants[563])*constants[566])/(constants[566]+states[16]) algebraic[263] = constants[558]*(power((constants[576]*algebraic[257]*states[71])/(constants[576]*algebraic[257]+algebraic[260]), 4.00000))*(states[69]-states[16]) algebraic[264] = (constants[561]*constants[566]+constants[564]*states[16])/(constants[566]+states[16]) algebraic[266] = (constants[575]*states[72])/((1.00000+power(constants[571]/states[16], constants[573]))*(1.00000+power(states[16]/constants[572], constants[574]))) algebraic[267] = ((constants[592]*constants[591])/constants[593])*log(constants[596]/states[73]) algebraic[268] = constants[595]*(states[1]-algebraic[267]) algebraic[275] = constants[604]*((power(states[73], constants[606]))/(power(constants[605], constants[606])+power(states[73], constants[606])))*(constants[607]-states[1]) algebraic[271] = constants[601]*((power(constants[600], constants[598]))/(power(constants[600], constants[598])+power(states[27], constants[598]))) algebraic[272] = algebraic[271]*(states[1]-constants[599]) algebraic[277] = constants[729]*(algebraic[272]+algebraic[268]+algebraic[275]) algebraic[22] = constants[15]*states[5]*states[6]*(states[1]-constants[649]) algebraic[23] = constants[24]*states[7]*(states[1]-constants[649]) algebraic[24] = constants[29]*(states[1]-constants[30]) algebraic[25] = constants[31]*(states[1]-constants[32]) algebraic[26] = algebraic[5]+algebraic[21]+algebraic[22]+algebraic[23]+algebraic[24]+algebraic[25] algebraic[64] = algebraic[58]+algebraic[62] algebraic[88] = algebraic[82]+algebraic[86] algebraic[112] = algebraic[106]+algebraic[110] algebraic[136] = algebraic[130]+algebraic[134] algebraic[160] = algebraic[154]+algebraic[158] algebraic[184] = algebraic[178]+algebraic[182] algebraic[208] = algebraic[202]+algebraic[206] algebraic[232] = algebraic[226]+algebraic[230] algebraic[256] = algebraic[250]+algebraic[254] algebraic[259] = constants[589]*exp(constants[590]*states[1]) algebraic[262] = ((constants[592]*constants[591])/(2.00000*constants[593]))*log(constants[594]/states[27]) algebraic[265] = algebraic[259]*(states[1]-algebraic[262]) algebraic[269] = constants[597]*((power(constants[600], constants[598]))/(power(constants[600], constants[598])+power(states[27], constants[598]))) algebraic[270] = algebraic[269]*(states[1]-constants[599]) algebraic[273] = constants[602]*((power(states[27], 2.00000))/(power(constants[603], 2.00000)+power(states[27], 2.00000))) algebraic[274] = constants[729]*(algebraic[265]+algebraic[270]+algebraic[273]) algebraic[279] = algebraic[274]+algebraic[277] algebraic[47] = constants[660]*(states[17]-states[27]) algebraic[276] = (constants[608]*(states[27]-constants[609]*states[74]))/(1.00000+constants[610]*states[27]+constants[611]*states[74]+constants[612]*states[27]*states[74]) algebraic[282] = constants[636]*(states[75]/(states[75]+constants[637])) algebraic[278] = (power(constants[615], constants[616]))/(power(constants[615], constants[616])+power(states[75], constants[616])) algebraic[280] = constants[613]*((power(states[17], 2.00000))/(power(constants[614], 2.00000)+power(states[17], 2.00000)))*algebraic[278] algebraic[284] = 1.00000/(1.00000+(constants[643]*constants[644])/(power(constants[643]+states[75], 2.00000))) algebraic[281] = (constants[619]*constants[624]+constants[621]*states[17])/(constants[624]+states[17]) algebraic[283] = ((constants[618]+constants[622])*constants[625])/(constants[625]+states[17]) algebraic[285] = constants[617]*(power((constants[635]*algebraic[281]*states[76])/(constants[635]*algebraic[281]+algebraic[283]), 4.00000))*(states[74]-states[17]) algebraic[286] = (constants[620]*constants[625]+constants[623]*states[17])/(constants[625]+states[17]) algebraic[287] = (constants[634]*states[77])/((1.00000+power(constants[630]/states[17], constants[632]))*(1.00000+power(states[17]/constants[631], constants[633]))) algebraic[6] = ((constants[35]*constants[34])/(2.00000*constants[36]))*log(constants[37]/states[4]) 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)