In this study, a thermodynamic model has been used to determine the solubility of carbon dioxide in an aqueous solution which is the combination of methyldiethanolamine (MDEA) and aminoethylethanolamine (AEEA). The physical equilibriums have been considered between the liquid and vapor phases and chemical equilibrium in the liquid phase. The SAFT-HR equation of state has been used to specify the fugacity coefficients of the components in the vapor phase. The liquid phase is considered as an electrolyte solution besides; the extended UNIQUAC has been applied to figure out the activity coefficients. The bubble point calculation has been used in this research. This method includes two main loops. Calculations related to chemical equilibrium are performed in the interior loop and the ones associated with phase equilibrium are done in the exterior loop. The solubility of carbon dioxide has been predicted by the optimized parameters of the model in the temperature range of 308.2–368.2 K. It has been calculated that the absolute average relative deviations of the model are 16.65, 19.33, 28.91 and 19.99 in the calculation of partial pressure of carbon dioxide in various loadings at the temperatures of 308.2, 328.2, 343.2 and 368.2 K.