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