% Suddividiamo l'intervallo dato in 1000 punti (per esempio)

%x_0 = input( 'Inserire estremo intervallo: ' ) ;
%N   = input( 'Inserire numero di elementi: ' ) ;
N = 6 ;
x = linspace( -pi/2 , pi/2 , N ) ;

% Calcolo funzione
f = sin(2*x) - x ;

% Segno della funzione
segno_f = sign(f) ;

% Cerchiamo dove il segno della funzione cambia
ix = find( segno_f(1:end-1) ~= segno_f(2:end) ) ;

% Approssimazione degli zeri (media di due punti adiacenti)
zeri = ( x(ix) + x(ix+1) ) / 2 ;

f_zeri = sin(2*zeri) - zeri 

%disp(zeri) ;