Fix some math typos
Fix a severe bug inside (old style) quadratures

This bug was responsible for the fact that some inner products
did not yield the expected norm of a braket anymore:

  <phi_k | phi_k> = 1  for all epsilon
Implement handling of higher order state phi_k and global phase difference