  The Non-Symmetrized Hyperspherical Harmonics method (NSHH) is introduced in
the hypernuclear sector and benchmarked with three different ab-initio methods,
namely the Auxiliary Field Diffusion Monte Carlo method, the Faddeev-Yakubovsky
approach and the Gaussian Expansion Method. Binding energies and hyperon
separation energies of three- to five-body hypernuclei are calculated by
employing the two-body $\Lambda$N component of the phenomenological
Bodmer-Usmani potential, and a hyperon-nucleon interaction simulating the
scattering phase shifts given by NSC97f. The range of applicability of the NSHH
method is briefly discussed.
