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