%例35-23：QR分解法求解下列方程組。
A=[3,1,-4,1;1,-5,0,7;0,2,1,-1;1,6,-1,-4]; %系數矩陣A
b=[13,-9,6,0]';  %矩陣B
[Q,R]=qr(A);  %QR分解
x=R\(Q\b)     % 求解
[Q,R,E]=qr(A); %  采用QR的第2種格式分解
x1=E*(R\(Q\b))   % 求解