gsxqf.m

来自「利用高斯消元法、雅克比消元法解决求解线性方程组的解的程序」· M 代码 · 共 22 行

function x=gausse(A,b)
%本程序用Guass消去法（无交换的消去法）解线性代数方程。
[n m]=size(b);
B = [A b];
for k = 1:n-1
   for i = (k+1):n
      for j = (k+1):(n+1)
         B(i,j) = B(i,j)-B(i,k)*B(k,j)/B(k,k);
      end
   end
end
%%此时B仅上三角元素改变，以下解方程也只用到了上三角部分。
%%以下解上三角方程。
x = ones(size(b));
for i = 1:n  %%MALAB5不认识i=n:1,但MATLAB6认识。
   i = n+1-i;
   ss = 0;
   for k = i+1:n
      ss = ss+B(i,k)*x(k);
   end
   x(i)=(B(i,n+1)-ss)/B(i,i);
end