Generated Code
The following is c_ida code generated by the CellML API from this CellML file. (Back to language selection)
The raw code is available.
/*
There are a total of 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)
{
}