dsolve(diff(y(x),x) = -1/2*x^2*(a*x-2*(a*(a*x^4+8*y(x)))^(1/2)),y(x))