Some code cleanup
Add symbolic BHZ
Small bugfixes
Add observables
Minor modifications since last version
Various updates
E1 shift; updates to plot
New method for 2nd order derivative plus few minor updates
More observables
To Dropbox
Test on cluster with magnetic field