We carefully study the energy eigen-values and splitting of heavy quarkonia
as there exist l/r(3) and delta(3)((r) over right arrow) singular terms in
the potential which make a direct numerical solution of the Schrodinger eq
uation impossible. We compare the results obtained in terms of perturbation
and variational methods with various treatments.