% Main C: discontinuità derivata, derivata approssimata
f = @(x) abs( cos(x) ) ;

% Definizione derivata approssimata per determinare gli intervalli di partenza
h_approx     = 1e-5 ; % Attenzione! h dev'essere sufficientemente piccolo
df_dx_approx = @(x) ( f(x+h_approx)-f(x-h_approx) ) / (2*h_approx) ;

x_ab = trova_discontinuita( df_dx_approx , [0 2*pi] ) ;
n_discontinuita = size( x_ab , 2 ) ;
x_discontinuita = zeros(1, n_discontinuita) ;

% Derivata mediante limite, g = df/dx
g = @(x0) limite( @(h) ( f(x0+h)-f(x0-h) ) / (2*h) , 0 , .1 , 1e-6 , true ) ;
for i = 1 : n_discontinuita
  ab    = x_ab(:,i) ;
  g_ab  = g( ab ) ;
  g_mid = sum( g_ab ) / 2 ;
  p = @(x) g(x) - g_mid ;
  tol_x = 1e-6 ;
  x_discontinuita(i) = metodo_bisezione(p, ab, tol_x) ;
end