clear
close

fig=figure('Units','centimeters','Position',[0.0 0.0 21.0 26.0]...
,'PaperUnits','centimeters','PaperPosition',[0.0 0.0 21.0 26.0]);

load wn_dat;
size(wn_dat)
load qxrr_dat;
load qxri_dat;

subplot(2,1,1)
plot(wn_dat,qxrr_dat,'-')
xlabel('Wavenumber')
ylabel('Real(Qx,r)')
title('Real(Qx,r)')

subplot(2,1,2)
plot(wn_dat,qxri_dat,'-')
xlabel('Wavenumber')
ylabel('Imag(Qx,r)')
title('Imag(Qx,r)')

print integ1.ps




