Generated Code
The following is c 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[58] is LB1_T in component environment (fmol).
* ALGEBRAIC[59] is RB1_T in component environment (fmol).
* ALGEBRAIC[43] is Gs_T in component environment (fmol).
* ALGEBRAIC[16] 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[13] is q_PDE_cAMP in component environment (fmol).
* ALGEBRAIC[14] is q_five_AMP in component environment (fmol).
* ALGEBRAIC[17] is q_IBMX in component environment (fmol).
* ALGEBRAIC[18] is q_PDEinh in component environment (fmol).
* ALGEBRAIC[19] is q_Gsa_GTP in component environment (fmol).
* ALGEBRAIC[20] is q_FSK in component environment (fmol).
* ALGEBRAIC[22] is q_Gia_GTP in component environment (fmol).
* ALGEBRAIC[23] is q_ACinh in component environment (fmol).
* ALGEBRAIC[25] is q_PPi in component environment (fmol).
* ALGEBRAIC[24] is q_RB1_inactive in component environment (fmol).
* ALGEBRAIC[26] is q_L_RB1_inactive in component environment (fmol).
* ALGEBRAIC[27] is q_LB1 in component environment (fmol).
* ALGEBRAIC[29] is q_RB1 in component environment (fmol).
* ALGEBRAIC[31] is q_Gs in component environment (fmol).
* ALGEBRAIC[33] is q_RB1_Gs in component environment (fmol).
* ALGEBRAIC[35] is q_L_RB1 in component environment (fmol).
* ALGEBRAIC[37] is q_L_RB1_Gs in component environment (fmol).
* ALGEBRAIC[39] is q_Gsbetagamma in component environment (fmol).
* ALGEBRAIC[41] is q_Gsa_GDP in component environment (fmol).
* ALGEBRAIC[44] is q_GTP in component environment (fmol).
* ALGEBRAIC[46] is q_GDP in component environment (fmol).
* ALGEBRAIC[48] is q_Pi in component environment (fmol).
* ALGEBRAIC[50] is q_RB1_tag in component environment (fmol).
* ALGEBRAIC[52] is q_L_RB1_tag in component environment (fmol).
* ALGEBRAIC[54] is q_RB1_GRKArr in component environment (fmol).
* ALGEBRAIC[56] is q_L_RB1_GRKArr in component environment (fmol).
* ALGEBRAIC[60] is q_GRKArr in component environment (fmol).
* ALGEBRAIC[62] is q_RM2_inactive in component environment (fmol).
* ALGEBRAIC[65] is q_L_RM2_inactive in component environment (fmol).
* ALGEBRAIC[68] is q_LM2 in component environment (fmol).
* ALGEBRAIC[71] is q_RM2 in component environment (fmol).
* ALGEBRAIC[74] is q_Gi in component environment (fmol).
* ALGEBRAIC[77] is q_RM2_Gi in component environment (fmol).
* ALGEBRAIC[80] is q_L_RM2 in component environment (fmol).
* ALGEBRAIC[85] is q_L_RM2_Gi in component environment (fmol).
* ALGEBRAIC[89] is q_Gibetagamma in component environment (fmol).
* ALGEBRAIC[93] is q_Gia_GDP in component environment (fmol).
* ALGEBRAIC[97] is q_RM2_tag in component environment (fmol).
* ALGEBRAIC[99] is q_L_RM2_tag in component environment (fmol).
* ALGEBRAIC[101] is q_RM2_GRKArr in component environment (fmol).
* ALGEBRAIC[103] 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[69] is v1a in component cAMP (fmol_per_sec).
* ALGEBRAIC[72] is v1b in component cAMP (fmol_per_sec).
* ALGEBRAIC[75] is v2a in component cAMP (fmol_per_sec).
* ALGEBRAIC[78] is v2b in component cAMP (fmol_per_sec).
* ALGEBRAIC[81] is v3a in component cAMP (fmol_per_sec).
* ALGEBRAIC[83] is v3b in component cAMP (fmol_per_sec).
* ALGEBRAIC[86] is v4a in component cAMP (fmol_per_sec).
* ALGEBRAIC[90] is v4b in component cAMP (fmol_per_sec).
* ALGEBRAIC[94] is v5 in component cAMP (fmol_per_sec).
* ALGEBRAIC[87] is v6 in component cAMP (fmol_per_sec).
* ALGEBRAIC[91] is v7 in component cAMP (fmol_per_sec).
* ALGEBRAIC[95] is vGiAC in component cAMP (fmol_per_sec).
* ALGEBRAIC[28] is mu_ATP in component cAMP (J_per_mol).
* ALGEBRAIC[32] is mu_AC in component cAMP (J_per_mol).
* ALGEBRAIC[30] is mu_cAMP in component cAMP (J_per_mol).
* ALGEBRAIC[34] is mu_AC_ATP in component cAMP (J_per_mol).
* ALGEBRAIC[57] is mu_FSK in component cAMP (J_per_mol).
* ALGEBRAIC[40] is mu_FSK_AC in component cAMP (J_per_mol).
* ALGEBRAIC[42] is mu_FSK_AC_ATP in component cAMP (J_per_mol).
* ALGEBRAIC[55] is mu_Gsa_GTP in component cAMP (J_per_mol).
* ALGEBRAIC[36] is mu_Gsa_GTP_AC in component cAMP (J_per_mol).
* ALGEBRAIC[38] is mu_Gsa_GTP_AC_ATP in component cAMP (J_per_mol).
* ALGEBRAIC[45] is mu_PDE in component cAMP (J_per_mol).
* ALGEBRAIC[53] is mu_PDEinh in component cAMP (J_per_mol).
* ALGEBRAIC[47] is mu_PDE_cAMP in component cAMP (J_per_mol).
* ALGEBRAIC[51] is mu_IBMX in component cAMP (J_per_mol).
* ALGEBRAIC[49] is mu_five_AMP in component cAMP (J_per_mol).
* ALGEBRAIC[61] is mu_Gia_GTP in component cAMP (J_per_mol).
* ALGEBRAIC[63] is mu_ACinh in component cAMP (J_per_mol).
* ALGEBRAIC[66] is mu_PPi in component cAMP (J_per_mol).
* CONSTANTS[147] is vol in component cAMP (pL).
* ALGEBRAIC[11] is ATP_T in component cAMP (fmol).
* ALGEBRAIC[12] is AC_T in component cAMP (fmol).
* ALGEBRAIC[21] is Gs_T in component cAMP (fmol).
* ALGEBRAIC[15] is cAMP_T in component cAMP (fmol).
* ALGEBRAIC[64] is mu_RB1_inactive in component GPCR_B1AR_reduced (J_per_mol).
* ALGEBRAIC[67] is mu_L_RB1_inactive in component GPCR_B1AR_reduced (J_per_mol).
* ALGEBRAIC[70] is mu_LB1 in component GPCR_B1AR_reduced (J_per_mol).
* ALGEBRAIC[73] is mu_RB1 in component GPCR_B1AR_reduced (J_per_mol).
* ALGEBRAIC[76] is mu_Gs in component GPCR_B1AR_reduced (J_per_mol).
* ALGEBRAIC[79] is mu_RB1_Gs in component GPCR_B1AR_reduced (J_per_mol).
* ALGEBRAIC[82] is mu_L_RB1 in component GPCR_B1AR_reduced (J_per_mol).
* ALGEBRAIC[84] is mu_L_RB1_Gs in component GPCR_B1AR_reduced (J_per_mol).
* ALGEBRAIC[88] is mu_Gsa_GTP in component GPCR_B1AR_reduced (J_per_mol).
* ALGEBRAIC[92] 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[98] is mu_GTP in component GPCR_B1AR_reduced (J_per_mol).
* ALGEBRAIC[100] is mu_GDP in component GPCR_B1AR_reduced (J_per_mol).
* ALGEBRAIC[102] is mu_Pi in component GPCR_B1AR_reduced (J_per_mol).
* ALGEBRAIC[104] is mu_RB1_tag in component GPCR_B1AR_reduced (J_per_mol).
* ALGEBRAIC[106] is mu_L_RB1_tag in component GPCR_B1AR_reduced (J_per_mol).
* ALGEBRAIC[108] is mu_RB1_GRKArr in component GPCR_B1AR_reduced (J_per_mol).
* ALGEBRAIC[110] is mu_L_RB1_GRKArr in component GPCR_B1AR_reduced (J_per_mol).
* ALGEBRAIC[112] is mu_GRKArr in component GPCR_B1AR_reduced (J_per_mol).
* ALGEBRAIC[114] is v_Rswitch_B1 in component GPCR_B1AR_reduced (fmol_per_sec).
* ALGEBRAIC[116] is v_LRswitch_B1 in component GPCR_B1AR_reduced (fmol_per_sec).
* ALGEBRAIC[118] is v_C_B1 in component GPCR_B1AR_reduced (fmol_per_sec).
* ALGEBRAIC[121] is v_R_B1 in component GPCR_B1AR_reduced (fmol_per_sec).
* ALGEBRAIC[119] is v_L_B1 in component GPCR_B1AR_reduced (fmol_per_sec).
* ALGEBRAIC[124] is v_Act1_Gs in component GPCR_B1AR_reduced (fmol_per_sec).
* ALGEBRAIC[126] is v_Act2_Gs in component GPCR_B1AR_reduced (fmol_per_sec).
* ALGEBRAIC[128] is v_Hyd_Gs in component GPCR_B1AR_reduced (fmol_per_sec).
* ALGEBRAIC[131] is v_Reassoc_Gs in component GPCR_B1AR_reduced (fmol_per_sec).
* ALGEBRAIC[129] is v_InternR_B1 in component GPCR_B1AR_reduced (fmol_per_sec).
* ALGEBRAIC[132] is v_InternLR_B1 in component GPCR_B1AR_reduced (fmol_per_sec).
* ALGEBRAIC[122] is v_Rsynthesis in component GPCR_B1AR_reduced (fmol_per_sec).
* ALGEBRAIC[105] is mu_RM2_inactive in component GPCR_M2_reduced (J_per_mol).
* ALGEBRAIC[107] is mu_L_RM2_inactive in component GPCR_M2_reduced (J_per_mol).
* ALGEBRAIC[109] is mu_LM2 in component GPCR_M2_reduced (J_per_mol).
* ALGEBRAIC[111] is mu_RM2 in component GPCR_M2_reduced (J_per_mol).
* ALGEBRAIC[113] is mu_Gi in component GPCR_M2_reduced (J_per_mol).
* ALGEBRAIC[115] is mu_RM2_Gi in component GPCR_M2_reduced (J_per_mol).
* ALGEBRAIC[117] is mu_L_RM2 in component GPCR_M2_reduced (J_per_mol).
* ALGEBRAIC[120] is mu_L_RM2_Gi in component GPCR_M2_reduced (J_per_mol).
* ALGEBRAIC[123] is mu_Gia_GTP in component GPCR_M2_reduced (J_per_mol).
* ALGEBRAIC[125] 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[130] is mu_GTP in component GPCR_M2_reduced (J_per_mol).
* ALGEBRAIC[133] is mu_GDP in component GPCR_M2_reduced (J_per_mol).
* ALGEBRAIC[134] is mu_Pi in component GPCR_M2_reduced (J_per_mol).
* ALGEBRAIC[135] is mu_RM2_tag in component GPCR_M2_reduced (J_per_mol).
* ALGEBRAIC[136] is mu_L_RM2_tag in component GPCR_M2_reduced (J_per_mol).
* ALGEBRAIC[137] is mu_RM2_GRKArr in component GPCR_M2_reduced (J_per_mol).
* ALGEBRAIC[138] is mu_L_RM2_GRKArr in component GPCR_M2_reduced (J_per_mol).
* ALGEBRAIC[139] is mu_GRKArr in component GPCR_M2_reduced (J_per_mol).
* ALGEBRAIC[140] is v_Rswitch_M2 in component GPCR_M2_reduced (fmol_per_sec).
* ALGEBRAIC[141] is v_LRswitch_M2 in component GPCR_M2_reduced (fmol_per_sec).
* ALGEBRAIC[142] is v_C_M2 in component GPCR_M2_reduced (fmol_per_sec).
* ALGEBRAIC[143] is v_R_M2 in component GPCR_M2_reduced (fmol_per_sec).
* ALGEBRAIC[144] is v_L_M2 in component GPCR_M2_reduced (fmol_per_sec).
* ALGEBRAIC[145] is v_Act1_Gi in component GPCR_M2_reduced (fmol_per_sec).
* ALGEBRAIC[146] is v_Act2_Gi in component GPCR_M2_reduced (fmol_per_sec).
* ALGEBRAIC[148] is v_Hyd_Gi in component GPCR_M2_reduced (fmol_per_sec).
* ALGEBRAIC[150] is v_Reassoc_Gi in component GPCR_M2_reduced (fmol_per_sec).
* ALGEBRAIC[149] is v_InternR_M2 in component GPCR_M2_reduced (fmol_per_sec).
* ALGEBRAIC[151] 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).
*/
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];
}
void
computeRates(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[2] = CONSTANTS[96]+STATES[0];
ALGEBRAIC[28] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[34]*ALGEBRAIC[2]);
ALGEBRAIC[4] = CONSTANTS[97]+STATES[2];
ALGEBRAIC[32] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[36]*ALGEBRAIC[4]);
ALGEBRAIC[5] = CONSTANTS[99]+STATES[3];
ALGEBRAIC[34] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[37]*ALGEBRAIC[5]);
ALGEBRAIC[69] = CONSTANTS[0]*(exp((ALGEBRAIC[32]+ALGEBRAIC[28])/( CONSTANTS[144]*CONSTANTS[145])) - exp(ALGEBRAIC[34]/( CONSTANTS[144]*CONSTANTS[145])));
ALGEBRAIC[3] = CONSTANTS[98]+STATES[1];
ALGEBRAIC[30] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[35]*ALGEBRAIC[3]);
ALGEBRAIC[25] = CONSTANTS[113]+STATES[17];
ALGEBRAIC[66] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[51]*ALGEBRAIC[25]);
ALGEBRAIC[72] = CONSTANTS[1]*(exp(ALGEBRAIC[34]/( CONSTANTS[144]*CONSTANTS[145])) - exp((ALGEBRAIC[32]+ALGEBRAIC[30]+ALGEBRAIC[66])/( CONSTANTS[144]*CONSTANTS[145])));
RATES[3] = ALGEBRAIC[69] - ALGEBRAIC[72];
ALGEBRAIC[6] = CONSTANTS[104]+STATES[4];
ALGEBRAIC[36] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[38]*ALGEBRAIC[6]);
ALGEBRAIC[7] = CONSTANTS[105]+STATES[5];
ALGEBRAIC[38] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[39]*ALGEBRAIC[7]);
ALGEBRAIC[75] = CONSTANTS[2]*(exp((ALGEBRAIC[36]+ALGEBRAIC[28])/( CONSTANTS[144]*CONSTANTS[145])) - exp(ALGEBRAIC[38]/( CONSTANTS[144]*CONSTANTS[145])));
ALGEBRAIC[78] = CONSTANTS[3]*(exp(ALGEBRAIC[38]/( CONSTANTS[144]*CONSTANTS[145])) - exp((ALGEBRAIC[36]+ALGEBRAIC[30]+ALGEBRAIC[66])/( CONSTANTS[144]*CONSTANTS[145])));
RATES[5] = ALGEBRAIC[75] - ALGEBRAIC[78];
ALGEBRAIC[8] = CONSTANTS[101]+STATES[6];
ALGEBRAIC[40] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[40]*ALGEBRAIC[8]);
ALGEBRAIC[9] = CONSTANTS[102]+STATES[7];
ALGEBRAIC[42] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[41]*ALGEBRAIC[9]);
ALGEBRAIC[81] = CONSTANTS[4]*(exp((ALGEBRAIC[40]+ALGEBRAIC[28])/( CONSTANTS[144]*CONSTANTS[145])) - exp(ALGEBRAIC[42]/( CONSTANTS[144]*CONSTANTS[145])));
RATES[0] = (- ALGEBRAIC[69] - ALGEBRAIC[81]) - ALGEBRAIC[75];
ALGEBRAIC[83] = CONSTANTS[5]*(exp(ALGEBRAIC[42]/( CONSTANTS[144]*CONSTANTS[145])) - exp((ALGEBRAIC[40]+ALGEBRAIC[30]+ALGEBRAIC[66])/( CONSTANTS[144]*CONSTANTS[145])));
RATES[7] = ALGEBRAIC[81] - ALGEBRAIC[83];
ALGEBRAIC[10] = CONSTANTS[106]+STATES[8];
ALGEBRAIC[45] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[42]*ALGEBRAIC[10]);
ALGEBRAIC[13] = CONSTANTS[108]+STATES[9];
ALGEBRAIC[47] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[43]*ALGEBRAIC[13]);
ALGEBRAIC[86] = CONSTANTS[6]*(exp((ALGEBRAIC[45]+ALGEBRAIC[30])/( CONSTANTS[144]*CONSTANTS[145])) - exp(ALGEBRAIC[47]/( CONSTANTS[144]*CONSTANTS[145])));
RATES[1] = (ALGEBRAIC[72]+ALGEBRAIC[83]+ALGEBRAIC[78]) - ALGEBRAIC[86];
ALGEBRAIC[19] = CONSTANTS[103]+STATES[13]+STATES[26];
ALGEBRAIC[55] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[47]*ALGEBRAIC[19]);
ALGEBRAIC[87] = CONSTANTS[9]*(exp((ALGEBRAIC[32]+ALGEBRAIC[55])/( CONSTANTS[144]*CONSTANTS[145])) - exp(ALGEBRAIC[36]/( CONSTANTS[144]*CONSTANTS[145])));
RATES[13] = - ALGEBRAIC[87];
RATES[4] = (ALGEBRAIC[87] - ALGEBRAIC[75])+ALGEBRAIC[78];
ALGEBRAIC[20] = CONSTANTS[100]+STATES[14];
ALGEBRAIC[57] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[48]*ALGEBRAIC[20]);
ALGEBRAIC[91] = CONSTANTS[10]*(exp((ALGEBRAIC[57]+ALGEBRAIC[32])/( CONSTANTS[144]*CONSTANTS[145])) - exp(ALGEBRAIC[40]/( CONSTANTS[144]*CONSTANTS[145])));
RATES[14] = - ALGEBRAIC[91];
RATES[6] = (ALGEBRAIC[91]+ALGEBRAIC[83]) - ALGEBRAIC[81];
ALGEBRAIC[14] = CONSTANTS[110]+STATES[10];
ALGEBRAIC[49] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[44]*ALGEBRAIC[14]);
ALGEBRAIC[90] = CONSTANTS[7]*(exp(ALGEBRAIC[47]/( CONSTANTS[144]*CONSTANTS[145])) - exp((ALGEBRAIC[45]+ALGEBRAIC[49])/( CONSTANTS[144]*CONSTANTS[145])));
RATES[9] = ALGEBRAIC[86] - ALGEBRAIC[90];
RATES[10] = ALGEBRAIC[90];
ALGEBRAIC[22] = CONSTANTS[111]+STATES[15]+STATES[45];
ALGEBRAIC[61] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[49]*ALGEBRAIC[22]);
ALGEBRAIC[23] = CONSTANTS[112]+STATES[16];
ALGEBRAIC[63] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[50]*ALGEBRAIC[23]);
ALGEBRAIC[95] = CONSTANTS[11]*(exp((ALGEBRAIC[32]+ALGEBRAIC[61])/( CONSTANTS[144]*CONSTANTS[145])) - exp(ALGEBRAIC[63]/( CONSTANTS[144]*CONSTANTS[145])));
RATES[2] = (((ALGEBRAIC[72] - ALGEBRAIC[69]) - ALGEBRAIC[87]) - ALGEBRAIC[91]) - ALGEBRAIC[95];
ALGEBRAIC[18] = CONSTANTS[107]+STATES[12];
ALGEBRAIC[53] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[46]*ALGEBRAIC[18]);
ALGEBRAIC[17] = CONSTANTS[109]+STATES[11];
ALGEBRAIC[51] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[45]*ALGEBRAIC[17]);
ALGEBRAIC[94] = CONSTANTS[8]*(exp((ALGEBRAIC[45]+ALGEBRAIC[51])/( CONSTANTS[144]*CONSTANTS[145])) - exp(ALGEBRAIC[53]/( CONSTANTS[144]*CONSTANTS[145])));
RATES[8] = (ALGEBRAIC[90] - ALGEBRAIC[86]) - ALGEBRAIC[94];
RATES[11] = - ALGEBRAIC[94];
RATES[12] = ALGEBRAIC[94];
RATES[15] = - ALGEBRAIC[95];
RATES[16] = ALGEBRAIC[95];
RATES[17] = ALGEBRAIC[95];
ALGEBRAIC[26] = CONSTANTS[115]+STATES[19];
ALGEBRAIC[67] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[53]*ALGEBRAIC[26]);
ALGEBRAIC[35] = CONSTANTS[119]+STATES[24];
ALGEBRAIC[82] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[58]*ALGEBRAIC[35]);
ALGEBRAIC[116] = CONSTANTS[13]*(exp(ALGEBRAIC[67]/( CONSTANTS[144]*CONSTANTS[145])) - exp(ALGEBRAIC[82]/( CONSTANTS[144]*CONSTANTS[145])));
ALGEBRAIC[24] = CONSTANTS[114]+STATES[18];
ALGEBRAIC[64] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[52]*ALGEBRAIC[24]);
ALGEBRAIC[0] = (VOI<CONSTANTS[86]&&VOI>CONSTANTS[86] - CONSTANTS[88] ? CONSTANTS[90]+ CONSTANTS[148]*((VOI - CONSTANTS[86])+CONSTANTS[88]) : VOI>=CONSTANTS[86]&&VOI<CONSTANTS[86]+CONSTANTS[87] ? CONSTANTS[89]+CONSTANTS[90] : VOI<CONSTANTS[86]+CONSTANTS[88]+CONSTANTS[87]&&VOI>=CONSTANTS[86]+CONSTANTS[87] ? CONSTANTS[90]+ - CONSTANTS[148]*(((VOI - CONSTANTS[86]) - CONSTANTS[88]) - CONSTANTS[87]) : CONSTANTS[90]);
ALGEBRAIC[27] = ALGEBRAIC[0]+STATES[20];
ALGEBRAIC[70] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[54]*ALGEBRAIC[27]);
ALGEBRAIC[119] = CONSTANTS[16]*(exp((ALGEBRAIC[64]+ALGEBRAIC[70])/( CONSTANTS[144]*CONSTANTS[145])) - exp(ALGEBRAIC[67]/( CONSTANTS[144]*CONSTANTS[145])));
RATES[19] = - ALGEBRAIC[116]+ALGEBRAIC[119];
RATES[20] = - ALGEBRAIC[119];
ALGEBRAIC[29] = CONSTANTS[116]+STATES[21];
ALGEBRAIC[73] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[55]*ALGEBRAIC[29]);
ALGEBRAIC[114] = CONSTANTS[12]*(exp(ALGEBRAIC[64]/( CONSTANTS[144]*CONSTANTS[145])) - exp(ALGEBRAIC[73]/( CONSTANTS[144]*CONSTANTS[145])));
ALGEBRAIC[31] = CONSTANTS[117]+STATES[22];
ALGEBRAIC[76] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[56]*ALGEBRAIC[31]);
ALGEBRAIC[33] = CONSTANTS[118]+STATES[23];
ALGEBRAIC[79] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[57]*ALGEBRAIC[33]);
ALGEBRAIC[118] = CONSTANTS[14]*(exp((ALGEBRAIC[73]+ALGEBRAIC[76])/( CONSTANTS[144]*CONSTANTS[145])) - exp(ALGEBRAIC[79]/( CONSTANTS[144]*CONSTANTS[145])));
RATES[21] = ALGEBRAIC[114] - ALGEBRAIC[118];
ALGEBRAIC[122] = ALGEBRAIC[114]+ALGEBRAIC[119];
RATES[18] = (- ALGEBRAIC[114] - ALGEBRAIC[119])+ 1.00000*ALGEBRAIC[122];
ALGEBRAIC[37] = CONSTANTS[120]+STATES[25];
ALGEBRAIC[84] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[59]*ALGEBRAIC[37]);
ALGEBRAIC[121] = CONSTANTS[15]*(exp((ALGEBRAIC[82]+ALGEBRAIC[76])/( CONSTANTS[144]*CONSTANTS[145])) - exp(ALGEBRAIC[84]/( CONSTANTS[144]*CONSTANTS[145])));
RATES[24] = ALGEBRAIC[116] - ALGEBRAIC[121];
ALGEBRAIC[88] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[47]*ALGEBRAIC[19]);
ALGEBRAIC[39] = CONSTANTS[121]+STATES[27];
ALGEBRAIC[92] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[60]*ALGEBRAIC[39]);
ALGEBRAIC[44] = CONSTANTS[123]+STATES[29]+STATES[48];
ALGEBRAIC[98] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[62]*ALGEBRAIC[44]);
ALGEBRAIC[46] = CONSTANTS[124]+STATES[30]+STATES[49];
ALGEBRAIC[100] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[63]*ALGEBRAIC[46]);
ALGEBRAIC[50] = CONSTANTS[126]+STATES[32];
ALGEBRAIC[104] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[65]*ALGEBRAIC[50]);
ALGEBRAIC[124] = CONSTANTS[17]*(exp((ALGEBRAIC[79]+ALGEBRAIC[98])/( CONSTANTS[144]*CONSTANTS[145])) - exp((ALGEBRAIC[88]+ALGEBRAIC[92]+ALGEBRAIC[104]+ALGEBRAIC[100])/( CONSTANTS[144]*CONSTANTS[145])));
RATES[23] = ALGEBRAIC[118] - ALGEBRAIC[124];
ALGEBRAIC[52] = CONSTANTS[127]+STATES[33];
ALGEBRAIC[106] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[66]*ALGEBRAIC[52]);
ALGEBRAIC[126] = CONSTANTS[18]*(exp((ALGEBRAIC[84]+ALGEBRAIC[98])/( CONSTANTS[144]*CONSTANTS[145])) - exp((ALGEBRAIC[88]+ALGEBRAIC[92]+ALGEBRAIC[106]+ALGEBRAIC[100])/( CONSTANTS[144]*CONSTANTS[145])));
RATES[25] = ALGEBRAIC[121] - ALGEBRAIC[126];
RATES[29] = - ALGEBRAIC[124] - ALGEBRAIC[126];
RATES[30] = ALGEBRAIC[124]+ALGEBRAIC[126];
ALGEBRAIC[41] = CONSTANTS[122]+STATES[28];
ALGEBRAIC[96] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[61]*ALGEBRAIC[41]);
ALGEBRAIC[48] = CONSTANTS[125]+STATES[31]+STATES[50];
ALGEBRAIC[102] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[64]*ALGEBRAIC[48]);
ALGEBRAIC[128] = CONSTANTS[19]*(exp(ALGEBRAIC[88]/( CONSTANTS[144]*CONSTANTS[145])) - exp((ALGEBRAIC[96]+ALGEBRAIC[102])/( CONSTANTS[144]*CONSTANTS[145])));
RATES[26] = (ALGEBRAIC[124]+ALGEBRAIC[126]) - ALGEBRAIC[128];
RATES[31] = ALGEBRAIC[128];
ALGEBRAIC[54] = CONSTANTS[128]+STATES[34];
ALGEBRAIC[108] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[67]*ALGEBRAIC[54]);
ALGEBRAIC[60] = CONSTANTS[130]+STATES[36]+STATES[55];
ALGEBRAIC[112] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[69]*ALGEBRAIC[60]);
ALGEBRAIC[129] = CONSTANTS[21]*(exp((ALGEBRAIC[104]+ALGEBRAIC[112])/( CONSTANTS[144]*CONSTANTS[145])) - exp(ALGEBRAIC[108]/( CONSTANTS[144]*CONSTANTS[145])));
RATES[32] = ALGEBRAIC[124] - ALGEBRAIC[129];
RATES[34] = ALGEBRAIC[129];
ALGEBRAIC[131] = CONSTANTS[20]*(exp((ALGEBRAIC[96]+ALGEBRAIC[92])/( CONSTANTS[144]*CONSTANTS[145])) - exp(ALGEBRAIC[76]/( CONSTANTS[144]*CONSTANTS[145])));
RATES[22] = (- ALGEBRAIC[118] - ALGEBRAIC[121])+ALGEBRAIC[131];
RATES[27] = (ALGEBRAIC[124]+ALGEBRAIC[126]) - ALGEBRAIC[131];
RATES[28] = ALGEBRAIC[128] - ALGEBRAIC[131];
ALGEBRAIC[56] = CONSTANTS[129]+STATES[35];
ALGEBRAIC[110] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[68]*ALGEBRAIC[56]);
ALGEBRAIC[132] = CONSTANTS[22]*(exp((ALGEBRAIC[106]+ALGEBRAIC[112])/( CONSTANTS[144]*CONSTANTS[145])) - exp(ALGEBRAIC[110]/( CONSTANTS[144]*CONSTANTS[145])));
RATES[33] = ALGEBRAIC[126] - ALGEBRAIC[132];
RATES[35] = ALGEBRAIC[132];
RATES[36] = - ALGEBRAIC[129] - ALGEBRAIC[132];
ALGEBRAIC[65] = CONSTANTS[132]+STATES[38];
ALGEBRAIC[107] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[71]*ALGEBRAIC[65]);
ALGEBRAIC[80] = CONSTANTS[136]+STATES[43];
ALGEBRAIC[117] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[76]*ALGEBRAIC[80]);
ALGEBRAIC[141] = CONSTANTS[24]*(exp(ALGEBRAIC[107]/( CONSTANTS[144]*CONSTANTS[145])) - exp(ALGEBRAIC[117]/( CONSTANTS[144]*CONSTANTS[145])));
RATES[38] = - ALGEBRAIC[141];
ALGEBRAIC[62] = CONSTANTS[131]+STATES[37];
ALGEBRAIC[105] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[70]*ALGEBRAIC[62]);
ALGEBRAIC[71] = CONSTANTS[133]+STATES[40];
ALGEBRAIC[111] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[73]*ALGEBRAIC[71]);
ALGEBRAIC[140] = CONSTANTS[23]*(exp(ALGEBRAIC[105]/( CONSTANTS[144]*CONSTANTS[145])) - exp(ALGEBRAIC[111]/( CONSTANTS[144]*CONSTANTS[145])));
ALGEBRAIC[1] = (VOI<CONSTANTS[91]&&VOI>CONSTANTS[91] - CONSTANTS[93] ? CONSTANTS[95]+ CONSTANTS[149]*((VOI - CONSTANTS[91])+CONSTANTS[93]) : VOI>=CONSTANTS[91]&&VOI<CONSTANTS[91]+CONSTANTS[92] ? CONSTANTS[94]+CONSTANTS[95] : VOI<CONSTANTS[91]+CONSTANTS[93]+CONSTANTS[92]&&VOI>=CONSTANTS[91]+CONSTANTS[92] ? CONSTANTS[95]+ - CONSTANTS[149]*(((VOI - CONSTANTS[91]) - CONSTANTS[93]) - CONSTANTS[92]) : CONSTANTS[95]);
ALGEBRAIC[68] = ALGEBRAIC[1]+STATES[39];
ALGEBRAIC[109] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[72]*ALGEBRAIC[68]);
ALGEBRAIC[144] = CONSTANTS[27]*(exp((ALGEBRAIC[105]+ALGEBRAIC[109])/( CONSTANTS[144]*CONSTANTS[145])) - exp(ALGEBRAIC[107]/( CONSTANTS[144]*CONSTANTS[145])));
RATES[37] = - ALGEBRAIC[140] - ALGEBRAIC[144];
RATES[39] = - ALGEBRAIC[144];
ALGEBRAIC[74] = CONSTANTS[134]+STATES[41];
ALGEBRAIC[113] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[74]*ALGEBRAIC[74]);
ALGEBRAIC[77] = CONSTANTS[135]+STATES[42];
ALGEBRAIC[115] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[75]*ALGEBRAIC[77]);
ALGEBRAIC[142] = CONSTANTS[25]*(exp((ALGEBRAIC[111]+ALGEBRAIC[113])/( CONSTANTS[144]*CONSTANTS[145])) - exp(ALGEBRAIC[115]/( CONSTANTS[144]*CONSTANTS[145])));
ALGEBRAIC[123] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[49]*ALGEBRAIC[22]);
ALGEBRAIC[89] = CONSTANTS[138]+STATES[46];
ALGEBRAIC[125] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[78]*ALGEBRAIC[89]);
ALGEBRAIC[130] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[62]*ALGEBRAIC[44]);
ALGEBRAIC[133] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[63]*ALGEBRAIC[46]);
ALGEBRAIC[97] = CONSTANTS[140]+STATES[51];
ALGEBRAIC[135] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[80]*ALGEBRAIC[97]);
ALGEBRAIC[145] = CONSTANTS[28]*(exp((ALGEBRAIC[115]+ALGEBRAIC[130])/( CONSTANTS[144]*CONSTANTS[145])) - exp((ALGEBRAIC[123]+ALGEBRAIC[125]+ALGEBRAIC[135]+ALGEBRAIC[133])/( CONSTANTS[144]*CONSTANTS[145])));
RATES[42] = ALGEBRAIC[142] - ALGEBRAIC[145];
ALGEBRAIC[85] = CONSTANTS[137]+STATES[44];
ALGEBRAIC[120] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[77]*ALGEBRAIC[85]);
ALGEBRAIC[143] = CONSTANTS[26]*(exp((ALGEBRAIC[117]+ALGEBRAIC[113])/( CONSTANTS[144]*CONSTANTS[145])) - exp(ALGEBRAIC[120]/( CONSTANTS[144]*CONSTANTS[145])));
RATES[43] = (ALGEBRAIC[141] - ALGEBRAIC[143])+ALGEBRAIC[144];
ALGEBRAIC[147] = (ALGEBRAIC[142])+ALGEBRAIC[144];
RATES[40] = (ALGEBRAIC[140] - ALGEBRAIC[142])+ 1.00000*ALGEBRAIC[147];
ALGEBRAIC[99] = CONSTANTS[141]+STATES[52];
ALGEBRAIC[136] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[81]*ALGEBRAIC[99]);
ALGEBRAIC[146] = CONSTANTS[29]*(exp((ALGEBRAIC[120]+ALGEBRAIC[130])/( CONSTANTS[144]*CONSTANTS[145])) - exp((ALGEBRAIC[123]+ALGEBRAIC[125]+ALGEBRAIC[136]+ALGEBRAIC[133])/( CONSTANTS[144]*CONSTANTS[145])));
RATES[44] = ALGEBRAIC[143] - ALGEBRAIC[146];
RATES[48] = - ALGEBRAIC[145] - ALGEBRAIC[146];
RATES[49] = ALGEBRAIC[145]+ALGEBRAIC[146];
ALGEBRAIC[93] = CONSTANTS[139]+STATES[47];
ALGEBRAIC[127] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[79]*ALGEBRAIC[93]);
ALGEBRAIC[134] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[64]*ALGEBRAIC[48]);
ALGEBRAIC[148] = CONSTANTS[30]*(exp(ALGEBRAIC[123]/( CONSTANTS[144]*CONSTANTS[145])) - exp((ALGEBRAIC[127]+ALGEBRAIC[134])/( CONSTANTS[144]*CONSTANTS[145])));
RATES[45] = (ALGEBRAIC[145]+ALGEBRAIC[146]) - ALGEBRAIC[148];
RATES[50] = ALGEBRAIC[148];
ALGEBRAIC[101] = CONSTANTS[142]+STATES[53];
ALGEBRAIC[137] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[82]*ALGEBRAIC[101]);
ALGEBRAIC[139] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[69]*ALGEBRAIC[60]);
ALGEBRAIC[149] = CONSTANTS[32]*(exp((ALGEBRAIC[135]+ALGEBRAIC[139])/( CONSTANTS[144]*CONSTANTS[145])) - exp(ALGEBRAIC[137]/( CONSTANTS[144]*CONSTANTS[145])));
RATES[51] = ALGEBRAIC[145] - ALGEBRAIC[149];
RATES[53] = ALGEBRAIC[149];
ALGEBRAIC[150] = CONSTANTS[31]*(exp((ALGEBRAIC[127]+ALGEBRAIC[125])/( CONSTANTS[144]*CONSTANTS[145])) - exp(ALGEBRAIC[113]/( CONSTANTS[144]*CONSTANTS[145])));
RATES[41] = (- ALGEBRAIC[142] - ALGEBRAIC[143])+ALGEBRAIC[150];
RATES[46] = (ALGEBRAIC[145]+ALGEBRAIC[146]) - ALGEBRAIC[150];
RATES[47] = ALGEBRAIC[148] - ALGEBRAIC[150];
ALGEBRAIC[103] = CONSTANTS[143]+STATES[54];
ALGEBRAIC[138] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[83]*ALGEBRAIC[103]);
ALGEBRAIC[151] = CONSTANTS[33]*(exp((ALGEBRAIC[136]+ALGEBRAIC[139])/( CONSTANTS[144]*CONSTANTS[145])) - exp(ALGEBRAIC[138]/( CONSTANTS[144]*CONSTANTS[145])));
RATES[52] = ALGEBRAIC[146] - ALGEBRAIC[151];
RATES[54] = ALGEBRAIC[151];
RATES[55] = - ALGEBRAIC[149] - ALGEBRAIC[151];
}
void
computeVariables(double VOI, double* CONSTANTS, double* RATES, double* STATES, double* ALGEBRAIC)
{
ALGEBRAIC[2] = CONSTANTS[96]+STATES[0];
ALGEBRAIC[28] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[34]*ALGEBRAIC[2]);
ALGEBRAIC[4] = CONSTANTS[97]+STATES[2];
ALGEBRAIC[32] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[36]*ALGEBRAIC[4]);
ALGEBRAIC[5] = CONSTANTS[99]+STATES[3];
ALGEBRAIC[34] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[37]*ALGEBRAIC[5]);
ALGEBRAIC[69] = CONSTANTS[0]*(exp((ALGEBRAIC[32]+ALGEBRAIC[28])/( CONSTANTS[144]*CONSTANTS[145])) - exp(ALGEBRAIC[34]/( CONSTANTS[144]*CONSTANTS[145])));
ALGEBRAIC[3] = CONSTANTS[98]+STATES[1];
ALGEBRAIC[30] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[35]*ALGEBRAIC[3]);
ALGEBRAIC[25] = CONSTANTS[113]+STATES[17];
ALGEBRAIC[66] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[51]*ALGEBRAIC[25]);
ALGEBRAIC[72] = CONSTANTS[1]*(exp(ALGEBRAIC[34]/( CONSTANTS[144]*CONSTANTS[145])) - exp((ALGEBRAIC[32]+ALGEBRAIC[30]+ALGEBRAIC[66])/( CONSTANTS[144]*CONSTANTS[145])));
ALGEBRAIC[6] = CONSTANTS[104]+STATES[4];
ALGEBRAIC[36] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[38]*ALGEBRAIC[6]);
ALGEBRAIC[7] = CONSTANTS[105]+STATES[5];
ALGEBRAIC[38] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[39]*ALGEBRAIC[7]);
ALGEBRAIC[75] = CONSTANTS[2]*(exp((ALGEBRAIC[36]+ALGEBRAIC[28])/( CONSTANTS[144]*CONSTANTS[145])) - exp(ALGEBRAIC[38]/( CONSTANTS[144]*CONSTANTS[145])));
ALGEBRAIC[78] = CONSTANTS[3]*(exp(ALGEBRAIC[38]/( CONSTANTS[144]*CONSTANTS[145])) - exp((ALGEBRAIC[36]+ALGEBRAIC[30]+ALGEBRAIC[66])/( CONSTANTS[144]*CONSTANTS[145])));
ALGEBRAIC[8] = CONSTANTS[101]+STATES[6];
ALGEBRAIC[40] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[40]*ALGEBRAIC[8]);
ALGEBRAIC[9] = CONSTANTS[102]+STATES[7];
ALGEBRAIC[42] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[41]*ALGEBRAIC[9]);
ALGEBRAIC[81] = CONSTANTS[4]*(exp((ALGEBRAIC[40]+ALGEBRAIC[28])/( CONSTANTS[144]*CONSTANTS[145])) - exp(ALGEBRAIC[42]/( CONSTANTS[144]*CONSTANTS[145])));
ALGEBRAIC[83] = CONSTANTS[5]*(exp(ALGEBRAIC[42]/( CONSTANTS[144]*CONSTANTS[145])) - exp((ALGEBRAIC[40]+ALGEBRAIC[30]+ALGEBRAIC[66])/( CONSTANTS[144]*CONSTANTS[145])));
ALGEBRAIC[10] = CONSTANTS[106]+STATES[8];
ALGEBRAIC[45] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[42]*ALGEBRAIC[10]);
ALGEBRAIC[13] = CONSTANTS[108]+STATES[9];
ALGEBRAIC[47] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[43]*ALGEBRAIC[13]);
ALGEBRAIC[86] = CONSTANTS[6]*(exp((ALGEBRAIC[45]+ALGEBRAIC[30])/( CONSTANTS[144]*CONSTANTS[145])) - exp(ALGEBRAIC[47]/( CONSTANTS[144]*CONSTANTS[145])));
ALGEBRAIC[19] = CONSTANTS[103]+STATES[13]+STATES[26];
ALGEBRAIC[55] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[47]*ALGEBRAIC[19]);
ALGEBRAIC[87] = CONSTANTS[9]*(exp((ALGEBRAIC[32]+ALGEBRAIC[55])/( CONSTANTS[144]*CONSTANTS[145])) - exp(ALGEBRAIC[36]/( CONSTANTS[144]*CONSTANTS[145])));
ALGEBRAIC[20] = CONSTANTS[100]+STATES[14];
ALGEBRAIC[57] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[48]*ALGEBRAIC[20]);
ALGEBRAIC[91] = CONSTANTS[10]*(exp((ALGEBRAIC[57]+ALGEBRAIC[32])/( CONSTANTS[144]*CONSTANTS[145])) - exp(ALGEBRAIC[40]/( CONSTANTS[144]*CONSTANTS[145])));
ALGEBRAIC[14] = CONSTANTS[110]+STATES[10];
ALGEBRAIC[49] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[44]*ALGEBRAIC[14]);
ALGEBRAIC[90] = CONSTANTS[7]*(exp(ALGEBRAIC[47]/( CONSTANTS[144]*CONSTANTS[145])) - exp((ALGEBRAIC[45]+ALGEBRAIC[49])/( CONSTANTS[144]*CONSTANTS[145])));
ALGEBRAIC[22] = CONSTANTS[111]+STATES[15]+STATES[45];
ALGEBRAIC[61] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[49]*ALGEBRAIC[22]);
ALGEBRAIC[23] = CONSTANTS[112]+STATES[16];
ALGEBRAIC[63] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[50]*ALGEBRAIC[23]);
ALGEBRAIC[95] = CONSTANTS[11]*(exp((ALGEBRAIC[32]+ALGEBRAIC[61])/( CONSTANTS[144]*CONSTANTS[145])) - exp(ALGEBRAIC[63]/( CONSTANTS[144]*CONSTANTS[145])));
ALGEBRAIC[18] = CONSTANTS[107]+STATES[12];
ALGEBRAIC[53] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[46]*ALGEBRAIC[18]);
ALGEBRAIC[17] = CONSTANTS[109]+STATES[11];
ALGEBRAIC[51] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[45]*ALGEBRAIC[17]);
ALGEBRAIC[94] = CONSTANTS[8]*(exp((ALGEBRAIC[45]+ALGEBRAIC[51])/( CONSTANTS[144]*CONSTANTS[145])) - exp(ALGEBRAIC[53]/( CONSTANTS[144]*CONSTANTS[145])));
ALGEBRAIC[26] = CONSTANTS[115]+STATES[19];
ALGEBRAIC[67] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[53]*ALGEBRAIC[26]);
ALGEBRAIC[35] = CONSTANTS[119]+STATES[24];
ALGEBRAIC[82] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[58]*ALGEBRAIC[35]);
ALGEBRAIC[116] = CONSTANTS[13]*(exp(ALGEBRAIC[67]/( CONSTANTS[144]*CONSTANTS[145])) - exp(ALGEBRAIC[82]/( CONSTANTS[144]*CONSTANTS[145])));
ALGEBRAIC[24] = CONSTANTS[114]+STATES[18];
ALGEBRAIC[64] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[52]*ALGEBRAIC[24]);
ALGEBRAIC[0] = (VOI<CONSTANTS[86]&&VOI>CONSTANTS[86] - CONSTANTS[88] ? CONSTANTS[90]+ CONSTANTS[148]*((VOI - CONSTANTS[86])+CONSTANTS[88]) : VOI>=CONSTANTS[86]&&VOI<CONSTANTS[86]+CONSTANTS[87] ? CONSTANTS[89]+CONSTANTS[90] : VOI<CONSTANTS[86]+CONSTANTS[88]+CONSTANTS[87]&&VOI>=CONSTANTS[86]+CONSTANTS[87] ? CONSTANTS[90]+ - CONSTANTS[148]*(((VOI - CONSTANTS[86]) - CONSTANTS[88]) - CONSTANTS[87]) : CONSTANTS[90]);
ALGEBRAIC[27] = ALGEBRAIC[0]+STATES[20];
ALGEBRAIC[70] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[54]*ALGEBRAIC[27]);
ALGEBRAIC[119] = CONSTANTS[16]*(exp((ALGEBRAIC[64]+ALGEBRAIC[70])/( CONSTANTS[144]*CONSTANTS[145])) - exp(ALGEBRAIC[67]/( CONSTANTS[144]*CONSTANTS[145])));
ALGEBRAIC[29] = CONSTANTS[116]+STATES[21];
ALGEBRAIC[73] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[55]*ALGEBRAIC[29]);
ALGEBRAIC[114] = CONSTANTS[12]*(exp(ALGEBRAIC[64]/( CONSTANTS[144]*CONSTANTS[145])) - exp(ALGEBRAIC[73]/( CONSTANTS[144]*CONSTANTS[145])));
ALGEBRAIC[31] = CONSTANTS[117]+STATES[22];
ALGEBRAIC[76] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[56]*ALGEBRAIC[31]);
ALGEBRAIC[33] = CONSTANTS[118]+STATES[23];
ALGEBRAIC[79] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[57]*ALGEBRAIC[33]);
ALGEBRAIC[118] = CONSTANTS[14]*(exp((ALGEBRAIC[73]+ALGEBRAIC[76])/( CONSTANTS[144]*CONSTANTS[145])) - exp(ALGEBRAIC[79]/( CONSTANTS[144]*CONSTANTS[145])));
ALGEBRAIC[122] = ALGEBRAIC[114]+ALGEBRAIC[119];
ALGEBRAIC[37] = CONSTANTS[120]+STATES[25];
ALGEBRAIC[84] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[59]*ALGEBRAIC[37]);
ALGEBRAIC[121] = CONSTANTS[15]*(exp((ALGEBRAIC[82]+ALGEBRAIC[76])/( CONSTANTS[144]*CONSTANTS[145])) - exp(ALGEBRAIC[84]/( CONSTANTS[144]*CONSTANTS[145])));
ALGEBRAIC[88] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[47]*ALGEBRAIC[19]);
ALGEBRAIC[39] = CONSTANTS[121]+STATES[27];
ALGEBRAIC[92] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[60]*ALGEBRAIC[39]);
ALGEBRAIC[44] = CONSTANTS[123]+STATES[29]+STATES[48];
ALGEBRAIC[98] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[62]*ALGEBRAIC[44]);
ALGEBRAIC[46] = CONSTANTS[124]+STATES[30]+STATES[49];
ALGEBRAIC[100] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[63]*ALGEBRAIC[46]);
ALGEBRAIC[50] = CONSTANTS[126]+STATES[32];
ALGEBRAIC[104] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[65]*ALGEBRAIC[50]);
ALGEBRAIC[124] = CONSTANTS[17]*(exp((ALGEBRAIC[79]+ALGEBRAIC[98])/( CONSTANTS[144]*CONSTANTS[145])) - exp((ALGEBRAIC[88]+ALGEBRAIC[92]+ALGEBRAIC[104]+ALGEBRAIC[100])/( CONSTANTS[144]*CONSTANTS[145])));
ALGEBRAIC[52] = CONSTANTS[127]+STATES[33];
ALGEBRAIC[106] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[66]*ALGEBRAIC[52]);
ALGEBRAIC[126] = CONSTANTS[18]*(exp((ALGEBRAIC[84]+ALGEBRAIC[98])/( CONSTANTS[144]*CONSTANTS[145])) - exp((ALGEBRAIC[88]+ALGEBRAIC[92]+ALGEBRAIC[106]+ALGEBRAIC[100])/( CONSTANTS[144]*CONSTANTS[145])));
ALGEBRAIC[41] = CONSTANTS[122]+STATES[28];
ALGEBRAIC[96] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[61]*ALGEBRAIC[41]);
ALGEBRAIC[48] = CONSTANTS[125]+STATES[31]+STATES[50];
ALGEBRAIC[102] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[64]*ALGEBRAIC[48]);
ALGEBRAIC[128] = CONSTANTS[19]*(exp(ALGEBRAIC[88]/( CONSTANTS[144]*CONSTANTS[145])) - exp((ALGEBRAIC[96]+ALGEBRAIC[102])/( CONSTANTS[144]*CONSTANTS[145])));
ALGEBRAIC[54] = CONSTANTS[128]+STATES[34];
ALGEBRAIC[108] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[67]*ALGEBRAIC[54]);
ALGEBRAIC[60] = CONSTANTS[130]+STATES[36]+STATES[55];
ALGEBRAIC[112] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[69]*ALGEBRAIC[60]);
ALGEBRAIC[129] = CONSTANTS[21]*(exp((ALGEBRAIC[104]+ALGEBRAIC[112])/( CONSTANTS[144]*CONSTANTS[145])) - exp(ALGEBRAIC[108]/( CONSTANTS[144]*CONSTANTS[145])));
ALGEBRAIC[131] = CONSTANTS[20]*(exp((ALGEBRAIC[96]+ALGEBRAIC[92])/( CONSTANTS[144]*CONSTANTS[145])) - exp(ALGEBRAIC[76]/( CONSTANTS[144]*CONSTANTS[145])));
ALGEBRAIC[56] = CONSTANTS[129]+STATES[35];
ALGEBRAIC[110] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[68]*ALGEBRAIC[56]);
ALGEBRAIC[132] = CONSTANTS[22]*(exp((ALGEBRAIC[106]+ALGEBRAIC[112])/( CONSTANTS[144]*CONSTANTS[145])) - exp(ALGEBRAIC[110]/( CONSTANTS[144]*CONSTANTS[145])));
ALGEBRAIC[65] = CONSTANTS[132]+STATES[38];
ALGEBRAIC[107] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[71]*ALGEBRAIC[65]);
ALGEBRAIC[80] = CONSTANTS[136]+STATES[43];
ALGEBRAIC[117] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[76]*ALGEBRAIC[80]);
ALGEBRAIC[141] = CONSTANTS[24]*(exp(ALGEBRAIC[107]/( CONSTANTS[144]*CONSTANTS[145])) - exp(ALGEBRAIC[117]/( CONSTANTS[144]*CONSTANTS[145])));
ALGEBRAIC[62] = CONSTANTS[131]+STATES[37];
ALGEBRAIC[105] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[70]*ALGEBRAIC[62]);
ALGEBRAIC[71] = CONSTANTS[133]+STATES[40];
ALGEBRAIC[111] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[73]*ALGEBRAIC[71]);
ALGEBRAIC[140] = CONSTANTS[23]*(exp(ALGEBRAIC[105]/( CONSTANTS[144]*CONSTANTS[145])) - exp(ALGEBRAIC[111]/( CONSTANTS[144]*CONSTANTS[145])));
ALGEBRAIC[1] = (VOI<CONSTANTS[91]&&VOI>CONSTANTS[91] - CONSTANTS[93] ? CONSTANTS[95]+ CONSTANTS[149]*((VOI - CONSTANTS[91])+CONSTANTS[93]) : VOI>=CONSTANTS[91]&&VOI<CONSTANTS[91]+CONSTANTS[92] ? CONSTANTS[94]+CONSTANTS[95] : VOI<CONSTANTS[91]+CONSTANTS[93]+CONSTANTS[92]&&VOI>=CONSTANTS[91]+CONSTANTS[92] ? CONSTANTS[95]+ - CONSTANTS[149]*(((VOI - CONSTANTS[91]) - CONSTANTS[93]) - CONSTANTS[92]) : CONSTANTS[95]);
ALGEBRAIC[68] = ALGEBRAIC[1]+STATES[39];
ALGEBRAIC[109] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[72]*ALGEBRAIC[68]);
ALGEBRAIC[144] = CONSTANTS[27]*(exp((ALGEBRAIC[105]+ALGEBRAIC[109])/( CONSTANTS[144]*CONSTANTS[145])) - exp(ALGEBRAIC[107]/( CONSTANTS[144]*CONSTANTS[145])));
ALGEBRAIC[74] = CONSTANTS[134]+STATES[41];
ALGEBRAIC[113] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[74]*ALGEBRAIC[74]);
ALGEBRAIC[77] = CONSTANTS[135]+STATES[42];
ALGEBRAIC[115] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[75]*ALGEBRAIC[77]);
ALGEBRAIC[142] = CONSTANTS[25]*(exp((ALGEBRAIC[111]+ALGEBRAIC[113])/( CONSTANTS[144]*CONSTANTS[145])) - exp(ALGEBRAIC[115]/( CONSTANTS[144]*CONSTANTS[145])));
ALGEBRAIC[123] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[49]*ALGEBRAIC[22]);
ALGEBRAIC[89] = CONSTANTS[138]+STATES[46];
ALGEBRAIC[125] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[78]*ALGEBRAIC[89]);
ALGEBRAIC[130] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[62]*ALGEBRAIC[44]);
ALGEBRAIC[133] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[63]*ALGEBRAIC[46]);
ALGEBRAIC[97] = CONSTANTS[140]+STATES[51];
ALGEBRAIC[135] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[80]*ALGEBRAIC[97]);
ALGEBRAIC[145] = CONSTANTS[28]*(exp((ALGEBRAIC[115]+ALGEBRAIC[130])/( CONSTANTS[144]*CONSTANTS[145])) - exp((ALGEBRAIC[123]+ALGEBRAIC[125]+ALGEBRAIC[135]+ALGEBRAIC[133])/( CONSTANTS[144]*CONSTANTS[145])));
ALGEBRAIC[85] = CONSTANTS[137]+STATES[44];
ALGEBRAIC[120] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[77]*ALGEBRAIC[85]);
ALGEBRAIC[143] = CONSTANTS[26]*(exp((ALGEBRAIC[117]+ALGEBRAIC[113])/( CONSTANTS[144]*CONSTANTS[145])) - exp(ALGEBRAIC[120]/( CONSTANTS[144]*CONSTANTS[145])));
ALGEBRAIC[147] = (ALGEBRAIC[142])+ALGEBRAIC[144];
ALGEBRAIC[99] = CONSTANTS[141]+STATES[52];
ALGEBRAIC[136] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[81]*ALGEBRAIC[99]);
ALGEBRAIC[146] = CONSTANTS[29]*(exp((ALGEBRAIC[120]+ALGEBRAIC[130])/( CONSTANTS[144]*CONSTANTS[145])) - exp((ALGEBRAIC[123]+ALGEBRAIC[125]+ALGEBRAIC[136]+ALGEBRAIC[133])/( CONSTANTS[144]*CONSTANTS[145])));
ALGEBRAIC[93] = CONSTANTS[139]+STATES[47];
ALGEBRAIC[127] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[79]*ALGEBRAIC[93]);
ALGEBRAIC[134] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[64]*ALGEBRAIC[48]);
ALGEBRAIC[148] = CONSTANTS[30]*(exp(ALGEBRAIC[123]/( CONSTANTS[144]*CONSTANTS[145])) - exp((ALGEBRAIC[127]+ALGEBRAIC[134])/( CONSTANTS[144]*CONSTANTS[145])));
ALGEBRAIC[101] = CONSTANTS[142]+STATES[53];
ALGEBRAIC[137] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[82]*ALGEBRAIC[101]);
ALGEBRAIC[139] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[69]*ALGEBRAIC[60]);
ALGEBRAIC[149] = CONSTANTS[32]*(exp((ALGEBRAIC[135]+ALGEBRAIC[139])/( CONSTANTS[144]*CONSTANTS[145])) - exp(ALGEBRAIC[137]/( CONSTANTS[144]*CONSTANTS[145])));
ALGEBRAIC[150] = CONSTANTS[31]*(exp((ALGEBRAIC[127]+ALGEBRAIC[125])/( CONSTANTS[144]*CONSTANTS[145])) - exp(ALGEBRAIC[113]/( CONSTANTS[144]*CONSTANTS[145])));
ALGEBRAIC[103] = CONSTANTS[143]+STATES[54];
ALGEBRAIC[138] = CONSTANTS[144]*CONSTANTS[145]*log( CONSTANTS[83]*ALGEBRAIC[103]);
ALGEBRAIC[151] = CONSTANTS[33]*(exp((ALGEBRAIC[136]+ALGEBRAIC[139])/( CONSTANTS[144]*CONSTANTS[145])) - exp(ALGEBRAIC[138]/( CONSTANTS[144]*CONSTANTS[145])));
ALGEBRAIC[11] = ALGEBRAIC[2]+ALGEBRAIC[5]+ALGEBRAIC[9]+ALGEBRAIC[7];
ALGEBRAIC[12] = ALGEBRAIC[4]+ALGEBRAIC[5]+ALGEBRAIC[8]+ALGEBRAIC[9]+ALGEBRAIC[6]+ALGEBRAIC[7];
ALGEBRAIC[15] = ALGEBRAIC[3]+ALGEBRAIC[13]+STATES[10];
ALGEBRAIC[16] = ALGEBRAIC[3]+ALGEBRAIC[13]+ALGEBRAIC[14]+ALGEBRAIC[2]+ALGEBRAIC[5]+ALGEBRAIC[7]+ALGEBRAIC[9];
ALGEBRAIC[21] = ALGEBRAIC[19]+ALGEBRAIC[6]+ALGEBRAIC[7];
ALGEBRAIC[43] = ALGEBRAIC[31]+ALGEBRAIC[33]+ALGEBRAIC[37]+ALGEBRAIC[19]+ALGEBRAIC[41];
ALGEBRAIC[58] = ALGEBRAIC[26]+ALGEBRAIC[27]+ALGEBRAIC[37]+ALGEBRAIC[35]+ALGEBRAIC[52]+ALGEBRAIC[56];
ALGEBRAIC[59] = ALGEBRAIC[24]+ALGEBRAIC[26]+ALGEBRAIC[29]+ALGEBRAIC[33]+ALGEBRAIC[35]+ALGEBRAIC[37]+ALGEBRAIC[50]+ALGEBRAIC[52]+ALGEBRAIC[54]+ALGEBRAIC[56];
}