Solving system of linear equation by Gaussian elimination

How I can get the solution of system of linear equation by Gaussian elimination in MATLAB ?

Answer 1

[Andrei Bobrov]

variant solution of system A*x = B, eg:

A= randi([-23 907],5)
B = randi([-23 907],5,1)
c = [A B];
[m,n] = size(c);
rezc = cell(m,1);
for j1 = 1:m
    [Y,ind] = sortrows(abs(c(j1:end,j1:end)),1);
    c(j1:end,j1:end) = c(j1-1+flipud(ind),j1:end);
    c1 = c(:,j1)*ones(1,n);
    c1(j1,:) = zeros(1,n);
    c2 = (ones(m,1)*c(j1,:)/c(j1,j1));
    c(j1,:) = c2(j1,:);
    c = c - c2.*c1;
    rezc{j1} = c;
end
x = c(:,end)
rezc{:} % so rrefmovie

OR

from book Numerical Methods Using MATLAB, 4e by John H. Mathews and Kurtis D. Fink files backsub.m and uptrbk.m

Answer 2

[shitao You]

just use rref(a), and a is a matrix.