% Create a smith chart. With hold on/hold off your data can be plotted on top of it.
% Author: Marcin Mleczko <address@hidden>

%retain hold status
hold_status=ishold();
hold on;
linecolor=[0.8 0.8 0.8];
plot(exp(1i.*[1:1:359]./360.*2.*pi),'color', linecolor)
z=0.2+[-10:0.1:10].*i;
plot((z-1)./(z+1),'color', linecolor)
z=0.5+[-10:0.1:10].*i;
plot((z-1)./(z+1),'color', linecolor)
z=1+[-10:0.1:10].*i;
plot((z-1)./(z+1),'color', linecolor)
z=5+[-10:0.1:10].*i;
plot((z-1)./(z+1),'color', linecolor)
z=2+[-10:0.1:10].*i;
plot((z-1)./(z+1),'color', linecolor)
z=10+[-10:0.1:10].*i;
plot((z-1)./(z+1),'color', linecolor)

z=[0:0.1:10]-0.2i;
plot((z-1)./(z+1),'color', linecolor)
z=[0:0.1:10]-0.5i;
plot((z-1)./(z+1),'color', linecolor)
z=[0:0.1:10]-1i;
plot((z-1)./(z+1),'color', linecolor)
z=[0:0.1:10]-5i;
plot((z-1)./(z+1),'color', linecolor)
z=[0:0.1:10]+0.0000001i;
plot((z-1)./(z+1),'color', linecolor)
z=[0:0.1:10]+0.2i;
plot((z-1)./(z+1),'color', linecolor)
z=[0:0.1:10]+0.5i;
plot((z-1)./(z+1),'color', linecolor)
z=[0:0.1:10]+1i;
plot((z-1)./(z+1),'color', linecolor)
z=[0:0.1:10]+5i;
plot((z-1)./(z+1),'color', linecolor)
%ax=gca;
%ax.Box='off'
%legend(' ',' ',' ',' ',' ',' ',' ',' ',' ',' ',' ',' ',' ',' ',' ',' ')
%legend( off);
axis off
if hold_status
  hold on;
else
  hold off;
end