% Suddividiamo l'intervallo dato in 1000 punti (per esempio)
x = linspace( -pi/2 , pi/2 , 1000 ) ;

% 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