QR分解法求解矩阵的特征值与特征向量

设计思想：用带双步位移的QR分解法求矩阵A（10*10）的全部特征值。在计算出A的基础上，先利用Householder矩阵对矩阵A作相似变换，把A化为拟上三角矩阵A(n-1)，然后进行带双步位移的QR分解（其中Mk的QR分解可调用子程序），通过调用一元二次方程的根解二阶块矩阵的特征值，最后计算出A(n-1)的特征值，即为A的特征值，然后对实特征值利用列主元高斯消元法求解其对应的特征向量。