function shiyan(r)
r=r;
t=0:0.01pi:2pi;
z=rexp(it);
w=z+1./z;
plot(w);
title('w=z+1/z');
axis equal
shiyan.m
function shiyan(r)
t=0:0.01*pi:2*pi;
z=r*exp(i*t);
w=z+1./z;
plot(w);
title('w=z+1/z');
axis equal
end
主函数
shiyan(3);