A series of computational experiments performed with various methods belonging to wave-function and density functional theories approaches the issue of bonding regime and exchange coupling in the title compounds. Gd 2 @C 80 is computed with a very weak exchange coupling, the sign depending on the method, while Gd 2 @C 79 N has resulted with a strong coupling and ferromagnetic ground state, irrespective of the computational approach. The multi-configuration calculation and broken symmetry estimation are yielding closely coincident coupling constants, of about J ∼ 400 cm −1 . No experimental estimation exists, but the ferromagnetic ground state of Gd 2 @C 79 N is confirmed from paramagnetic resonance data. The different behaviour is due to particularities of electron accommodation in the orbital scheme. The exchange effects localised on atom lead to preference for parallel alignment of the electrons placed in the 4f and 5d lanthanide shells, determining also a ferromagnetic inter-centre coupling. The structural insight is completed with a ligand field analysis of the density functional theory results in the context of frozen density embedding. The energy decomposition analysis of bonding effects is also discussed. Finally, with the help of home-made codes (named Xatom+Xsphere), a model for the atom encapsulated in a cage is designed, the exemplified numeric experiments showing relevance for the considered endohedral metallo-fullerene issues.