Preservation of stability type of difference schemes when solving stiff differential algebraic equations

Implicit methods applied to the numerical solution of systems of ordinary differential equations (ODEs) with an identically singular matrix multiplying the derivative of the sought-for vector-function are considered. The effects produced by losing $L$-stability of a classical implicit Euler scheme when solving such stiff systems are discussed.