A method is described to calculate molar solution energies from solubility data of binary low-volatility organic compound + compressed carbon dioxide system. It has been applied to equilibrium measurements of the two anthraquinone derivatives 1,4-bis-(methylamino)-9,10-anthraquinone and 1,4-bis-(1-methylethylamino) -9,10-anthraquinone in CO2 within p,T-range of 6.79-18.80 MPa and 299 345 K, respectively. Since solubility data do not exceed 0.135 g dm(-3), the fluid phase can be regarded as infinitely dilute. The data have been analyzed in way that from the two derivative plots logarithm of the solubility vs. the reciprocal temperature at constant pressure as well as constant density, two expressions for the molar solution energy, Delta(sol)E(p) and Delta(sol)E(rho), have been derived and interpreted. We summarize our findings in four main statements: Delta(sol)E(rho) remains almost constant within significant density range, here approximately between 400 and 800 kg m(-3); Delta(sol)E(p) is drastically pressure (and density) dependent; both expressions are connected with each other via the isobaric expansion coefficient alpha(p) of the solvent; and finally for these two systems the Delta(sol)E(p) and Delta(sol)E(p) values, respectively, are quite similar, whereas the solubilities differ by about factor of ten.