%
% file:      	complex_view_3.m, (c) Matthew Roughan, Sun Aug 22 2004
%
x = -1.5:0.01:1.5;
y = -1.5:0.01:1.5;
[X,Y] = meshgrid(x,y);
Z = X + sqrt(-1) * Y;
a = 0.9;
H = abs((1-a)./(1 - a*Z.^(-1)));
circ_z = circ_x + sqrt(-1)*circ_y;
circ_h = abs((1-a)./(1 - a*circ_z.^(-1)));

figure(3)
hold off
surf(X,Y,log10(H));
shading interp
hold on
plot3(circ_x, circ_y, log10(circ_h)+0.01, 'r', 'linewidth', 3);

set(gca,'xlim',[-1.5 1.5],'ylim',[-1.5 1.5]);
set(gca, 'linewidth', 3, 'fontsize', 18);
xlabel('real');
ylabel('imaginary');
zlabel('log_{10} Y(z)');
az = -10; el = 30;
view(az,el);

set(gcf, 'PaperUnits', 'centimeters', 'PaperPosition', [0 0 18 10])
print('-depsc', 'Plots/complex_view_2.eps');
