Numerical solution of differential algebraic equations and computation of consistent initial/boundary conditions

Boundary value methods for the solution of differential-algebraic equations are described. We consider both initial and boundary value problems and derive an algorithm that does not require additional information from the user, but only the initial or boundary conditions needed in theory to obtain a unique solution.