Implement error-based adaptive timestepping

[dwillcox]

Adapt the timestep to the current error estimate, perhaps taking a min of the adaptive timestep and the timestep determined from the neutrino potential as in #57

We could use, e.g. RK45 or Richardson extrapolation.