\[(\ell+1)rf'\left(\epsilon,\ell;r\right)=\left((\ell+1)^{2}-r\right)f\left(%
\epsilon,\ell;r\right)-\left(1+(\ell+1)^{2}\epsilon\right)rf\left(\epsilon,%
\ell+1;r\right),\]