function [L, U]=my_lutrid(A) % trova i fattori LU della fattorizzazione LU di A % A e' tridiagonale senza pivoting % [n, m]=size(A); L=eye(n); U=A; for k=1:n-1 pivot=U(k, k); %sono gli alpha(i) %L(k+1:n, k)=U(k+1:n, k)/pivot; L(k+1, k)=U(k+1, k)/pivot; % elementi sottodiag L %colonnaL= L(k+1:n, k); %rigaU=U(k, k+1:n); U(k+1, k+1)= U(k+1, k+1)-L(k+1, k)*U(k, k+1); U(k+1:n, k)=0; end end