Lots of debugging print statements i.e. a big mess
Cleaned up numerical Gaussian notebook a bit
Confirmed that numerical version basically works, although hard to tell for under-resolved or over-resolved versions