% 《工程線性代數（MATLAB版）》第一章例題1.4程序ea104a
% 陳懷琛，高淑萍，楊威合編，電子工業出版社，2007年6月
% 用消元子程序消元法舉例
% 
A=[2,-2,0,6;2,-1,2,4;3,-1,4,4;1,1,1,8] ,pause
b=[-2;-2;-3;2] ,pause			% 給系數矩陣A,b賦值
B=[A,b] ,pause					% 賦值B為增廣矩陣
B1=gauss(B,1,2,1),pause			% 消去B(2,1)
B1=gauss(B1,1,3,1),pause		% 消去B1(3,1)
B1=gauss(B1,1,4,1),pause		% 消去B1(4,1)
B1=gauss(B1,2,3,2),pause		% 消去B1(3,2)
B1=gauss(B1,2,4,2),pause		% 消去B1(4,2)
B1([3,4],:)= B1([4,3],:) 		% 由于B1(3,3)=0,需交換B1中的第3,4行
B1=gauss(B1,3,4,3),pause		% 消去B1(4,3)
% 至此完成了行階梯變換,最后的B1就是行階梯形式的等價增廣矩陣
B2= gauss(B1,4,3,4),pause		% 消去B1(3,4)
B2= gauss(B2,4,2,4),pause		% 消去B2(2,4)
B2= gauss(B2,4,1,4),pause		% 消去B2(1,4)
B2= gauss(B2,3,2,3),pause		% 消去B2(2,3)
B2= gauss(B2,3,1,3),pause		% 消去B2(1,3)
B2= gauss(B2,2,1,2),pause		% 消去B2(1,2)
% 至此完成了對角形的行階梯變換 
B2(1,:)=B2(1,:)/B2(1,1),pause		% B2第一行對角元素歸一
B2(2,:)=B2(2,:)/B2(2,2),pause		% B2第二行對角元素歸一
B2(3,:)=B2(3,:)/B2(3,3),pause		% B2第三行對角元素歸一
B2(4,:)=B2(4,:)/B2(4,4) 			% B2第四行對角元素歸一
% 至此完成了最簡形行階梯變換，B2是所需的結果。