/* 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]); }