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. 2