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