Suszeptibility fixed

First version of what seems to be a properly working programm to get
x,xb,y,z1 and z2, to calculate the susceptibility and its poles.
plot_chi.py adapted to show the result.

main fix: Q set to (pi,pi).