clear all
close all
fig1=figure(1)
grid on
hold on
axis([-1 1 -1 1])

plot(0,0)
axis('square')
fig2=figure(2)
grid on
hold on
axis([0 20 -20 20])

plot(0,0)
axis('square')
r=(0:.08:80)';
x=[-5 -2 -1 -.5 -.2 0 .2 .5 1 2 5];
for u=1:11
   z=r+j*x(u);
   Gamma=(z-1)./(z+1);
   figure(1)
   plot(real(Gamma),imag(Gamma),'b')
   figure(2)
   plot(real(z),imag(z),'b')
end
x=[0 .2 .5 1 2 5];
r=(-40:.08:40)';
for u=1:6
  z=x(u)+j*r;
  figure(1)
  Gamma=(z-1)./(z+1);
  plot(real(Gamma),imag(Gamma),'r')
   figure(2)
   plot(real(z),imag(z),'r')
end
% Finalize the plots with legends and titles
figure(1)
title('Gamma Mapping')
xlabel('Real Part')
ylabel('Imaginary Part')
%legend('Initial', 'Modified')