%
%  code10.m -- stationary phase (djm - 10 oct 01)
%

k = -6:0.01:6;
z = -9;

subplot(3,1,1)
plot(k,real(exp(i*(z*k + k.^3/3))),'r')
hold on;
plot(k,imag(exp(i*(z*k + k.^3/3))),'b--')
set(gca,'xtick',[-6:3:6])
title('\bf Airy exponent (z=-9)')
ylabel('\bf Re (red) & Im (blue)')
xlabel('\bf k-axis')