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