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