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.