% Input: f, function handle della funzione da integrare
%        x0, primo estremo dell'integrazione, scalare
%        x, vettore del secondo estremo d'integrazione
% Output: F, funzione integrale valutata sui valori di x
function F = F_integrale( f , x0 , x )
  I = 0*x ;
  if size(x,1) > 1
    x = [ x0 ; sort( x ) ] ;
  else
    x = [ x0   sort( x ) ] ;
  end
  for i = 1 : length(x)-1
    I(i) = integral( f , x(i) , x(i+1) ) ;
  end
  F = cumsum( I ) ;
end