Numerical Differentiation Formula with piecewise constant step for efficient time-domain simulation of stiff systems.

The paper is dedicated to an idea of a variable step and order numerical differentiation method where a step prediction is performed with a considerably improved accuracy thanks to a piecewise constant method of step selection.