Numerical solution of dynamical systems with equality state-space constraints

Numerical simulation of constrained dynamical systems is known to exhibit stability problems even when the unconstrained system can be simulated in a stable manner. We show that not the constraints themselves, but the transformation of a continuous set of equations to a discrete set of equations is the true source of the stability problem....