dsolve(diff(diff(y(x),x),x) = -1/sin(x)*cos(x)*diff(y(x),x)-(v*(v+1)*sin(x)^2-n^2)/sin(x)^2*y(x),y(x))