% Input: f, function handle alla funzione f
%        x0, valore di x dove si vuole controllare la continuità
% Output: b, booleano: vero se f è continua in x0, falso altrimenti
function b = is_continua( f , x0 )
  h = 0.1 ;
  salto_tollerato_y = 1e-8 ;               % salto funzione
  salto_tollerato_x = abs(x0) * 10 * eps ; % salto incremento in x = h
  while abs( f(x0+h)-f(x0-h) ) > salto_tollerato_y && h > salto_tollerato_x
    h = h / 2 ;
  end
  b = h > salto_tollerato_x ;
end