Hoi Andreas,
Es scheint tatsächlich ein Problem der Rechengenauigkeit zu sein. Durch
Verändern der Präzision erhält man das richtige Ergebnis (siehe Anhang).
In[20]:= Clear[w]
w = Eigensystem[SetPrecision[M, 4]][[2]]
Out[21]= {{0.05714, 0.05650, 0.05921, -0.9950},
{0.5884, 0.5838, 0.5505, 0.09975},
{0.7071, -0.7071, -0.00001116, -2.544*10^-10},
{0.3917, 0.3917, -0.8326, -0.004860}}
die Angabe der Matrixelemente als Festkommazahlen mit vier
Nachkommastellen reicht schlicht und einfach, mit anderen Worten
In[44]:= Clear[c1, c2, c3, c4, e1, e2, e3]
c1 = 34.8669`4;
c2 = 36.6194`4;
c3 = -3.7212`4;
c4 = -4.5547`4;
e1 = 1958.56`4;
e2 = 1950.11`4;
e3 = 2096.11`4;
und es funktioniert.
Gruss
Udo.