x=[-20:.1:20];
r=.5;
%x=-1;
%r=[0:.1:20];
z=r+j*x;
gamma=(z-1)./(z+1);
plot(gamma);