Fix FSAL bug in Bogacki-Shampine method

Previously, taking a step from a pervious step with a different stepsize
would result in an incorrect evolution due to a bug in the FSAL
implementation.
Add OdeIntegrator interface
Implement ConstantDtIntegrator