z+xz'=(x+xz)/(x-xz)
xz'=(1+z)/(1-z)-z
xz'=(1+z-z+z^2)/(1-z)
x*dz/dx=(1+z^2)/(1-z)
(1-z)dz / (1+z^2)=dx/x
Дальше уже просто надо интегралы посчитать... там получается
arctgz-1/2 ln(1+z^2)=ln|x|+C
z=y/x
arctg(y/x)-1/2 ln(1+(y/x)^2)=ln|x|+C
arctg(y/x)-1/2 ln(x^2+y^2) + 1/2ln x^2 = ln|x|+C
arctg(y/x)-1/2 ln(x^2+y^2) = ln|x|+C-ln|x|
arctg(y/x)-1/2 ln(x^2+y^2) = C
Вроде так получается