Hello there,

I am trying to solve a system of non-linear differential-algebraic equations using NDSolve, but keep getting the following error :

NDSolve::icfail: Unable to find initial conditions that satisfy the residual function within specified tolerances. Try giving initial conditions for both values and derivatives of the functions.

The code looks like following :

Sol = NDSolve[{DGL, (xp[t] /. 
      t -> 0) == {-0.05, -0.052, -0.054, -0.06, -0.057, 0, 0, 0, 0, 0,
      1, 1.5, 2, 1.1, 1.2}, (D[xp[t], t] /. t -> 0) == {0, 0, 0, 0, 0,
      0, 0, 0, 0, 0, 0, 0, 0, 0, 0}, 
   Table[\[Omega][t][[k]] == Subscript[L, k][t]/M[[k, k]], {k, 1, 
     Dimensions[AdjazenzMatrix][[1]], 1}]}, xp[t], {t, 0, 100}, 
  MaxSteps -> 100, AccuracyGoal -> 8]

The system of differential-algebraic equations is defined as :

DGL = Chop[
  Table[(D[xp[t], t])[[k]] == (JR.DeltaHpNeu [t])[[k]] + 
     ColG[t][[k]] + (G.up)[[k]], {k, 1, 
    3*Dimensions[AdjazenzMatrix][[1]], 1}]]

Does the error mean that there exists no solution for the system or is there something wrong with my code?

Thanks a lot! P.S. - I am attaching the Mathematica source code.

Attachments:

Beispielnetz - P...nb