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 164 entries in the algebraic variable array.
There are a total of 60 entries in each of the rate and state variable arrays.
There are a total of 139 entries in the constant variable array.
*/
/*
* CONSTANTS[0] is kappa_1a in component BG_parameters (fmol_per_sec).
* CONSTANTS[1] is kappa_1b in component BG_parameters (fmol_per_sec).
* CONSTANTS[2] is kappa_2a in component BG_parameters (fmol_per_sec).
* CONSTANTS[3] is kappa_2b in component BG_parameters (fmol_per_sec).
* CONSTANTS[4] is kappa_3a in component BG_parameters (fmol_per_sec).
* CONSTANTS[5] is kappa_3b in component BG_parameters (fmol_per_sec).
* CONSTANTS[6] is kappa_4a in component BG_parameters (fmol_per_sec).
* CONSTANTS[7] is kappa_4b in component BG_parameters (fmol_per_sec).
* CONSTANTS[8] is kappa_5 in component BG_parameters (fmol_per_sec).
* CONSTANTS[9] is kappa_6 in component BG_parameters (fmol_per_sec).
* CONSTANTS[10] is kappa_7 in component BG_parameters (fmol_per_sec).
* CONSTANTS[11] is kappa_GiAC in component BG_parameters (fmol_per_sec).
* CONSTANTS[12] is kappa_sig1_B1 in component BG_parameters (fmol_per_sec).
* CONSTANTS[13] is kappa_sig2_B1 in component BG_parameters (fmol_per_sec).
* CONSTANTS[14] is kappa_sig3_B1 in component BG_parameters (fmol_per_sec).
* CONSTANTS[15] is kappa_sig4_B1 in component BG_parameters (fmol_per_sec).
* CONSTANTS[16] is kappa_sig1_M2 in component BG_parameters (fmol_per_sec).
* CONSTANTS[17] is kappa_sig2_M2 in component BG_parameters (fmol_per_sec).
* CONSTANTS[18] is kappa_sig3_M2 in component BG_parameters (fmol_per_sec).
* CONSTANTS[19] is kappa_sig4_M2 in component BG_parameters (fmol_per_sec).
* CONSTANTS[20] is kappa_Doff1_Gs in component BG_parameters (fmol_per_sec).
* CONSTANTS[21] is kappa_Ton1_Gs in component BG_parameters (fmol_per_sec).
* CONSTANTS[22] is kappa_Act1_Gs in component BG_parameters (fmol_per_sec).
* CONSTANTS[23] is kappa_Doff2_Gs in component BG_parameters (fmol_per_sec).
* CONSTANTS[24] is kappa_Ton2_Gs in component BG_parameters (fmol_per_sec).
* CONSTANTS[25] is kappa_Act2_Gs in component BG_parameters (fmol_per_sec).
* CONSTANTS[26] is kappa_hyd_Gs in component BG_parameters (fmol_per_sec).
* CONSTANTS[27] is kappa_reassoc_Gs in component BG_parameters (fmol_per_sec).
* CONSTANTS[28] is kappa_Doff1_Gi in component BG_parameters (fmol_per_sec).
* CONSTANTS[29] is kappa_Ton1_Gi in component BG_parameters (fmol_per_sec).
* CONSTANTS[30] is kappa_Act1_Gi in component BG_parameters (fmol_per_sec).
* CONSTANTS[31] is kappa_Doff2_Gi in component BG_parameters (fmol_per_sec).
* CONSTANTS[32] is kappa_Ton2_Gi in component BG_parameters (fmol_per_sec).
* CONSTANTS[33] is kappa_Act2_Gi in component BG_parameters (fmol_per_sec).
* CONSTANTS[34] is kappa_hyd_Gi in component BG_parameters (fmol_per_sec).
* CONSTANTS[35] is kappa_reassoc_Gi in component BG_parameters (fmol_per_sec).
* CONSTANTS[36] is K_ATP in component BG_parameters (per_fmol).
* CONSTANTS[37] is K_cAMP in component BG_parameters (per_fmol).
* CONSTANTS[38] is K_AC in component BG_parameters (per_fmol).
* CONSTANTS[39] is K_AC_ATP in component BG_parameters (per_fmol).
* CONSTANTS[40] is K_a_Gs_GTP_AC in component BG_parameters (per_fmol).
* CONSTANTS[41] is K_a_Gs_GTP_AC_ATP in component BG_parameters (per_fmol).
* CONSTANTS[42] is K_FSK_AC in component BG_parameters (per_fmol).
* CONSTANTS[43] is K_FSK_AC_ATP in component BG_parameters (per_fmol).
* CONSTANTS[44] is K_PDE in component BG_parameters (per_fmol).
* CONSTANTS[45] is K_PDE_cAMP in component BG_parameters (per_fmol).
* CONSTANTS[46] is K_five_AMP in component BG_parameters (per_fmol).
* CONSTANTS[47] is K_IBMX in component BG_parameters (per_fmol).
* CONSTANTS[48] is K_PDEinh in component BG_parameters (per_fmol).
* CONSTANTS[49] is K_a_Gs_GTP in component BG_parameters (per_fmol).
* CONSTANTS[50] is K_FSK in component BG_parameters (per_fmol).
* CONSTANTS[51] is K_a_Gi_GTP in component BG_parameters (per_fmol).
* CONSTANTS[52] is K_ACinh in component BG_parameters (per_fmol).
* CONSTANTS[53] is K_PPi in component BG_parameters (per_fmol).
* CONSTANTS[54] is K_L_B1 in component BG_parameters (per_fmol).
* CONSTANTS[55] is K_R_B1 in component BG_parameters (per_fmol).
* CONSTANTS[56] is K_Gs in component BG_parameters (per_fmol).
* CONSTANTS[57] is K_LR_B1 in component BG_parameters (per_fmol).
* CONSTANTS[58] is K_R_B1Gs in component BG_parameters (per_fmol).
* CONSTANTS[59] is K_LR_B1Gs in component BG_parameters (per_fmol).
* CONSTANTS[60] is K_L_M2 in component BG_parameters (per_fmol).
* CONSTANTS[61] is K_R_M2 in component BG_parameters (per_fmol).
* CONSTANTS[62] is K_Gi in component BG_parameters (per_fmol).
* CONSTANTS[63] is K_LR_M2 in component BG_parameters (per_fmol).
* CONSTANTS[64] is K_R_M2Gi in component BG_parameters (per_fmol).
* CONSTANTS[65] is K_LR_M2Gi in component BG_parameters (per_fmol).
* CONSTANTS[66] is K_beta_gamma_Gs in component BG_parameters (per_fmol).
* CONSTANTS[67] is K_a_Gs_GDP in component BG_parameters (per_fmol).
* CONSTANTS[68] is K_Pi in component BG_parameters (per_fmol).
* CONSTANTS[69] is K_R_B1_aby in component BG_parameters (per_fmol).
* CONSTANTS[70] is K_R_B1_aby_T in component BG_parameters (per_fmol).
* CONSTANTS[71] is K_LR_B1_aby in component BG_parameters (per_fmol).
* CONSTANTS[72] is K_LR_B1_aby_T in component BG_parameters (per_fmol).
* CONSTANTS[73] is K_GTP in component BG_parameters (per_fmol).
* CONSTANTS[74] is K_GDP in component BG_parameters (per_fmol).
* CONSTANTS[75] is K_beta_gamma_Gi in component BG_parameters (per_fmol).
* CONSTANTS[76] is K_a_Gi_GDP in component BG_parameters (per_fmol).
* CONSTANTS[77] is K_R_M2_aby in component BG_parameters (per_fmol).
* CONSTANTS[78] is K_R_M2_aby_T in component BG_parameters (per_fmol).
* CONSTANTS[79] is K_LR_M2_aby in component BG_parameters (per_fmol).
* CONSTANTS[80] is K_LR_M2_aby_T in component BG_parameters (per_fmol).
* VOI is time in component environment (second).
* CONSTANTS[81] is vol_myo in component environment (pL).
* CONSTANTS[82] is freq in component environment (dimensionless).
* CONSTANTS[83] is stimSt in component environment (second).
* CONSTANTS[84] is stimSt2 in component environment (second).
* CONSTANTS[85] is stimDur in component environment (second).
* CONSTANTS[86] is tRamp in component environment (second).
* CONSTANTS[87] is stimMag in component environment (fmol).
* CONSTANTS[88] is stimHolding in component environment (fmol).
* CONSTANTS[138] is m in component environment (fmol_per_sec).
* CONSTANTS[89] is q_ATP_init in component environment (fmol).
* CONSTANTS[90] is q_AC_init in component environment (fmol).
* CONSTANTS[91] is q_cAMP_init in component environment (fmol).
* CONSTANTS[92] is q_AC_ATP_init in component environment (fmol).
* CONSTANTS[93] is q_FSK_init in component environment (fmol).
* CONSTANTS[94] is q_FSK_AC_init in component environment (fmol).
* CONSTANTS[95] is q_FSK_AC_ATP_init in component environment (fmol).
* CONSTANTS[96] is q_a_Gs_GTP_init in component environment (fmol).
* CONSTANTS[97] is q_a_Gs_GTP_AC_init in component environment (fmol).
* CONSTANTS[98] is q_a_Gs_GTP_AC_ATP_init in component environment (fmol).
* CONSTANTS[99] is q_PDE_init in component environment (fmol).
* CONSTANTS[100] is q_PDEinh_init in component environment (fmol).
* CONSTANTS[101] is q_PDE_cAMP_init in component environment (fmol).
* CONSTANTS[102] is q_IBMX_init in component environment (fmol).
* CONSTANTS[103] is q_five_AMP_init in component environment (fmol).
* CONSTANTS[104] is q_a_Gi_GTP_init in component environment (fmol).
* CONSTANTS[105] is q_ACinh_init in component environment (fmol).
* CONSTANTS[106] is q_PPi_init in component environment (fmol).
* CONSTANTS[107] is q_R_B1_init in component environment (fmol).
* CONSTANTS[108] is q_Gs_init in component environment (fmol).
* ALGEBRAIC[0] is q_L_B1_init in component environment (fmol).
* CONSTANTS[109] is q_LR_B1_init in component environment (fmol).
* CONSTANTS[110] is q_R_B1Gs_init in component environment (fmol).
* CONSTANTS[111] is q_LR_B1Gs_init in component environment (fmol).
* ALGEBRAIC[1] is q_L_M2_init in component environment (fmol).
* CONSTANTS[112] is q_R_M2_init in component environment (fmol).
* CONSTANTS[113] is q_Gi_init in component environment (fmol).
* CONSTANTS[114] is q_LR_M2_init in component environment (fmol).
* CONSTANTS[115] is q_R_M2Gi_init in component environment (fmol).
* CONSTANTS[116] is q_LR_M2Gi_init in component environment (fmol).
* CONSTANTS[117] is q_beta_gamma_Gs_init in component environment (fmol).
* CONSTANTS[118] is q_a_Gs_GDP_init in component environment (fmol).
* CONSTANTS[119] is q_Pi_init in component environment (fmol).
* CONSTANTS[120] is q_R_B1_aby_init in component environment (fmol).
* CONSTANTS[121] is q_R_B1_aby_T_init in component environment (fmol).
* CONSTANTS[122] is q_LR_B1_aby_init in component environment (fmol).
* CONSTANTS[123] is q_LR_B1_aby_T_init in component environment (fmol).
* CONSTANTS[124] is q_GTP_init in component environment (fmol).
* CONSTANTS[125] is q_GDP_init in component environment (fmol).
* CONSTANTS[126] is q_beta_gamma_Gi_init in component environment (fmol).
* CONSTANTS[127] is q_a_Gi_GDP_init in component environment (fmol).
* CONSTANTS[128] is q_R_M2_aby_init in component environment (fmol).
* CONSTANTS[129] is q_R_M2_aby_T_init in component environment (fmol).
* CONSTANTS[130] is q_LR_M2_aby_init in component environment (fmol).
* CONSTANTS[131] is q_LR_M2_aby_T_init in component environment (fmol).
* ALGEBRAIC[157] is L_B1_T in component environment (fmol).
* ALGEBRAIC[158] is L_M2_T in component environment (fmol).
* ALGEBRAIC[159] is R_B1_T in component environment (fmol).
* ALGEBRAIC[160] is R_M2_T in component environment (fmol).
* ALGEBRAIC[161] is Gs_T in component environment (fmol).
* ALGEBRAIC[162] is Gi_T in component environment (fmol).
* ALGEBRAIC[163] is adenosine_T in component environment (fmol).
* ALGEBRAIC[2] is q_ATP in component environment (fmol).
* ALGEBRAIC[3] is q_cAMP in component environment (fmol).
* ALGEBRAIC[4] is q_AC in component environment (fmol).
* ALGEBRAIC[5] is q_AC_ATP in component environment (fmol).
* ALGEBRAIC[6] is q_a_Gs_GTP_AC in component environment (fmol).
* ALGEBRAIC[7] is q_a_Gs_GTP_AC_ATP in component environment (fmol).
* ALGEBRAIC[8] is q_FSK_AC in component environment (fmol).
* ALGEBRAIC[9] is q_FSK_AC_ATP in component environment (fmol).
* ALGEBRAIC[10] is q_PDE in component environment (fmol).
* ALGEBRAIC[11] is q_PDE_cAMP in component environment (fmol).
* ALGEBRAIC[12] is q_five_AMP in component environment (fmol).
* ALGEBRAIC[13] is q_IBMX in component environment (fmol).
* ALGEBRAIC[14] is q_PDEinh in component environment (fmol).
* ALGEBRAIC[15] is q_a_Gs_GTP in component environment (fmol).
* ALGEBRAIC[16] is q_FSK in component environment (fmol).
* ALGEBRAIC[17] is q_a_Gi_GTP in component environment (fmol).
* ALGEBRAIC[18] is q_ACinh in component environment (fmol).
* ALGEBRAIC[19] is q_PPi in component environment (fmol).
* ALGEBRAIC[20] is q_L_B1 in component environment (fmol).
* ALGEBRAIC[21] is q_R_B1 in component environment (fmol).
* ALGEBRAIC[22] is q_Gs in component environment (fmol).
* ALGEBRAIC[23] is q_LR_B1 in component environment (fmol).
* ALGEBRAIC[24] is q_R_B1Gs in component environment (fmol).
* ALGEBRAIC[25] is q_LR_B1Gs in component environment (fmol).
* ALGEBRAIC[26] is q_L_M2 in component environment (fmol).
* ALGEBRAIC[27] is q_R_M2 in component environment (fmol).
* ALGEBRAIC[28] is q_Gi in component environment (fmol).
* ALGEBRAIC[29] is q_LR_M2 in component environment (fmol).
* ALGEBRAIC[30] is q_R_M2Gi in component environment (fmol).
* ALGEBRAIC[31] is q_LR_M2Gi in component environment (fmol).
* ALGEBRAIC[32] is q_beta_gamma_Gs in component environment (fmol).
* ALGEBRAIC[33] is q_a_Gs_GDP in component environment (fmol).
* ALGEBRAIC[34] is q_Pi in component environment (fmol).
* ALGEBRAIC[35] is q_R_B1_aby in component environment (fmol).
* ALGEBRAIC[36] is q_R_B1_aby_T in component environment (fmol).
* ALGEBRAIC[37] is q_LR_B1_aby in component environment (fmol).
* ALGEBRAIC[38] is q_LR_B1_aby_T in component environment (fmol).
* ALGEBRAIC[39] is q_GTP in component environment (fmol).
* ALGEBRAIC[40] is q_GDP in component environment (fmol).
* ALGEBRAIC[41] is q_beta_gamma_Gi in component environment (fmol).
* ALGEBRAIC[42] is q_a_Gi_GDP in component environment (fmol).
* ALGEBRAIC[43] is q_R_M2_aby in component environment (fmol).
* ALGEBRAIC[44] is q_R_M2_aby_T in component environment (fmol).
* ALGEBRAIC[45] is q_LR_M2_aby in component environment (fmol).
* ALGEBRAIC[46] is q_LR_M2_aby_T in component environment (fmol).
* STATES[0] is q_ATP in component cAMP (fmol).
* STATES[1] is q_cAMP in component cAMP (fmol).
* STATES[2] is q_AC in component cAMP (fmol).
* STATES[3] is q_AC_ATP in component cAMP (fmol).
* STATES[4] is q_a_Gs_GTP_AC in component cAMP (fmol).
* STATES[5] is q_a_Gs_GTP_AC_ATP in component cAMP (fmol).
* STATES[6] is q_FSK_AC in component cAMP (fmol).
* STATES[7] is q_FSK_AC_ATP in component cAMP (fmol).
* STATES[8] is q_PDE in component cAMP (fmol).
* STATES[9] is q_PDE_cAMP in component cAMP (fmol).
* STATES[10] is q_five_AMP in component cAMP (fmol).
* STATES[11] is q_IBMX in component cAMP (fmol).
* STATES[12] is q_PDEinh in component cAMP (fmol).
* STATES[13] is q_a_Gs_GTP in component cAMP (fmol).
* STATES[14] is q_FSK in component cAMP (fmol).
* STATES[15] is q_a_Gi_GTP in component cAMP (fmol).
* STATES[16] is q_ACinh in component cAMP (fmol).
* STATES[17] is q_PPi in component cAMP (fmol).
* STATES[18] is q_L_B1 in component LRGbinding_B1AR (fmol).
* STATES[19] is q_R_B1 in component LRGbinding_B1AR (fmol).
* STATES[20] is q_Gs in component LRGbinding_B1AR (fmol).
* STATES[21] is q_LR_B1 in component LRGbinding_B1AR (fmol).
* STATES[22] is q_R_B1Gs in component LRGbinding_B1AR (fmol).
* STATES[23] is q_LR_B1Gs in component LRGbinding_B1AR (fmol).
* STATES[24] is q_L_M2 in component LRGbinding_M2 (fmol).
* STATES[25] is q_R_M2 in component LRGbinding_M2 (fmol).
* STATES[26] is q_Gi in component LRGbinding_M2 (fmol).
* STATES[27] is q_LR_M2 in component LRGbinding_M2 (fmol).
* STATES[28] is q_R_M2Gi in component LRGbinding_M2 (fmol).
* STATES[29] is q_LR_M2Gi in component LRGbinding_M2 (fmol).
* STATES[30] is q_R_B1 in component GsProtein (fmol).
* STATES[31] is q_Gs in component GsProtein (fmol).
* STATES[32] is q_R_B1Gs in component GsProtein (fmol).
* STATES[33] is q_LR_B1 in component GsProtein (fmol).
* STATES[34] is q_LR_B1Gs in component GsProtein (fmol).
* STATES[35] is q_a_Gs_GTP in component GsProtein (fmol).
* STATES[36] is q_beta_gamma_Gs in component GsProtein (fmol).
* STATES[37] is q_a_Gs_GDP in component GsProtein (fmol).
* STATES[38] is q_Pi in component GsProtein (fmol).
* STATES[39] is q_R_B1_aby in component GsProtein (fmol).
* STATES[40] is q_R_B1_aby_T in component GsProtein (fmol).
* STATES[41] is q_LR_B1_aby in component GsProtein (fmol).
* STATES[42] is q_LR_B1_aby_T in component GsProtein (fmol).
* STATES[43] is q_GTP in component GsProtein (fmol).
* STATES[44] is q_GDP in component GsProtein (fmol).
* STATES[45] is q_R_M2 in component GiProtein (fmol).
* STATES[46] is q_Gi in component GiProtein (fmol).
* STATES[47] is q_R_M2Gi in component GiProtein (fmol).
* STATES[48] is q_LR_M2 in component GiProtein (fmol).
* STATES[49] is q_LR_M2Gi in component GiProtein (fmol).
* STATES[50] is q_a_Gi_GTP in component GiProtein (fmol).
* STATES[51] is q_beta_gamma_Gi in component GiProtein (fmol).
* STATES[52] is q_a_Gi_GDP in component GiProtein (fmol).
* STATES[53] is q_Pi in component GiProtein (fmol).
* STATES[54] is q_R_M2_aby in component GiProtein (fmol).
* STATES[55] is q_R_M2_aby_T in component GiProtein (fmol).
* STATES[56] is q_LR_M2_aby in component GiProtein (fmol).
* STATES[57] is q_LR_M2_aby_T in component GiProtein (fmol).
* STATES[58] is q_GTP in component GiProtein (fmol).
* STATES[59] is q_GDP in component GiProtein (fmol).
* CONSTANTS[132] is R in component constants (J_per_K_per_mol).
* CONSTANTS[133] is T in component constants (kelvin).
* CONSTANTS[134] is F in component constants (C_per_mol).
* ALGEBRAIC[65] is v1a in component cAMP (fmol_per_sec).
* ALGEBRAIC[66] is v1b in component cAMP (fmol_per_sec).
* ALGEBRAIC[67] is v2a in component cAMP (fmol_per_sec).
* ALGEBRAIC[68] is v2b in component cAMP (fmol_per_sec).
* ALGEBRAIC[69] is v3a in component cAMP (fmol_per_sec).
* ALGEBRAIC[70] is v3b in component cAMP (fmol_per_sec).
* ALGEBRAIC[71] is v4a in component cAMP (fmol_per_sec).
* ALGEBRAIC[72] is v4b in component cAMP (fmol_per_sec).
* ALGEBRAIC[73] is v5 in component cAMP (fmol_per_sec).
* ALGEBRAIC[74] is v6 in component cAMP (fmol_per_sec).
* ALGEBRAIC[75] is v7 in component cAMP (fmol_per_sec).
* ALGEBRAIC[76] is vGiAC in component cAMP (fmol_per_sec).
* ALGEBRAIC[47] is mu_ATP in component cAMP (J_per_mol).
* ALGEBRAIC[49] is mu_AC in component cAMP (J_per_mol).
* ALGEBRAIC[48] is mu_cAMP in component cAMP (J_per_mol).
* ALGEBRAIC[50] is mu_AC_ATP in component cAMP (J_per_mol).
* ALGEBRAIC[61] is mu_FSK in component cAMP (J_per_mol).
* ALGEBRAIC[53] is mu_FSK_AC in component cAMP (J_per_mol).
* ALGEBRAIC[54] is mu_FSK_AC_ATP in component cAMP (J_per_mol).
* ALGEBRAIC[60] is mu_a_Gs_GTP in component cAMP (J_per_mol).
* ALGEBRAIC[51] is mu_a_Gs_GTP_AC in component cAMP (J_per_mol).
* ALGEBRAIC[52] is mu_a_Gs_GTP_AC_ATP in component cAMP (J_per_mol).
* ALGEBRAIC[55] is mu_PDE in component cAMP (J_per_mol).
* ALGEBRAIC[59] is mu_PDEinh in component cAMP (J_per_mol).
* ALGEBRAIC[56] is mu_PDE_cAMP in component cAMP (J_per_mol).
* ALGEBRAIC[58] is mu_IBMX in component cAMP (J_per_mol).
* ALGEBRAIC[57] is mu_five_AMP in component cAMP (J_per_mol).
* ALGEBRAIC[62] is mu_a_Gi_GTP in component cAMP (J_per_mol).
* ALGEBRAIC[63] is mu_ACinh in component cAMP (J_per_mol).
* ALGEBRAIC[64] is mu_PPi in component cAMP (J_per_mol).
* CONSTANTS[135] is vol in component cAMP (pL).
* ALGEBRAIC[77] is ATP_T in component cAMP (fmol).
* ALGEBRAIC[78] is AC_T in component cAMP (fmol).
* ALGEBRAIC[79] is Gs_T in component cAMP (fmol).
* ALGEBRAIC[80] is cAMP_T in component cAMP (fmol).
* ALGEBRAIC[81] is mu_L_B1 in component LRGbinding_B1AR (J_per_mol).
* ALGEBRAIC[82] is mu_R_B1 in component LRGbinding_B1AR (J_per_mol).
* ALGEBRAIC[83] is mu_Gs in component LRGbinding_B1AR (J_per_mol).
* ALGEBRAIC[84] is mu_LR_B1 in component LRGbinding_B1AR (J_per_mol).
* ALGEBRAIC[85] is mu_R_B1Gs in component LRGbinding_B1AR (J_per_mol).
* ALGEBRAIC[86] is mu_LR_B1Gs in component LRGbinding_B1AR (J_per_mol).
* ALGEBRAIC[87] is vsig1_B1 in component LRGbinding_B1AR (fmol_per_sec).
* ALGEBRAIC[88] is vsig2_B1 in component LRGbinding_B1AR (fmol_per_sec).
* ALGEBRAIC[89] is vsig3_B1 in component LRGbinding_B1AR (fmol_per_sec).
* ALGEBRAIC[90] is vsig4_B1 in component LRGbinding_B1AR (fmol_per_sec).
* CONSTANTS[136] is vol in component LRGbinding_B1AR (pL).
* ALGEBRAIC[91] is L_T in component LRGbinding_B1AR (fmol).
* ALGEBRAIC[92] is R_T in component LRGbinding_B1AR (fmol).
* ALGEBRAIC[93] is Gs_T in component LRGbinding_B1AR (fmol).
* ALGEBRAIC[94] is mu_L_M2 in component LRGbinding_M2 (J_per_mol).
* ALGEBRAIC[95] is mu_R_M2 in component LRGbinding_M2 (J_per_mol).
* ALGEBRAIC[96] is mu_Gi in component LRGbinding_M2 (J_per_mol).
* ALGEBRAIC[97] is mu_LR_M2 in component LRGbinding_M2 (J_per_mol).
* ALGEBRAIC[98] is mu_R_M2Gi in component LRGbinding_M2 (J_per_mol).
* ALGEBRAIC[99] is mu_LR_M2Gi in component LRGbinding_M2 (J_per_mol).
* ALGEBRAIC[100] is vsig1_M2 in component LRGbinding_M2 (fmol_per_sec).
* ALGEBRAIC[101] is vsig2_M2 in component LRGbinding_M2 (fmol_per_sec).
* ALGEBRAIC[102] is vsig3_M2 in component LRGbinding_M2 (fmol_per_sec).
* ALGEBRAIC[103] is vsig4_M2 in component LRGbinding_M2 (fmol_per_sec).
* CONSTANTS[137] is vol in component LRGbinding_M2 (pL).
* ALGEBRAIC[104] is L_T in component LRGbinding_M2 (fmol).
* ALGEBRAIC[105] is R_T in component LRGbinding_M2 (fmol).
* ALGEBRAIC[106] is Gi_T in component LRGbinding_M2 (fmol).
* ALGEBRAIC[122] is vDoff1_Gs in component GsProtein (fmol_per_sec).
* ALGEBRAIC[123] is vTon1_Gs in component GsProtein (fmol_per_sec).
* ALGEBRAIC[124] is vAct1_Gs in component GsProtein (fmol_per_sec).
* ALGEBRAIC[125] is vDoff2_Gs in component GsProtein (fmol_per_sec).
* ALGEBRAIC[126] is vTon2_Gs in component GsProtein (fmol_per_sec).
* ALGEBRAIC[127] is vAct2_Gs in component GsProtein (fmol_per_sec).
* ALGEBRAIC[128] is vhyd_Gs in component GsProtein (fmol_per_sec).
* ALGEBRAIC[129] is vreassoc_Gs in component GsProtein (fmol_per_sec).
* ALGEBRAIC[107] is mu_R_B1 in component GsProtein (J_per_mol).
* ALGEBRAIC[108] is mu_Gs in component GsProtein (J_per_mol).
* ALGEBRAIC[109] is mu_R_B1Gs in component GsProtein (J_per_mol).
* ALGEBRAIC[110] is mu_LR_B1 in component GsProtein (J_per_mol).
* ALGEBRAIC[111] is mu_LR_B1Gs in component GsProtein (J_per_mol).
* ALGEBRAIC[112] is mu_a_Gs_GTP in component GsProtein (J_per_mol).
* ALGEBRAIC[113] is mu_beta_gamma_Gs in component GsProtein (J_per_mol).
* ALGEBRAIC[114] is mu_a_Gs_GDP in component GsProtein (J_per_mol).
* ALGEBRAIC[115] is mu_Pi in component GsProtein (J_per_mol).
* ALGEBRAIC[116] is mu_R_B1_aby in component GsProtein (J_per_mol).
* ALGEBRAIC[117] is mu_R_B1_aby_T in component GsProtein (J_per_mol).
* ALGEBRAIC[118] is mu_LR_B1_aby in component GsProtein (J_per_mol).
* ALGEBRAIC[119] is mu_LR_B1_aby_T in component GsProtein (J_per_mol).
* ALGEBRAIC[120] is mu_GTP in component GsProtein (J_per_mol).
* ALGEBRAIC[121] is mu_GDP in component GsProtein (J_per_mol).
* ALGEBRAIC[130] is R_T in component GsProtein (fmol).
* ALGEBRAIC[131] is Gs_T in component GsProtein (fmol).
* ALGEBRAIC[147] is vDoff1_Gi in component GiProtein (fmol_per_sec).
* ALGEBRAIC[148] is vTon1_Gi in component GiProtein (fmol_per_sec).
* ALGEBRAIC[149] is vAct1_Gi in component GiProtein (fmol_per_sec).
* ALGEBRAIC[150] is vDoff2_Gi in component GiProtein (fmol_per_sec).
* ALGEBRAIC[151] is vTon2_Gi in component GiProtein (fmol_per_sec).
* ALGEBRAIC[152] is vAct2_Gi in component GiProtein (fmol_per_sec).
* ALGEBRAIC[153] is vhyd_Gi in component GiProtein (fmol_per_sec).
* ALGEBRAIC[154] is vreassoc_Gi in component GiProtein (fmol_per_sec).
* ALGEBRAIC[132] is mu_R_M2 in component GiProtein (J_per_mol).
* ALGEBRAIC[133] is mu_Gi in component GiProtein (J_per_mol).
* ALGEBRAIC[134] is mu_R_M2Gi in component GiProtein (J_per_mol).
* ALGEBRAIC[135] is mu_LR_M2 in component GiProtein (J_per_mol).
* ALGEBRAIC[136] is mu_LR_M2Gi in component GiProtein (J_per_mol).
* ALGEBRAIC[137] is mu_a_Gi_GTP in component GiProtein (J_per_mol).
* ALGEBRAIC[138] is mu_beta_gamma_Gi in component GiProtein (J_per_mol).
* ALGEBRAIC[139] is mu_a_Gi_GDP in component GiProtein (J_per_mol).
* ALGEBRAIC[140] is mu_Pi in component GiProtein (J_per_mol).
* ALGEBRAIC[141] is mu_R_M2_aby in component GiProtein (J_per_mol).
* ALGEBRAIC[142] is mu_R_M2_aby_T in component GiProtein (J_per_mol).
* ALGEBRAIC[143] is mu_LR_M2_aby in component GiProtein (J_per_mol).
* ALGEBRAIC[144] is mu_LR_M2_aby_T in component GiProtein (J_per_mol).
* ALGEBRAIC[145] is mu_GTP in component GiProtein (J_per_mol).
* ALGEBRAIC[146] is mu_GDP in component GiProtein (J_per_mol).
* ALGEBRAIC[155] is R_T in component GiProtein (fmol).
* ALGEBRAIC[156] is Gi_T in component GiProtein (fmol).
* RATES[0] is d/dt q_ATP in component cAMP (fmol).
* RATES[2] is d/dt q_AC in component cAMP (fmol).
* RATES[3] is d/dt q_AC_ATP in component cAMP (fmol).
* RATES[1] is d/dt q_cAMP in component cAMP (fmol).
* RATES[14] is d/dt q_FSK in component cAMP (fmol).
* RATES[6] is d/dt q_FSK_AC in component cAMP (fmol).
* RATES[7] is d/dt q_FSK_AC_ATP in component cAMP (fmol).
* RATES[13] is d/dt q_a_Gs_GTP in component cAMP (fmol).
* RATES[4] is d/dt q_a_Gs_GTP_AC in component cAMP (fmol).
* RATES[5] is d/dt q_a_Gs_GTP_AC_ATP in component cAMP (fmol).
* RATES[9] is d/dt q_PDE_cAMP in component cAMP (fmol).
* RATES[8] is d/dt q_PDE in component cAMP (fmol).
* RATES[11] is d/dt q_IBMX in component cAMP (fmol).
* RATES[12] is d/dt q_PDEinh in component cAMP (fmol).
* RATES[10] is d/dt q_five_AMP in component cAMP (fmol).
* RATES[15] is d/dt q_a_Gi_GTP in component cAMP (fmol).
* RATES[16] is d/dt q_ACinh in component cAMP (fmol).
* RATES[17] is d/dt q_PPi in component cAMP (fmol).
* RATES[18] is d/dt q_L_B1 in component LRGbinding_B1AR (fmol).
* RATES[19] is d/dt q_R_B1 in component LRGbinding_B1AR (fmol).
* RATES[20] is d/dt q_Gs in component LRGbinding_B1AR (fmol).
* RATES[21] is d/dt q_LR_B1 in component LRGbinding_B1AR (fmol).
* RATES[22] is d/dt q_R_B1Gs in component LRGbinding_B1AR (fmol).
* RATES[23] is d/dt q_LR_B1Gs in component LRGbinding_B1AR (fmol).
* RATES[24] is d/dt q_L_M2 in component LRGbinding_M2 (fmol).
* RATES[25] is d/dt q_R_M2 in component LRGbinding_M2 (fmol).
* RATES[26] is d/dt q_Gi in component LRGbinding_M2 (fmol).
* RATES[27] is d/dt q_LR_M2 in component LRGbinding_M2 (fmol).
* RATES[28] is d/dt q_R_M2Gi in component LRGbinding_M2 (fmol).
* RATES[29] is d/dt q_LR_M2Gi in component LRGbinding_M2 (fmol).
* RATES[30] is d/dt q_R_B1 in component GsProtein (fmol).
* RATES[32] is d/dt q_R_B1Gs in component GsProtein (fmol).
* RATES[31] is d/dt q_Gs in component GsProtein (fmol).
* RATES[33] is d/dt q_LR_B1 in component GsProtein (fmol).
* RATES[34] is d/dt q_LR_B1Gs in component GsProtein (fmol).
* RATES[35] is d/dt q_a_Gs_GTP in component GsProtein (fmol).
* RATES[36] is d/dt q_beta_gamma_Gs in component GsProtein (fmol).
* RATES[37] is d/dt q_a_Gs_GDP in component GsProtein (fmol).
* RATES[38] is d/dt q_Pi in component GsProtein (fmol).
* RATES[39] is d/dt q_R_B1_aby in component GsProtein (fmol).
* RATES[40] is d/dt q_R_B1_aby_T in component GsProtein (fmol).
* RATES[41] is d/dt q_LR_B1_aby in component GsProtein (fmol).
* RATES[42] is d/dt q_LR_B1_aby_T in component GsProtein (fmol).
* RATES[43] is d/dt q_GTP in component GsProtein (fmol).
* RATES[44] is d/dt q_GDP in component GsProtein (fmol).
* RATES[45] is d/dt q_R_M2 in component GiProtein (fmol).
* RATES[47] is d/dt q_R_M2Gi in component GiProtein (fmol).
* RATES[46] is d/dt q_Gi in component GiProtein (fmol).
* RATES[48] is d/dt q_LR_M2 in component GiProtein (fmol).
* RATES[49] is d/dt q_LR_M2Gi in component GiProtein (fmol).
* RATES[50] is d/dt q_a_Gi_GTP in component GiProtein (fmol).
* RATES[51] is d/dt q_beta_gamma_Gi in component GiProtein (fmol).
* RATES[52] is d/dt q_a_Gi_GDP in component GiProtein (fmol).
* RATES[53] is d/dt q_Pi in component GiProtein (fmol).
* RATES[54] is d/dt q_R_M2_aby in component GiProtein (fmol).
* RATES[55] is d/dt q_R_M2_aby_T in component GiProtein (fmol).
* RATES[56] is d/dt q_LR_M2_aby in component GiProtein (fmol).
* RATES[57] is d/dt q_LR_M2_aby_T in component GiProtein (fmol).
* RATES[58] is d/dt q_GTP in component GiProtein (fmol).
* RATES[59] is d/dt q_GDP in component GiProtein (fmol).
* There are a total of 12 condition variables.
*/
void
initConsts(double* CONSTANTS, double* RATES, double *STATES)
{
CONSTANTS[0] = 1.8998e+06;
CONSTANTS[1] = 0.00039666;
CONSTANTS[2] = 585859;
CONSTANTS[3] = 0.158089;
CONSTANTS[4] = 2.31489e+08;
CONSTANTS[5] = 2.69173e-17;
CONSTANTS[6] = 31490.2;
CONSTANTS[7] = 0.121117;
CONSTANTS[8] = 815.554;
CONSTANTS[9] = 8751.63;
CONSTANTS[10] = 345800;
CONSTANTS[11] = 2878.6;
CONSTANTS[12] = 2687.66;
CONSTANTS[13] = 63.5391;
CONSTANTS[14] = 14.2293;
CONSTANTS[15] = 1.13636e+07;
CONSTANTS[16] = 4069.25;
CONSTANTS[17] = 51.1466;
CONSTANTS[18] = 135.504;
CONSTANTS[19] = 1.78875e+06;
CONSTANTS[20] = 0.150576;
CONSTANTS[21] = 0.00442248;
CONSTANTS[22] = 163.698;
CONSTANTS[23] = 0.110539;
CONSTANTS[24] = 0.000540176;
CONSTANTS[25] = 3.32675;
CONSTANTS[26] = 0.86602;
CONSTANTS[27] = 3.20932e+06;
CONSTANTS[28] = 0.226169;
CONSTANTS[29] = 0.00556858;
CONSTANTS[30] = 26.9986;
CONSTANTS[31] = 0.209142;
CONSTANTS[32] = 0.00153751;
CONSTANTS[33] = 0.0445153;
CONSTANTS[34] = 7.99723;
CONSTANTS[35] = 1.84852e+06;
CONSTANTS[36] = 4.34247e-05;
CONSTANTS[37] = 0.0150226;
CONSTANTS[38] = 9.52621;
CONSTANTS[39] = 14.6573;
CONSTANTS[40] = 3.32164;
CONSTANTS[41] = 1.563;
CONSTANTS[42] = 0.0840652;
CONSTANTS[43] = 0.107997;
CONSTANTS[44] = 1.78633;
CONSTANTS[45] = 1.20007;
CONSTANTS[46] = 0.0150226;
CONSTANTS[47] = 0.0193351;
CONSTANTS[48] = 35.6442;
CONSTANTS[49] = 0.0253404;
CONSTANTS[50] = 5.83022e-06;
CONSTANTS[51] = 0.00308164;
CONSTANTS[52] = 10.0986;
CONSTANTS[53] = 8.40298e-05;
CONSTANTS[54] = 0.0361965;
CONSTANTS[55] = 0.00257076;
CONSTANTS[56] = 148.667;
CONSTANTS[57] = 0.310121;
CONSTANTS[58] = 293.639;
CONSTANTS[59] = 163.157;
CONSTANTS[60] = 0.079362;
CONSTANTS[61] = 0.00692368;
CONSTANTS[62] = 35.9123;
CONSTANTS[63] = 0.153698;
CONSTANTS[64] = 178.993;
CONSTANTS[65] = 105.77;
CONSTANTS[66] = 1.51453e-06;
CONSTANTS[67] = 0.198513;
CONSTANTS[68] = 5.20902e-06;
CONSTANTS[69] = 0.0798558;
CONSTANTS[70] = 0.00432163;
CONSTANTS[71] = 0.266681;
CONSTANTS[72] = 0.0867409;
CONSTANTS[73] = 3639.47;
CONSTANTS[74] = 4.62054e-05;
CONSTANTS[75] = 3.06222e-05;
CONSTANTS[76] = 0.0191425;
CONSTANTS[77] = 0.0580669;
CONSTANTS[78] = 0.0037486;
CONSTANTS[79] = 0.114918;
CONSTANTS[80] = 0.0248464;
CONSTANTS[81] = 34.4;
CONSTANTS[82] = 500;
CONSTANTS[83] = 3.5e-4;
CONSTANTS[84] = 3e-4;
CONSTANTS[85] = 0.25e-4;
CONSTANTS[86] = 1.8e-4;
CONSTANTS[87] = 1e1;
CONSTANTS[88] = 1e-5;
CONSTANTS[89] = 190;
CONSTANTS[90] = 1.889E-03;
CONSTANTS[91] = 1e-18;
CONSTANTS[92] = 1e-18;
CONSTANTS[93] = 3.800E-05;
CONSTANTS[94] = 1e-18;
CONSTANTS[95] = 1e-18;
CONSTANTS[96] = 1e-18;
CONSTANTS[97] = 1e-18;
CONSTANTS[98] = 1e-18;
CONSTANTS[99] = 1.482E-03;
CONSTANTS[100] = 1e-18;
CONSTANTS[101] = 1e-18;
CONSTANTS[102] = 3.80E-02;
CONSTANTS[103] = 1e-18;
CONSTANTS[104] = 1e-18;
CONSTANTS[105] = 1e-18;
CONSTANTS[106] = 1e-18;
CONSTANTS[107] = 0.0004579000;
CONSTANTS[108] = 0.1455400000;
CONSTANTS[109] = 1e-18;
CONSTANTS[110] = 1e-18;
CONSTANTS[111] = 1e-18;
CONSTANTS[112] = 0.00072;
CONSTANTS[113] = 0.00836;
CONSTANTS[114] = 1e-18;
CONSTANTS[115] = 1e-18;
CONSTANTS[116] = 1e-18;
CONSTANTS[117] = 1e-18;
CONSTANTS[118] = 1e-18;
CONSTANTS[119] = 570;
CONSTANTS[120] = 1e-18;
CONSTANTS[121] = 1e-18;
CONSTANTS[122] = 1e-18;
CONSTANTS[123] = 1e-18;
CONSTANTS[124] = 2.2;
CONSTANTS[125] = 1.1;
CONSTANTS[126] = 1e-18;
CONSTANTS[127] = 1e-18;
CONSTANTS[128] = 1e-18;
CONSTANTS[129] = 1e-18;
CONSTANTS[130] = 1e-18;
CONSTANTS[131] = 1e-18;
STATES[0] = 1e-16;
STATES[1] = 1e-16;
STATES[2] = 1e-16;
STATES[3] = 1e-16;
STATES[4] = 1e-16;
STATES[5] = 1e-16;
STATES[6] = 1e-16;
STATES[7] = 1e-16;
STATES[8] = 1e-16;
STATES[9] = 1e-16;
STATES[10] = 1e-16;
STATES[11] = 1e-16;
STATES[12] = 1e-16;
STATES[13] = 1e-16;
STATES[14] = 1e-16;
STATES[15] = 1e-16;
STATES[16] = 1e-16;
STATES[17] = 1e-16;
STATES[18] = 1e-18;
STATES[19] = 1e-18;
STATES[20] = 1e-18;
STATES[21] = 1e-18;
STATES[22] = 1e-18;
STATES[23] = 1e-18;
STATES[24] = 1e-16;
STATES[25] = 1e-16;
STATES[26] = 1e-16;
STATES[27] = 1e-16;
STATES[28] = 1e-16;
STATES[29] = 1e-16;
STATES[30] = 1e-16;
STATES[31] = 1e-16;
STATES[32] = 1e-16;
STATES[33] = 1e-16;
STATES[34] = 1e-16;
STATES[35] = 1e-16;
STATES[36] = 1e-16;
STATES[37] = 1e-16;
STATES[38] = 1e-16;
STATES[39] = 1e-16;
STATES[40] = 1e-16;
STATES[41] = 1e-16;
STATES[42] = 1e-16;
STATES[43] = 1e-16;
STATES[44] = 1e-16;
STATES[45] = 1e-16;
STATES[46] = 1e-16;
STATES[47] = 1e-16;
STATES[48] = 1e-16;
STATES[49] = 1e-16;
STATES[50] = 1e-16;
STATES[51] = 1e-16;
STATES[52] = 1e-16;
STATES[53] = 1e-16;
STATES[54] = 1e-16;
STATES[55] = 1e-16;
STATES[56] = 1e-16;
STATES[57] = 1e-16;
STATES[58] = 1e-16;
STATES[59] = 1e-16;
CONSTANTS[132] = 8.31;
CONSTANTS[133] = 310;
CONSTANTS[134] = 96485;
CONSTANTS[135] = 38.0;
CONSTANTS[136] = 34.4;
CONSTANTS[137] = 34.4;
CONSTANTS[138] = CONSTANTS[87]/CONSTANTS[86];
RATES[0] = 0.1001;
RATES[2] = 0.1001;
RATES[3] = 0.1001;
RATES[1] = 0.1001;
RATES[14] = 0.1001;
RATES[6] = 0.1001;
RATES[7] = 0.1001;
RATES[13] = 0.1001;
RATES[4] = 0.1001;
RATES[5] = 0.1001;
RATES[9] = 0.1001;
RATES[8] = 0.1001;
RATES[11] = 0.1001;
RATES[12] = 0.1001;
RATES[10] = 0.1001;
RATES[15] = 0.1001;
RATES[16] = 0.1001;
RATES[17] = 0.1001;
RATES[18] = 0.1001;
RATES[19] = 0.1001;
RATES[20] = 0.1001;
RATES[21] = 0.1001;
RATES[22] = 0.1001;
RATES[23] = 0.1001;
RATES[24] = 0.1001;
RATES[25] = 0.1001;
RATES[26] = 0.1001;
RATES[27] = 0.1001;
RATES[28] = 0.1001;
RATES[29] = 0.1001;
RATES[30] = 0.1001;
RATES[32] = 0.1001;
RATES[31] = 0.1001;
RATES[33] = 0.1001;
RATES[34] = 0.1001;
RATES[35] = 0.1001;
RATES[36] = 0.1001;
RATES[37] = 0.1001;
RATES[38] = 0.1001;
RATES[39] = 0.1001;
RATES[40] = 0.1001;
RATES[41] = 0.1001;
RATES[42] = 0.1001;
RATES[43] = 0.1001;
RATES[44] = 0.1001;
RATES[45] = 0.1001;
RATES[47] = 0.1001;
RATES[46] = 0.1001;
RATES[48] = 0.1001;
RATES[49] = 0.1001;
RATES[50] = 0.1001;
RATES[51] = 0.1001;
RATES[52] = 0.1001;
RATES[53] = 0.1001;
RATES[54] = 0.1001;
RATES[55] = 0.1001;
RATES[56] = 0.1001;
RATES[57] = 0.1001;
RATES[58] = 0.1001;
RATES[59] = 0.1001;
}
void
computeResiduals(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES,
double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS)
{
resid[0] = RATES[0] - (- ALGEBRAIC[65] - ALGEBRAIC[69]) - ALGEBRAIC[67];
resid[1] = RATES[2] - (((ALGEBRAIC[66] - ALGEBRAIC[65]) - ALGEBRAIC[74]) - ALGEBRAIC[75]) - ALGEBRAIC[76];
resid[2] = RATES[3] - ALGEBRAIC[65] - ALGEBRAIC[66];
resid[3] = RATES[1] - (ALGEBRAIC[66]+ALGEBRAIC[70]+ALGEBRAIC[68]) - ALGEBRAIC[71];
resid[4] = RATES[14] - - ALGEBRAIC[75];
resid[5] = RATES[6] - (ALGEBRAIC[75]+ALGEBRAIC[70]) - ALGEBRAIC[69];
resid[6] = RATES[7] - ALGEBRAIC[69] - ALGEBRAIC[70];
resid[7] = RATES[13] - - ALGEBRAIC[74];
resid[8] = RATES[4] - (ALGEBRAIC[74] - ALGEBRAIC[67])+ALGEBRAIC[68];
resid[9] = RATES[5] - ALGEBRAIC[67] - ALGEBRAIC[68];
resid[10] = RATES[9] - ALGEBRAIC[71] - ALGEBRAIC[72];
resid[11] = RATES[8] - (ALGEBRAIC[72] - ALGEBRAIC[71]) - ALGEBRAIC[73];
resid[12] = RATES[11] - - ALGEBRAIC[73];
resid[13] = RATES[12] - ALGEBRAIC[73];
resid[14] = RATES[10] - ALGEBRAIC[72];
resid[15] = RATES[15] - - ALGEBRAIC[76];
resid[16] = RATES[16] - ALGEBRAIC[76];
resid[17] = RATES[17] - ALGEBRAIC[76];
resid[18] = RATES[18] - - ALGEBRAIC[88] - ALGEBRAIC[90];
resid[19] = RATES[19] - - ALGEBRAIC[87] - ALGEBRAIC[90];
resid[20] = RATES[20] - - ALGEBRAIC[87] - ALGEBRAIC[89];
resid[21] = RATES[21] - - ALGEBRAIC[89]+ALGEBRAIC[90];
resid[22] = RATES[22] - ALGEBRAIC[87] - ALGEBRAIC[88];
resid[23] = RATES[23] - ALGEBRAIC[88]+ALGEBRAIC[89];
resid[24] = RATES[24] - - ALGEBRAIC[101] - ALGEBRAIC[103];
resid[25] = RATES[25] - - ALGEBRAIC[100] - ALGEBRAIC[103];
resid[26] = RATES[26] - - ALGEBRAIC[100] - ALGEBRAIC[102];
resid[27] = RATES[27] - - ALGEBRAIC[102]+ALGEBRAIC[103];
resid[28] = RATES[28] - ALGEBRAIC[100] - ALGEBRAIC[101];
resid[29] = RATES[29] - ALGEBRAIC[101]+ALGEBRAIC[102];
resid[30] = RATES[30] - ALGEBRAIC[124];
resid[31] = RATES[32] - - ALGEBRAIC[122];
resid[32] = RATES[31] - ALGEBRAIC[129];
resid[33] = RATES[33] - ALGEBRAIC[127];
resid[34] = RATES[34] - - ALGEBRAIC[125];
resid[35] = RATES[35] - (ALGEBRAIC[124]+ALGEBRAIC[127]) - ALGEBRAIC[128];
resid[36] = RATES[36] - (ALGEBRAIC[124]+ALGEBRAIC[127]) - ALGEBRAIC[129];
resid[37] = RATES[37] - ALGEBRAIC[128] - ALGEBRAIC[129];
resid[38] = RATES[38] - ALGEBRAIC[128];
resid[39] = RATES[39] - ALGEBRAIC[122] - ALGEBRAIC[123];
resid[40] = RATES[40] - ALGEBRAIC[123] - ALGEBRAIC[124];
resid[41] = RATES[41] - ALGEBRAIC[125] - ALGEBRAIC[126];
resid[42] = RATES[42] - ALGEBRAIC[126] - ALGEBRAIC[127];
resid[43] = RATES[43] - - ALGEBRAIC[123] - ALGEBRAIC[126];
resid[44] = RATES[44] - ALGEBRAIC[122]+ALGEBRAIC[125];
resid[45] = RATES[45] - ALGEBRAIC[149];
resid[46] = RATES[47] - - ALGEBRAIC[147];
resid[47] = RATES[46] - ALGEBRAIC[154];
resid[48] = RATES[48] - ALGEBRAIC[152];
resid[49] = RATES[49] - - ALGEBRAIC[150];
resid[50] = RATES[50] - (ALGEBRAIC[149]+ALGEBRAIC[152]) - ALGEBRAIC[153];
resid[51] = RATES[51] - (ALGEBRAIC[149]+ALGEBRAIC[152]) - ALGEBRAIC[154];
resid[52] = RATES[52] - ALGEBRAIC[153] - ALGEBRAIC[154];
resid[53] = RATES[53] - ALGEBRAIC[153];
resid[54] = RATES[54] - ALGEBRAIC[147] - ALGEBRAIC[148];
resid[55] = RATES[55] - ALGEBRAIC[148] - ALGEBRAIC[149];
resid[56] = RATES[56] - ALGEBRAIC[150] - ALGEBRAIC[151];
resid[57] = RATES[57] - ALGEBRAIC[151] - ALGEBRAIC[152];
resid[58] = RATES[58] - - ALGEBRAIC[148] - ALGEBRAIC[151];
resid[59] = RATES[59] - ALGEBRAIC[147]+ALGEBRAIC[150];
}
void
computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[77] = ALGEBRAIC[2]+ALGEBRAIC[5]+ALGEBRAIC[9]+ALGEBRAIC[7];
ALGEBRAIC[78] = ALGEBRAIC[4]+ALGEBRAIC[5]+ALGEBRAIC[8]+ALGEBRAIC[9]+ALGEBRAIC[6]+ALGEBRAIC[7];
ALGEBRAIC[79] = ALGEBRAIC[15]+ALGEBRAIC[6]+ALGEBRAIC[7];
ALGEBRAIC[80] = ALGEBRAIC[3]+ALGEBRAIC[11]+STATES[10];
ALGEBRAIC[91] = ALGEBRAIC[20]+ALGEBRAIC[23]+ALGEBRAIC[25];
ALGEBRAIC[92] = ALGEBRAIC[21]+ALGEBRAIC[23]+ALGEBRAIC[24]+ALGEBRAIC[25];
ALGEBRAIC[93] = ALGEBRAIC[22]+ALGEBRAIC[24]+ALGEBRAIC[25];
ALGEBRAIC[104] = ALGEBRAIC[26]+ALGEBRAIC[29]+ALGEBRAIC[31];
ALGEBRAIC[105] = ALGEBRAIC[27]+ALGEBRAIC[29]+ALGEBRAIC[30]+ALGEBRAIC[31];
ALGEBRAIC[106] = ALGEBRAIC[28]+ALGEBRAIC[30]+ALGEBRAIC[31];
ALGEBRAIC[130] = ALGEBRAIC[21]+ALGEBRAIC[24]+ALGEBRAIC[23]+ALGEBRAIC[25]+ALGEBRAIC[35]+ALGEBRAIC[36]+ALGEBRAIC[37]+ALGEBRAIC[38];
ALGEBRAIC[131] = ALGEBRAIC[22]+ALGEBRAIC[24]+ALGEBRAIC[25]+ALGEBRAIC[15]+ALGEBRAIC[33]+ALGEBRAIC[35]+ALGEBRAIC[36]+ALGEBRAIC[37]+ALGEBRAIC[38];
ALGEBRAIC[155] = ALGEBRAIC[27]+ALGEBRAIC[30]+ALGEBRAIC[29]+ALGEBRAIC[31]+ALGEBRAIC[43]+ALGEBRAIC[44]+ALGEBRAIC[45]+ALGEBRAIC[46];
ALGEBRAIC[156] = ALGEBRAIC[28]+ALGEBRAIC[30]+ALGEBRAIC[31]+ALGEBRAIC[17]+ALGEBRAIC[42]+ALGEBRAIC[43]+ALGEBRAIC[44]+ALGEBRAIC[45]+ALGEBRAIC[46];
ALGEBRAIC[157] = ALGEBRAIC[20]+ALGEBRAIC[25]+ALGEBRAIC[23]+ALGEBRAIC[37]+ALGEBRAIC[38];
ALGEBRAIC[158] = ALGEBRAIC[26]+ALGEBRAIC[31]+ALGEBRAIC[29]+ALGEBRAIC[45]+ALGEBRAIC[46];
ALGEBRAIC[159] = ALGEBRAIC[21]+ALGEBRAIC[24]+ALGEBRAIC[23]+ALGEBRAIC[25]+ALGEBRAIC[35]+ALGEBRAIC[36]+ALGEBRAIC[37]+ALGEBRAIC[38];
ALGEBRAIC[160] = ALGEBRAIC[27]+ALGEBRAIC[30]+ALGEBRAIC[29]+ALGEBRAIC[31]+ALGEBRAIC[43]+ALGEBRAIC[44]+ALGEBRAIC[45]+ALGEBRAIC[46];
ALGEBRAIC[161] = ALGEBRAIC[22]+ALGEBRAIC[24]+ALGEBRAIC[25]+ALGEBRAIC[15]+ALGEBRAIC[33]+ALGEBRAIC[35]+ALGEBRAIC[36]+ALGEBRAIC[37]+ALGEBRAIC[38];
ALGEBRAIC[162] = ALGEBRAIC[28]+ALGEBRAIC[30]+ALGEBRAIC[31]+ALGEBRAIC[17]+ALGEBRAIC[42]+ALGEBRAIC[43]+ALGEBRAIC[44]+ALGEBRAIC[45]+ALGEBRAIC[46];
ALGEBRAIC[163] = ALGEBRAIC[3]+ALGEBRAIC[11]+ALGEBRAIC[12]+ALGEBRAIC[2]+ALGEBRAIC[5]+ALGEBRAIC[7]+ALGEBRAIC[9];
}
void
computeEssentialVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[2] = CONSTANTS[89]+STATES[0];
ALGEBRAIC[47] = CONSTANTS[132]*CONSTANTS[133]*log( CONSTANTS[36]*ALGEBRAIC[2]);
ALGEBRAIC[4] = CONSTANTS[90]+STATES[2];
ALGEBRAIC[49] = CONSTANTS[132]*CONSTANTS[133]*log( CONSTANTS[38]*ALGEBRAIC[4]);
ALGEBRAIC[5] = CONSTANTS[92]+STATES[3];
ALGEBRAIC[50] = CONSTANTS[132]*CONSTANTS[133]*log( CONSTANTS[39]*ALGEBRAIC[5]);
ALGEBRAIC[65] = CONSTANTS[0]*(exp((ALGEBRAIC[49]+ALGEBRAIC[47])/( CONSTANTS[132]*CONSTANTS[133])) - exp(ALGEBRAIC[50]/( CONSTANTS[132]*CONSTANTS[133])));
ALGEBRAIC[3] = CONSTANTS[91]+STATES[1];
ALGEBRAIC[48] = CONSTANTS[132]*CONSTANTS[133]*log( CONSTANTS[37]*ALGEBRAIC[3]);
ALGEBRAIC[19] = CONSTANTS[106]+STATES[17];
ALGEBRAIC[64] = CONSTANTS[132]*CONSTANTS[133]*log( CONSTANTS[53]*ALGEBRAIC[19]);
ALGEBRAIC[66] = CONSTANTS[1]*(exp(ALGEBRAIC[50]/( CONSTANTS[132]*CONSTANTS[133])) - exp((ALGEBRAIC[49]+ALGEBRAIC[48]+ALGEBRAIC[64])/( CONSTANTS[132]*CONSTANTS[133])));
ALGEBRAIC[6] = CONSTANTS[97]+STATES[4];
ALGEBRAIC[51] = CONSTANTS[132]*CONSTANTS[133]*log( CONSTANTS[40]*ALGEBRAIC[6]);
ALGEBRAIC[7] = CONSTANTS[98]+STATES[5];
ALGEBRAIC[52] = CONSTANTS[132]*CONSTANTS[133]*log( CONSTANTS[41]*ALGEBRAIC[7]);
ALGEBRAIC[67] = CONSTANTS[2]*(exp((ALGEBRAIC[51]+ALGEBRAIC[47])/( CONSTANTS[132]*CONSTANTS[133])) - exp(ALGEBRAIC[52]/( CONSTANTS[132]*CONSTANTS[133])));
ALGEBRAIC[68] = CONSTANTS[3]*(exp(ALGEBRAIC[52]/( CONSTANTS[132]*CONSTANTS[133])) - exp((ALGEBRAIC[51]+ALGEBRAIC[48]+ALGEBRAIC[64])/( CONSTANTS[132]*CONSTANTS[133])));
ALGEBRAIC[8] = CONSTANTS[94]+STATES[6];
ALGEBRAIC[53] = CONSTANTS[132]*CONSTANTS[133]*log( CONSTANTS[42]*ALGEBRAIC[8]);
ALGEBRAIC[9] = CONSTANTS[95]+STATES[7];
ALGEBRAIC[54] = CONSTANTS[132]*CONSTANTS[133]*log( CONSTANTS[43]*ALGEBRAIC[9]);
ALGEBRAIC[69] = CONSTANTS[4]*(exp((ALGEBRAIC[53]+ALGEBRAIC[47])/( CONSTANTS[132]*CONSTANTS[133])) - exp(ALGEBRAIC[54]/( CONSTANTS[132]*CONSTANTS[133])));
ALGEBRAIC[70] = CONSTANTS[5]*(exp(ALGEBRAIC[54]/( CONSTANTS[132]*CONSTANTS[133])) - exp((ALGEBRAIC[53]+ALGEBRAIC[48]+ALGEBRAIC[64])/( CONSTANTS[132]*CONSTANTS[133])));
ALGEBRAIC[10] = CONSTANTS[99]+STATES[8];
ALGEBRAIC[55] = CONSTANTS[132]*CONSTANTS[133]*log( CONSTANTS[44]*ALGEBRAIC[10]);
ALGEBRAIC[11] = CONSTANTS[101]+STATES[9];
ALGEBRAIC[56] = CONSTANTS[132]*CONSTANTS[133]*log( CONSTANTS[45]*ALGEBRAIC[11]);
ALGEBRAIC[71] = CONSTANTS[6]*(exp((ALGEBRAIC[55]+ALGEBRAIC[48])/( CONSTANTS[132]*CONSTANTS[133])) - exp(ALGEBRAIC[56]/( CONSTANTS[132]*CONSTANTS[133])));
ALGEBRAIC[12] = CONSTANTS[103]+STATES[10];
ALGEBRAIC[57] = CONSTANTS[132]*CONSTANTS[133]*log( CONSTANTS[46]*ALGEBRAIC[12]);
ALGEBRAIC[72] = CONSTANTS[7]*(exp(ALGEBRAIC[56]/( CONSTANTS[132]*CONSTANTS[133])) - exp((ALGEBRAIC[55]+ALGEBRAIC[57])/( CONSTANTS[132]*CONSTANTS[133])));
ALGEBRAIC[14] = CONSTANTS[100]+STATES[12];
ALGEBRAIC[59] = CONSTANTS[132]*CONSTANTS[133]*log( CONSTANTS[48]*ALGEBRAIC[14]);
ALGEBRAIC[13] = CONSTANTS[102]+STATES[11];
ALGEBRAIC[58] = CONSTANTS[132]*CONSTANTS[133]*log( CONSTANTS[47]*ALGEBRAIC[13]);
ALGEBRAIC[73] = CONSTANTS[8]*(exp((ALGEBRAIC[55]+ALGEBRAIC[58])/( CONSTANTS[132]*CONSTANTS[133])) - exp(ALGEBRAIC[59]/( CONSTANTS[132]*CONSTANTS[133])));
ALGEBRAIC[15] = CONSTANTS[96]+STATES[13]+STATES[35];
ALGEBRAIC[60] = CONSTANTS[132]*CONSTANTS[133]*log( CONSTANTS[49]*ALGEBRAIC[15]);
ALGEBRAIC[74] = CONSTANTS[9]*(exp((ALGEBRAIC[49]+ALGEBRAIC[60])/( CONSTANTS[132]*CONSTANTS[133])) - exp(ALGEBRAIC[51]/( CONSTANTS[132]*CONSTANTS[133])));
ALGEBRAIC[16] = CONSTANTS[93]+STATES[14];
ALGEBRAIC[61] = CONSTANTS[132]*CONSTANTS[133]*log( CONSTANTS[50]*ALGEBRAIC[16]);
ALGEBRAIC[75] = CONSTANTS[10]*(exp((ALGEBRAIC[61]+ALGEBRAIC[49])/( CONSTANTS[132]*CONSTANTS[133])) - exp(ALGEBRAIC[53]/( CONSTANTS[132]*CONSTANTS[133])));
ALGEBRAIC[17] = CONSTANTS[104]+STATES[15]+STATES[50];
ALGEBRAIC[62] = CONSTANTS[132]*CONSTANTS[133]*log( CONSTANTS[51]*ALGEBRAIC[17]);
ALGEBRAIC[18] = CONSTANTS[105]+STATES[16];
ALGEBRAIC[63] = CONSTANTS[132]*CONSTANTS[133]*log( CONSTANTS[52]*ALGEBRAIC[18]);
ALGEBRAIC[76] = CONSTANTS[11]*(exp((ALGEBRAIC[49]+ALGEBRAIC[62])/( CONSTANTS[132]*CONSTANTS[133])) - exp(ALGEBRAIC[63]/( CONSTANTS[132]*CONSTANTS[133])));
ALGEBRAIC[21] = CONSTANTS[107]+STATES[19]+STATES[30];
ALGEBRAIC[82] = CONSTANTS[132]*CONSTANTS[133]*log( CONSTANTS[55]*ALGEBRAIC[21]);
ALGEBRAIC[22] = CONSTANTS[108]+STATES[20]+STATES[31];
ALGEBRAIC[83] = CONSTANTS[132]*CONSTANTS[133]*log( CONSTANTS[56]*ALGEBRAIC[22]);
ALGEBRAIC[24] = CONSTANTS[110]+STATES[22]+STATES[32];
ALGEBRAIC[85] = CONSTANTS[132]*CONSTANTS[133]*log( CONSTANTS[58]*ALGEBRAIC[24]);
ALGEBRAIC[87] = CONSTANTS[12]*exp((ALGEBRAIC[82]+ALGEBRAIC[83])/( CONSTANTS[132]*CONSTANTS[133])) - exp(ALGEBRAIC[85]/( CONSTANTS[132]*CONSTANTS[133]));
ALGEBRAIC[0] = (CONDVAR[0]<0.00000&&CONDVAR[1]>0.00000 ? CONSTANTS[88]+ CONSTANTS[138]*((VOI - CONSTANTS[83])+CONSTANTS[86]) : CONDVAR[2]>=0.00000&&CONDVAR[3]<0.00000 ? CONSTANTS[87]+CONSTANTS[88] : CONDVAR[4]<0.00000&&CONDVAR[5]>=0.00000 ? CONSTANTS[88]+ - CONSTANTS[138]*(((VOI - CONSTANTS[83]) - CONSTANTS[86]) - CONSTANTS[85]) : CONSTANTS[88]);
ALGEBRAIC[20] = ALGEBRAIC[0]+STATES[18];
ALGEBRAIC[81] = CONSTANTS[132]*CONSTANTS[133]*log( CONSTANTS[54]*ALGEBRAIC[20]);
ALGEBRAIC[25] = CONSTANTS[111]+STATES[23]+STATES[34];
ALGEBRAIC[86] = CONSTANTS[132]*CONSTANTS[133]*log( CONSTANTS[59]*ALGEBRAIC[25]);
ALGEBRAIC[88] = CONSTANTS[13]*exp((ALGEBRAIC[85]+ALGEBRAIC[81])/( CONSTANTS[132]*CONSTANTS[133])) - exp(ALGEBRAIC[86]/( CONSTANTS[132]*CONSTANTS[133]));
ALGEBRAIC[23] = CONSTANTS[109]+STATES[21]+STATES[33];
ALGEBRAIC[84] = CONSTANTS[132]*CONSTANTS[133]*log( CONSTANTS[57]*ALGEBRAIC[23]);
ALGEBRAIC[89] = CONSTANTS[14]*exp((ALGEBRAIC[84]+ALGEBRAIC[83])/( CONSTANTS[132]*CONSTANTS[133])) - exp(ALGEBRAIC[86]/( CONSTANTS[132]*CONSTANTS[133]));
ALGEBRAIC[90] = CONSTANTS[15]*exp((ALGEBRAIC[82]+ALGEBRAIC[81])/( CONSTANTS[132]*CONSTANTS[133])) - exp(ALGEBRAIC[84]/( CONSTANTS[132]*CONSTANTS[133]));
ALGEBRAIC[27] = CONSTANTS[112]+STATES[25]+STATES[45];
ALGEBRAIC[95] = CONSTANTS[132]*CONSTANTS[133]*log( CONSTANTS[61]*ALGEBRAIC[27]);
ALGEBRAIC[28] = CONSTANTS[113]+STATES[26]+STATES[46];
ALGEBRAIC[96] = CONSTANTS[132]*CONSTANTS[133]*log( CONSTANTS[62]*ALGEBRAIC[28]);
ALGEBRAIC[30] = CONSTANTS[115]+STATES[28]+STATES[47];
ALGEBRAIC[98] = CONSTANTS[132]*CONSTANTS[133]*log( CONSTANTS[64]*ALGEBRAIC[30]);
ALGEBRAIC[100] = CONSTANTS[16]*exp((ALGEBRAIC[95]+ALGEBRAIC[96])/( CONSTANTS[132]*CONSTANTS[133])) - exp(ALGEBRAIC[98]/( CONSTANTS[132]*CONSTANTS[133]));
ALGEBRAIC[1] = (CONDVAR[6]<0.00000&&CONDVAR[7]>0.00000 ? CONSTANTS[88]+ CONSTANTS[138]*((VOI - CONSTANTS[84])+CONSTANTS[86]) : CONDVAR[8]>=0.00000&&CONDVAR[9]<0.00000 ? CONSTANTS[87]+CONSTANTS[88] : CONDVAR[10]<0.00000&&CONDVAR[11]>=0.00000 ? CONSTANTS[88]+ - CONSTANTS[138]*(((VOI - CONSTANTS[84]) - CONSTANTS[86]) - CONSTANTS[85]) : CONSTANTS[88]);
ALGEBRAIC[26] = ALGEBRAIC[1]+STATES[24];
ALGEBRAIC[94] = CONSTANTS[132]*CONSTANTS[133]*log( CONSTANTS[60]*ALGEBRAIC[26]);
ALGEBRAIC[31] = CONSTANTS[116]+STATES[29]+STATES[49];
ALGEBRAIC[99] = CONSTANTS[132]*CONSTANTS[133]*log( CONSTANTS[65]*ALGEBRAIC[31]);
ALGEBRAIC[101] = CONSTANTS[17]*exp((ALGEBRAIC[98]+ALGEBRAIC[94])/( CONSTANTS[132]*CONSTANTS[133])) - exp(ALGEBRAIC[99]/( CONSTANTS[132]*CONSTANTS[133]));
ALGEBRAIC[29] = CONSTANTS[114]+STATES[27]+STATES[48];
ALGEBRAIC[97] = CONSTANTS[132]*CONSTANTS[133]*log( CONSTANTS[63]*ALGEBRAIC[29]);
ALGEBRAIC[102] = CONSTANTS[18]*exp((ALGEBRAIC[97]+ALGEBRAIC[96])/( CONSTANTS[132]*CONSTANTS[133])) - exp(ALGEBRAIC[99]/( CONSTANTS[132]*CONSTANTS[133]));
ALGEBRAIC[103] = CONSTANTS[19]*exp((ALGEBRAIC[95]+ALGEBRAIC[94])/( CONSTANTS[132]*CONSTANTS[133])) - exp(ALGEBRAIC[97]/( CONSTANTS[132]*CONSTANTS[133]));
ALGEBRAIC[109] = CONSTANTS[132]*CONSTANTS[133]*log( CONSTANTS[58]*ALGEBRAIC[24]);
ALGEBRAIC[35] = CONSTANTS[120]+STATES[39];
ALGEBRAIC[116] = CONSTANTS[132]*CONSTANTS[133]*log( CONSTANTS[69]*ALGEBRAIC[35]);
ALGEBRAIC[40] = CONSTANTS[125]+STATES[44]+STATES[59];
ALGEBRAIC[121] = CONSTANTS[132]*CONSTANTS[133]*log( CONSTANTS[74]*ALGEBRAIC[40]);
ALGEBRAIC[122] = CONSTANTS[20]*(exp(ALGEBRAIC[109]/( CONSTANTS[132]*CONSTANTS[133])) - exp((ALGEBRAIC[116]+ALGEBRAIC[121])/( CONSTANTS[132]*CONSTANTS[133])));
ALGEBRAIC[36] = CONSTANTS[121]+STATES[40];
ALGEBRAIC[117] = CONSTANTS[132]*CONSTANTS[133]*log( CONSTANTS[70]*ALGEBRAIC[36]);
ALGEBRAIC[39] = CONSTANTS[124]+STATES[43]+STATES[58];
ALGEBRAIC[120] = CONSTANTS[132]*CONSTANTS[133]*log( CONSTANTS[73]*ALGEBRAIC[39]);
ALGEBRAIC[123] = CONSTANTS[21]*(exp((ALGEBRAIC[116]+ALGEBRAIC[120])/( CONSTANTS[132]*CONSTANTS[133])) - exp(ALGEBRAIC[117]/( CONSTANTS[132]*CONSTANTS[133])));
ALGEBRAIC[107] = CONSTANTS[132]*CONSTANTS[133]*log( CONSTANTS[55]*ALGEBRAIC[21]);
ALGEBRAIC[112] = CONSTANTS[132]*CONSTANTS[133]*log( CONSTANTS[49]*ALGEBRAIC[15]);
ALGEBRAIC[32] = CONSTANTS[117]+STATES[36];
ALGEBRAIC[113] = CONSTANTS[132]*CONSTANTS[133]*log( CONSTANTS[66]*ALGEBRAIC[32]);
ALGEBRAIC[124] = CONSTANTS[22]*(exp(ALGEBRAIC[117]/( CONSTANTS[132]*CONSTANTS[133])) - exp((ALGEBRAIC[112]+ALGEBRAIC[113]+ALGEBRAIC[107])/( CONSTANTS[132]*CONSTANTS[133])));
ALGEBRAIC[111] = CONSTANTS[132]*CONSTANTS[133]*log( CONSTANTS[59]*ALGEBRAIC[25]);
ALGEBRAIC[37] = CONSTANTS[122]+STATES[41];
ALGEBRAIC[118] = CONSTANTS[132]*CONSTANTS[133]*log( CONSTANTS[71]*ALGEBRAIC[37]);
ALGEBRAIC[125] = CONSTANTS[23]*(exp(ALGEBRAIC[111]/( CONSTANTS[132]*CONSTANTS[133])) - exp((ALGEBRAIC[118]+ALGEBRAIC[121])/( CONSTANTS[132]*CONSTANTS[133])));
ALGEBRAIC[38] = CONSTANTS[123]+STATES[42];
ALGEBRAIC[119] = CONSTANTS[132]*CONSTANTS[133]*log( CONSTANTS[72]*ALGEBRAIC[38]);
ALGEBRAIC[126] = CONSTANTS[24]*(exp((ALGEBRAIC[118]+ALGEBRAIC[120])/( CONSTANTS[132]*CONSTANTS[133])) - exp(ALGEBRAIC[119]/( CONSTANTS[132]*CONSTANTS[133])));
ALGEBRAIC[110] = CONSTANTS[132]*CONSTANTS[133]*log( CONSTANTS[57]*ALGEBRAIC[23]);
ALGEBRAIC[127] = CONSTANTS[25]*(exp(ALGEBRAIC[119]/( CONSTANTS[132]*CONSTANTS[133])) - exp((ALGEBRAIC[112]+ALGEBRAIC[113]+ALGEBRAIC[110])/( CONSTANTS[132]*CONSTANTS[133])));
ALGEBRAIC[33] = CONSTANTS[118]+STATES[37];
ALGEBRAIC[114] = CONSTANTS[132]*CONSTANTS[133]*log( CONSTANTS[67]*ALGEBRAIC[33]);
ALGEBRAIC[34] = CONSTANTS[119]+STATES[38]+STATES[53];
ALGEBRAIC[115] = CONSTANTS[132]*CONSTANTS[133]*log( CONSTANTS[68]*ALGEBRAIC[34]);
ALGEBRAIC[128] = CONSTANTS[26]*(exp(ALGEBRAIC[112]/( CONSTANTS[132]*CONSTANTS[133])) - exp((ALGEBRAIC[114]+ALGEBRAIC[115])/( CONSTANTS[132]*CONSTANTS[133])));
ALGEBRAIC[108] = CONSTANTS[132]*CONSTANTS[133]*log( CONSTANTS[56]*ALGEBRAIC[22]);
ALGEBRAIC[129] = CONSTANTS[27]*(exp((ALGEBRAIC[114]+ALGEBRAIC[113])/( CONSTANTS[132]*CONSTANTS[133])) - exp(ALGEBRAIC[108]/( CONSTANTS[132]*CONSTANTS[133])));
ALGEBRAIC[134] = CONSTANTS[132]*CONSTANTS[133]*log( CONSTANTS[64]*ALGEBRAIC[30]);
ALGEBRAIC[43] = CONSTANTS[128]+STATES[54];
ALGEBRAIC[141] = CONSTANTS[132]*CONSTANTS[133]*log( CONSTANTS[77]*ALGEBRAIC[43]);
ALGEBRAIC[146] = CONSTANTS[132]*CONSTANTS[133]*log( CONSTANTS[74]*ALGEBRAIC[40]);
ALGEBRAIC[147] = CONSTANTS[28]*(exp(ALGEBRAIC[134]/( CONSTANTS[132]*CONSTANTS[133])) - exp((ALGEBRAIC[141]+ALGEBRAIC[146])/( CONSTANTS[132]*CONSTANTS[133])));
ALGEBRAIC[44] = CONSTANTS[129]+STATES[55];
ALGEBRAIC[142] = CONSTANTS[132]*CONSTANTS[133]*log( CONSTANTS[78]*ALGEBRAIC[44]);
ALGEBRAIC[145] = CONSTANTS[132]*CONSTANTS[133]*log( CONSTANTS[73]*ALGEBRAIC[39]);
ALGEBRAIC[148] = CONSTANTS[29]*(exp((ALGEBRAIC[141]+ALGEBRAIC[145])/( CONSTANTS[132]*CONSTANTS[133])) - exp(ALGEBRAIC[142]/( CONSTANTS[132]*CONSTANTS[133])));
ALGEBRAIC[132] = CONSTANTS[132]*CONSTANTS[133]*log( CONSTANTS[61]*ALGEBRAIC[27]);
ALGEBRAIC[137] = CONSTANTS[132]*CONSTANTS[133]*log( CONSTANTS[51]*ALGEBRAIC[17]);
ALGEBRAIC[41] = CONSTANTS[126]+STATES[51];
ALGEBRAIC[138] = CONSTANTS[132]*CONSTANTS[133]*log( CONSTANTS[75]*ALGEBRAIC[41]);
ALGEBRAIC[149] = CONSTANTS[30]*(exp(ALGEBRAIC[142]/( CONSTANTS[132]*CONSTANTS[133])) - exp((ALGEBRAIC[137]+ALGEBRAIC[138]+ALGEBRAIC[132])/( CONSTANTS[132]*CONSTANTS[133])));
ALGEBRAIC[136] = CONSTANTS[132]*CONSTANTS[133]*log( CONSTANTS[65]*ALGEBRAIC[31]);
ALGEBRAIC[45] = CONSTANTS[130]+STATES[56];
ALGEBRAIC[143] = CONSTANTS[132]*CONSTANTS[133]*log( CONSTANTS[79]*ALGEBRAIC[45]);
ALGEBRAIC[150] = CONSTANTS[31]*(exp(ALGEBRAIC[136]/( CONSTANTS[132]*CONSTANTS[133])) - exp((ALGEBRAIC[143]+ALGEBRAIC[146])/( CONSTANTS[132]*CONSTANTS[133])));
ALGEBRAIC[46] = CONSTANTS[131]+STATES[57];
ALGEBRAIC[144] = CONSTANTS[132]*CONSTANTS[133]*log( CONSTANTS[80]*ALGEBRAIC[46]);
ALGEBRAIC[151] = CONSTANTS[32]*(exp((ALGEBRAIC[143]+ALGEBRAIC[145])/( CONSTANTS[132]*CONSTANTS[133])) - exp(ALGEBRAIC[144]/( CONSTANTS[132]*CONSTANTS[133])));
ALGEBRAIC[135] = CONSTANTS[132]*CONSTANTS[133]*log( CONSTANTS[63]*ALGEBRAIC[29]);
ALGEBRAIC[152] = CONSTANTS[33]*(exp(ALGEBRAIC[144]/( CONSTANTS[132]*CONSTANTS[133])) - exp((ALGEBRAIC[137]+ALGEBRAIC[138]+ALGEBRAIC[135])/( CONSTANTS[132]*CONSTANTS[133])));
ALGEBRAIC[42] = CONSTANTS[127]+STATES[52];
ALGEBRAIC[139] = CONSTANTS[132]*CONSTANTS[133]*log( CONSTANTS[76]*ALGEBRAIC[42]);
ALGEBRAIC[140] = CONSTANTS[132]*CONSTANTS[133]*log( CONSTANTS[68]*ALGEBRAIC[34]);
ALGEBRAIC[153] = CONSTANTS[34]*(exp(ALGEBRAIC[137]/( CONSTANTS[132]*CONSTANTS[133])) - exp((ALGEBRAIC[139]+ALGEBRAIC[140])/( CONSTANTS[132]*CONSTANTS[133])));
ALGEBRAIC[133] = CONSTANTS[132]*CONSTANTS[133]*log( CONSTANTS[62]*ALGEBRAIC[28]);
ALGEBRAIC[154] = CONSTANTS[35]*(exp((ALGEBRAIC[139]+ALGEBRAIC[138])/( CONSTANTS[132]*CONSTANTS[133])) - exp(ALGEBRAIC[133]/( CONSTANTS[132]*CONSTANTS[133])));
}
void
getStateInformation(double* SI)
{
SI[0] = 1.0;
SI[1] = 1.0;
SI[2] = 1.0;
SI[3] = 1.0;
SI[4] = 1.0;
SI[5] = 1.0;
SI[6] = 1.0;
SI[7] = 1.0;
SI[8] = 1.0;
SI[9] = 1.0;
SI[10] = 1.0;
SI[11] = 1.0;
SI[12] = 1.0;
SI[13] = 1.0;
SI[14] = 1.0;
SI[15] = 1.0;
SI[16] = 1.0;
SI[17] = 1.0;
SI[18] = 1.0;
SI[19] = 1.0;
SI[20] = 1.0;
SI[21] = 1.0;
SI[22] = 1.0;
SI[23] = 1.0;
SI[24] = 1.0;
SI[25] = 1.0;
SI[26] = 1.0;
SI[27] = 1.0;
SI[28] = 1.0;
SI[29] = 1.0;
SI[30] = 1.0;
SI[31] = 1.0;
SI[32] = 1.0;
SI[33] = 1.0;
SI[34] = 1.0;
SI[35] = 1.0;
SI[36] = 1.0;
SI[37] = 1.0;
SI[38] = 1.0;
SI[39] = 1.0;
SI[40] = 1.0;
SI[41] = 1.0;
SI[42] = 1.0;
SI[43] = 1.0;
SI[44] = 1.0;
SI[45] = 1.0;
SI[46] = 1.0;
SI[47] = 1.0;
SI[48] = 1.0;
SI[49] = 1.0;
SI[50] = 1.0;
SI[51] = 1.0;
SI[52] = 1.0;
SI[53] = 1.0;
SI[54] = 1.0;
SI[55] = 1.0;
SI[56] = 1.0;
SI[57] = 1.0;
SI[58] = 1.0;
SI[59] = 1.0;
}
void
computeRoots(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES,
double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS)
{
CONDVAR[0] = VOI - CONSTANTS[83];
CONDVAR[1] = VOI - (CONSTANTS[83] - CONSTANTS[86]);
CONDVAR[2] = VOI - CONSTANTS[83];
CONDVAR[3] = VOI - (CONSTANTS[83]+CONSTANTS[85]);
CONDVAR[4] = VOI - (CONSTANTS[83]+CONSTANTS[86]+CONSTANTS[85]);
CONDVAR[5] = VOI - (CONSTANTS[83]+CONSTANTS[85]);
CONDVAR[6] = VOI - CONSTANTS[84];
CONDVAR[7] = VOI - (CONSTANTS[84] - CONSTANTS[86]);
CONDVAR[8] = VOI - CONSTANTS[84];
CONDVAR[9] = VOI - (CONSTANTS[84]+CONSTANTS[85]);
CONDVAR[10] = VOI - (CONSTANTS[84]+CONSTANTS[86]+CONSTANTS[85]);
CONDVAR[11] = VOI - (CONSTANTS[84]+CONSTANTS[85]);
}