dsolve(diff(y(x),x) = (a*y(x)^2+b*x^2)^2*x/a^(5/2)/y(x),y(x))