The Fermi surface (FS) in paramagnetic ErGa 3 is estimated by various band structure methods: linear muffin-tin orbital (LMTO) in the atomic sphere approximation (ASA) and three new codes contained full potential (FP) instead of ASA, namely FP-LMTO, FP linear augmented plane wave (FLAPW), and FP local orbitals (FPLO) methods. Three dimensional (3D) electron-positron (e-p) momentum densities ρ(p) are reconstructed from two dimensional (2D) angular correlation of annihilation radiation (ACAR) spectra measured for paramagnetic ErGa 3 . Reconstructed densities in the extended p-space, folded into the reduced k-space, are compared with the theoretical results. Unfortunately, none of these modern FP codes is able to give satisfying description of the experimental data that are in perfect agreement with previous LMTO-ASA results.The rare-earth (RE) compounds constitute an interesting subject of investigations due to their magnetic properties [1]. Lanthanides model assumes that the magnetic coupling between the 4f electrons is carried out indirectly via polarization of conducting electrons, according to the RudermanKittel-Kasuya-Yosida (RKKY) theory [2] and that magnetic properties are controlled by the FS topology (so-called nesting -flat FS sheets) [3]. The first direct observation of such an FS feature was possible thanks to the study of 2D ACAR spectra [4][5][6][7][8][9][10]. In this unique technique one measures line projections of the e-p pairs in the extended momentum space ρ(p). From such spectra for different directions of projections one can reconstruct ρ(p) in 3D momentum space applying tomography techniques [11]. Next, using the LCW folding [12], i.e., a conversion from the extended p into reduced k space, one may reproduce the FS.In the case of ErGa 3 2D-ACAR experiment was carried out in the paramagnetic state at temperature of 60 K and in the magnetic field 1.8 T (for more details see [9]). Four projections were measured and next deconvoluted employing the Van Citter-Gerhardt algorithm [13]. The e-p momentum densities ρ(p) were reconstructed by Cormack's method (CM) [14] and the lattice harmonics expansion (LHE) algorithm [15]. Finally, the 3D-LCW transformation [12] was applied to get ρ(k) (occupation numbers "seen" by positrons). These results were compared with band structure calculations performed via four different methods: LMTO-ASA [16]