clear
close

load wn_dat;
size(wn_dat)
load qxzr_dat;
load qxzi_dat;

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

print qxz.ps




