WIP: Polynomial propagators with variable time step

Polynomial propagators (except for the two-step symmetric) allow variable time step that can be adapted so that it is as large as possible while the local truncation error is below a pre-specified tolerance.