% Main A: discontinuità funzione
% funzione Heaviside
f = @heaviside_f ;

% Approssimazione iniziale
x_ab = trova_discontinuita( f , [-1 1] ) ;
n_discontinuita = size( x_ab , 2 ) ;
x_discontinuita = zeros(1, n_discontinuita) ;

% Calcolo più accurato punto per punto
for i = 1 : n_discontinuita
  ab    = x_ab(:,i) ;
  f_ab  = f( ab ) ;
  f_mid = sum( f_ab ) / 2 ;
  g = @(x) f(x) - f_mid ;
  tol_x = 1e-10 ;
  x_discontinuita(i) = metodo_bisezione(g, ab, tol_x) ;
end