Ich mach das jetzt mit einem Matrizenrechner.
https://matrixcalc.org/de/#{{1,1,0},{1,a,0},{1,1,a}}^2
Der Grundgedanke; jede matrix löst ihre eigene Säkulardeterminante:
A ³ + b2 A ² + b1 A + b0 * 1| = 0 | * A ^ -1 ( 1a )
A ² + b2 A + b1 * 1| + b0 A ^ -1 = 0 ( 1b )
1 1 0
1 a 0 =: A ( 2a )
1 1 a
2 a+1 0
a+1 a^2+1 0 = A ² ( 2b )
a+2 2*a+1 a^2
a+3 a^2+a+2 0
a^2+a+2 a^3+2*a+1 0 = A ³ ( 2c )
a^2+3*a+3 3*a^2+2*a+2 a^3
Bereits in ( 2a ) wird klar: a ist Eigenwert von A ; für a = 0 wird A singulär . Wir berechnen jetzt ( 1a ) für Matrixelement ( 1 ; 2 ) und anschließend für ( 3 ; 1 )
( a + 1 ) b2 + b1 = - ( a ² + a + 2 ) ( 3a )
( a + 2 ) b2 + b1 = - ( a ² + 3 a + 3 ) ( 3b )
b2 = - ( 2 a + 1 ) ( 3c )
b1 = 2 a ² - ( a - 1 ) ² ( 3d )
( ( 3d ) mit Unterstützung von Wolfram ) Für ( 3c ) hast du eine Probe aus dem Satz von Vieta
b2 = - ( E1 + E2 + E3 ) = - Sp ( A ) ( 4 )
Jetzt zwecks Ermittlung von a0 Element ( 1 ; 1 )
2 b2 + b1 + b0 = a ² - 2 a - 3 + b0 = - ( a + 3 ) ( 5a )
b0 = - a ( a - 1 ) ( 5b )
Abermals hast du die Vietaprobe; denn
b0 = - E1 E2 E3 = - det ( A ) ( 5c )
Wir hatten gesagt, E3 = a ist Eigenwert. Welche gibt es noch? Ich mach das jetzt nicht mit Polynomdivision, sondern mittels der von mir entwickelten ersten und zweiten ===> Alfonsinischen pq-Formel. Mit Vieta sind die ( fast ) selbst erklärend:
b2 = - ( p + E3 ) = - ( 2 a + 1 ) ===> p = a + 1 ( 6a )
b0 = - q E3 = - a ( a - 1 ) ===> q = a - 1 ( 6b )
x ² - ( a + 1 ) x + a - 1 = 0 ( 6c )
Warum hab ich das jetzt eigens gemacht? Weil du in ( 6c ) die zweite Singularität kriegst, wenn a = 1 . Das haben wir jetzt auch verstanden.