dsolve(diff(diff(y(x),x),x) = -sin(x)/cos(x)*diff(y(x),x)-1/4*(2*x^2+x^2*sin(x)^2-24*cos(x)^2)/x^2/cos(x)^2*y(x)+cos(x)^(1/2),y(x))