matlab实现高斯消元法

% gauss.m
function x = gauss（AB）
%The sizes of matrices AB are supposed to be NA x NA and NA x NB.
%This function solves Ax = B by Gauss elimination algorithm.
NA = size（A2）; [NB1NB] = size（B）;
if NB1 ~= NA error（‘A and B must have compatible dimensions‘）; end
N = NA + NB; AB = [A（1:NA1:NA） B（1:NA1:NB）]; % Augmented matrix
epss = eps*ones（NA1）;
for k = 1:NA
%Scaled Partial Pivoting at AB（kk） by Eq.（2.2.20）
[akxkx] = max（abs（AB（k:NAk））./max（abs（[AB（k:NAk + 1:NA） epss（1:NA - k + 1）]‘））‘）;
if akx < eps error（‘Singular matrix and No unique solution‘）; end