高斯消元 matlab

clear
n=input('请输入矩阵阶数:\n');
for i = 1:n
    fprintf('请输入矩阵第%d行',i);
    A(i,:) = input('\n');
end
fprintf('请输入矩阵B')
B(:,1)=input('');

for i = 1:n-1
   for j = i+1:n
       L=A(j,i)/A(i,i);
       for k = i:n
           A(j,k) = A(j,k)-A(i,k)*L;
       end
       B(j) = B(j) - B(i)*L;
   end
end
A
for i = n:-1:1
    m=n;
    while i<m
        B(i) = B(i)-A(i,m)*X(m);
        m=m-1;
    end
    X(i)=B(i)/A(i,i);
end