I have been attempting to solve the following central potential using schrodinger equation.$$V(r)=a_{1}r^{2}+\frac{a_{2}}{r^{2}}+a_{3}$$
The radial equation becomes, $$\Bigg[\frac{-\hbar^{2}}{2m}\dfrac{d^{2}}{dr^{2}}+V(r)+\frac{l(l+1)}{r^{2}}\Bigg]u(r)=Eu(r)$$ I have used laplace transform to solve this problem to get the eigen values. Defining $$u(y)=y^{-\frac{\nu}{2}}\phi(y)$$ I get the following differential equation,$$\Bigg[y\dfrac{d^{2}}{dy^{2}}-\bigg(\nu-\dfrac{1}{2}\bigg)\dfrac{d}{dy}-\dfrac{1}{4}\bigg(\mu^{2}y-\epsilon^{2}\bigg)\Bigg]\phi(y)=0$$ where, $$\nu(\nu+1)=\dfrac{2ma_{2}}{\hbar^{2}}+l(l+1)$$ $$\mu^{2}=\dfrac{2ma_{1}}{\hbar^{2}}$$ $$\epsilon^{2}=\dfrac{2mE}{\hbar^{2}}$$ Now after doing laplace transform I have the following result $$\bigg(t^{2}-\dfrac{\mu^{2}}{4}\bigg)\dfrac{df(t)}{dt}+\bigg[\bigg(\nu+\dfrac{3}{2}\bigg)t-\dfrac{\epsilon^{2}}{4}\bigg]f(t)$$ The soluion in laplace space is, $$f(t)=N\bigg(t+\dfrac{\mu}{2}\bigg)^{-\dfrac{\epsilon^{2}}{4\mu} -\dfrac{1}{2}\bigg(\nu+\dfrac{3}{2}\bigg)}\bigg(t-\frac{\mu}{2}\bigg)^{\dfrac{\epsilon^{2}}{4\mu}+\dfrac{1}{2}\bigg(\nu+\dfrac{3}{2}\Bigg)}$$ To get a finite wave functions the solutions have to be a finite polynomial $$\dfrac{\epsilon^{2}}{4\mu}-\dfrac{1}{2}\bigg(\nu + \dfrac{3}{2}\bigg)=n$$ where, $$n=0,1,2,3...$$ The energy eighen values are ,$$E_{nl}=\sqrt{\dfrac{8\hbar^{2}a_{1}}{m}}\bigg(\dfrac{1}{2}+\dfrac{1}{4}\sqrt{1+4l(l+1)+\dfrac{8ma_{2}}{\hbar^{2}}}\bigg)$$
I have understood everything except the condition to be imposed on the angular quantum number $l$ for a given value of radial quantum number $n$. For instance, suppose $n$=3 the minimum value of $l$ = 0, but the what is the maximum value of $l$ ? What is the condition on $l$ for a given $n$ ? I have followed the following research article
*Exact solutions of the Schrödinger equation via Laplace transform approach: pseudoharmonic potential and Mie-type potentialshttps://link.springer.com/journal/10910 *
