decomLU

Esta é uma pré-visualização de arquivo. Entre para ver o arquivo original

function [L, U] = decomLU(A)
 [lin col] = size(A);
 L = eye(lin, col);
 
 for i = 1:lin-1
 pivo = A(i, i)
 for j = i+1:lin
 m = A(j,i)/pivo;
 A(j,:) = A(j,:) - m*A(i,:);
 L(j,i) =m; 
 end
 end
 U = A;
endfunction