% Plot the scattering transforms computed 
% (1) from Lippmann-Schwinger type equation by tLS_comp.m, and 
% (2) from the boundary integral equation by tBIE_comp.m
%
% Samuli Siltanen May 2012

% Plot parameters
lwidth = .5;
fsize  = 16;

% Create plot objects
load data/tBIE tBIE 
load data/BIEindex BIEindex
load data/tLS tLS
MAXk = max(kvec);

% Plot
figure(1)
clf
plot(kvec(1:BIEindex),real(tBIE(1:BIEindex)),'m','linewidth',4*lwidth)
hold on
plot(kvec,real(tLS),'k','linewidth',2*lwidth)
axis([0 round(MAXk) -3 4])
set(gca,'xtick',[0:5:round(MAXk)],'fontsize',fsize)
box off
legend('tBIE','tLS')

% Save image to file
%print -dpng scatLSBIE.png
