  An efficient procedure of orthonormalisation of the Bargmann-Moshinsky (BM)
basis is examined using analytical formulas of the overlap integrals of the BM
basis. Calculations of components of the quadrupole operator between the both
BM and the orthonormalised bases needed for construction of the nuclear models
are tested. The proposed procedure is also implemented as the Fortran program.
