Решить СЛАУ методом Гауса. Matlab

% А - матрица коэффициентов левой части
% b - вектор-столбец правой части
% На выходе - вектор-столбец значений неизвестных
C=[A,b];
[m,n]=size(C);
for k=1:m-1
p=C(k,k)./C(k+1:m,k);
P=repmat(p,1,n-k+1);
C(k+1:m,k:n)=repmat(C(k,k:n),m-k,1)-P.*C(k+1:m,k:n);
end
X(m)=C(m,n)/C(m,n-1);
for k=m-1:-1:1
C(1:k,k+1)=C(1:k,k+2)-X(k+1)*C(1:k,k+1);
X(k)=C(k,k+1)/C(k,k);
end
X'