A note on numerically consistent initial values for high index differential-algebraic equations
When differential-algebraic equations of index 3 or higher are solved with backward differentiation formulas, the solution in the first few steps can have gross errors, the solution can have gross errors in the first few steps, even if the initial values are equal to the exact solution and even if the step size is kept constant. This raises the question of what are consistent initial values f
