  Phenomenological potentials describe the quarkonium systems like Charmonia,
Bottomonia and $B_c$ Meson. They give a good accuracy for the mass spectra. In
the present work we extend one of our previous works in the central case by
adding spin dependent terms to allow for relativistic corrections. By using non
central terms, we get better accuracy than other previous theoretical
calculation. In the present work, the mass spectra of the bound states of heavy
quarks are studied within the framework of the non relativistic Schrodinger s
equation. First, we solve Schrodinger's equation by Nikiforov Uvarov (NU)
method. The energy eigenvalues are presented using a non central potential. The
results obtained are in good agreement with the experiment and are better than
previous theoretical estimates.
