DSolve[Derivative[2][y][x] == Sqrt[Cos[x]] - (Sec[x]^2*(2*x^2 - 24*Cos[x]^2 + x^2*Sin[x]^2)*y[x])/(4*x^2) - Tan[x]*Derivative[1][y][x],y[x],x]