dsolve(diff(y(x),x) = (b*x^3+c^2*a^(1/2)-2*c*b*x^2*a^(1/2)+2*c*y(x)^2*a^(3/2)+b^2*x^4*a^(1/2)-2*y(x)^2*a^(3/2)*b*x^2+a^(5/2)*y(x)^4)/a/x^2/y(x),y(x))