Ajee, mir fällt gerade auf, dass Deine Drehmatrix nicht passt:
Du hast S A S^-1 =D üblich ist S^-1 A S = D
Du müsstes also drehen mit: (x,y) S A S^-1 (x,y) + a S^-1 (x,y) = a_0
kommst Du überhaupt auf die gedrehte Form q_D (oben)?
hast Du dir den Link angesehen?
von q_D mit der quadratischen Ergänzung T zu q_N
! je nach Eigenvektoren können x,y,Vorzeichen getauscht sein!
z.B.