BDF Methods for DAEs in Multi-body Dynamics: Shortcomings and Improvements

If backward differentiation formula (BDF) methods are used to solve sets of differential-algebraic equations (DAEs), problems may arise if the step size (h), becomes very small. The step size may decrease for a variety of reasons, but the most difficult problem to solve involves a sudden change, or (near) discontinuity in any of the system variables.