  Using the variational method and supersymmetric quantum mechanic we calculate
in a approximate way eigenvalues, eigenfunctions and wave functions at origin
of Cornell potential. We compare results with numerical solutions for heavy
quarkonia $c\bar{c}$, $b\bar{b}$ y $b\bar{c}$.
