% Code for plot of WKB wave in string of exponentially-growing density
%
% Barnett 2/10/04

clear;
a = .1;           % set the exponential growth rate (small)
x = 0:.02:50;     % vector of x values to plot
% do the plot
figure; plot(x, real(exp(-a*x/2 +i*(exp(a*x)-1)/a)),'-');
xlabel('x'); ylabel('Re u(x)'); title('WKB soln, exp-growing density');
