Solution Numerical of Stiff ODEs and DAEs in Mechanical and Other Systems
Abstract Presented is a formal solution of the initial-value problem of the system of general implicit differential-algebraic equations (DAEs) F(x, y, y’) = 0 of index zero or higher, based on perturbations of the polynomial coefficients of the vector y(x). The equation is linearized with respect to the coefficients and brought into a form suitable for implementation of the weighted residual methods. The solution is advanced by a single-step multi-stage collocation qadrature formula which is stiffly accurate and suitable for solving stiff differential equations and DAEs that arise in many mechanical and other systems. The algorithm is illustrated by two index-2 and index-3 examples — one of which is the well known pendulum problem.