NDSolve::indexss
NDSolveValue::indexss
ParametricNDSolve::indexss
ParametricNDSolveValue::indexss
The mass matrix associated with certain differential algebraic systems is singular even after index reduction at certain places:
NDSolve[{x'[t] + t y[t] == Sin[t], x[t] + t y[t] == Cos[t], x[0] == 1}, {x, y}, {t, 0, 1}]The integration is successful if the singularity can be avoided:
NDSolve[{x'[t] + t y[t] == Sin[t], x[t] + t y[t] == Cos[t], x[0.001] == 1}, {x, y}, {t, 0.001, 1}]