function [x, n, Xs, om]=filterDT(om1,om2,N)
%function [x, n, Xs, om]=filterDT(om1,om2,N)
%om1, om2 frequenze di taglio
%om1>=0, om2<=pi greco, om1<om2 
%campioni risposta = 2N+1
%x segnale
%n asse dei tempi
%Xs trasformata di Fourier
%om pulsazione
n=-N:N;
if om1<0 || om1 > pi || om2<0 || om2>pi || om2<=om1
    disp('error')
    return
end
if om1>=0 && om2 <pi
    x=(sin(om2*n)-sin(om1*n))./(n*pi);
    x(n==0)=(om2-om1)/pi;
else
    x=-sin(om1*n)./(n*pi);
    x(n==0)=1-om1/pi;
end
Xt=fft(x);
Xs=fftshift(Xt);
om=-pi+pi/length(n):2*pi/length(n):pi-pi/length(n);