The B-spline basis-set method is applied to determining the rovibrational eigenspectrum of diatomic molecules. Particular attention is paid to a challenging numerical task of an accurate and efficient description of the vibrational levels near the dissociation limit (halo-state and Feshbach molecules). Advantages of using B-splines are highlighted by comparing the performance of the method with...