  In this paper, we calculate the mass spectra of heavy quarkonium by using
matrix Numerov's method to make the predictions of F and G states for further
experiments. The method gives a very reasonable result which is in a good
agreement with other methods and with recently published theoretical data. From
the yielded wave functions we calculate the root mean square radius r_ms and
\b{eta} coefficients of heavy quarkonium
