There exists several algorithms to hadle adaptive time stepping, We shall start with BDF.

The methodology will be :

1. support BDF variable time stepping using automated backward differences
2. implement variable time stepping 1st and 2nd order possibly higher order as well since we will automate coefficient computations
3. implement 1st 2nd variable order variable time stepping strategies
4. implement variable order variable time stepping strategies up to order 4