  In this work, we used a tool of conventional Nikiforov-Uvarov method to
determine bound state solution of Schrodinger equation with quantum interaction
potential called Hulthen-Yukawa inversely quadratic potential (HYIQP). We
obtained the energy eigen values and the total wave function . We employed
Hellmann-Feynmann Theorem (HFT) to compute expectation values for four
different diatomic molecules: Hydrogen molecule (H2), Lithium hydride molecule
(LiH), Hydrogen Chloride molecule (HCl) and Carbon(II)Oxide molecule. The
resulting energy equation reduces to three well known potentials which are:
Hulthen potential, Yukawa potential and inversely quadratic potential. We
obtained the numerical bound state energies of the expectation values by
implementing Matlab algorithm using experimentally determined spectroscopic
constant for the different diatomic molecules. We developed a mathematica
programming to obtain wave function and probability density plots for different
orbital angular quantum number.
