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