function circle(z, r)
close all 
t = 0:.001:2*pi ;
plot(real(z)+r*cos(t), imag(z)+r*sin(t), "m");
axis equal
end

