Algorithme de la Factorisation de Gauss en Matlab

Algorithme de la Factorisation de Gauss en Matlab :

function A=lugauss (A) %LUGAUSS Factorisation LU sans pivot. % A = LUGAUSS(A) stocke une matrice triangulaire % supérieure dans la partie triangulaire supérieure de % A et une matrice triangulaire inférieure dans la % partie strictement triangulaire inférieure A (les % termes diagonaux de L valant 1). [n,m]=size(A); if n ~= m error(’A n’’est pas une matrice carrée’); else for k = 1:n-1 for i = k+1:n A(i,k) = A(i,k)/A(k,k); if A(k,k) == 0, error(’Elément diagonal nul’); end j = [k+1:n]; A(i,j) = A(i,j) - A(i,k)*A(k,j); end end end return