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