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