clear
close

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

plot(wn_dat,qxrr_dat,'-',wn_dat,qxri_dat,'--')
xlabel('Wavenumber')
ylabel('Integrands of Qx,r')
title('Solid Lines: Real Parts,  Dashed Lines: Imaginary Parts')

print qxr.ps




