[nm]=?size（A）;??%取A的行數n和列數m
Q?=?zeros（nm）;?%Q為n階單位陣
Q1?=?zeros（nm）;
R?=?zeros（mm）;
Q（:1）?=?A（:1）;
Q1（:1）?=?Q（:1）/norm（Q（:1））;
R（11）?=?norm（Q（:1））;
for?i?=?2:m
?%???Q（:i）?=?A（:i）;
????B?=?zeros（n1）;
????for?j?=?1:i-1
????????B?=?B?+?A（:i）