clear;augm=input（‘Please?input?the?augumental?matric:‘）;
[nt]=size（augm）;
for?i=1:n?????????????%d控制列三角主元的層數
????b=augm（i:）;p=i;??%4到10行是選列主元并交換
????for?w=i+1:n
????????if?b（11）<=augm（wi）
????????????b=augm（w:）;p=w;
????????end
????end
????c=augm（i:）;?augm（i:）=b;augm（p:）=c;
????if?p~=i???????????????%只有發生了換行才將這種效果輸出
????????augm
????end

????for?j=i:t?????????%首先變換與augm（ii）同一行的元素，其列指標從i到t
????????s=0;
????????for?k=1:i-1
????????????s=s+