function?[xQR]=qrxq（Ab）?
N?=?size（A）;
n?=?N（1）;
B?=?A;????????????????????????%保存系數矩陣

A（1:n1）=A（1:n1）/norm（A（1:n1））;?%將A的第一列正規化
for?i=2:n
????for?j=1:（i-1）
????????A（1:ni）=?A（1:ni）-dot（A（1:ni）A（1:nj））*