NU method is based on solving hypergeometric second order differential equations with special orthogonal functions [1]}. Any second order differential equations which have three singular points can be converted to a hypergeometric type differential equation by a suitable change of variables. In spherical coordinates, for a given potential Schrödinger equation is reduced to a generalized hypergeometric-type equation with appropriate coordinate transformation. Then this reduced equation can be solved systematically in order to obtain eigenvalues and eigenfunctions. The main equation of NU method is given in the following form [2]};
\(\psi ^{\prime \prime }(z)+\frac{\widetilde{\tau }(z)}{\sigma (z)}\psi ^{\prime }(z)+\frac{\widetilde{\sigma }(z)}{\sigma ^{2}(z)}\psi (z)=0\) 
