Change Integrator interface so that it takes time step as an argumant
Currently time step is a property of an integrator. In case of varying time step, this parameter needs to be changed before every call to integrator, which is ugly. Let's make it another argument of integrator call.