Untitled
function [L U]=LUFactorization(A) [nr nc]=size(A); L=eye(nr); if nr==nc %allora possiamo effettuare i calcoli for i=1:nr for j=i+1:nr if A(i,i)==0 disp(['errore: pivot nullo!']) else L(j,i)=A(j,i)/A(i,i); for k=i:nr A(j,k)=A(j,k)-L(j,i)*A(i,k) end end end end else disp(['errore: la matrice dei coefficienti non è quadrata']) end U=A; end
Leave a Comment