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