clear
close

load wn_dat;
size(wn_dat)
load qxr_dat;
load qxi_dat;

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

print qx.ps




