Work on Squeezed PFO spectra for Mariano.
Before I optimise the code a bit...
Optimisations of Sturm sequence; now uses temporary variable, not vector.
Hamiltonian extracted from global scope to its own function.
Loops over count for lots of data...
Edits; preparing to integrate me Z, woo!
Numerical integration for partition function Z! Woo! Numerically odd
Wo hoo! Working numeric integration of partition function Z. I love Physics.
Basic sampling of thetas a fn of pressure. :^) Probs bugs galore.
Produces DoSes as a fn(P). They seem fairly sensible!
More physical form of the pressure potential: +P.sin(theta)^2