Should provide a user-facing interface for adapting step size

For methods that support adapting the step size, we should let the user pass the next value of <dt> to support the use case of obtaining the solution at a set of grid points.

Fortran allows this already (though only because <dt> is exposed in the state data structure), but Python and the interpreter don't.