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