%% Plot Airy Functions
% Plot the Airy Functions, Ai(x) and Bi(x), over the interval [-10 2] using fplot.

syms x
fplot(airy(x), [-10 2])
hold on
fplot(airy(2,x), [-10 2])
legend('Ai(x)','Bi(x)','Location','Best')
title('Airy functions Ai(x) and Bi(x)')
grid on

%
% Plot the absolute value of Ai(z) over the complex plane.

syms y
z = x + 1i*y;
figure(2)
fsurf(abs(airy(z)))
title('|Ai(z)|')
a = gca;
a.ZLim = [0 10];
caxis([0 10])

% Copyright 2012 The MathWorks, Inc.