计算方法Gauss-Jordan消去法求线性方程组的解

function x=GaussJordan(A,b)
% 用Gauss-Jordan消去法求线性方程组的解
% A为系数矩阵(要求解惟一，则必须为满秩方阵 n*n),b为右值 n*1
% x为所求的根

%作者：野渡无人
% 最后修改日期：2008.3.27

% 例如教材45页例3.1：
% A=[1 4 -5;1 3 -2;6 -1 18], b=[3 2 2]'
% x=GaussJordan(A,b)

% 结果为：
%     1     4    -5     3
%     1     3    -2     2
%     6    -1    18     2
%
%     1     4    -5     3
%     0    -1     3    -1
%     0     0   -27     9
%  x =
%
%     1.3333
%          0
%