/* There are a total of 152 entries in the algebraic variable array. There are a total of 56 entries in each of the rate and state variable arrays. There are a total of 150 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_Rswitch_B1 in component BG_parameters (fmol_per_sec). * CONSTANTS[13] is kappa_LRswitch_B1 in component BG_parameters (fmol_per_sec). * CONSTANTS[14] is kappa_C_B1 in component BG_parameters (fmol_per_sec). * CONSTANTS[15] is kappa_R_B1 in component BG_parameters (fmol_per_sec). * CONSTANTS[16] is kappa_L_B1 in component BG_parameters (fmol_per_sec). * CONSTANTS[17] is kappa_Act1_Gs in component BG_parameters (fmol_per_sec). * CONSTANTS[18] is kappa_Act2_Gs in component BG_parameters (fmol_per_sec). * CONSTANTS[19] is kappa_Hyd_Gs in component BG_parameters (fmol_per_sec). * CONSTANTS[20] is kappa_Reassoc_Gs in component BG_parameters (fmol_per_sec). * CONSTANTS[21] is kappa_InternR_B1 in component BG_parameters (fmol_per_sec). * CONSTANTS[22] is kappa_InternLR_B1 in component BG_parameters (fmol_per_sec). * CONSTANTS[23] is kappa_Rswitch_M2 in component BG_parameters (fmol_per_sec). * CONSTANTS[24] is kappa_LRswitch_M2 in component BG_parameters (fmol_per_sec). * CONSTANTS[25] is kappa_C_M2 in component BG_parameters (fmol_per_sec). * CONSTANTS[26] is kappa_R_M2 in component BG_parameters (fmol_per_sec). * CONSTANTS[27] is kappa_L_M2 in component BG_parameters (fmol_per_sec). * CONSTANTS[28] is kappa_Act1_Gi in component BG_parameters (fmol_per_sec). * CONSTANTS[29] is kappa_Act2_Gi in component BG_parameters (fmol_per_sec). * CONSTANTS[30] is kappa_Hyd_Gi in component BG_parameters (fmol_per_sec). * CONSTANTS[31] is kappa_Reassoc_Gi in component BG_parameters (fmol_per_sec). * CONSTANTS[32] is kappa_InternR_M2 in component BG_parameters (fmol_per_sec). * CONSTANTS[33] is kappa_InternLR_M2 in component BG_parameters (fmol_per_sec). * CONSTANTS[34] is K_ATP in component BG_parameters (per_fmol). * CONSTANTS[35] is K_cAMP in component BG_parameters (per_fmol). * CONSTANTS[36] is K_AC in component BG_parameters (per_fmol). * CONSTANTS[37] is K_AC_ATP in component BG_parameters (per_fmol). * CONSTANTS[38] is K_Gsa_GTP_AC in component BG_parameters (per_fmol). * CONSTANTS[39] is K_Gsa_GTP_AC_ATP in component BG_parameters (per_fmol). * CONSTANTS[40] is K_FSK_AC in component BG_parameters (per_fmol). * CONSTANTS[41] is K_FSK_AC_ATP in component BG_parameters (per_fmol). * CONSTANTS[42] is K_PDE in component BG_parameters (per_fmol). * CONSTANTS[43] is K_PDE_cAMP in component BG_parameters (per_fmol). * CONSTANTS[44] is K_five_AMP in component BG_parameters (per_fmol). * CONSTANTS[45] is K_IBMX in component BG_parameters (per_fmol). * CONSTANTS[46] is K_PDEinh in component BG_parameters (per_fmol). * CONSTANTS[47] is K_Gsa_GTP in component BG_parameters (per_fmol). * CONSTANTS[48] is K_FSK in component BG_parameters (per_fmol). * CONSTANTS[49] is K_Gia_GTP in component BG_parameters (per_fmol). * CONSTANTS[50] is K_ACinh in component BG_parameters (per_fmol). * CONSTANTS[51] is K_PPi in component BG_parameters (per_fmol). * CONSTANTS[52] is K_RB1_inactive in component BG_parameters (per_fmol). * CONSTANTS[53] is K_L_RB1_inactive in component BG_parameters (per_fmol). * CONSTANTS[54] is K_LB1 in component BG_parameters (per_fmol). * CONSTANTS[55] is K_RB1 in component BG_parameters (per_fmol). * CONSTANTS[56] is K_Gs in component BG_parameters (per_fmol). * CONSTANTS[57] is K_RB1_Gs in component BG_parameters (per_fmol). * CONSTANTS[58] is K_L_RB1 in component BG_parameters (per_fmol). * CONSTANTS[59] is K_L_RB1_Gs in component BG_parameters (per_fmol). * CONSTANTS[60] is K_Gsbetagamma in component BG_parameters (per_fmol). * CONSTANTS[61] is K_Gsa_GDP in component BG_parameters (per_fmol). * CONSTANTS[62] is K_GTP in component BG_parameters (per_fmol). * CONSTANTS[63] is K_GDP in component BG_parameters (per_fmol). * CONSTANTS[64] is K_Pi in component BG_parameters (per_fmol). * CONSTANTS[65] is K_RB1_tag in component BG_parameters (per_fmol). * CONSTANTS[66] is K_L_RB1_tag in component BG_parameters (per_fmol). * CONSTANTS[67] is K_RB1_GRKArr in component BG_parameters (per_fmol). * CONSTANTS[68] is K_L_RB1_GRKArr in component BG_parameters (per_fmol). * CONSTANTS[69] is K_GRKArr in component BG_parameters (per_fmol). * CONSTANTS[70] is K_RM2_inactive in component BG_parameters (per_fmol). * CONSTANTS[71] is K_L_RM2_inactive in component BG_parameters (per_fmol). * CONSTANTS[72] is K_LM2 in component BG_parameters (per_fmol). * CONSTANTS[73] is K_RM2 in component BG_parameters (per_fmol). * CONSTANTS[74] is K_Gi in component BG_parameters (per_fmol). * CONSTANTS[75] is K_RM2_Gi in component BG_parameters (per_fmol). * CONSTANTS[76] is K_L_RM2 in component BG_parameters (per_fmol). * CONSTANTS[77] is K_L_RM2_Gi in component BG_parameters (per_fmol). * CONSTANTS[78] is K_Gibetagamma in component BG_parameters (per_fmol). * CONSTANTS[79] is K_Gia_GDP in component BG_parameters (per_fmol). * CONSTANTS[80] is K_RM2_tag in component BG_parameters (per_fmol). * CONSTANTS[81] is K_L_RM2_tag in component BG_parameters (per_fmol). * CONSTANTS[82] is K_RM2_GRKArr in component BG_parameters (per_fmol). * CONSTANTS[83] is K_L_RM2_GRKArr in component BG_parameters (per_fmol). * VOI is time in component environment (second). * CONSTANTS[84] is vol_myo in component environment (pL). * CONSTANTS[85] is freq in component environment (dimensionless). * CONSTANTS[86] is stimSt1 in component environment (second). * CONSTANTS[87] is stimDur1 in component environment (second). * CONSTANTS[88] is tRamp1 in component environment (second). * CONSTANTS[89] is stimMag1 in component environment (fmol). * CONSTANTS[90] is stimHolding1 in component environment (fmol). * CONSTANTS[148] is m1 in component environment (fmol_per_sec). * CONSTANTS[91] is stimSt2 in component environment (second). * CONSTANTS[92] is stimDur2 in component environment (second). * CONSTANTS[93] is tRamp2 in component environment (second). * CONSTANTS[94] is stimMag2 in component environment (fmol). * CONSTANTS[95] is stimHolding2 in component environment (fmol). * CONSTANTS[149] is m2 in component environment (fmol_per_sec). * CONSTANTS[96] is q_ATP_init in component environment (fmol). * CONSTANTS[97] is q_AC_init in component environment (fmol). * CONSTANTS[98] is q_cAMP_init in component environment (fmol). * CONSTANTS[99] is q_AC_ATP_init in component environment (fmol). * CONSTANTS[100] is q_FSK_init in component environment (fmol). * CONSTANTS[101] is q_FSK_AC_init in component environment (fmol). * CONSTANTS[102] is q_FSK_AC_ATP_init in component environment (fmol). * CONSTANTS[103] is q_Gsa_GTP_init in component environment (fmol). * CONSTANTS[104] is q_Gsa_GTP_AC_init in component environment (fmol). * CONSTANTS[105] is q_Gsa_GTP_AC_ATP_init in component environment (fmol). * CONSTANTS[106] is q_PDE_init in component environment (fmol). * CONSTANTS[107] is q_PDEinh_init in component environment (fmol). * CONSTANTS[108] is q_PDE_cAMP_init in component environment (fmol). * CONSTANTS[109] is q_IBMX_init in component environment (fmol). * CONSTANTS[110] is q_five_AMP_init in component environment (fmol). * CONSTANTS[111] is q_Gia_GTP_init in component environment (fmol). * CONSTANTS[112] is q_ACinh_init in component environment (fmol). * CONSTANTS[113] is q_PPi_init in component environment (fmol). * CONSTANTS[114] is q_RB1_inactive_init in component environment (fmol). * CONSTANTS[115] is q_L_RB1_inactive_init in component environment (fmol). * ALGEBRAIC[0] is q_LB1_init in component environment (fmol). * CONSTANTS[116] is q_RB1_init in component environment (fmol). * CONSTANTS[117] is q_Gs_init in component environment (fmol). * CONSTANTS[118] is q_RB1_Gs_init in component environment (fmol). * CONSTANTS[119] is q_L_RB1_init in component environment (fmol). * CONSTANTS[120] is q_L_RB1_Gs_init in component environment (fmol). * CONSTANTS[121] is q_Gsbetagamma_init in component environment (fmol). * CONSTANTS[122] is q_Gsa_GDP_init in component environment (fmol). * CONSTANTS[123] is q_GTP_init in component environment (fmol). * CONSTANTS[124] is q_GDP_init in component environment (fmol). * CONSTANTS[125] is q_Pi_init in component environment (fmol). * CONSTANTS[126] is q_RB1_tag_init in component environment (fmol). * CONSTANTS[127] is q_L_RB1_tag_init in component environment (fmol). * CONSTANTS[128] is q_RB1_GRKArr_init in component environment (fmol). * CONSTANTS[129] is q_L_RB1_GRKArr_init in component environment (fmol). * CONSTANTS[130] is q_GRKArr_init in component environment (fmol). * CONSTANTS[131] is q_RM2_inactive_init in component environment (fmol). * CONSTANTS[132] is q_L_RM2_inactive_init in component environment (fmol). * ALGEBRAIC[1] is q_LM2_init in component environment (fmol). * CONSTANTS[133] is q_RM2_init in component environment (fmol). * CONSTANTS[134] is q_Gi_init in component environment (fmol). * CONSTANTS[135] is q_RM2_Gi_init in component environment (fmol). * CONSTANTS[136] is q_L_RM2_init in component environment (fmol). * CONSTANTS[137] is q_L_RM2_Gi_init in component environment (fmol). * CONSTANTS[138] is q_Gibetagamma_init in component environment (fmol). * CONSTANTS[139] is q_Gia_GDP_init in component environment (fmol). * CONSTANTS[140] is q_RM2_tag_init in component environment (fmol). * CONSTANTS[141] is q_L_RM2_tag_init in component environment (fmol). * CONSTANTS[142] is q_RM2_GRKArr_init in component environment (fmol). * CONSTANTS[143] is q_L_RM2_GRKArr_init in component environment (fmol). * ALGEBRAIC[148] is LB1_T in component environment (fmol). * ALGEBRAIC[149] is RB1_T in component environment (fmol). * ALGEBRAIC[150] is Gs_T in component environment (fmol). * ALGEBRAIC[151] 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_Gsa_GTP_AC in component environment (fmol). * ALGEBRAIC[7] is q_Gsa_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_Gsa_GTP in component environment (fmol). * ALGEBRAIC[16] is q_FSK in component environment (fmol). * ALGEBRAIC[17] is q_Gia_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_RB1_inactive in component environment (fmol). * ALGEBRAIC[21] is q_L_RB1_inactive in component environment (fmol). * ALGEBRAIC[22] is q_LB1 in component environment (fmol). * ALGEBRAIC[23] is q_RB1 in component environment (fmol). * ALGEBRAIC[24] is q_Gs in component environment (fmol). * ALGEBRAIC[25] is q_RB1_Gs in component environment (fmol). * ALGEBRAIC[26] is q_L_RB1 in component environment (fmol). * ALGEBRAIC[27] is q_L_RB1_Gs in component environment (fmol). * ALGEBRAIC[28] is q_Gsbetagamma in component environment (fmol). * ALGEBRAIC[29] is q_Gsa_GDP in component environment (fmol). * ALGEBRAIC[30] is q_GTP in component environment (fmol). * ALGEBRAIC[31] is q_GDP in component environment (fmol). * ALGEBRAIC[32] is q_Pi in component environment (fmol). * ALGEBRAIC[33] is q_RB1_tag in component environment (fmol). * ALGEBRAIC[34] is q_L_RB1_tag in component environment (fmol). * ALGEBRAIC[35] is q_RB1_GRKArr in component environment (fmol). * ALGEBRAIC[36] is q_L_RB1_GRKArr in component environment (fmol). * ALGEBRAIC[37] is q_GRKArr in component environment (fmol). * ALGEBRAIC[38] is q_RM2_inactive in component environment (fmol). * ALGEBRAIC[39] is q_L_RM2_inactive in component environment (fmol). * ALGEBRAIC[40] is q_LM2 in component environment (fmol). * ALGEBRAIC[41] is q_RM2 in component environment (fmol). * ALGEBRAIC[42] is q_Gi in component environment (fmol). * ALGEBRAIC[43] is q_RM2_Gi in component environment (fmol). * ALGEBRAIC[44] is q_L_RM2 in component environment (fmol). * ALGEBRAIC[45] is q_L_RM2_Gi in component environment (fmol). * ALGEBRAIC[46] is q_Gibetagamma in component environment (fmol). * ALGEBRAIC[47] is q_Gia_GDP in component environment (fmol). * ALGEBRAIC[48] is q_RM2_tag in component environment (fmol). * ALGEBRAIC[49] is q_L_RM2_tag in component environment (fmol). * ALGEBRAIC[50] is q_RM2_GRKArr in component environment (fmol). * ALGEBRAIC[51] is q_L_RM2_GRKArr 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_Gsa_GTP_AC in component cAMP (fmol). * STATES[5] is q_Gsa_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_Gsa_GTP in component cAMP (fmol). * STATES[14] is q_FSK in component cAMP (fmol). * STATES[15] is q_Gia_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_RB1_inactive in component GPCR_B1AR_reduced (fmol). * STATES[19] is q_L_RB1_inactive in component GPCR_B1AR_reduced (fmol). * STATES[20] is q_LB1 in component GPCR_B1AR_reduced (fmol). * STATES[21] is q_RB1 in component GPCR_B1AR_reduced (fmol). * STATES[22] is q_Gs in component GPCR_B1AR_reduced (fmol). * STATES[23] is q_RB1_Gs in component GPCR_B1AR_reduced (fmol). * STATES[24] is q_L_RB1 in component GPCR_B1AR_reduced (fmol). * STATES[25] is q_L_RB1_Gs in component GPCR_B1AR_reduced (fmol). * STATES[26] is q_Gsa_GTP in component GPCR_B1AR_reduced (fmol). * STATES[27] is q_Gsbetagamma in component GPCR_B1AR_reduced (fmol). * STATES[28] is q_Gsa_GDP in component GPCR_B1AR_reduced (fmol). * STATES[29] is q_GTP in component GPCR_B1AR_reduced (fmol). * STATES[30] is q_GDP in component GPCR_B1AR_reduced (fmol). * STATES[31] is q_Pi in component GPCR_B1AR_reduced (fmol). * STATES[32] is q_RB1_tag in component GPCR_B1AR_reduced (fmol). * STATES[33] is q_L_RB1_tag in component GPCR_B1AR_reduced (fmol). * STATES[34] is q_RB1_GRKArr in component GPCR_B1AR_reduced (fmol). * STATES[35] is q_L_RB1_GRKArr in component GPCR_B1AR_reduced (fmol). * STATES[36] is q_GRKArr in component GPCR_B1AR_reduced (fmol). * STATES[37] is q_RM2_inactive in component GPCR_M2_reduced (fmol). * STATES[38] is q_L_RM2_inactive in component GPCR_M2_reduced (fmol). * STATES[39] is q_LM2 in component GPCR_M2_reduced (fmol). * STATES[40] is q_RM2 in component GPCR_M2_reduced (fmol). * STATES[41] is q_Gi in component GPCR_M2_reduced (fmol). * STATES[42] is q_RM2_Gi in component GPCR_M2_reduced (fmol). * STATES[43] is q_L_RM2 in component GPCR_M2_reduced (fmol). * STATES[44] is q_L_RM2_Gi in component GPCR_M2_reduced (fmol). * STATES[45] is q_Gia_GTP in component GPCR_M2_reduced (fmol). * STATES[46] is q_Gibetagamma in component GPCR_M2_reduced (fmol). * STATES[47] is q_Gia_GDP in component GPCR_M2_reduced (fmol). * STATES[48] is q_GTP in component GPCR_M2_reduced (fmol). * STATES[49] is q_GDP in component GPCR_M2_reduced (fmol). * STATES[50] is q_Pi in component GPCR_M2_reduced (fmol). * STATES[51] is q_RM2_tag in component GPCR_M2_reduced (fmol). * STATES[52] is q_L_RM2_tag in component GPCR_M2_reduced (fmol). * STATES[53] is q_RM2_GRKArr in component GPCR_M2_reduced (fmol). * STATES[54] is q_L_RM2_GRKArr in component GPCR_M2_reduced (fmol). * STATES[55] is q_GRKArr in component GPCR_M2_reduced (fmol). * CONSTANTS[144] is R in component constants (J_per_K_per_mol). * CONSTANTS[145] is T in component constants (kelvin). * CONSTANTS[146] is F in component constants (C_per_mol). * ALGEBRAIC[70] is v1a in component cAMP (fmol_per_sec). * ALGEBRAIC[71] is v1b in component cAMP (fmol_per_sec). * ALGEBRAIC[72] is v2a in component cAMP (fmol_per_sec). * ALGEBRAIC[73] is v2b in component cAMP (fmol_per_sec). * ALGEBRAIC[74] is v3a in component cAMP (fmol_per_sec). * ALGEBRAIC[75] is v3b in component cAMP (fmol_per_sec). * ALGEBRAIC[76] is v4a in component cAMP (fmol_per_sec). * ALGEBRAIC[77] is v4b in component cAMP (fmol_per_sec). * ALGEBRAIC[78] is v5 in component cAMP (fmol_per_sec). * ALGEBRAIC[79] is v6 in component cAMP (fmol_per_sec). * ALGEBRAIC[80] is v7 in component cAMP (fmol_per_sec). * ALGEBRAIC[81] is vGiAC in component cAMP (fmol_per_sec). * ALGEBRAIC[52] is mu_ATP in component cAMP (J_per_mol). * ALGEBRAIC[54] is mu_AC in component cAMP (J_per_mol). * ALGEBRAIC[53] is mu_cAMP in component cAMP (J_per_mol). * ALGEBRAIC[55] is mu_AC_ATP in component cAMP (J_per_mol). * ALGEBRAIC[66] is mu_FSK in component cAMP (J_per_mol). * ALGEBRAIC[58] is mu_FSK_AC in component cAMP (J_per_mol). * ALGEBRAIC[59] is mu_FSK_AC_ATP in component cAMP (J_per_mol). * ALGEBRAIC[65] is mu_Gsa_GTP in component cAMP (J_per_mol). * ALGEBRAIC[56] is mu_Gsa_GTP_AC in component cAMP (J_per_mol). * ALGEBRAIC[57] is mu_Gsa_GTP_AC_ATP in component cAMP (J_per_mol). * ALGEBRAIC[60] is mu_PDE in component cAMP (J_per_mol). * ALGEBRAIC[64] is mu_PDEinh in component cAMP (J_per_mol). * ALGEBRAIC[61] is mu_PDE_cAMP in component cAMP (J_per_mol). * ALGEBRAIC[63] is mu_IBMX in component cAMP (J_per_mol). * ALGEBRAIC[62] is mu_five_AMP in component cAMP (J_per_mol). * ALGEBRAIC[67] is mu_Gia_GTP in component cAMP (J_per_mol). * ALGEBRAIC[68] is mu_ACinh in component cAMP (J_per_mol). * ALGEBRAIC[69] is mu_PPi in component cAMP (J_per_mol). * CONSTANTS[147] is vol in component cAMP (pL). * ALGEBRAIC[82] is ATP_T in component cAMP (fmol). * ALGEBRAIC[83] is AC_T in component cAMP (fmol). * ALGEBRAIC[84] is Gs_T in component cAMP (fmol). * ALGEBRAIC[85] is cAMP_T in component cAMP (fmol). * ALGEBRAIC[86] is mu_RB1_inactive in component GPCR_B1AR_reduced (J_per_mol). * ALGEBRAIC[87] is mu_L_RB1_inactive in component GPCR_B1AR_reduced (J_per_mol). * ALGEBRAIC[88] is mu_LB1 in component GPCR_B1AR_reduced (J_per_mol). * ALGEBRAIC[89] is mu_RB1 in component GPCR_B1AR_reduced (J_per_mol). * ALGEBRAIC[90] is mu_Gs in component GPCR_B1AR_reduced (J_per_mol). * ALGEBRAIC[91] is mu_RB1_Gs in component GPCR_B1AR_reduced (J_per_mol). * ALGEBRAIC[92] is mu_L_RB1 in component GPCR_B1AR_reduced (J_per_mol). * ALGEBRAIC[93] is mu_L_RB1_Gs in component GPCR_B1AR_reduced (J_per_mol). * ALGEBRAIC[94] is mu_Gsa_GTP in component GPCR_B1AR_reduced (J_per_mol). * ALGEBRAIC[95] is mu_Gsbetagamma in component GPCR_B1AR_reduced (J_per_mol). * ALGEBRAIC[96] is mu_Gsa_GDP in component GPCR_B1AR_reduced (J_per_mol). * ALGEBRAIC[97] is mu_GTP in component GPCR_B1AR_reduced (J_per_mol). * ALGEBRAIC[98] is mu_GDP in component GPCR_B1AR_reduced (J_per_mol). * ALGEBRAIC[99] is mu_Pi in component GPCR_B1AR_reduced (J_per_mol). * ALGEBRAIC[100] is mu_RB1_tag in component GPCR_B1AR_reduced (J_per_mol). * ALGEBRAIC[101] is mu_L_RB1_tag in component GPCR_B1AR_reduced (J_per_mol). * ALGEBRAIC[102] is mu_RB1_GRKArr in component GPCR_B1AR_reduced (J_per_mol). * ALGEBRAIC[103] is mu_L_RB1_GRKArr in component GPCR_B1AR_reduced (J_per_mol). * ALGEBRAIC[104] is mu_GRKArr in component GPCR_B1AR_reduced (J_per_mol). * ALGEBRAIC[105] is v_Rswitch_B1 in component GPCR_B1AR_reduced (fmol_per_sec). * ALGEBRAIC[106] is v_LRswitch_B1 in component GPCR_B1AR_reduced (fmol_per_sec). * ALGEBRAIC[107] is v_C_B1 in component GPCR_B1AR_reduced (fmol_per_sec). * ALGEBRAIC[108] is v_R_B1 in component GPCR_B1AR_reduced (fmol_per_sec). * ALGEBRAIC[109] is v_L_B1 in component GPCR_B1AR_reduced (fmol_per_sec). * ALGEBRAIC[110] is v_Act1_Gs in component GPCR_B1AR_reduced (fmol_per_sec). * ALGEBRAIC[111] is v_Act2_Gs in component GPCR_B1AR_reduced (fmol_per_sec). * ALGEBRAIC[112] is v_Hyd_Gs in component GPCR_B1AR_reduced (fmol_per_sec). * ALGEBRAIC[113] is v_Reassoc_Gs in component GPCR_B1AR_reduced (fmol_per_sec). * ALGEBRAIC[114] is v_InternR_B1 in component GPCR_B1AR_reduced (fmol_per_sec). * ALGEBRAIC[115] is v_InternLR_B1 in component GPCR_B1AR_reduced (fmol_per_sec). * ALGEBRAIC[116] is v_Rsynthesis in component GPCR_B1AR_reduced (fmol_per_sec). * ALGEBRAIC[117] is mu_RM2_inactive in component GPCR_M2_reduced (J_per_mol). * ALGEBRAIC[118] is mu_L_RM2_inactive in component GPCR_M2_reduced (J_per_mol). * ALGEBRAIC[119] is mu_LM2 in component GPCR_M2_reduced (J_per_mol). * ALGEBRAIC[120] is mu_RM2 in component GPCR_M2_reduced (J_per_mol). * ALGEBRAIC[121] is mu_Gi in component GPCR_M2_reduced (J_per_mol). * ALGEBRAIC[122] is mu_RM2_Gi in component GPCR_M2_reduced (J_per_mol). * ALGEBRAIC[123] is mu_L_RM2 in component GPCR_M2_reduced (J_per_mol). * ALGEBRAIC[124] is mu_L_RM2_Gi in component GPCR_M2_reduced (J_per_mol). * ALGEBRAIC[125] is mu_Gia_GTP in component GPCR_M2_reduced (J_per_mol). * ALGEBRAIC[126] is mu_Gibetagamma in component GPCR_M2_reduced (J_per_mol). * ALGEBRAIC[127] is mu_Gia_GDP in component GPCR_M2_reduced (J_per_mol). * ALGEBRAIC[128] is mu_GTP in component GPCR_M2_reduced (J_per_mol). * ALGEBRAIC[129] is mu_GDP in component GPCR_M2_reduced (J_per_mol). * ALGEBRAIC[130] is mu_Pi in component GPCR_M2_reduced (J_per_mol). * ALGEBRAIC[131] is mu_RM2_tag in component GPCR_M2_reduced (J_per_mol). * ALGEBRAIC[132] is mu_L_RM2_tag in component GPCR_M2_reduced (J_per_mol). * ALGEBRAIC[133] is mu_RM2_GRKArr in component GPCR_M2_reduced (J_per_mol). * ALGEBRAIC[134] is mu_L_RM2_GRKArr in component GPCR_M2_reduced (J_per_mol). * ALGEBRAIC[135] is mu_GRKArr in component GPCR_M2_reduced (J_per_mol). * ALGEBRAIC[136] is v_Rswitch_M2 in component GPCR_M2_reduced (fmol_per_sec). * ALGEBRAIC[137] is v_LRswitch_M2 in component GPCR_M2_reduced (fmol_per_sec). * ALGEBRAIC[138] is v_C_M2 in component GPCR_M2_reduced (fmol_per_sec). * ALGEBRAIC[139] is v_R_M2 in component GPCR_M2_reduced (fmol_per_sec). * ALGEBRAIC[140] is v_L_M2 in component GPCR_M2_reduced (fmol_per_sec). * ALGEBRAIC[141] is v_Act1_Gi in component GPCR_M2_reduced (fmol_per_sec). * ALGEBRAIC[142] is v_Act2_Gi in component GPCR_M2_reduced (fmol_per_sec). * ALGEBRAIC[143] is v_Hyd_Gi in component GPCR_M2_reduced (fmol_per_sec). * ALGEBRAIC[144] is v_Reassoc_Gi in component GPCR_M2_reduced (fmol_per_sec). * ALGEBRAIC[145] is v_InternR_M2 in component GPCR_M2_reduced (fmol_per_sec). * ALGEBRAIC[146] is v_InternLR_M2 in component GPCR_M2_reduced (fmol_per_sec). * ALGEBRAIC[147] is v_Rsynthesis in component GPCR_M2_reduced (fmol_per_sec). * 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_Gsa_GTP in component cAMP (fmol). * RATES[4] is d/dt q_Gsa_GTP_AC in component cAMP (fmol). * RATES[5] is d/dt q_Gsa_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_Gia_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_RB1_inactive in component GPCR_B1AR_reduced (fmol). * RATES[19] is d/dt q_L_RB1_inactive in component GPCR_B1AR_reduced (fmol). * RATES[20] is d/dt q_LB1 in component GPCR_B1AR_reduced (fmol). * RATES[21] is d/dt q_RB1 in component GPCR_B1AR_reduced (fmol). * RATES[22] is d/dt q_Gs in component GPCR_B1AR_reduced (fmol). * RATES[23] is d/dt q_RB1_Gs in component GPCR_B1AR_reduced (fmol). * RATES[24] is d/dt q_L_RB1 in component GPCR_B1AR_reduced (fmol). * RATES[25] is d/dt q_L_RB1_Gs in component GPCR_B1AR_reduced (fmol). * RATES[26] is d/dt q_Gsa_GTP in component GPCR_B1AR_reduced (fmol). * RATES[27] is d/dt q_Gsbetagamma in component GPCR_B1AR_reduced (fmol). * RATES[28] is d/dt q_Gsa_GDP in component GPCR_B1AR_reduced (fmol). * RATES[29] is d/dt q_GTP in component GPCR_B1AR_reduced (fmol). * RATES[30] is d/dt q_GDP in component GPCR_B1AR_reduced (fmol). * RATES[31] is d/dt q_Pi in component GPCR_B1AR_reduced (fmol). * RATES[32] is d/dt q_RB1_tag in component GPCR_B1AR_reduced (fmol). * RATES[33] is d/dt q_L_RB1_tag in component GPCR_B1AR_reduced (fmol). * RATES[34] is d/dt q_RB1_GRKArr in component GPCR_B1AR_reduced (fmol). * RATES[35] is d/dt q_L_RB1_GRKArr in component GPCR_B1AR_reduced (fmol). * RATES[36] is d/dt q_GRKArr in component GPCR_B1AR_reduced (fmol). * RATES[37] is d/dt q_RM2_inactive in component GPCR_M2_reduced (fmol). * RATES[38] is d/dt q_L_RM2_inactive in component GPCR_M2_reduced (fmol). * RATES[39] is d/dt q_LM2 in component GPCR_M2_reduced (fmol). * RATES[40] is d/dt q_RM2 in component GPCR_M2_reduced (fmol). * RATES[41] is d/dt q_Gi in component GPCR_M2_reduced (fmol). * RATES[42] is d/dt q_RM2_Gi in component GPCR_M2_reduced (fmol). * RATES[43] is d/dt q_L_RM2 in component GPCR_M2_reduced (fmol). * RATES[44] is d/dt q_L_RM2_Gi in component GPCR_M2_reduced (fmol). * RATES[45] is d/dt q_Gia_GTP in component GPCR_M2_reduced (fmol). * RATES[46] is d/dt q_Gibetagamma in component GPCR_M2_reduced (fmol). * RATES[47] is d/dt q_Gia_GDP in component GPCR_M2_reduced (fmol). * RATES[48] is d/dt q_GTP in component GPCR_M2_reduced (fmol). * RATES[49] is d/dt q_GDP in component GPCR_M2_reduced (fmol). * RATES[50] is d/dt q_Pi in component GPCR_M2_reduced (fmol). * RATES[51] is d/dt q_RM2_tag in component GPCR_M2_reduced (fmol). * RATES[52] is d/dt q_L_RM2_tag in component GPCR_M2_reduced (fmol). * RATES[53] is d/dt q_RM2_GRKArr in component GPCR_M2_reduced (fmol). * RATES[54] is d/dt q_L_RM2_GRKArr in component GPCR_M2_reduced (fmol). * RATES[55] is d/dt q_GRKArr in component GPCR_M2_reduced (fmol). * There are a total of 12 condition variables. */ void initConsts(double* CONSTANTS, double* RATES, double *STATES) { CONSTANTS[0] = 7.05076e+06; CONSTANTS[1] = 0.00147213; CONSTANTS[2] = 186465; CONSTANTS[3] = 0.0503158; CONSTANTS[4] = 2.50455e+08; CONSTANTS[5] = 2.91227e-17; CONSTANTS[6] = 30352.7; CONSTANTS[7] = 0.116741; CONSTANTS[8] = 833.759; CONSTANTS[9] = 4236.98; CONSTANTS[10] = 569102; CONSTANTS[11] = 1052.47; CONSTANTS[12] = 93.1496; CONSTANTS[13] = 0.00202063; CONSTANTS[14] = 1459.71; CONSTANTS[15] = 31664.5; CONSTANTS[16] = 0.00057588; CONSTANTS[17] = 0.00326844; CONSTANTS[18] = 37.7371; CONSTANTS[19] = 0.0787034; CONSTANTS[20] = 0.00156706; CONSTANTS[21] = 0.00130092; CONSTANTS[22] = 15.0203; CONSTANTS[23] = 271.765; CONSTANTS[24] = 0.0101296; CONSTANTS[25] = 7036.75; CONSTANTS[26] = 262283; CONSTANTS[27] = 0.111426; CONSTANTS[28] = 0.00270806; CONSTANTS[29] = 1.32138; CONSTANTS[30] = 0.48875; CONSTANTS[31] = 0.00258927; CONSTANTS[32] = 0.000652343; CONSTANTS[33] = 0.318306; CONSTANTS[34] = 6.60544e-05; CONSTANTS[35] = 0.0166528; CONSTANTS[36] = 1.68744; CONSTANTS[37] = 3.94935; CONSTANTS[38] = 6.86096; CONSTANTS[39] = 4.91084; CONSTANTS[40] = 0.0510801; CONSTANTS[41] = 0.0998184; CONSTANTS[42] = 1.67186; CONSTANTS[43] = 1.24505; CONSTANTS[44] = 0.0166528; CONSTANTS[45] = 0.020208; CONSTANTS[46] = 34.8659; CONSTANTS[47] = 0.295487; CONSTANTS[48] = 1.99992e-05; CONSTANTS[49] = 0.0475823; CONSTANTS[50] = 27.6205; CONSTANTS[51] = 0.000115307; CONSTANTS[52] = 3.12076; CONSTANTS[53] = 143865; CONSTANTS[54] = 4702.09; CONSTANTS[55] = 0.0312076; CONSTANTS[56] = 0.185505; CONSTANTS[57] = 6.57187; CONSTANTS[58] = 0.00143865; CONSTANTS[59] = 0.000569194; CONSTANTS[60] = 9.19376; CONSTANTS[61] = 70.9724; CONSTANTS[62] = 0.629469; CONSTANTS[63] = 5.92518e-07; CONSTANTS[64] = 1.51286e-06; CONSTANTS[65] = 1.35735; CONSTANTS[66] = 0.000117561; CONSTANTS[67] = 0.00223456; CONSTANTS[68] = 1.93537e-07; CONSTANTS[69] = 4785.64; CONSTANTS[70] = 1.06966; CONSTANTS[71] = 28697.8; CONSTANTS[72] = 70.9007; CONSTANTS[73] = 0.0106966; CONSTANTS[74] = 0.11227; CONSTANTS[75] = 1.23934; CONSTANTS[76] = 0.000286978; CONSTANTS[77] = 0.00253994; CONSTANTS[78] = 34.5537; CONSTANTS[79] = 11.4287; CONSTANTS[80] = 2.70687; CONSTANTS[81] = 0.00554752; CONSTANTS[82] = 0.00445621; CONSTANTS[83] = 9.13266e-06; CONSTANTS[84] = 34.4; CONSTANTS[85] = 500; CONSTANTS[86] = 3.5e-4; CONSTANTS[87] = 0.25e-4; CONSTANTS[88] = 1.8e-4; CONSTANTS[89] = 1e1; CONSTANTS[90] = 1e-5; CONSTANTS[91] = 7e-4; CONSTANTS[92] = 0.25e-4; CONSTANTS[93] = 1.8e-4; CONSTANTS[94] = 1e1; CONSTANTS[95] = 1e-5; CONSTANTS[96] = 190; CONSTANTS[97] = 1.889E-03; CONSTANTS[98] = 1e-18; CONSTANTS[99] = 1e-18; CONSTANTS[100] = 3.800E-05; CONSTANTS[101] = 1e-18; CONSTANTS[102] = 1e-18; CONSTANTS[103] = 1e-18; CONSTANTS[104] = 1e-18; CONSTANTS[105] = 1e-18; CONSTANTS[106] = 1.482E-03; CONSTANTS[107] = 1e-18; CONSTANTS[108] = 1e-18; CONSTANTS[109] = 3.80E-02; CONSTANTS[110] = 1e-18; CONSTANTS[111] = 1e-18; CONSTANTS[112] = 1e-18; CONSTANTS[113] = 1e-18; CONSTANTS[114] = 0.0004579000; CONSTANTS[115] = 1e-18; CONSTANTS[116] = 1e-18; CONSTANTS[117] = 0.1455400000; CONSTANTS[118] = 1e-18; CONSTANTS[119] = 1e-18; CONSTANTS[120] = 1e-18; CONSTANTS[121] = 1e-18; CONSTANTS[122] = 1e-18; CONSTANTS[123] = 2.2; CONSTANTS[124] = 1.1; CONSTANTS[125] = 570; CONSTANTS[126] = 1e-18; CONSTANTS[127] = 1e-18; CONSTANTS[128] = 1e-18; CONSTANTS[129] = 1e-18; CONSTANTS[130] = 1e-8; CONSTANTS[131] = 0.00072; CONSTANTS[132] = 1e-18; CONSTANTS[133] = 1e-18; CONSTANTS[134] = 0.00836; CONSTANTS[135] = 1e-18; CONSTANTS[136] = 1e-18; CONSTANTS[137] = 1e-18; CONSTANTS[138] = 1e-18; CONSTANTS[139] = 1e-18; CONSTANTS[140] = 1e-18; CONSTANTS[141] = 1e-18; CONSTANTS[142] = 1e-18; CONSTANTS[143] = 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-16; STATES[19] = 1e-16; STATES[20] = 1e-16; STATES[21] = 1e-16; STATES[22] = 1e-16; STATES[23] = 1e-16; 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; CONSTANTS[144] = 8.31; CONSTANTS[145] = 310; CONSTANTS[146] = 96485; CONSTANTS[147] = 38.0; CONSTANTS[148] = CONSTANTS[89]/CONSTANTS[88]; CONSTANTS[149] = CONSTANTS[94]/CONSTANTS[93]; 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[31] = 0.1001; RATES[32] = 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[46] = 0.1001; RATES[47] = 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; } void computeResiduals(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES, double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS) { resid[0] = RATES[0] - (- ALGEBRAIC[70] - ALGEBRAIC[74]) - ALGEBRAIC[72]; resid[1] = RATES[2] - (((ALGEBRAIC[71] - ALGEBRAIC[70]) - ALGEBRAIC[79]) - ALGEBRAIC[80]) - ALGEBRAIC[81]; resid[2] = RATES[3] - ALGEBRAIC[70] - ALGEBRAIC[71]; resid[3] = RATES[1] - (ALGEBRAIC[71]+ALGEBRAIC[75]+ALGEBRAIC[73]) - ALGEBRAIC[76]; resid[4] = RATES[14] - - ALGEBRAIC[80]; resid[5] = RATES[6] - (ALGEBRAIC[80]+ALGEBRAIC[75]) - ALGEBRAIC[74]; resid[6] = RATES[7] - ALGEBRAIC[74] - ALGEBRAIC[75]; resid[7] = RATES[13] - - ALGEBRAIC[79]; resid[8] = RATES[4] - (ALGEBRAIC[79] - ALGEBRAIC[72])+ALGEBRAIC[73]; resid[9] = RATES[5] - ALGEBRAIC[72] - ALGEBRAIC[73]; resid[10] = RATES[9] - ALGEBRAIC[76] - ALGEBRAIC[77]; resid[11] = RATES[8] - (ALGEBRAIC[77] - ALGEBRAIC[76]) - ALGEBRAIC[78]; resid[12] = RATES[11] - - ALGEBRAIC[78]; resid[13] = RATES[12] - ALGEBRAIC[78]; resid[14] = RATES[10] - ALGEBRAIC[77]; resid[15] = RATES[15] - - ALGEBRAIC[81]; resid[16] = RATES[16] - ALGEBRAIC[81]; resid[17] = RATES[17] - ALGEBRAIC[81]; resid[18] = RATES[18] - (- ALGEBRAIC[105] - ALGEBRAIC[109])+ 1.00000*ALGEBRAIC[116]; resid[19] = RATES[19] - - ALGEBRAIC[106]+ALGEBRAIC[109]; resid[20] = RATES[20] - - ALGEBRAIC[109]; resid[21] = RATES[21] - ALGEBRAIC[105] - ALGEBRAIC[107]; resid[22] = RATES[22] - (- ALGEBRAIC[107] - ALGEBRAIC[108])+ALGEBRAIC[113]; resid[23] = RATES[23] - ALGEBRAIC[107] - ALGEBRAIC[110]; resid[24] = RATES[24] - ALGEBRAIC[106] - ALGEBRAIC[108]; resid[25] = RATES[25] - ALGEBRAIC[108] - ALGEBRAIC[111]; resid[26] = RATES[26] - (ALGEBRAIC[110]+ALGEBRAIC[111]) - ALGEBRAIC[112]; resid[27] = RATES[27] - (ALGEBRAIC[110]+ALGEBRAIC[111]) - ALGEBRAIC[113]; resid[28] = RATES[28] - ALGEBRAIC[112] - ALGEBRAIC[113]; resid[29] = RATES[29] - - ALGEBRAIC[110] - ALGEBRAIC[111]; resid[30] = RATES[30] - ALGEBRAIC[110]+ALGEBRAIC[111]; resid[31] = RATES[31] - ALGEBRAIC[112]; resid[32] = RATES[32] - ALGEBRAIC[110] - ALGEBRAIC[114]; resid[33] = RATES[33] - ALGEBRAIC[111] - ALGEBRAIC[115]; resid[34] = RATES[34] - ALGEBRAIC[114]; resid[35] = RATES[35] - ALGEBRAIC[115]; resid[36] = RATES[36] - - ALGEBRAIC[114] - ALGEBRAIC[115]; resid[37] = RATES[37] - - ALGEBRAIC[136] - ALGEBRAIC[140]; resid[38] = RATES[38] - - ALGEBRAIC[137]; resid[39] = RATES[39] - - ALGEBRAIC[140]; resid[40] = RATES[40] - (ALGEBRAIC[136] - ALGEBRAIC[138])+ 1.00000*ALGEBRAIC[147]; resid[41] = RATES[41] - (- ALGEBRAIC[138] - ALGEBRAIC[139])+ALGEBRAIC[144]; resid[42] = RATES[42] - ALGEBRAIC[138] - ALGEBRAIC[141]; resid[43] = RATES[43] - (ALGEBRAIC[137] - ALGEBRAIC[139])+ALGEBRAIC[140]; resid[44] = RATES[44] - ALGEBRAIC[139] - ALGEBRAIC[142]; resid[45] = RATES[45] - (ALGEBRAIC[141]+ALGEBRAIC[142]) - ALGEBRAIC[143]; resid[46] = RATES[46] - (ALGEBRAIC[141]+ALGEBRAIC[142]) - ALGEBRAIC[144]; resid[47] = RATES[47] - ALGEBRAIC[143] - ALGEBRAIC[144]; resid[48] = RATES[48] - - ALGEBRAIC[141] - ALGEBRAIC[142]; resid[49] = RATES[49] - ALGEBRAIC[141]+ALGEBRAIC[142]; resid[50] = RATES[50] - ALGEBRAIC[143]; resid[51] = RATES[51] - ALGEBRAIC[141] - ALGEBRAIC[145]; resid[52] = RATES[52] - ALGEBRAIC[142] - ALGEBRAIC[146]; resid[53] = RATES[53] - ALGEBRAIC[145]; resid[54] = RATES[54] - ALGEBRAIC[146]; resid[55] = RATES[55] - - ALGEBRAIC[145] - ALGEBRAIC[146]; } void computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC) { ALGEBRAIC[82] = ALGEBRAIC[2]+ALGEBRAIC[5]+ALGEBRAIC[9]+ALGEBRAIC[7]; ALGEBRAIC[83] = ALGEBRAIC[4]+ALGEBRAIC[5]+ALGEBRAIC[8]+ALGEBRAIC[9]+ALGEBRAIC[6]+ALGEBRAIC[7]; ALGEBRAIC[84] = ALGEBRAIC[15]+ALGEBRAIC[6]+ALGEBRAIC[7]; ALGEBRAIC[85] = ALGEBRAIC[3]+ALGEBRAIC[11]+STATES[10]; ALGEBRAIC[148] = ALGEBRAIC[21]+ALGEBRAIC[22]+ALGEBRAIC[27]+ALGEBRAIC[26]+ALGEBRAIC[34]+ALGEBRAIC[36]; ALGEBRAIC[149] = ALGEBRAIC[20]+ALGEBRAIC[21]+ALGEBRAIC[23]+ALGEBRAIC[25]+ALGEBRAIC[26]+ALGEBRAIC[27]+ALGEBRAIC[33]+ALGEBRAIC[34]+ALGEBRAIC[35]+ALGEBRAIC[36]; ALGEBRAIC[150] = ALGEBRAIC[24]+ALGEBRAIC[25]+ALGEBRAIC[27]+ALGEBRAIC[15]+ALGEBRAIC[29]; ALGEBRAIC[151] = 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[96]+STATES[0]; ALGEBRAIC[52] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[34]*ALGEBRAIC[2]); ALGEBRAIC[4] = CONSTANTS[97]+STATES[2]; ALGEBRAIC[54] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[36]*ALGEBRAIC[4]); ALGEBRAIC[5] = CONSTANTS[99]+STATES[3]; ALGEBRAIC[55] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[37]*ALGEBRAIC[5]); ALGEBRAIC[70] = CONSTANTS[0]*(exp((ALGEBRAIC[54]+ALGEBRAIC[52])/( CONSTANTS[144]*CONSTANTS[145])) - exp(ALGEBRAIC[55]/( CONSTANTS[144]*CONSTANTS[145]))); ALGEBRAIC[3] = CONSTANTS[98]+STATES[1]; ALGEBRAIC[53] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[35]*ALGEBRAIC[3]); ALGEBRAIC[19] = CONSTANTS[113]+STATES[17]; ALGEBRAIC[69] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[51]*ALGEBRAIC[19]); ALGEBRAIC[71] = CONSTANTS[1]*(exp(ALGEBRAIC[55]/( CONSTANTS[144]*CONSTANTS[145])) - exp((ALGEBRAIC[54]+ALGEBRAIC[53]+ALGEBRAIC[69])/( CONSTANTS[144]*CONSTANTS[145]))); ALGEBRAIC[6] = CONSTANTS[104]+STATES[4]; ALGEBRAIC[56] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[38]*ALGEBRAIC[6]); ALGEBRAIC[7] = CONSTANTS[105]+STATES[5]; ALGEBRAIC[57] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[39]*ALGEBRAIC[7]); ALGEBRAIC[72] = CONSTANTS[2]*(exp((ALGEBRAIC[56]+ALGEBRAIC[52])/( CONSTANTS[144]*CONSTANTS[145])) - exp(ALGEBRAIC[57]/( CONSTANTS[144]*CONSTANTS[145]))); ALGEBRAIC[73] = CONSTANTS[3]*(exp(ALGEBRAIC[57]/( CONSTANTS[144]*CONSTANTS[145])) - exp((ALGEBRAIC[56]+ALGEBRAIC[53]+ALGEBRAIC[69])/( CONSTANTS[144]*CONSTANTS[145]))); ALGEBRAIC[8] = CONSTANTS[101]+STATES[6]; ALGEBRAIC[58] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[40]*ALGEBRAIC[8]); ALGEBRAIC[9] = CONSTANTS[102]+STATES[7]; ALGEBRAIC[59] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[41]*ALGEBRAIC[9]); ALGEBRAIC[74] = CONSTANTS[4]*(exp((ALGEBRAIC[58]+ALGEBRAIC[52])/( CONSTANTS[144]*CONSTANTS[145])) - exp(ALGEBRAIC[59]/( CONSTANTS[144]*CONSTANTS[145]))); ALGEBRAIC[75] = CONSTANTS[5]*(exp(ALGEBRAIC[59]/( CONSTANTS[144]*CONSTANTS[145])) - exp((ALGEBRAIC[58]+ALGEBRAIC[53]+ALGEBRAIC[69])/( CONSTANTS[144]*CONSTANTS[145]))); ALGEBRAIC[10] = CONSTANTS[106]+STATES[8]; ALGEBRAIC[60] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[42]*ALGEBRAIC[10]); ALGEBRAIC[11] = CONSTANTS[108]+STATES[9]; ALGEBRAIC[61] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[43]*ALGEBRAIC[11]); ALGEBRAIC[76] = CONSTANTS[6]*(exp((ALGEBRAIC[60]+ALGEBRAIC[53])/( CONSTANTS[144]*CONSTANTS[145])) - exp(ALGEBRAIC[61]/( CONSTANTS[144]*CONSTANTS[145]))); ALGEBRAIC[12] = CONSTANTS[110]+STATES[10]; ALGEBRAIC[62] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[44]*ALGEBRAIC[12]); ALGEBRAIC[77] = CONSTANTS[7]*(exp(ALGEBRAIC[61]/( CONSTANTS[144]*CONSTANTS[145])) - exp((ALGEBRAIC[60]+ALGEBRAIC[62])/( CONSTANTS[144]*CONSTANTS[145]))); ALGEBRAIC[14] = CONSTANTS[107]+STATES[12]; ALGEBRAIC[64] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[46]*ALGEBRAIC[14]); ALGEBRAIC[13] = CONSTANTS[109]+STATES[11]; ALGEBRAIC[63] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[45]*ALGEBRAIC[13]); ALGEBRAIC[78] = CONSTANTS[8]*(exp((ALGEBRAIC[60]+ALGEBRAIC[63])/( CONSTANTS[144]*CONSTANTS[145])) - exp(ALGEBRAIC[64]/( CONSTANTS[144]*CONSTANTS[145]))); ALGEBRAIC[15] = CONSTANTS[103]+STATES[13]+STATES[26]; ALGEBRAIC[65] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[47]*ALGEBRAIC[15]); ALGEBRAIC[79] = CONSTANTS[9]*(exp((ALGEBRAIC[54]+ALGEBRAIC[65])/( CONSTANTS[144]*CONSTANTS[145])) - exp(ALGEBRAIC[56]/( CONSTANTS[144]*CONSTANTS[145]))); ALGEBRAIC[16] = CONSTANTS[100]+STATES[14]; ALGEBRAIC[66] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[48]*ALGEBRAIC[16]); ALGEBRAIC[80] = CONSTANTS[10]*(exp((ALGEBRAIC[66]+ALGEBRAIC[54])/( CONSTANTS[144]*CONSTANTS[145])) - exp(ALGEBRAIC[58]/( CONSTANTS[144]*CONSTANTS[145]))); ALGEBRAIC[17] = CONSTANTS[111]+STATES[15]+STATES[45]; ALGEBRAIC[67] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[49]*ALGEBRAIC[17]); ALGEBRAIC[18] = CONSTANTS[112]+STATES[16]; ALGEBRAIC[68] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[50]*ALGEBRAIC[18]); ALGEBRAIC[81] = CONSTANTS[11]*(exp((ALGEBRAIC[54]+ALGEBRAIC[67])/( CONSTANTS[144]*CONSTANTS[145])) - exp(ALGEBRAIC[68]/( CONSTANTS[144]*CONSTANTS[145]))); ALGEBRAIC[20] = CONSTANTS[114]+STATES[18]; ALGEBRAIC[86] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[52]*ALGEBRAIC[20]); ALGEBRAIC[23] = CONSTANTS[116]+STATES[21]; ALGEBRAIC[89] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[55]*ALGEBRAIC[23]); ALGEBRAIC[105] = CONSTANTS[12]*(exp(ALGEBRAIC[86]/( CONSTANTS[144]*CONSTANTS[145])) - exp(ALGEBRAIC[89]/( CONSTANTS[144]*CONSTANTS[145]))); ALGEBRAIC[21] = CONSTANTS[115]+STATES[19]; ALGEBRAIC[87] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[53]*ALGEBRAIC[21]); ALGEBRAIC[26] = CONSTANTS[119]+STATES[24]; ALGEBRAIC[92] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[58]*ALGEBRAIC[26]); ALGEBRAIC[106] = CONSTANTS[13]*(exp(ALGEBRAIC[87]/( CONSTANTS[144]*CONSTANTS[145])) - exp(ALGEBRAIC[92]/( CONSTANTS[144]*CONSTANTS[145]))); ALGEBRAIC[24] = CONSTANTS[117]+STATES[22]; ALGEBRAIC[90] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[56]*ALGEBRAIC[24]); ALGEBRAIC[25] = CONSTANTS[118]+STATES[23]; ALGEBRAIC[91] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[57]*ALGEBRAIC[25]); ALGEBRAIC[107] = CONSTANTS[14]*(exp((ALGEBRAIC[89]+ALGEBRAIC[90])/( CONSTANTS[144]*CONSTANTS[145])) - exp(ALGEBRAIC[91]/( CONSTANTS[144]*CONSTANTS[145]))); ALGEBRAIC[27] = CONSTANTS[120]+STATES[25]; ALGEBRAIC[93] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[59]*ALGEBRAIC[27]); ALGEBRAIC[108] = CONSTANTS[15]*(exp((ALGEBRAIC[92]+ALGEBRAIC[90])/( CONSTANTS[144]*CONSTANTS[145])) - exp(ALGEBRAIC[93]/( CONSTANTS[144]*CONSTANTS[145]))); ALGEBRAIC[0] = (CONDVAR[0]<0.00000&&CONDVAR[1]>0.00000 ? CONSTANTS[90]+ CONSTANTS[148]*((VOI - CONSTANTS[86])+CONSTANTS[88]) : CONDVAR[2]>=0.00000&&CONDVAR[3]<0.00000 ? CONSTANTS[89]+CONSTANTS[90] : CONDVAR[4]<0.00000&&CONDVAR[5]>=0.00000 ? CONSTANTS[90]+ - CONSTANTS[148]*(((VOI - CONSTANTS[86]) - CONSTANTS[88]) - CONSTANTS[87]) : CONSTANTS[90]); ALGEBRAIC[22] = ALGEBRAIC[0]+STATES[20]; ALGEBRAIC[88] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[54]*ALGEBRAIC[22]); ALGEBRAIC[109] = CONSTANTS[16]*(exp((ALGEBRAIC[86]+ALGEBRAIC[88])/( CONSTANTS[144]*CONSTANTS[145])) - exp(ALGEBRAIC[87]/( CONSTANTS[144]*CONSTANTS[145]))); ALGEBRAIC[94] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[47]*ALGEBRAIC[15]); ALGEBRAIC[28] = CONSTANTS[121]+STATES[27]; ALGEBRAIC[95] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[60]*ALGEBRAIC[28]); ALGEBRAIC[30] = CONSTANTS[123]+STATES[29]+STATES[48]; ALGEBRAIC[97] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[62]*ALGEBRAIC[30]); ALGEBRAIC[31] = CONSTANTS[124]+STATES[30]+STATES[49]; ALGEBRAIC[98] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[63]*ALGEBRAIC[31]); ALGEBRAIC[33] = CONSTANTS[126]+STATES[32]; ALGEBRAIC[100] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[65]*ALGEBRAIC[33]); ALGEBRAIC[110] = CONSTANTS[17]*(exp((ALGEBRAIC[91]+ALGEBRAIC[97])/( CONSTANTS[144]*CONSTANTS[145])) - exp((ALGEBRAIC[94]+ALGEBRAIC[95]+ALGEBRAIC[100]+ALGEBRAIC[98])/( CONSTANTS[144]*CONSTANTS[145]))); ALGEBRAIC[34] = CONSTANTS[127]+STATES[33]; ALGEBRAIC[101] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[66]*ALGEBRAIC[34]); ALGEBRAIC[111] = CONSTANTS[18]*(exp((ALGEBRAIC[93]+ALGEBRAIC[97])/( CONSTANTS[144]*CONSTANTS[145])) - exp((ALGEBRAIC[94]+ALGEBRAIC[95]+ALGEBRAIC[101]+ALGEBRAIC[98])/( CONSTANTS[144]*CONSTANTS[145]))); ALGEBRAIC[29] = CONSTANTS[122]+STATES[28]; ALGEBRAIC[96] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[61]*ALGEBRAIC[29]); ALGEBRAIC[32] = CONSTANTS[125]+STATES[31]+STATES[50]; ALGEBRAIC[99] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[64]*ALGEBRAIC[32]); ALGEBRAIC[112] = CONSTANTS[19]*(exp(ALGEBRAIC[94]/( CONSTANTS[144]*CONSTANTS[145])) - exp((ALGEBRAIC[96]+ALGEBRAIC[99])/( CONSTANTS[144]*CONSTANTS[145]))); ALGEBRAIC[113] = CONSTANTS[20]*(exp((ALGEBRAIC[96]+ALGEBRAIC[95])/( CONSTANTS[144]*CONSTANTS[145])) - exp(ALGEBRAIC[90]/( CONSTANTS[144]*CONSTANTS[145]))); ALGEBRAIC[35] = CONSTANTS[128]+STATES[34]; ALGEBRAIC[102] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[67]*ALGEBRAIC[35]); ALGEBRAIC[37] = CONSTANTS[130]+STATES[36]+STATES[55]; ALGEBRAIC[104] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[69]*ALGEBRAIC[37]); ALGEBRAIC[114] = CONSTANTS[21]*(exp((ALGEBRAIC[100]+ALGEBRAIC[104])/( CONSTANTS[144]*CONSTANTS[145])) - exp(ALGEBRAIC[102]/( CONSTANTS[144]*CONSTANTS[145]))); ALGEBRAIC[36] = CONSTANTS[129]+STATES[35]; ALGEBRAIC[103] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[68]*ALGEBRAIC[36]); ALGEBRAIC[115] = CONSTANTS[22]*(exp((ALGEBRAIC[101]+ALGEBRAIC[104])/( CONSTANTS[144]*CONSTANTS[145])) - exp(ALGEBRAIC[103]/( CONSTANTS[144]*CONSTANTS[145]))); ALGEBRAIC[116] = ALGEBRAIC[105]+ALGEBRAIC[109]; ALGEBRAIC[38] = CONSTANTS[131]+STATES[37]; ALGEBRAIC[117] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[70]*ALGEBRAIC[38]); ALGEBRAIC[41] = CONSTANTS[133]+STATES[40]; ALGEBRAIC[120] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[73]*ALGEBRAIC[41]); ALGEBRAIC[136] = CONSTANTS[23]*(exp(ALGEBRAIC[117]/( CONSTANTS[144]*CONSTANTS[145])) - exp(ALGEBRAIC[120]/( CONSTANTS[144]*CONSTANTS[145]))); ALGEBRAIC[39] = CONSTANTS[132]+STATES[38]; ALGEBRAIC[118] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[71]*ALGEBRAIC[39]); ALGEBRAIC[44] = CONSTANTS[136]+STATES[43]; ALGEBRAIC[123] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[76]*ALGEBRAIC[44]); ALGEBRAIC[137] = CONSTANTS[24]*(exp(ALGEBRAIC[118]/( CONSTANTS[144]*CONSTANTS[145])) - exp(ALGEBRAIC[123]/( CONSTANTS[144]*CONSTANTS[145]))); ALGEBRAIC[42] = CONSTANTS[134]+STATES[41]; ALGEBRAIC[121] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[74]*ALGEBRAIC[42]); ALGEBRAIC[43] = CONSTANTS[135]+STATES[42]; ALGEBRAIC[122] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[75]*ALGEBRAIC[43]); ALGEBRAIC[138] = CONSTANTS[25]*(exp((ALGEBRAIC[120]+ALGEBRAIC[121])/( CONSTANTS[144]*CONSTANTS[145])) - exp(ALGEBRAIC[122]/( CONSTANTS[144]*CONSTANTS[145]))); ALGEBRAIC[45] = CONSTANTS[137]+STATES[44]; ALGEBRAIC[124] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[77]*ALGEBRAIC[45]); ALGEBRAIC[139] = CONSTANTS[26]*(exp((ALGEBRAIC[123]+ALGEBRAIC[121])/( CONSTANTS[144]*CONSTANTS[145])) - exp(ALGEBRAIC[124]/( CONSTANTS[144]*CONSTANTS[145]))); ALGEBRAIC[1] = (CONDVAR[6]<0.00000&&CONDVAR[7]>0.00000 ? CONSTANTS[95]+ CONSTANTS[149]*((VOI - CONSTANTS[91])+CONSTANTS[93]) : CONDVAR[8]>=0.00000&&CONDVAR[9]<0.00000 ? CONSTANTS[94]+CONSTANTS[95] : CONDVAR[10]<0.00000&&CONDVAR[11]>=0.00000 ? CONSTANTS[95]+ - CONSTANTS[149]*(((VOI - CONSTANTS[91]) - CONSTANTS[93]) - CONSTANTS[92]) : CONSTANTS[95]); ALGEBRAIC[40] = ALGEBRAIC[1]+STATES[39]; ALGEBRAIC[119] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[72]*ALGEBRAIC[40]); ALGEBRAIC[140] = CONSTANTS[27]*(exp((ALGEBRAIC[117]+ALGEBRAIC[119])/( CONSTANTS[144]*CONSTANTS[145])) - exp(ALGEBRAIC[118]/( CONSTANTS[144]*CONSTANTS[145]))); ALGEBRAIC[125] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[49]*ALGEBRAIC[17]); ALGEBRAIC[46] = CONSTANTS[138]+STATES[46]; ALGEBRAIC[126] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[78]*ALGEBRAIC[46]); ALGEBRAIC[128] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[62]*ALGEBRAIC[30]); ALGEBRAIC[129] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[63]*ALGEBRAIC[31]); ALGEBRAIC[48] = CONSTANTS[140]+STATES[51]; ALGEBRAIC[131] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[80]*ALGEBRAIC[48]); ALGEBRAIC[141] = CONSTANTS[28]*(exp((ALGEBRAIC[122]+ALGEBRAIC[128])/( CONSTANTS[144]*CONSTANTS[145])) - exp((ALGEBRAIC[125]+ALGEBRAIC[126]+ALGEBRAIC[131]+ALGEBRAIC[129])/( CONSTANTS[144]*CONSTANTS[145]))); ALGEBRAIC[49] = CONSTANTS[141]+STATES[52]; ALGEBRAIC[132] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[81]*ALGEBRAIC[49]); ALGEBRAIC[142] = CONSTANTS[29]*(exp((ALGEBRAIC[124]+ALGEBRAIC[128])/( CONSTANTS[144]*CONSTANTS[145])) - exp((ALGEBRAIC[125]+ALGEBRAIC[126]+ALGEBRAIC[132]+ALGEBRAIC[129])/( CONSTANTS[144]*CONSTANTS[145]))); ALGEBRAIC[47] = CONSTANTS[139]+STATES[47]; ALGEBRAIC[127] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[79]*ALGEBRAIC[47]); ALGEBRAIC[130] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[64]*ALGEBRAIC[32]); ALGEBRAIC[143] = CONSTANTS[30]*(exp(ALGEBRAIC[125]/( CONSTANTS[144]*CONSTANTS[145])) - exp((ALGEBRAIC[127]+ALGEBRAIC[130])/( CONSTANTS[144]*CONSTANTS[145]))); ALGEBRAIC[144] = CONSTANTS[31]*(exp((ALGEBRAIC[127]+ALGEBRAIC[126])/( CONSTANTS[144]*CONSTANTS[145])) - exp(ALGEBRAIC[121]/( CONSTANTS[144]*CONSTANTS[145]))); ALGEBRAIC[50] = CONSTANTS[142]+STATES[53]; ALGEBRAIC[133] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[82]*ALGEBRAIC[50]); ALGEBRAIC[135] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[69]*ALGEBRAIC[37]); ALGEBRAIC[145] = CONSTANTS[32]*(exp((ALGEBRAIC[131]+ALGEBRAIC[135])/( CONSTANTS[144]*CONSTANTS[145])) - exp(ALGEBRAIC[133]/( CONSTANTS[144]*CONSTANTS[145]))); ALGEBRAIC[51] = CONSTANTS[143]+STATES[54]; ALGEBRAIC[134] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[83]*ALGEBRAIC[51]); ALGEBRAIC[146] = CONSTANTS[33]*(exp((ALGEBRAIC[132]+ALGEBRAIC[135])/( CONSTANTS[144]*CONSTANTS[145])) - exp(ALGEBRAIC[134]/( CONSTANTS[144]*CONSTANTS[145]))); ALGEBRAIC[147] = (ALGEBRAIC[138])+ALGEBRAIC[140]; } 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; } void computeRoots(double VOI, double* CONSTANTS, double* RATES, double* OLDRATES, double* STATES, double* OLDSTATES, double* ALGEBRAIC, double* CONDVARS) { CONDVAR[0] = VOI - CONSTANTS[86]; CONDVAR[1] = VOI - (CONSTANTS[86] - CONSTANTS[88]); CONDVAR[2] = VOI - CONSTANTS[86]; CONDVAR[3] = VOI - (CONSTANTS[86]+CONSTANTS[87]); CONDVAR[4] = VOI - (CONSTANTS[86]+CONSTANTS[88]+CONSTANTS[87]); CONDVAR[5] = VOI - (CONSTANTS[86]+CONSTANTS[87]); CONDVAR[6] = VOI - CONSTANTS[91]; CONDVAR[7] = VOI - (CONSTANTS[91] - CONSTANTS[93]); CONDVAR[8] = VOI - CONSTANTS[91]; CONDVAR[9] = VOI - (CONSTANTS[91]+CONSTANTS[92]); CONDVAR[10] = VOI - (CONSTANTS[91]+CONSTANTS[93]+CONSTANTS[92]); CONDVAR[11] = VOI - (CONSTANTS[91]+CONSTANTS[92]); }