GitHub

lu_decomposition_doolittle.m

U(k, k:n) = M(k, k:n) - L(k, 1:k-1) * U(1:k-1, k:n); L((k+1):n, k) = (M((k+1):n, k) - L((k+1):n, 1:k-1) * U(1:k-1, k)) / U(k, k); ...