An algorithm for a numerical solution of differential-algebraic equations

Presented is the solution algorithm of initial value problem of the general implicit system of differential-algebraic equations (DAE) f(x,y,y´)=0. System is linearized with respect to polynomial coefficients in y and the solution is advanced by a single-step multi-stage collocation method. The algorithm turns out to be robust and stable, as well as a convenient tool for derivation of all possible collocation quadrature formulae and for designing their desired properties. The method is suitable for solving stiff differential equations and DAE that arise in many mechanical and control systems