NDSolve::index
NDSolveValue::index
ParametricNDSolve::index
ParametricNDSolveValue::index
The differential algebraic equation system
is an index-3 system that cannot be solved without performing index reduction:
NDSolve[{x'[t] + y[t] == 0, y'[t] + z[t] == 0, z'[t] + w[t] == 0, y[t] == Sin[t], x[0] == 1}, {x, y, z}, {t, 0, 1}]Differentiating the last equation and substituting into the second equation gives
. Differentiating again and substituting the third equation gives
. Differentiating again gives the final system. Since the DAE was differentiated three times, the index is said to be 3:
NDSolve[{x'[t] + y[t] == 0, y'[t] + z[t] == 0, z'[t] + w[t] == 0, y[t] == Sin[t], x[0] == 1}, {x, y, z}, {t, 0, 1}, Method -> {"IndexReduction" -> True}]