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