(p1, q1), (p2, q2)を固有ベクトルとして
(x1, x2) = a1 * (p1, q1) + b1 * (p2, q2)
と分解すれば、(an, bn)が求まるというやり方は分かります。
でも、>>56のやり方はそれとは違います。