((r^2-r) d^2f/dr^2+ df/dr - (rho^2 r^3 /(r-1)+(l+1)l-3/r)f)/r^2, where f= (r^2/(r-1) e^(r-1))^rho g(r)