lu.m

来自「数值分析, 同济大学教材<<现代数值数学和计算>>数值分」· M 代码 · 共 13 行

function [L U] = LU(A)
n = size(A, 1);
L = eye(n);
U = zeros(n);
for k = 1 : n
  U(k, k:n) = A(k, k:n) - L(k, 1:k-1) * U(1:k-1, k:n);
  if U(k, k) == 0
    error('LU Failed.');
  end
  L(k+1:n, k) = (A(k+1:n, k) - L(k+1:n, 1:k-1) * U(1:k-1, k)) / U(k, k);
end