function r=eig_qr(A)
%qr方法計算矩陣特征值
A=hessenberg(A);
%將矩陣轉化為上Hessenberg矩陣
[q1,r1]=qr_decomposition(A);
%對上Hessenberg矩陣進行qr分解，直到得到的新矩陣的對角線元素趨近于其特征值
A1=r1*q1;
[q2,r2]=qr_decomposition(A1);
A2=r2*q2;
while norm(diag(A1-A2))>1e-5
%設定允許誤差
[q1,r1]=qr_decomposition(A2);
A1=r1*q1;
[q2,r2]=qr_decomposition(A1);
A2=r2*q2;
end
r=diag(A2);