A central ΛN potential containing spin and charge dependences is constructed on the basis of an analysis of fundamental hypernuclear data (Λp scattering, the binding energies of the ground and excited states of hypernuclei). The lightest hypernuclei ( Λ 3 H, Λ 4 H, Λ 4 He) are calculated as a system of three or four particles on the basis of the highly efficient variational method. Heavier hypernuclei are considered in terms of a twobody model (Λ + undeformed nucleon core) corrected by comparison with direct 3–5-particle calculations. By comparing the results of the calculations of single and double Λ hypernuclei and taking the considerations of meson theory into account, it is possible to estimate the binding energies of charmed and other (b- and t-quark) supernuclei.