/* There are a total of 179 entries in the algebraic variable array. There are a total of 17 entries in each of the rate and state variable arrays. There are a total of 230 entries in the constant variable array. */ /* * VOI is time in component environment (day). * ALGEBRAIC[57] is U_Pd_PdPs in component Q_Pd (mole_per_day). * STATES[0] is Q_Pd in component Q_Pd (mole). * ALGEBRAIC[0] is C_Pd in component Q_Pd (molar). * ALGEBRAIC[1] is U_Pd_PdEx in component Q_Pd (mole_per_day). * CONSTANTS[155] is P_Pd_InPd in component Q_Pd (mole_per_day). * ALGEBRAIC[2] is v_PdPs in component Q_Pd (mole_per_day). * CONSTANTS[0] is v_PdPs_star in component Q_Pd (mole_per_gram_day). * CONSTANTS[174] is M_Pd_PdPs in component Q_Pd (molar). * CONSTANTS[1] is M_Pd_PdPs_star in component Q_Pd (molar). * CONSTANTS[2] is Y_Pd_InPd in component Q_Pd (molar). * CONSTANTS[3] is D_Pd in component Q_Pd (litre_per_day). * CONSTANTS[154] is k_PdEx in component Q_Pd (first_order_rate_constant). * CONSTANTS[4] is T_Pd in component Q_Pd (day). * CONSTANTS[5] is T_Pd_star in component Q_Pd (day). * STATES[1] is Q_Ma in component Q_Ma (gram). * STATES[2] is Q_Mc in component Q_Mc (gram). * CONSTANTS[6] is V_Ru in component model_parameters (litre). * CONSTANTS[7] is k_SoEx in component model_parameters (first_order_rate_constant). * ALGEBRAIC[3] is C_Ps in component Q_Ps (molar). * ALGEBRAIC[61] is U_McPs_PsAm in component Q_Ps (mole_per_day). * ALGEBRAIC[62] is U_Ps_PsMa in component Q_Ps (mole_per_day). * ALGEBRAIC[63] is U_Ps_PsMc in component Q_Ps (mole_per_day). * ALGEBRAIC[60] is U_MaPs_PsAm in component Q_Ps (mole_per_day). * CONSTANTS[8] is M_Ha_McMa in component Q_Ps (molar). * STATES[3] is Q_Ps in component Q_Ps (mole). * CONSTANTS[176] is P_Ps_InPs in component Q_Ps (mole_per_day). * CONSTANTS[179] is P_Ps_SaPs in component Q_Ps (mole_per_day). * ALGEBRAIC[58] is P_Ps_PdPs in component Q_Ps (mole_per_day). * ALGEBRAIC[59] is P_Ps_MaMd in component Q_Ps (mole_per_day). * ALGEBRAIC[167] is P_Ps_McPs in component Q_Ps (mole_per_day). * ALGEBRAIC[4] is U_Ps_PsEx in component Q_Ps (mole_per_day). * ALGEBRAIC[136] is U_Mc_McPs in component Q_Ps (mole_per_day). * ALGEBRAIC[5] is v_Ma_PsAm in component Q_Ps (mole_per_day). * ALGEBRAIC[6] is v_Mc_PsAm in component Q_Ps (mole_per_day). * ALGEBRAIC[7] is v_PsMa in component Q_Ps (mole_per_day). * ALGEBRAIC[8] is v_PsMc in component Q_Ps (mole_per_day). * CONSTANTS[9] is v_PsAm_star in component Q_Ps (mole_per_gram_day). * CONSTANTS[10] is v_PsMa_star in component Q_Ps (mole_per_gram_day). * CONSTANTS[11] is v_PsMc_star in component Q_Ps (mole_per_gram_day). * CONSTANTS[156] is k_PsEx in component Q_Ps (first_order_rate_constant). * CONSTANTS[12] is Y_Ps_InPs in component Q_Ps (molar). * CONSTANTS[13] is Y_Ps_SaPs in component Q_Ps (molar). * CONSTANTS[14] is Y_Ps_PdPs in component Q_Ps (dimensionless). * CONSTANTS[15] is Y_Ps_MaMd in component Q_Ps (mole_per_gram). * CONSTANTS[16] is Y_Ps_McPs in component Q_Ps (dimensionless). * CONSTANTS[17] is D_Ps in component Q_Ps (litre_per_day). * CONSTANTS[18] is M_Ps_PsAm in component Q_Ps (molar). * CONSTANTS[19] is M_Ps_PsMa in component Q_Ps (molar). * CONSTANTS[20] is M_Ps_PsMc in component Q_Ps (molar). * CONSTANTS[21] is M_Ha_PsMa in component Q_Ps (molar). * CONSTANTS[22] is M_Hc_PsMc in component Q_Ps (molar). * CONSTANTS[23] is J_Ha_PsAm in component Q_Ps (molar). * CONSTANTS[24] is J_Hc_PsAm in component Q_Ps (molar). * ALGEBRAIC[87] is U_Mc_McEg in component Q_Mc (gram_per_day). * ALGEBRAIC[35] is U_Ma_MaMd in component Q_Ma (gram_per_day). * ALGEBRAIC[80] is U_Hc_PsMc in component Q_Hc (mole_per_day). * ALGEBRAIC[169] is U_Ha_McMa in component Q_Ha (mole_per_day). * ALGEBRAIC[27] is C_Ha in component Q_Ha (molar). * ALGEBRAIC[31] is C_Hc in component Q_Hc (molar). * CONSTANTS[25] is k_FlEx in component model_parameters (first_order_rate_constant). * CONSTANTS[26] is J_Ha_McAm in component model_parameters (molar). * STATES[4] is Q_Pu in component Q_Pu (mole). * ALGEBRAIC[9] is C_Pu in component Q_Pu (molar). * CONSTANTS[157] is P_Pu_InPu in component Q_Pu (mole_per_day). * ALGEBRAIC[10] is U_Pu_PuEx in component Q_Pu (mole_per_day). * CONSTANTS[27] is Y_Pu_InPu in component Q_Pu (molar). * CONSTANTS[28] is D_Pu in component Q_Pu (litre_per_day). * CONSTANTS[158] is k_PuEx in component Q_Pu (first_order_rate_constant). * ALGEBRAIC[11] is C_Am in component Q_Am (molar). * ALGEBRAIC[68] is U_Am_AmMc in component Q_Am (mole_per_day). * ALGEBRAIC[67] is U_Am_AmMa in component Q_Am (mole_per_day). * ALGEBRAIC[137] is U_Mc_McAm in component Q_Am (mole_per_day). * STATES[5] is Q_Am in component Q_Am (mole). * CONSTANTS[180] is P_Am_InAm in component Q_Am (mole_per_day). * ALGEBRAIC[64] is P_Am_UeAm in component Q_Am (mole_per_day). * ALGEBRAIC[65] is P_MaAm_PsAm in component Q_Am (mole_per_day). * ALGEBRAIC[66] is P_McAm_PsAm in component Q_Am (mole_per_day). * ALGEBRAIC[168] is P_Am_McAm in component Q_Am (mole_per_day). * ALGEBRAIC[12] is U_Am_AmAb in component Q_Am (mole_per_day). * ALGEBRAIC[13] is U_Am_AmEx in component Q_Am (mole_per_day). * ALGEBRAIC[14] is v_UeAm in component Q_Am (mole_per_day). * CONSTANTS[181] is v_AmAb in component Q_Am (mole_per_day). * ALGEBRAIC[15] is v_AmMa in component Q_Am (mole_per_day). * ALGEBRAIC[16] is v_AmMc in component Q_Am (mole_per_day). * CONSTANTS[29] is v_UeAm_star in component Q_Am (mole_per_litre_day). * CONSTANTS[30] is v_AmAb_star in component Q_Am (mole_per_litre_day). * CONSTANTS[31] is v_AmMa_star in component Q_Am (mole_per_gram_day). * CONSTANTS[32] is v_AmMc_star in component Q_Am (mole_per_gram_day). * CONSTANTS[159] is k_AmEx in component Q_Am (first_order_rate_constant). * CONSTANTS[33] is Y_Am_InAm in component Q_Am (molar). * CONSTANTS[34] is Y_Am_UeAm in component Q_Am (dimensionless). * CONSTANTS[35] is Y_Am_PsAm in component Q_Am (dimensionless). * CONSTANTS[36] is Y_Am_McAm in component Q_Am (dimensionless). * CONSTANTS[37] is M_Am_AmMa in component Q_Am (molar). * CONSTANTS[38] is M_Am_AmMc in component Q_Am (molar). * CONSTANTS[39] is M_Am_AmAb in component Q_Am (molar). * CONSTANTS[40] is M_pH_AmAb in component Q_Am (dimensionless). * CONSTANTS[41] is phi_pH_AmAb in component Q_Am (dimensionless). * CONSTANTS[42] is M_Ha_AmMa in component Q_Am (molar). * CONSTANTS[43] is M_Hc_AmMc in component Q_Am (molar). * CONSTANTS[44] is J_Am_UeAm in component Q_Am (molar). * CONSTANTS[45] is D_Am in component Q_Am (litre_per_day). * CONSTANTS[46] is pH in component model_parameters (dimensionless). * CONSTANTS[47] is f_Po_Ma in component model_parameters (dimensionless). * CONSTANTS[48] is f_Ni_Fe in component model_parameters (dimensionless). * STATES[6] is Q_Li in component Q_Li (mole). * ALGEBRAIC[17] is C_Li in component Q_Li (molar). * CONSTANTS[182] is P_Li_InLi in component Q_Li (mole_per_day). * ALGEBRAIC[69] is P_Li_MaMd in component Q_Li (mole_per_day). * ALGEBRAIC[18] is U_Li_LiEx in component Q_Li (mole_per_day). * CONSTANTS[49] is Y_Li_InLi in component Q_Li (molar). * CONSTANTS[50] is Y_Li_MaMd in component Q_Li (mole_per_gram). * CONSTANTS[51] is D_Li in component Q_Li (litre_per_day). * CONSTANTS[160] is k_LiEx in component Q_Li (first_order_rate_constant). * ALGEBRAIC[70] is U_Sd_SdHa in component Q_Sd (mole_per_day). * STATES[7] is Q_Sd in component Q_Sd (mole). * ALGEBRAIC[19] is C_Sd in component Q_Sd (molar). * CONSTANTS[183] is P_Sd_InSd in component Q_Sd (mole_per_day). * ALGEBRAIC[20] is U_Sd_SdEx in component Q_Sd (mole_per_day). * CONSTANTS[52] is Y_Sd_InSd in component Q_Sd (molar). * CONSTANTS[53] is D_Sd in component Q_Sd (litre_per_day). * CONSTANTS[161] is k_SdEx in component Q_Sd (first_order_rate_constant). * ALGEBRAIC[21] is v_SdHa in component Q_Sd (mole_per_day). * CONSTANTS[54] is v_SdHa_star in component Q_Sd (mole_per_gram_day). * CONSTANTS[184] is M_Sd_SdHa in component Q_Sd (molar). * CONSTANTS[55] is M_Sd_SdHa_star in component Q_Sd (molar). * CONSTANTS[56] is T_Sd in component Q_Sd (day). * CONSTANTS[57] is T_Sd_star in component Q_Sd (day). * ALGEBRAIC[71] is U_Fd_FdHc in component Q_Fd (mole_per_day). * STATES[8] is Q_Fd in component Q_Fd (mole). * ALGEBRAIC[22] is C_Fd in component Q_Fd (molar). * CONSTANTS[185] is P_Fd_InFd in component Q_Fd (mole_per_day). * ALGEBRAIC[23] is U_Fd_FdEx in component Q_Fd (mole_per_day). * CONSTANTS[58] is Y_Fd_InFd in component Q_Fd (molar). * CONSTANTS[59] is D_Fd in component Q_Fd (litre_per_day). * CONSTANTS[162] is k_FdEx in component Q_Fd (first_order_rate_constant). * ALGEBRAIC[24] is v_FdHc in component Q_Fd (mole_per_day). * CONSTANTS[60] is v_FdHc_star in component Q_Fd (mole_per_gram_day). * CONSTANTS[186] is M_Fd_FdHc in component Q_Fd (molar). * CONSTANTS[61] is M_Fd_FdHc_star in component Q_Fd (molar). * CONSTANTS[62] is T_Fd in component Q_Fd (day). * CONSTANTS[63] is T_Fd_star in component Q_Fd (day). * CONSTANTS[64] is phi_pH_FdHc in component Q_Fd (dimensionless). * CONSTANTS[65] is M_pH_FdHc in component Q_Fd (dimensionless). * CONSTANTS[66] is pm in component model_parameters (dimensionless). * CONSTANTS[67] is f in component model_parameters (dimensionless). * STATES[9] is Q_Fu in component Q_Fu (mole). * ALGEBRAIC[25] is C_Fu in component Q_Fu (molar). * CONSTANTS[163] is P_Fu_InFu in component Q_Fu (mole_per_day). * ALGEBRAIC[26] is U_Fu_FuEx in component Q_Fu (mole_per_day). * CONSTANTS[68] is Y_Fu_InFu in component Q_Fu (molar). * CONSTANTS[69] is D_Fu in component Q_Fu (litre_per_day). * CONSTANTS[164] is k_FuEx in component Q_Fu (first_order_rate_constant). * ALGEBRAIC[76] is U_Ha_HaAs in component Q_Ha (mole_per_day). * ALGEBRAIC[74] is U_Ha_AmMa in component Q_Ha (mole_per_day). * ALGEBRAIC[75] is U_Ha_PsMa in component Q_Ha (mole_per_day). * ALGEBRAIC[77] is U_Ha_HaVa in component Q_Ha (mole_per_day). * CONSTANTS[187] is P_Ha_WrHa in component Q_Ha (mole_per_day). * CONSTANTS[189] is P_Ha_LaHa in component Q_Ha (mole_per_day). * CONSTANTS[188] is P_Ha_SrHa in component Q_Ha (mole_per_day). * ALGEBRAIC[72] is P_Ha_SdHa in component Q_Ha (mole_per_day). * ALGEBRAIC[73] is P_Ha_MaMd in component Q_Ha (mole_per_day). * STATES[10] is Q_Ha in component Q_Ha (mole). * ALGEBRAIC[28] is U_Ha_HaEx in component Q_Ha (mole_per_day). * CONSTANTS[70] is R_Ha_AmMa in component Q_Ha (dimensionless). * CONSTANTS[71] is R_Ha_PsMa in component Q_Ha (dimensionless). * CONSTANTS[72] is R_Ha_McMa in component Q_Ha (mole_per_gram). * CONSTANTS[73] is M_Ha_HaAs in component Q_Ha (molar). * CONSTANTS[74] is M_Ha_HaVa in component Q_Ha (molar). * ALGEBRAIC[29] is v_HaAs in component Q_Ha (mole_per_day). * ALGEBRAIC[30] is v_HaVa in component Q_Ha (mole_per_day). * CONSTANTS[75] is v_HaAs_star in component Q_Ha (mole_per_gram_day). * CONSTANTS[76] is v_HaVa_star in component Q_Ha (mole_per_gram_day). * CONSTANTS[165] is k_HaEx in component Q_Ha (first_order_rate_constant). * CONSTANTS[77] is Y_Ha_WrHa in component Q_Ha (molar). * CONSTANTS[78] is Y_Ha_SrHa in component Q_Ha (molar). * CONSTANTS[79] is Y_Ha_LaHa in component Q_Ha (molar). * CONSTANTS[80] is Y_Ha_SdHa in component Q_Ha (dimensionless). * CONSTANTS[81] is Y_Ha_MaMd in component Q_Ha (mole_per_gram). * CONSTANTS[82] is D_Wr in component Q_Ha (litre_per_day). * CONSTANTS[83] is D_Sr in component Q_Ha (litre_per_day). * CONSTANTS[84] is D_La in component Q_Ha (litre_per_day). * CONSTANTS[85] is J_Am_HaVa in component Q_Ha (molar). * CONSTANTS[86] is J_Ps_HaVa in component Q_Ha (molar). * ALGEBRAIC[39] is U_As_MaMd in component Q_As (gram_per_day). * ALGEBRAIC[138] is U_Mc_McMa in component Q_Ma (gram_per_day). * ALGEBRAIC[79] is U_Hc_AmMc in component Q_Hc (mole_per_day). * ALGEBRAIC[81] is U_Hc_HcVa in component Q_Hc (mole_per_day). * STATES[11] is Q_Hc in component Q_Hc (mole). * ALGEBRAIC[78] is P_Hc_FdHc in component Q_Hc (mole_per_day). * ALGEBRAIC[32] is U_Hc_HcEx in component Q_Hc (mole_per_day). * CONSTANTS[87] is R_Hc_AmMc in component Q_Hc (dimensionless). * CONSTANTS[88] is R_Hc_PsMc in component Q_Hc (dimensionless). * ALGEBRAIC[33] is v_HcVa in component Q_Hc (mole_per_day). * CONSTANTS[89] is v_HcVa_star in component Q_Hc (mole_per_gram_day). * CONSTANTS[166] is k_HcEx in component Q_Hc (first_order_rate_constant). * CONSTANTS[90] is Y_Hc_FdHc in component Q_Hc (dimensionless). * CONSTANTS[91] is M_Hc_HcVa in component Q_Hc (molar). * CONSTANTS[92] is J_Am_HcVa in component Q_Hc (molar). * CONSTANTS[93] is J_Ps_HcVa in component Q_Hc (molar). * ALGEBRAIC[34] is C_Ma in component Q_Ma (gram_per_litre). * ALGEBRAIC[82] is P_Ma_AmMa in component Q_Ma (gram_per_day). * ALGEBRAIC[83] is P_Ma_PsMa in component Q_Ma (gram_per_day). * ALGEBRAIC[170] is P_Ma_McMa in component Q_Ma (gram_per_day). * ALGEBRAIC[37] is U_Ma_PoEx in component Q_Ma (gram_per_day). * ALGEBRAIC[36] is U_Ma_MaEx in component Q_Ma (gram_per_day). * CONSTANTS[177] is k_MaMd in component Q_Ma (first_order_rate_constant). * CONSTANTS[94] is Y_Ma_AmMa in component Q_Ma (gram_per_mole). * CONSTANTS[95] is Y_Ma_PsMa in component Q_Ma (gram_per_mole). * CONSTANTS[96] is Y_Ma_McMa in component Q_Ma (dimensionless). * CONSTANTS[167] is k_MaEx in component model_parameters (first_order_rate_constant). * CONSTANTS[175] is k_PoEx in component model_parameters (first_order_rate_constant). * STATES[12] is Q_As in component Q_As (gram). * ALGEBRAIC[38] is C_As in component Q_As (gram_per_litre). * ALGEBRAIC[84] is P_As_HaAs in component Q_As (gram_per_day). * ALGEBRAIC[41] is U_As_PoEx in component Q_As (gram_per_day). * ALGEBRAIC[40] is U_As_AsEx in component Q_As (gram_per_day). * CONSTANTS[168] is k_AsEx in component Q_As (first_order_rate_constant). * CONSTANTS[178] is k_AsMd in component Q_As (first_order_rate_constant). * CONSTANTS[97] is Y_As_HaAs in component Q_As (gram_per_mole). * ALGEBRAIC[42] is C_Mc in component Q_Mc (gram_per_litre). * ALGEBRAIC[85] is P_Mc_AmMc in component Q_Mc (gram_per_day). * ALGEBRAIC[86] is P_Mc_PsMc in component Q_Mc (gram_per_day). * ALGEBRAIC[43] is U_Mc_McEx in component Q_Mc (gram_per_day). * CONSTANTS[169] is k_McEx in component Q_Mc (first_order_rate_constant). * CONSTANTS[98] is Y_Mc_AmMc in component Q_Mc (gram_per_mole). * CONSTANTS[99] is Y_Mc_PsMc in component Q_Mc (gram_per_mole). * ALGEBRAIC[44] is v_McEg in component Q_Mc (gram_per_day). * CONSTANTS[100] is v_McEg_star in component Q_Mc (first_order_rate_constant). * CONSTANTS[101] is M_Mc_McEg in component Q_Mc (gram_per_litre). * STATES[13] is Q_Ac in component Q_Ac (mole). * ALGEBRAIC[45] is C_Ac in component Q_Ac (molar). * CONSTANTS[190] is P_Ac_InAc in component Q_Ac (mole_per_day). * ALGEBRAIC[139] is P_Ac_AmMa in component Q_Ac (mole_per_day). * ALGEBRAIC[140] is P_Ac_AmMc in component Q_Ac (mole_per_day). * ALGEBRAIC[141] is P_Ac_PsMa in component Q_Ac (mole_per_day). * ALGEBRAIC[142] is P_Ac_PsMc in component Q_Ac (mole_per_day). * ALGEBRAIC[143] is P_Ac_HaAs in component Q_Ac (mole_per_day). * ALGEBRAIC[144] is P_Ac_HaVa in component Q_Ac (mole_per_day). * ALGEBRAIC[88] is P_Ac_HcVa in component Q_Ac (mole_per_day). * ALGEBRAIC[89] is P_MaAc_PsAm in component Q_Ac (mole_per_day). * ALGEBRAIC[90] is P_McAc_PsAm in component Q_Ac (mole_per_day). * ALGEBRAIC[175] is P_Ac_McMa in component Q_Ac (mole_per_day). * ALGEBRAIC[145] is P_Ac_McAm in component Q_Ac (mole_per_day). * ALGEBRAIC[46] is U_Ac_AcAb in component Q_Ac (mole_per_day). * ALGEBRAIC[47] is U_Ac_AcEx in component Q_Ac (mole_per_day). * ALGEBRAIC[91] is U_Hf_AmMa in component Q_Ac (mole_per_day). * ALGEBRAIC[92] is U_Hf_AmMc in component Q_Ac (mole_per_day). * ALGEBRAIC[93] is U_Hf_PsMa in component Q_Ac (mole_per_day). * ALGEBRAIC[94] is U_Hf_PsMc in component Q_Ac (mole_per_day). * ALGEBRAIC[95] is U_Hf_HaAs in component Q_Ac (mole_per_day). * ALGEBRAIC[171] is U_Hf_McMa in component Q_Ac (mole_per_day). * CONSTANTS[102] is M_Ac_AcAb in component Q_Ac (molar). * CONSTANTS[199] is v_AcAb in component Q_Ac (mole_per_day). * CONSTANTS[103] is v_AcAb_star in component Q_Ac (mole_per_litre_day). * CONSTANTS[191] is Y_Ac_LaAc in component Q_Ac (dimensionless). * ALGEBRAIC[96] is Y_Ac_HaVa in component Q_Ac (dimensionless). * CONSTANTS[193] is Y_Ac_PsVa in component Q_Ac (dimensionless). * ALGEBRAIC[97] is Y_Ac_AmMa in component Q_Ac (dimensionless). * CONSTANTS[194] is Y_Ac_AmMc in component Q_Ac (dimensionless). * ALGEBRAIC[98] is Y_Ac_PsMa in component Q_Ac (dimensionless). * CONSTANTS[195] is Y_Ac_PsMc in component Q_Ac (dimensionless). * ALGEBRAIC[99] is Y_Ac_HaAs in component Q_Ac (dimensionless). * CONSTANTS[196] is Y_Ac_McMa in component Q_Ac (dimensionless). * CONSTANTS[197] is Y_Ac_PsAm in component Q_Ac (dimensionless). * CONSTANTS[198] is Y_Ac_McAm in component Q_Ac (dimensionless). * CONSTANTS[104] is Y_Ac_WrAc in component Q_Ac (dimensionless). * CONSTANTS[105] is Y_Ac_InAc in component Q_Ac (molar). * CONSTANTS[106] is Y_Ac_StAc in component Q_Ac (dimensionless). * CONSTANTS[107] is Y_Ac_CeAc in component Q_Ac (dimensionless). * CONSTANTS[108] is Y_Ac_HeAc in component Q_Ac (dimensionless). * CONSTANTS[109] is Y_Ac_PsAc in component Q_Ac (dimensionless). * CONSTANTS[192] is Y_Ac_HcVa in component Q_Ac (dimensionless). * CONSTANTS[170] is k_AcEx in component Q_Ac (first_order_rate_constant). * CONSTANTS[110] is D_Ac in component Q_Ac (litre_per_day). * CONSTANTS[111] is J_pH_AcAb in component Q_Ac (dimensionless). * CONSTANTS[112] is phi_pH_AcAb in component Q_Ac (dimensionless). * CONSTANTS[113] is f_Hf_AmMa in component model_parameters (dimensionless). * CONSTANTS[114] is f_Hf_AmMc in component model_parameters (dimensionless). * CONSTANTS[115] is f_Hf_PsMa in component model_parameters (dimensionless). * CONSTANTS[116] is f_Hf_PsMc in component model_parameters (dimensionless). * CONSTANTS[117] is f_Hf_HaAs in component model_parameters (dimensionless). * CONSTANTS[118] is f_Hf_McMa in component model_parameters (dimensionless). * CONSTANTS[119] is f_Lc_Le in component model_parameters (dimensionless). * CONSTANTS[120] is f_Ce_Fd in component model_parameters (dimensionless). * STATES[14] is Q_Pr in component Q_Pr (mole). * ALGEBRAIC[48] is C_Pr in component Q_Pr (molar). * CONSTANTS[200] is P_Pr_InPr in component Q_Pr (mole_per_day). * ALGEBRAIC[146] is P_Pr_AmMa in component Q_Pr (mole_per_day). * ALGEBRAIC[147] is P_Pr_AmMc in component Q_Pr (mole_per_day). * ALGEBRAIC[148] is P_Pr_PsMa in component Q_Pr (mole_per_day). * ALGEBRAIC[149] is P_Pr_PsMc in component Q_Pr (mole_per_day). * ALGEBRAIC[150] is P_Pr_HaAs in component Q_Pr (mole_per_day). * ALGEBRAIC[151] is P_Pr_HaVa in component Q_Pr (mole_per_day). * ALGEBRAIC[100] is P_Pr_HcVa in component Q_Pr (mole_per_day). * ALGEBRAIC[101] is P_MaPr_PsAm in component Q_Pr (mole_per_day). * ALGEBRAIC[102] is P_McPr_PsAm in component Q_Pr (mole_per_day). * ALGEBRAIC[176] is P_Pr_McMa in component Q_Pr (mole_per_day). * ALGEBRAIC[152] is P_Pr_McAm in component Q_Pr (mole_per_day). * ALGEBRAIC[49] is U_Pr_PrAb in component Q_Pr (mole_per_day). * ALGEBRAIC[50] is U_Pr_PrEx in component Q_Pr (mole_per_day). * ALGEBRAIC[103] is U_Hf_AmMa in component Q_Pr (mole_per_day). * ALGEBRAIC[104] is U_Hf_AmMc in component Q_Pr (mole_per_day). * ALGEBRAIC[105] is U_Hf_PsMa in component Q_Pr (mole_per_day). * ALGEBRAIC[106] is U_Hf_PsMc in component Q_Pr (mole_per_day). * ALGEBRAIC[107] is U_Hf_HaAs in component Q_Pr (mole_per_day). * ALGEBRAIC[172] is U_Hf_McMa in component Q_Pr (mole_per_day). * CONSTANTS[121] is M_Pr_PrAb in component Q_Pr (molar). * CONSTANTS[209] is v_PrAb in component Q_Pr (mole_per_day). * CONSTANTS[122] is v_PrAb_star in component Q_Pr (mole_per_litre_day). * CONSTANTS[201] is Y_Pr_LaPr in component Q_Pr (dimensionless). * ALGEBRAIC[108] is Y_Pr_HaVa in component Q_Pr (dimensionless). * CONSTANTS[203] is Y_Pr_PsVa in component Q_Pr (dimensionless). * ALGEBRAIC[109] is Y_Pr_AmMa in component Q_Pr (dimensionless). * CONSTANTS[204] is Y_Pr_AmMc in component Q_Pr (dimensionless). * ALGEBRAIC[110] is Y_Pr_PsMa in component Q_Pr (dimensionless). * CONSTANTS[205] is Y_Pr_PsMc in component Q_Pr (dimensionless). * ALGEBRAIC[111] is Y_Pr_HaAs in component Q_Pr (dimensionless). * CONSTANTS[206] is Y_Pr_McMa in component Q_Pr (dimensionless). * CONSTANTS[207] is Y_Pr_PsAm in component Q_Pr (dimensionless). * CONSTANTS[208] is Y_Pr_McAm in component Q_Pr (dimensionless). * CONSTANTS[123] is Y_Pr_WrPr in component Q_Pr (dimensionless). * CONSTANTS[124] is Y_Pr_InPr in component Q_Pr (molar). * CONSTANTS[202] is Y_Pr_HcVa in component Q_Pr (dimensionless). * CONSTANTS[125] is Y_Pr_StPr in component Q_Pr (dimensionless). * CONSTANTS[126] is Y_Pr_CePr in component Q_Pr (dimensionless). * CONSTANTS[127] is Y_Pr_HePr in component Q_Pr (dimensionless). * CONSTANTS[128] is Y_Pr_PsPr in component Q_Pr (dimensionless). * CONSTANTS[171] is k_PrEx in component Q_Pr (first_order_rate_constant). * CONSTANTS[129] is D_Pr in component Q_Pr (litre_per_day). * CONSTANTS[130] is phi_pH_PrAb in component Q_Pr (dimensionless). * CONSTANTS[131] is J_pH_PrAb in component Q_Pr (dimensionless). * STATES[15] is Q_Bu in component Q_Bu (mole). * ALGEBRAIC[51] is C_Bu in component Q_Bu (molar). * CONSTANTS[210] is P_Bu_InBu in component Q_Bu (mole_per_day). * ALGEBRAIC[153] is P_Bu_AmMa in component Q_Bu (mole_per_day). * ALGEBRAIC[154] is P_Bu_AmMc in component Q_Bu (mole_per_day). * ALGEBRAIC[155] is P_Bu_PsMa in component Q_Bu (mole_per_day). * ALGEBRAIC[156] is P_Bu_PsMc in component Q_Bu (mole_per_day). * ALGEBRAIC[157] is P_Bu_HaAs in component Q_Bu (mole_per_day). * ALGEBRAIC[158] is P_Bu_HaVa in component Q_Bu (mole_per_day). * ALGEBRAIC[112] is P_Bu_HcVa in component Q_Bu (mole_per_day). * ALGEBRAIC[113] is P_MaBu_PsAm in component Q_Bu (mole_per_day). * ALGEBRAIC[114] is P_McBu_PsAm in component Q_Bu (mole_per_day). * ALGEBRAIC[177] is P_Bu_McMa in component Q_Bu (mole_per_day). * ALGEBRAIC[159] is P_Bu_McAm in component Q_Bu (mole_per_day). * ALGEBRAIC[52] is U_Bu_BuAb in component Q_Bu (mole_per_day). * ALGEBRAIC[53] is U_Bu_BuEx in component Q_Bu (mole_per_day). * ALGEBRAIC[115] is U_Hf_AmMa in component Q_Bu (mole_per_day). * ALGEBRAIC[116] is U_Hf_AmMc in component Q_Bu (mole_per_day). * ALGEBRAIC[117] is U_Hf_PsMa in component Q_Bu (mole_per_day). * ALGEBRAIC[118] is U_Hf_PsMc in component Q_Bu (mole_per_day). * ALGEBRAIC[119] is U_Hf_HaAs in component Q_Bu (mole_per_day). * ALGEBRAIC[173] is U_Hf_McMa in component Q_Bu (mole_per_day). * CONSTANTS[132] is M_Bu_BuAb in component Q_Bu (molar). * CONSTANTS[219] is v_BuAb in component Q_Bu (mole_per_day). * CONSTANTS[133] is v_BuAb_star in component Q_Bu (mole_per_litre_day). * CONSTANTS[211] is Y_Bu_LaBu in component Q_Bu (dimensionless). * ALGEBRAIC[120] is Y_Bu_HaVa in component Q_Bu (dimensionless). * CONSTANTS[213] is Y_Bu_PsVa in component Q_Bu (dimensionless). * ALGEBRAIC[121] is Y_Bu_AmMa in component Q_Bu (dimensionless). * CONSTANTS[214] is Y_Bu_AmMc in component Q_Bu (dimensionless). * ALGEBRAIC[122] is Y_Bu_PsMa in component Q_Bu (dimensionless). * CONSTANTS[215] is Y_Bu_PsMc in component Q_Bu (dimensionless). * ALGEBRAIC[123] is Y_Bu_HaAs in component Q_Bu (dimensionless). * CONSTANTS[216] is Y_Bu_McMa in component Q_Bu (dimensionless). * CONSTANTS[217] is Y_Bu_PsAm in component Q_Bu (dimensionless). * CONSTANTS[218] is Y_Bu_McAm in component Q_Bu (dimensionless). * CONSTANTS[134] is Y_Bu_WrBu in component Q_Bu (dimensionless). * CONSTANTS[135] is Y_Bu_InBu in component Q_Bu (molar). * CONSTANTS[212] is Y_Bu_HcVa in component Q_Bu (dimensionless). * CONSTANTS[136] is Y_Bu_StBu in component Q_Bu (dimensionless). * CONSTANTS[137] is Y_Bu_CeBu in component Q_Bu (dimensionless). * CONSTANTS[138] is Y_Bu_HeBu in component Q_Bu (dimensionless). * CONSTANTS[139] is Y_Bu_PsBu in component Q_Bu (dimensionless). * CONSTANTS[172] is k_BuEx in component Q_Bu (first_order_rate_constant). * CONSTANTS[140] is D_Bu in component Q_Bu (litre_per_day). * CONSTANTS[141] is phi_pH_BuAb in component Q_Bu (dimensionless). * CONSTANTS[142] is J_pH_BuAb in component Q_Bu (dimensionless). * STATES[16] is Q_Vl in component Q_Vl (mole). * ALGEBRAIC[54] is C_Vl in component Q_Vl (molar). * CONSTANTS[220] is P_Vl_InVl in component Q_Vl (mole_per_day). * ALGEBRAIC[160] is P_Vl_AmMa in component Q_Vl (mole_per_day). * ALGEBRAIC[161] is P_Vl_AmMc in component Q_Vl (mole_per_day). * ALGEBRAIC[162] is P_Vl_PsMa in component Q_Vl (mole_per_day). * ALGEBRAIC[163] is P_Vl_PsMc in component Q_Vl (mole_per_day). * ALGEBRAIC[164] is P_Vl_HaAs in component Q_Vl (mole_per_day). * ALGEBRAIC[165] is P_Vl_HaVa in component Q_Vl (mole_per_day). * ALGEBRAIC[124] is P_Vl_HcVa in component Q_Vl (mole_per_day). * ALGEBRAIC[125] is P_MaVl_PsAm in component Q_Vl (mole_per_day). * ALGEBRAIC[126] is P_McVl_PsAm in component Q_Vl (mole_per_day). * ALGEBRAIC[178] is P_Vl_McMa in component Q_Vl (mole_per_day). * ALGEBRAIC[166] is P_Vl_McAm in component Q_Vl (mole_per_day). * ALGEBRAIC[55] is U_Vl_VlAb in component Q_Vl (mole_per_day). * ALGEBRAIC[56] is U_Vl_VlEx in component Q_Vl (mole_per_day). * ALGEBRAIC[127] is U_Hf_AmMa in component Q_Vl (mole_per_day). * ALGEBRAIC[128] is U_Hf_AmMc in component Q_Vl (mole_per_day). * ALGEBRAIC[129] is U_Hf_PsMa in component Q_Vl (mole_per_day). * ALGEBRAIC[130] is U_Hf_PsMc in component Q_Vl (mole_per_day). * ALGEBRAIC[131] is U_Hf_HaAs in component Q_Vl (mole_per_day). * ALGEBRAIC[174] is U_Hf_McMa in component Q_Vl (mole_per_day). * CONSTANTS[143] is M_Vl_VlAb in component Q_Vl (molar). * CONSTANTS[229] is v_VlAb in component Q_Vl (mole_per_day). * CONSTANTS[144] is v_VlAb_star in component Q_Vl (mole_per_litre_day). * CONSTANTS[221] is Y_Vl_LaVl in component Q_Vl (dimensionless). * ALGEBRAIC[132] is Y_Vl_HaVa in component Q_Vl (dimensionless). * CONSTANTS[223] is Y_Vl_PsVa in component Q_Vl (dimensionless). * ALGEBRAIC[133] is Y_Vl_AmMa in component Q_Vl (dimensionless). * CONSTANTS[224] is Y_Vl_AmMc in component Q_Vl (dimensionless). * ALGEBRAIC[134] is Y_Vl_PsMa in component Q_Vl (dimensionless). * CONSTANTS[225] is Y_Vl_PsMc in component Q_Vl (dimensionless). * ALGEBRAIC[135] is Y_Vl_HaAs in component Q_Vl (dimensionless). * CONSTANTS[226] is Y_Vl_McMa in component Q_Vl (dimensionless). * CONSTANTS[227] is Y_Vl_PsAm in component Q_Vl (dimensionless). * CONSTANTS[228] is Y_Vl_McAm in component Q_Vl (dimensionless). * CONSTANTS[145] is Y_Vl_WrVl in component Q_Vl (dimensionless). * CONSTANTS[146] is Y_Vl_InVl in component Q_Vl (molar). * CONSTANTS[222] is Y_Vl_HcVa in component Q_Vl (dimensionless). * CONSTANTS[147] is Y_Vl_StVl in component Q_Vl (dimensionless). * CONSTANTS[148] is Y_Vl_CeVl in component Q_Vl (dimensionless). * CONSTANTS[149] is Y_Vl_HeVl in component Q_Vl (dimensionless). * CONSTANTS[150] is Y_Vl_PsVl in component Q_Vl (dimensionless). * CONSTANTS[173] is k_VlEx in component Q_Vl (first_order_rate_constant). * CONSTANTS[151] is D_Vl in component Q_Vl (litre_per_day). * CONSTANTS[152] is phi_pH_VlAb in component Q_Vl (dimensionless). * CONSTANTS[153] is J_pH_VlAb in component Q_Vl (dimensionless). * RATES[0] is d/dt Q_Pd in component Q_Pd (mole). * RATES[3] is d/dt Q_Ps in component Q_Ps (mole). * RATES[4] is d/dt Q_Pu in component Q_Pu (mole). * RATES[5] is d/dt Q_Am in component Q_Am (mole). * RATES[6] is d/dt Q_Li in component Q_Li (mole). * RATES[7] is d/dt Q_Sd in component Q_Sd (mole). * RATES[8] is d/dt Q_Fd in component Q_Fd (mole). * RATES[9] is d/dt Q_Fu in component Q_Fu (mole). * RATES[10] is d/dt Q_Ha in component Q_Ha (mole). * RATES[11] is d/dt Q_Hc in component Q_Hc (mole). * RATES[1] is d/dt Q_Ma in component Q_Ma (gram). * RATES[12] is d/dt Q_As in component Q_As (gram). * RATES[2] is d/dt Q_Mc in component Q_Mc (gram). * RATES[13] is d/dt Q_Ac in component Q_Ac (mole). * RATES[14] is d/dt Q_Pr in component Q_Pr (mole). * RATES[15] is d/dt Q_Bu in component Q_Bu (mole). * RATES[16] is d/dt Q_Vl in component Q_Vl (mole). * There are a total of 0 condition variables. */ void initConsts(double* CONSTANTS, double* RATES, double *STATES) { STATES[0] = 1.000; CONSTANTS[0] = 0.0576; CONSTANTS[1] = 0.264; CONSTANTS[2] = 0.0091; CONSTANTS[3] = 1.000; CONSTANTS[4] = 1.000; CONSTANTS[5] = 0.66; STATES[1] = 1.000; STATES[2] = 1.000; CONSTANTS[6] = 1.000; CONSTANTS[7] = 1.000; CONSTANTS[8] = 0.0248; STATES[3] = 1.000; CONSTANTS[9] = 0.0144; CONSTANTS[10] = 0.0576; CONSTANTS[11] = 0.0576; CONSTANTS[12] = 0.0091; CONSTANTS[13] = 0.0010; CONSTANTS[14] = 1.0; CONSTANTS[15] = 0.0067; CONSTANTS[16] = 0.0067; CONSTANTS[17] = 1.000; CONSTANTS[18] = 0.0289; CONSTANTS[19] = 0.0224; CONSTANTS[20] = 0.0224; CONSTANTS[21] = 0.0248; CONSTANTS[22] = 0.0248; CONSTANTS[23] = 0.0165; CONSTANTS[24] = 0.0165; CONSTANTS[25] = 1.000; CONSTANTS[26] = 0.0165; STATES[4] = 1.000; CONSTANTS[27] = 0.0091; CONSTANTS[28] = 1.000; STATES[5] = 1.000; CONSTANTS[29] = 0.00165; CONSTANTS[30] = 1.10; CONSTANTS[31] = 0.0528; CONSTANTS[32] = 0.0528; CONSTANTS[33] = 0.0588; CONSTANTS[34] = 2.0; CONSTANTS[35] = 1.257; CONSTANTS[36] = 0.0084; CONSTANTS[37] = 0.00135; CONSTANTS[38] = 0.00135; CONSTANTS[39] = 0.0132; CONSTANTS[40] = 7.5; CONSTANTS[41] = 7.85; CONSTANTS[42] = 0.0159; CONSTANTS[43] = 0.0159; CONSTANTS[44] = 0.00621; CONSTANTS[45] = 1.000; CONSTANTS[46] = 7.000; CONSTANTS[47] = 0.35; CONSTANTS[48] = 1.000; STATES[6] = 1.000; CONSTANTS[49] = 0.0015; CONSTANTS[50] = 0.00021; CONSTANTS[51] = 1.000; STATES[7] = 1.000; CONSTANTS[52] = 0.0062; CONSTANTS[53] = 1.000; CONSTANTS[54] = 0.2179; CONSTANTS[55] = 0.416; CONSTANTS[56] = 1.000; CONSTANTS[57] = 1.04; STATES[8] = 1.000; CONSTANTS[58] = 0.0062; CONSTANTS[59] = 1.000; CONSTANTS[60] = 0.1646; CONSTANTS[61] = 0.332; CONSTANTS[62] = 1.000; CONSTANTS[63] = 0.83; CONSTANTS[64] = 22.9; CONSTANTS[65] = 5.97; CONSTANTS[66] = 1.000; CONSTANTS[67] = 1.000; STATES[9] = 1.000; CONSTANTS[68] = 0.0062; CONSTANTS[69] = 1.000; STATES[10] = 1.000; CONSTANTS[70] = 1.793; CONSTANTS[71] = 1.291; CONSTANTS[72] = 0.0086; CONSTANTS[73] = 0.0268; CONSTANTS[74] = 0.055; CONSTANTS[75] = 0.053; CONSTANTS[76] = 0.1646; CONSTANTS[77] = 0.0062; CONSTANTS[78] = 0.0062; CONSTANTS[79] = 0.0025; CONSTANTS[80] = 1.0; CONSTANTS[81] = 0.0062; CONSTANTS[82] = 1.000; CONSTANTS[83] = 1.000; CONSTANTS[84] = 1.000; CONSTANTS[85] = 0.00861; CONSTANTS[86] = 0.01465; STATES[11] = 1.000; CONSTANTS[87] = 1.793; CONSTANTS[88] = 1.291; CONSTANTS[89] = 0.1646; CONSTANTS[90] = 1.0; CONSTANTS[91] = 0.055; CONSTANTS[92] = 0.00861; CONSTANTS[93] = 0.01465; CONSTANTS[94] = 118.91; CONSTANTS[95] = 149.48; CONSTANTS[96] = 149.48; STATES[12] = 1.000; CONSTANTS[97] = 112.5; CONSTANTS[98] = 118.91; CONSTANTS[99] = 149.48; CONSTANTS[100] = 15.439; CONSTANTS[101] = 34.694; STATES[13] = 1.000; CONSTANTS[102] = 0.338; CONSTANTS[103] = 7.86; CONSTANTS[104] = 1.000; CONSTANTS[105] = 0.0167; CONSTANTS[106] = 1.000; CONSTANTS[107] = 1.000; CONSTANTS[108] = 1.000; CONSTANTS[109] = 1.000; CONSTANTS[110] = 1.000; CONSTANTS[111] = 6.45; CONSTANTS[112] = 6.48; CONSTANTS[113] = 0.526; CONSTANTS[114] = 0.526; CONSTANTS[115] = 0.711; CONSTANTS[116] = 0.711; CONSTANTS[117] = 0.306; CONSTANTS[118] = 0.711; CONSTANTS[119] = 2.250; CONSTANTS[120] = 1.000; STATES[14] = 1.000; CONSTANTS[121] = 0.338; CONSTANTS[122] = 7.86; CONSTANTS[123] = 1.000; CONSTANTS[124] = 0.0135; CONSTANTS[125] = 1.000; CONSTANTS[126] = 1.000; CONSTANTS[127] = 1.000; CONSTANTS[128] = 1.000; CONSTANTS[129] = 1.000; CONSTANTS[130] = 6.48; CONSTANTS[131] = 6.45; STATES[15] = 1.000; CONSTANTS[132] = 0.338; CONSTANTS[133] = 7.86; CONSTANTS[134] = 1.000; CONSTANTS[135] = 0.0114; CONSTANTS[136] = 1.000; CONSTANTS[137] = 1.000; CONSTANTS[138] = 1.000; CONSTANTS[139] = 1.000; CONSTANTS[140] = 1.000; CONSTANTS[141] = 6.48; CONSTANTS[142] = 6.45; STATES[16] = 1.000; CONSTANTS[143] = 0.338; CONSTANTS[144] = 7.86; CONSTANTS[145] = 1.000; CONSTANTS[146] = 0.0098; CONSTANTS[147] = 1.000; CONSTANTS[148] = 1.000; CONSTANTS[149] = 1.000; CONSTANTS[150] = 1.000; CONSTANTS[151] = 1.000; CONSTANTS[152] = 6.48; CONSTANTS[153] = 6.45; CONSTANTS[154] = CONSTANTS[7]; CONSTANTS[155] = CONSTANTS[2]*CONSTANTS[3]; CONSTANTS[156] = CONSTANTS[25]; CONSTANTS[157] = CONSTANTS[27]*CONSTANTS[28]; CONSTANTS[158] = CONSTANTS[7]; CONSTANTS[159] = CONSTANTS[25]; CONSTANTS[160] = CONSTANTS[25]; CONSTANTS[161] = CONSTANTS[7]; CONSTANTS[162] = CONSTANTS[7]; CONSTANTS[163] = CONSTANTS[68]*CONSTANTS[69]; CONSTANTS[164] = CONSTANTS[7]; CONSTANTS[165] = CONSTANTS[25]; CONSTANTS[166] = CONSTANTS[25]; CONSTANTS[167] = CONSTANTS[25]; CONSTANTS[168] = CONSTANTS[25]; CONSTANTS[169] = CONSTANTS[7]; CONSTANTS[170] = CONSTANTS[25]; CONSTANTS[171] = CONSTANTS[25]; CONSTANTS[172] = CONSTANTS[25]; CONSTANTS[173] = CONSTANTS[25]; CONSTANTS[174] = CONSTANTS[1]*(CONSTANTS[4]/CONSTANTS[5]); CONSTANTS[175] = CONSTANTS[7]/2.00000; CONSTANTS[176] = CONSTANTS[12]*CONSTANTS[17]; CONSTANTS[177] = CONSTANTS[167] - CONSTANTS[175]; CONSTANTS[178] = CONSTANTS[167] - CONSTANTS[175]; CONSTANTS[179] = CONSTANTS[13]*CONSTANTS[17]; CONSTANTS[180] = CONSTANTS[33]*CONSTANTS[45]; CONSTANTS[181] = CONSTANTS[30]*pow(1.00000, 0.250000)*pow(CONSTANTS[6], 0.750000)*(1.00000+pow(CONSTANTS[40]/CONSTANTS[46], CONSTANTS[41])); CONSTANTS[182] = CONSTANTS[49]*CONSTANTS[51]; CONSTANTS[183] = CONSTANTS[52]*CONSTANTS[53]; CONSTANTS[184] = CONSTANTS[55]*(CONSTANTS[56]/CONSTANTS[57]); CONSTANTS[185] = CONSTANTS[58]*CONSTANTS[59]; CONSTANTS[186] = CONSTANTS[61]*(CONSTANTS[62]/CONSTANTS[63]); CONSTANTS[187] = CONSTANTS[77]*CONSTANTS[82]; CONSTANTS[188] = CONSTANTS[78]*CONSTANTS[83]; CONSTANTS[189] = CONSTANTS[79]*CONSTANTS[84]; CONSTANTS[190] = CONSTANTS[105]*CONSTANTS[110]; CONSTANTS[191] = CONSTANTS[119]*CONSTANTS[104]; CONSTANTS[192] = CONSTANTS[120]*CONSTANTS[107]+ (1.00000 - CONSTANTS[120])*CONSTANTS[108]; CONSTANTS[193] = CONSTANTS[109]; CONSTANTS[194] = CONSTANTS[192]; CONSTANTS[195] = CONSTANTS[192]; CONSTANTS[196] = CONSTANTS[192]; CONSTANTS[197] = CONSTANTS[193]; CONSTANTS[198] = CONSTANTS[193]; CONSTANTS[199] = CONSTANTS[103]*pow(1.00000, 0.250000)*(pow(CONSTANTS[6], 0.750000)/(1.00000+pow(CONSTANTS[46]/CONSTANTS[111], CONSTANTS[112]))); CONSTANTS[200] = CONSTANTS[124]*CONSTANTS[129]; CONSTANTS[201] = CONSTANTS[119]*CONSTANTS[123]; CONSTANTS[202] = CONSTANTS[120]*CONSTANTS[126]+ (1.00000 - CONSTANTS[120])*CONSTANTS[127]; CONSTANTS[203] = CONSTANTS[128]; CONSTANTS[204] = CONSTANTS[202]; CONSTANTS[205] = CONSTANTS[202]; CONSTANTS[206] = CONSTANTS[202]; CONSTANTS[207] = CONSTANTS[203]; CONSTANTS[208] = CONSTANTS[203]; CONSTANTS[209] = CONSTANTS[122]*pow(1.00000, 0.250000)*(pow(CONSTANTS[6], 0.750000)/(1.00000+pow(CONSTANTS[46]/CONSTANTS[131], CONSTANTS[130]))); CONSTANTS[210] = CONSTANTS[135]*CONSTANTS[140]; CONSTANTS[211] = CONSTANTS[119]*CONSTANTS[134]; CONSTANTS[212] = CONSTANTS[120]*CONSTANTS[137]+ (1.00000 - CONSTANTS[120])*CONSTANTS[138]; CONSTANTS[213] = CONSTANTS[139]; CONSTANTS[214] = CONSTANTS[212]; CONSTANTS[215] = CONSTANTS[212]; CONSTANTS[216] = CONSTANTS[212]; CONSTANTS[217] = CONSTANTS[213]; CONSTANTS[218] = CONSTANTS[213]; CONSTANTS[219] = CONSTANTS[133]*pow(1.00000, 0.250000)*(pow(CONSTANTS[6], 0.750000)/(1.00000+pow(CONSTANTS[46]/CONSTANTS[142], CONSTANTS[141]))); CONSTANTS[220] = CONSTANTS[146]*CONSTANTS[151]; CONSTANTS[221] = CONSTANTS[119]*CONSTANTS[145]; CONSTANTS[222] = CONSTANTS[120]*CONSTANTS[148]+ (1.00000 - CONSTANTS[120])*CONSTANTS[149]; CONSTANTS[223] = CONSTANTS[150]; CONSTANTS[224] = CONSTANTS[222]; CONSTANTS[225] = CONSTANTS[222]; CONSTANTS[226] = CONSTANTS[222]; CONSTANTS[227] = CONSTANTS[223]; CONSTANTS[228] = CONSTANTS[223]; CONSTANTS[229] = CONSTANTS[144]*pow(1.00000, 0.250000)*(pow(CONSTANTS[6], 0.750000)/(1.00000+pow(CONSTANTS[46]/CONSTANTS[153], CONSTANTS[152]))); RATES[0] = 0.1001; RATES[3] = 0.1001; RATES[4] = 0.1001; RATES[5] = 0.1001; RATES[6] = 0.1001; RATES[7] = 0.1001; RATES[8] = 0.1001; RATES[9] = 0.1001; RATES[10] = 0.1001; RATES[11] = 0.1001; RATES[1] = 0.1001; RATES[12] = 0.1001; RATES[2] = 0.1001; RATES[13] = 0.1001; RATES[14] = 0.1001; RATES[15] = 0.1001; RATES[16] = 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[155] - (ALGEBRAIC[57]+ALGEBRAIC[1]); resid[1] = RATES[3] - (CONSTANTS[176]+ALGEBRAIC[58]+ALGEBRAIC[59]+CONSTANTS[179]+ALGEBRAIC[167]) - (ALGEBRAIC[60]+ALGEBRAIC[61]+ALGEBRAIC[62]+ALGEBRAIC[63]+ALGEBRAIC[4]); resid[2] = RATES[4] - CONSTANTS[157] - ALGEBRAIC[10]; resid[3] = RATES[5] - (CONSTANTS[180]+ALGEBRAIC[65]+ALGEBRAIC[66]+ALGEBRAIC[64]+ALGEBRAIC[168]) - (ALGEBRAIC[12]+ALGEBRAIC[67]+ALGEBRAIC[68]+ALGEBRAIC[13]); resid[4] = RATES[6] - (CONSTANTS[182]+ALGEBRAIC[69]) - ALGEBRAIC[18]; resid[5] = RATES[7] - CONSTANTS[183] - (ALGEBRAIC[70]+ALGEBRAIC[20]); resid[6] = RATES[8] - CONSTANTS[185] - (ALGEBRAIC[71]+ALGEBRAIC[23]); resid[7] = RATES[9] - CONSTANTS[163] - ALGEBRAIC[26]; resid[8] = RATES[10] - (CONSTANTS[187]+CONSTANTS[188]+ALGEBRAIC[72]+ALGEBRAIC[73]+CONSTANTS[189]) - (ALGEBRAIC[74]+ALGEBRAIC[75]+ALGEBRAIC[169]+ALGEBRAIC[76]+ALGEBRAIC[77]+ALGEBRAIC[28]); resid[9] = RATES[11] - ALGEBRAIC[78] - (ALGEBRAIC[79]+ALGEBRAIC[80]+ALGEBRAIC[81]+ALGEBRAIC[32]); resid[10] = RATES[1] - (ALGEBRAIC[82]+ALGEBRAIC[83]+ALGEBRAIC[170]) - (ALGEBRAIC[35]+ALGEBRAIC[138]+ALGEBRAIC[37]); resid[11] = RATES[12] - ALGEBRAIC[84] - (ALGEBRAIC[39]+ALGEBRAIC[40]+ALGEBRAIC[41]); resid[12] = RATES[2] - (ALGEBRAIC[85]+ALGEBRAIC[86]) - (ALGEBRAIC[87]+ALGEBRAIC[43]); resid[13] = RATES[13] - (CONSTANTS[190]+ALGEBRAIC[139]+ALGEBRAIC[140]+ALGEBRAIC[141]+ALGEBRAIC[142]+ALGEBRAIC[143]+ALGEBRAIC[144]+ALGEBRAIC[88]+ALGEBRAIC[89]+ALGEBRAIC[90]+ALGEBRAIC[175]+ALGEBRAIC[145]) - (ALGEBRAIC[46]+ALGEBRAIC[47]); resid[14] = RATES[14] - (CONSTANTS[200]+ALGEBRAIC[146]+ALGEBRAIC[147]+ALGEBRAIC[148]+ALGEBRAIC[149]+ALGEBRAIC[150]+ALGEBRAIC[151]+ALGEBRAIC[100]+ALGEBRAIC[101]+ALGEBRAIC[102]+ALGEBRAIC[176]+ALGEBRAIC[152]) - (ALGEBRAIC[49]+ALGEBRAIC[50]); resid[15] = RATES[15] - (CONSTANTS[210]+ALGEBRAIC[153]+ALGEBRAIC[154]+ALGEBRAIC[155]+ALGEBRAIC[156]+ALGEBRAIC[157]+ALGEBRAIC[158]+ALGEBRAIC[112]+ALGEBRAIC[113]+ALGEBRAIC[114]+ALGEBRAIC[177]+ALGEBRAIC[159]) - (ALGEBRAIC[52]+ALGEBRAIC[53]); resid[16] = RATES[16] - (CONSTANTS[220]+ALGEBRAIC[160]+ALGEBRAIC[161]+ALGEBRAIC[162]+ALGEBRAIC[163]+ALGEBRAIC[164]+ALGEBRAIC[165]+ALGEBRAIC[124]+ALGEBRAIC[125]+ALGEBRAIC[126]+ALGEBRAIC[178]+ALGEBRAIC[166]) - (ALGEBRAIC[55]+ALGEBRAIC[56]); } void computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC) { ALGEBRAIC[9] = STATES[4]/CONSTANTS[6]; ALGEBRAIC[17] = STATES[6]/CONSTANTS[6]; ALGEBRAIC[25] = STATES[9]/CONSTANTS[6]; ALGEBRAIC[36] = CONSTANTS[167]*STATES[1]*(1.00000 - CONSTANTS[47]); ALGEBRAIC[38] = STATES[12]/CONSTANTS[6]; } void computeEssentialVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC) { ALGEBRAIC[1] = CONSTANTS[154]*STATES[0]; ALGEBRAIC[4] = CONSTANTS[156]*STATES[3]; ALGEBRAIC[10] = CONSTANTS[158]*STATES[4]; ALGEBRAIC[11] = STATES[5]/CONSTANTS[6]; ALGEBRAIC[12] = CONSTANTS[181]/(1.00000+CONSTANTS[39]/ALGEBRAIC[11]); ALGEBRAIC[13] = CONSTANTS[159]*STATES[5]; ALGEBRAIC[18] = CONSTANTS[160]*STATES[6]; ALGEBRAIC[20] = CONSTANTS[161]*STATES[7]; ALGEBRAIC[23] = CONSTANTS[162]*STATES[8]; ALGEBRAIC[26] = CONSTANTS[164]*STATES[9]; ALGEBRAIC[28] = CONSTANTS[165]*STATES[10]; ALGEBRAIC[32] = CONSTANTS[166]*STATES[11]; ALGEBRAIC[35] = CONSTANTS[177]*STATES[1]*CONSTANTS[47]; ALGEBRAIC[37] = CONSTANTS[175]*STATES[1]*CONSTANTS[47]; ALGEBRAIC[39] = CONSTANTS[178]*STATES[12]*CONSTANTS[47]; ALGEBRAIC[40] = CONSTANTS[168]*STATES[12]*(1.00000 - CONSTANTS[47]); ALGEBRAIC[41] = CONSTANTS[175]*STATES[12]*CONSTANTS[47]; ALGEBRAIC[43] = CONSTANTS[169]*STATES[2]; ALGEBRAIC[45] = STATES[13]/CONSTANTS[6]; ALGEBRAIC[46] = CONSTANTS[199]/(1.00000+CONSTANTS[102]/ALGEBRAIC[45]); ALGEBRAIC[47] = CONSTANTS[170]*STATES[13]; ALGEBRAIC[48] = STATES[14]/CONSTANTS[6]; ALGEBRAIC[49] = CONSTANTS[209]/(1.00000+CONSTANTS[121]/ALGEBRAIC[48]); ALGEBRAIC[50] = CONSTANTS[171]*STATES[14]; ALGEBRAIC[51] = STATES[15]/CONSTANTS[6]; ALGEBRAIC[52] = CONSTANTS[219]/(1.00000+CONSTANTS[132]/ALGEBRAIC[51]); ALGEBRAIC[53] = CONSTANTS[172]*STATES[15]; ALGEBRAIC[54] = STATES[16]/CONSTANTS[6]; ALGEBRAIC[55] = CONSTANTS[229]/(1.00000+CONSTANTS[143]/ALGEBRAIC[54]); ALGEBRAIC[56] = CONSTANTS[173]*STATES[16]; ALGEBRAIC[0] = STATES[0]/CONSTANTS[6]; ALGEBRAIC[2] = CONSTANTS[0]*(STATES[1]+STATES[2]); ALGEBRAIC[57] = ALGEBRAIC[2]/(1.00000+CONSTANTS[174]/ALGEBRAIC[0]); ALGEBRAIC[58] = CONSTANTS[14]*ALGEBRAIC[57]; ALGEBRAIC[59] = CONSTANTS[15]*ALGEBRAIC[35]; ALGEBRAIC[3] = STATES[3]/CONSTANTS[6]; ALGEBRAIC[5] = CONSTANTS[9]*STATES[1]; ALGEBRAIC[27] = STATES[10]/CONSTANTS[6]; ALGEBRAIC[60] = ALGEBRAIC[5]/(1.00000+CONSTANTS[18]/ALGEBRAIC[3]+ALGEBRAIC[27]/CONSTANTS[23]); ALGEBRAIC[6] = CONSTANTS[9]*STATES[2]; ALGEBRAIC[31] = STATES[11]/CONSTANTS[6]; ALGEBRAIC[61] = ALGEBRAIC[6]/(1.00000+CONSTANTS[18]/ALGEBRAIC[3]+ALGEBRAIC[31]/CONSTANTS[24]); ALGEBRAIC[7] = CONSTANTS[10]*STATES[1]; ALGEBRAIC[62] = ALGEBRAIC[7]/(1.00000+CONSTANTS[19]/ALGEBRAIC[3]+CONSTANTS[21]/ALGEBRAIC[27]); ALGEBRAIC[8] = CONSTANTS[11]*STATES[2]; ALGEBRAIC[63] = ALGEBRAIC[8]/(1.00000+CONSTANTS[20]/ALGEBRAIC[3]+CONSTANTS[22]/ALGEBRAIC[31]); ALGEBRAIC[14] = CONSTANTS[29]*CONSTANTS[6]*(CONSTANTS[48]/(1.00000+ALGEBRAIC[11]/CONSTANTS[44])); ALGEBRAIC[64] = CONSTANTS[34]*ALGEBRAIC[14]; ALGEBRAIC[65] = CONSTANTS[35]*ALGEBRAIC[60]; ALGEBRAIC[66] = CONSTANTS[35]*ALGEBRAIC[61]; ALGEBRAIC[15] = CONSTANTS[31]*STATES[1]*(1.00000 - CONSTANTS[47]); ALGEBRAIC[67] = ALGEBRAIC[15]/(1.00000+CONSTANTS[37]/ALGEBRAIC[11]+CONSTANTS[42]/ALGEBRAIC[27]); ALGEBRAIC[16] = CONSTANTS[32]*STATES[2]; ALGEBRAIC[68] = ALGEBRAIC[16]/(1.00000+CONSTANTS[38]/ALGEBRAIC[11]+CONSTANTS[43]/ALGEBRAIC[31]); ALGEBRAIC[69] = CONSTANTS[50]*ALGEBRAIC[35]; ALGEBRAIC[19] = STATES[7]/CONSTANTS[6]; ALGEBRAIC[21] = CONSTANTS[54]*STATES[1]; ALGEBRAIC[70] = ALGEBRAIC[21]/(1.00000+CONSTANTS[184]/ALGEBRAIC[19]); ALGEBRAIC[22] = STATES[8]/CONSTANTS[6]; ALGEBRAIC[24] = CONSTANTS[60]*STATES[2]*((1.00000 - ( VOI*CONSTANTS[67])/1.00000)+(( VOI*CONSTANTS[67])/1.00000)/(1.00000+pow(CONSTANTS[65]/CONSTANTS[66], CONSTANTS[64]))); ALGEBRAIC[71] = ALGEBRAIC[24]/(1.00000+CONSTANTS[186]/ALGEBRAIC[22]); ALGEBRAIC[72] = CONSTANTS[80]*ALGEBRAIC[70]; ALGEBRAIC[73] = CONSTANTS[81]*ALGEBRAIC[39]; ALGEBRAIC[74] = CONSTANTS[70]*ALGEBRAIC[67]; ALGEBRAIC[75] = CONSTANTS[71]*ALGEBRAIC[62]; ALGEBRAIC[29] = CONSTANTS[75]*STATES[1]; ALGEBRAIC[76] = ALGEBRAIC[29]/(1.00000+CONSTANTS[73]/ALGEBRAIC[27]); ALGEBRAIC[30] = CONSTANTS[76]*STATES[1]; ALGEBRAIC[77] = ALGEBRAIC[30]*((1.00000 - CONSTANTS[47])/( (1.00000+CONSTANTS[74]/ALGEBRAIC[27])*(1.00000+ALGEBRAIC[11]/CONSTANTS[85])*(ALGEBRAIC[3]/CONSTANTS[86])))+( ALGEBRAIC[30]*CONSTANTS[47])/(1.00000+CONSTANTS[74]/ALGEBRAIC[27]); ALGEBRAIC[78] = CONSTANTS[90]*ALGEBRAIC[71]; ALGEBRAIC[79] = CONSTANTS[87]*ALGEBRAIC[68]; ALGEBRAIC[80] = CONSTANTS[88]*ALGEBRAIC[63]; ALGEBRAIC[33] = CONSTANTS[89]*STATES[2]; ALGEBRAIC[81] = ALGEBRAIC[33]/( (1.00000+CONSTANTS[91]/ALGEBRAIC[31])*(1.00000+ALGEBRAIC[11]/CONSTANTS[92]+ALGEBRAIC[3]/CONSTANTS[93])); ALGEBRAIC[82] = CONSTANTS[94]*ALGEBRAIC[67]; ALGEBRAIC[83] = CONSTANTS[95]*ALGEBRAIC[62]; ALGEBRAIC[84] = CONSTANTS[97]*ALGEBRAIC[76]; ALGEBRAIC[85] = CONSTANTS[98]*ALGEBRAIC[68]; ALGEBRAIC[86] = CONSTANTS[99]*ALGEBRAIC[63]; ALGEBRAIC[34] = STATES[1]/CONSTANTS[6]; ALGEBRAIC[42] = STATES[2]/CONSTANTS[6]; ALGEBRAIC[44] = ( CONSTANTS[100]*STATES[1]*CONSTANTS[47]*ALGEBRAIC[42])/(ALGEBRAIC[42]+ ALGEBRAIC[34]*(1.00000 - CONSTANTS[47])); ALGEBRAIC[87] = ALGEBRAIC[44]/(1.00000+CONSTANTS[101]/(ALGEBRAIC[42]+ ALGEBRAIC[34]*(1.00000 - CONSTANTS[47]))); ALGEBRAIC[88] = CONSTANTS[192]*ALGEBRAIC[81]; ALGEBRAIC[89] = CONSTANTS[197]*ALGEBRAIC[60]; ALGEBRAIC[90] = CONSTANTS[197]*ALGEBRAIC[61]; ALGEBRAIC[100] = CONSTANTS[202]*ALGEBRAIC[81]; ALGEBRAIC[101] = CONSTANTS[207]*ALGEBRAIC[60]; ALGEBRAIC[102] = CONSTANTS[207]*ALGEBRAIC[61]; ALGEBRAIC[112] = CONSTANTS[212]*ALGEBRAIC[81]; ALGEBRAIC[113] = CONSTANTS[217]*ALGEBRAIC[60]; ALGEBRAIC[114] = CONSTANTS[217]*ALGEBRAIC[61]; ALGEBRAIC[124] = CONSTANTS[222]*ALGEBRAIC[81]; ALGEBRAIC[125] = CONSTANTS[227]*ALGEBRAIC[60]; ALGEBRAIC[126] = CONSTANTS[227]*ALGEBRAIC[61]; ALGEBRAIC[138] = ALGEBRAIC[87]/(1.00000+CONSTANTS[8]/ALGEBRAIC[27]); ALGEBRAIC[91] = CONSTANTS[113]*ALGEBRAIC[74]; ALGEBRAIC[96] = ( CONSTANTS[104]*CONSTANTS[187]+ CONSTANTS[191]*CONSTANTS[189]+ CONSTANTS[106]*(CONSTANTS[188]+ALGEBRAIC[72]+ALGEBRAIC[73]))/(CONSTANTS[187]+CONSTANTS[189]+CONSTANTS[188]+ALGEBRAIC[72]+ALGEBRAIC[73]); ALGEBRAIC[97] = ALGEBRAIC[96]; ALGEBRAIC[139] = ALGEBRAIC[97]*ALGEBRAIC[91]; ALGEBRAIC[92] = CONSTANTS[114]*ALGEBRAIC[79]; ALGEBRAIC[140] = CONSTANTS[194]*ALGEBRAIC[92]; ALGEBRAIC[93] = CONSTANTS[115]*ALGEBRAIC[75]; ALGEBRAIC[98] = ALGEBRAIC[96]; ALGEBRAIC[141] = ALGEBRAIC[98]*ALGEBRAIC[93]; ALGEBRAIC[94] = CONSTANTS[116]*ALGEBRAIC[80]; ALGEBRAIC[142] = CONSTANTS[195]*ALGEBRAIC[94]; ALGEBRAIC[95] = CONSTANTS[117]*ALGEBRAIC[76]; ALGEBRAIC[99] = ALGEBRAIC[96]; ALGEBRAIC[143] = ALGEBRAIC[99]*ALGEBRAIC[95]; ALGEBRAIC[144] = ALGEBRAIC[96]*ALGEBRAIC[77]; ALGEBRAIC[137] = ( 1.00000*ALGEBRAIC[87])/(1.00000+ALGEBRAIC[27]/CONSTANTS[26]); ALGEBRAIC[145] = CONSTANTS[198]*ALGEBRAIC[137]; ALGEBRAIC[103] = CONSTANTS[113]*ALGEBRAIC[74]; ALGEBRAIC[108] = ( CONSTANTS[123]*CONSTANTS[187]+ CONSTANTS[201]*CONSTANTS[189]+ CONSTANTS[125]*(CONSTANTS[188]+ALGEBRAIC[72]+ALGEBRAIC[73]))/(CONSTANTS[187]+CONSTANTS[189]+CONSTANTS[188]+ALGEBRAIC[72]+ALGEBRAIC[73]); ALGEBRAIC[109] = ALGEBRAIC[108]; ALGEBRAIC[146] = ALGEBRAIC[109]*ALGEBRAIC[103]; ALGEBRAIC[104] = CONSTANTS[114]*ALGEBRAIC[79]; ALGEBRAIC[147] = CONSTANTS[204]*ALGEBRAIC[104]; ALGEBRAIC[105] = CONSTANTS[115]*ALGEBRAIC[75]; ALGEBRAIC[110] = ALGEBRAIC[108]; ALGEBRAIC[148] = ALGEBRAIC[110]*ALGEBRAIC[105]; ALGEBRAIC[106] = CONSTANTS[116]*ALGEBRAIC[80]; ALGEBRAIC[149] = CONSTANTS[205]*ALGEBRAIC[106]; ALGEBRAIC[107] = CONSTANTS[117]*ALGEBRAIC[76]; ALGEBRAIC[111] = ALGEBRAIC[108]; ALGEBRAIC[150] = ALGEBRAIC[111]*ALGEBRAIC[107]; ALGEBRAIC[151] = ALGEBRAIC[108]*ALGEBRAIC[77]; ALGEBRAIC[152] = CONSTANTS[208]*ALGEBRAIC[137]; ALGEBRAIC[115] = CONSTANTS[113]*ALGEBRAIC[74]; ALGEBRAIC[120] = ( CONSTANTS[134]*CONSTANTS[187]+ CONSTANTS[211]*CONSTANTS[189]+ CONSTANTS[136]*(CONSTANTS[188]+ALGEBRAIC[72]+ALGEBRAIC[73]))/(CONSTANTS[187]+CONSTANTS[189]+CONSTANTS[188]+ALGEBRAIC[72]+ALGEBRAIC[73]); ALGEBRAIC[121] = ALGEBRAIC[120]; ALGEBRAIC[153] = ALGEBRAIC[121]*ALGEBRAIC[115]; ALGEBRAIC[116] = CONSTANTS[114]*ALGEBRAIC[79]; ALGEBRAIC[154] = CONSTANTS[214]*ALGEBRAIC[116]; ALGEBRAIC[117] = CONSTANTS[115]*ALGEBRAIC[75]; ALGEBRAIC[122] = ALGEBRAIC[120]; ALGEBRAIC[155] = ALGEBRAIC[122]*ALGEBRAIC[117]; ALGEBRAIC[118] = CONSTANTS[116]*ALGEBRAIC[80]; ALGEBRAIC[156] = CONSTANTS[215]*ALGEBRAIC[118]; ALGEBRAIC[119] = CONSTANTS[117]*ALGEBRAIC[76]; ALGEBRAIC[123] = ALGEBRAIC[120]; ALGEBRAIC[157] = ALGEBRAIC[123]*ALGEBRAIC[119]; ALGEBRAIC[158] = ALGEBRAIC[120]*ALGEBRAIC[77]; ALGEBRAIC[159] = CONSTANTS[218]*ALGEBRAIC[137]; ALGEBRAIC[127] = CONSTANTS[113]*ALGEBRAIC[74]; ALGEBRAIC[132] = ( CONSTANTS[145]*CONSTANTS[187]+ CONSTANTS[221]*CONSTANTS[189]+ CONSTANTS[147]*(CONSTANTS[188]+ALGEBRAIC[72]+ALGEBRAIC[73]))/(CONSTANTS[187]+CONSTANTS[189]+CONSTANTS[188]+ALGEBRAIC[72]+ALGEBRAIC[73]); ALGEBRAIC[133] = ALGEBRAIC[132]; ALGEBRAIC[160] = ALGEBRAIC[133]*ALGEBRAIC[127]; ALGEBRAIC[128] = CONSTANTS[114]*ALGEBRAIC[79]; ALGEBRAIC[161] = CONSTANTS[224]*ALGEBRAIC[128]; ALGEBRAIC[129] = CONSTANTS[115]*ALGEBRAIC[75]; ALGEBRAIC[134] = ALGEBRAIC[132]; ALGEBRAIC[162] = ALGEBRAIC[134]*ALGEBRAIC[129]; ALGEBRAIC[130] = CONSTANTS[116]*ALGEBRAIC[80]; ALGEBRAIC[163] = CONSTANTS[225]*ALGEBRAIC[130]; ALGEBRAIC[131] = CONSTANTS[117]*ALGEBRAIC[76]; ALGEBRAIC[135] = ALGEBRAIC[132]; ALGEBRAIC[164] = ALGEBRAIC[135]*ALGEBRAIC[131]; ALGEBRAIC[165] = ALGEBRAIC[132]*ALGEBRAIC[77]; ALGEBRAIC[166] = CONSTANTS[228]*ALGEBRAIC[137]; ALGEBRAIC[136] = ALGEBRAIC[87]*(1.00000 - (1.00000/(1.00000+ALGEBRAIC[27]/CONSTANTS[26])+1.00000/(1.00000+CONSTANTS[8]/ALGEBRAIC[27]))); ALGEBRAIC[167] = CONSTANTS[16]*ALGEBRAIC[136]; ALGEBRAIC[168] = CONSTANTS[36]*ALGEBRAIC[137]; ALGEBRAIC[169] = CONSTANTS[72]*ALGEBRAIC[138]; ALGEBRAIC[170] = CONSTANTS[96]*ALGEBRAIC[138]; ALGEBRAIC[171] = CONSTANTS[118]*ALGEBRAIC[169]; ALGEBRAIC[175] = CONSTANTS[196]*ALGEBRAIC[171]; ALGEBRAIC[172] = CONSTANTS[118]*ALGEBRAIC[169]; ALGEBRAIC[176] = CONSTANTS[206]*ALGEBRAIC[172]; ALGEBRAIC[173] = CONSTANTS[118]*ALGEBRAIC[169]; ALGEBRAIC[177] = CONSTANTS[216]*ALGEBRAIC[173]; ALGEBRAIC[174] = CONSTANTS[118]*ALGEBRAIC[169]; ALGEBRAIC[178] = CONSTANTS[226]*ALGEBRAIC[174]; } 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; } void computeRoots(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES, double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS) { }