Drazin inverse based numerical methods for singular linear differential systems

In this paper, numerical methods for the solution of linear singular differential system are analysed. The numerical solution of initial value problems by means of a finite difference approach and a possible implementation of the Drazin inverse vector product is discussed. Examples of index 1 and index 2 differential algebraic equations have been studied numerically.