function xi = SI(T,c)
//
// Funzione SI: calcola la colonna xi approssimazione della soluzione del
// sistema di equazioni lineari
//
//       T x = c
//
// con T matrice N x N triangolare superiore 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(N) = c(N) / T(N,N);
for k = N-1:-1:1,
    xi(k) = (c(k) - T(k,k+1:N)*xi(k+1:N)) / T(k,k);
end;
endfunction
