dsolve(diff(diff(y(x),x),x) = -(-(a^2*b^2-(a+1)^2)*sin(x)^2-a*(a+1)*b*sin(2*x)-a*(a-1))/sin(x)^2*y(x),y(x))