Run the solver and compare energy decay at the ten fixed times

Run the solver and compare energy decay at the ten fixed times

Integrate the frozen initial field and record total kinetic energy at each of the ten frozen times. Compute the relative difference from the analytic value at every time. Deliverable: A CSV with time, simulated_energy, analytic_energy, and relative_error, plus the solver and time step used.