Fullerenes are reactive as dienophiles in Diels–Alder reactions. Their distinctive molecular shape and properties result in interesting and sometimes elusive reaction patterns. Herein, to contribute to the understanding of fullerene reactivity, we evaluate the energies of reactions for Diels–Alder cycloadditions of C60, C70, and IC60MA with anthracene (Ant), by means of DFT computational analysis in vacuum and solution. The methods used showed little differentiation between the reactivity of the different fullerenes. The C70-Ant adducts where addition takes place near the edge of the fullerene were found to be the most stable regioisomers. For the IC60MA-Ant adducts, the calculated energies of reaction increase in the order: equatorial > trans-3 > trans-2 ≈ trans-4 ≈ trans-1 > cis-3 > cis-2. The change in the functional suggests the existence of stabilizing dispersive interactions between the surface of the fullerene and the addends. HOMA (harmonic oscillator model of aromaticity) analysis indicated an increase in aromaticity in the fullerene hexagons adjacent to the bonded addend. This increase is bigger in the rings of bisadduct isomers that are simultaneously adjacent to both addends, which helps explain the extra stability of the equatorial isomers. Solvation by m-xylene decreases the exothermicity of the reactions studied but has little distinguishing effect on the possible isomers. Thermal corrections reduce the exothermicity of the reactions by ~10 kJ∙mol−1.