% 《工程線性代數（MATLAB版）》第一章例題1.5程序ea105
% 《工程線性代數（MATLAB版）》引用
% 陳懷琛，高淑萍，楊威合編，電子工業出版社，2007年6月
% 消元法舉例
% 
clear
A=[0,1,-2,6,7;1,1,2,-1,2;2,3,2,4,3;-3,-4,-4,-3,-5] ,pause
A1=A;A1([1,2],:)=A1([2,1],:),pause    % 交換第1，2行
%A1(2,:)=A1(2,:)-A1(1,:)*A1(2,1)/A1(1,1),pause			% 消去A(2,1)
A1(3,:)=A1(3,:)-A1(1,:)*A1(3,1)/A1(1,1),pause			% 消去A(3,1)
A1(4,:)=A1(4,:)-A1(1,:)*A1(4,1)/A1(1,1),pause			% 消去A(4,1)
A1(3,:)=A1(3,:)-A1(2,:)*A1(3,2)/A1(2,2),pause			% 消去A(3,1)
A1(4,:)=A1(4,:)-A1(2,:)*A1(4,2)/A1(2,2),pause			% 消去A(4,1)
A1(4,:)=A1(3,:)-A1(3,:)*A1(3,5)/A1(3,5),pause			% 消去A(4,5)
A1(2,:)=A1(2,:)-A1(3,:)*A1(2,5)/A1(3,5),pause			% 消去A(2,5)
A1(1,:)=A1(1,:)-A1(3,:)*A1(1,5)/A1(3,5),pause			% 消去A(1,5)
A1(1,:)=A1(1,:)-A1(2,:)*A1(1,2)/A1(2,2),pause			% 消去A(1,2)