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