Add OFMD Energies notebook
XX
Reuse ne
Update a test
Add back the old Gaussian
Use the new gaussian
Expose energy_eps
Plot new energies
Add ofdft_fe
Use it from ofmd.f90
Converge SCF to 1e-9, correct final energies
Use the same free_energy_min()

And expose Nq, quad_type and energy_eps.