gauss1.m

来自「用matlab写的一些数值算法」· M 代码 · 共 14 行

function  [U, c] = gauss1(A, b)
% Gaussian elimination applied to Ax = b

% Version 11.12.2003.  INCBOX

n = size(A,1); 
b = b(:);  % ensure that b is a column vector
for  k = 1 : n-1
  ii = k+1 : n;
  A(ii,k) = A(ii,k)/A(k,k);  % multipliers
  A(ii,ii) = A(ii,ii) - A(ii,k)*A(k,ii);  % modify A
  b(ii) = b(ii) - b(k)*A(ii,k);  % modify b
end
U = triu(A);  c = b;  % return transformed problem