Abstract—Difference schemes of the Euler and trapezoidal types for the numerical solution of the ini tial value problem for linear differential algebraic equations are examined. These schemes are ana lyzed for model examples, and their superiority over