%?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?