% Input: f, cell array dei function handle delle funzioni f_i
%        ab, intervallo dentro il dominio
% Output: I, booleano che indica se le f_i sono linearmente indipendenti
function I = f_indipendenti( f , ab )
  n = length( f ) ; % numero di funzioni
  m = 2*n ;         % numero di punti
  x = linspace( ab(1) , ab(2) , m )' ;
  F = zeros( m , n ) ; % allocazione matrice F
  for i = 1 : n
    F(:,i) = f{i}(x) ; % valutazione f_i( x1 , ... , xm )
  end
  kF = null( F ) ; % ker( F )
  I  = size( kF , 2 ) == 0 ;
end