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