It is clear that the problem in the potential ##V(r)+W(r)*(l,s)##, but without the (l,s) term: ##V(r)+W(r)##, simply reduces to solving the radial Schrodinger equation. And I understand how to find energy levels in such a task. But what to do when ls interaction also appears in the delta layer?