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

x_0 = input( 'Inserire estremo intervallo: ' ) ;
N   = input( 'Inserire numero di elementi: ' ) ;

x = linspace( -x_0 , x_0 , 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 ;

disp(zeri) ;