%
%  w02backSp.m -- djm -- 15 jan 06
%

load BSdata2

figure(3);  clf ; hold on

plot((data1(:,1)),data1(:,1) .* (data1(:,3)),'bx')

title('\bf back-scattered wave amplitude')
xlabel('\bf \omega')
ylabel('\bf \omega |b_R|    (N = 2^{12} points)')

a = -0.5;  sqa = sqrt(1-a^2);  izs = acosh(-1/a);
dom = 0.01;  om = dom:dom:data1(end,1);
wave = (a/2)*i*exp(i*pi*om).*cos(2*(a+pi/2)*om);
wave = wave + i*(pi/2)*exp(i*pi*om).*exp(-2*(izs - sqa)*om).*om;
%     - (i/4)*sqa*sqrt(pi/2/sqa)* sqrt(om).*exp(-2*sinv*om);
plot(om,abs(wave),'g');