function xi = SA(T,c)
//
// Funzione SA: calcola la colonna xi approssimazione della soluzione del
// sistema di equazioni lineari
//
//       T x = c
//
// con T matrice N x N triangolare inferiore invertibile e c colonna di N
// componenti.
//
[NR, NC] = size(T);
N = NR; // // il sistema ha dimensione N
//
// *** Calcolo soluzione
//
xi = zeros(N,1);
//
xi(1) = c(1) / T(1,1);
for k = 2:N,
    xi(k) = (c(k) - T(k,1:k-1)*xi(1:k-1)) / T(k,k);
end;
endfunction
