\left(x\frac{d^2}{dx^2}+(k+1-x)\frac{d}{dx}+(n-k)\right)L_n^k(x)=0 \tag{1} \Psi(r, \theta, \phi)=N_{l,n} \left(\frac{\rho}{n}\right)^l e^{-\frac{\rho}{n}} L_{n+l}^{2l+1}\left(\frac{\rho}{n}\right) \times M_{l,m} P_l^{|m|}(\cos\theta)e^{im\phi}\\ \ \\ \rho=\frac{Z}{a_0}r,\; a_0=\frac{4\pi\varepsilon_0\hbar^2}{\mu e^2},\\ N_{l,n}=-\sqrt{\frac{4(n-l-1)!}{n^4[(n+l)!]^3}}\left(\frac{Z}{a_0}\right)^{\fr