Add adaptive and multistep methods

[SuperFluffy]

At the moment, we only support fixed timestep methods. Adaptive and multistep methods are useful to minimize numerical errors and solve stiff problems, respectively.

This will require thinking about whether the integrator or the stepper should control the timestep.