We apply matrix Numerov’s method to obtain the radial wave functions; from these wave functions we calculate the root mean square radius rms and β coefficients of bottomonium bb . The
obtained results have implications for decay constants, decay widths and differential cross sections of heavy mesons.