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