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