function x=solveLTcs(A,b) %SOLVELTCS Risolve il sistema lineare Ax=b con A matrice triangolare % inferiore utilizzando il metodo column sweep % % x=SOLVELTCS(A,b) % % I parametri della funzione sono: % A -> la matrice dei coefficienti del sistema lineare % b -> il vettore dei termini noti % % I valori di ritorno sono: % x -> il vettore soluzione del sistema lineare % % See Also SOLVELT n=length(b); x=b; for j=1:n if A(j,j)==0 disp('Matrice non risolvibile') break end x(j)=x(j)/A(j,j); for i=j+1:n x(i)=x(i)-A(i,j)*x(j); end end