x=linspace(-5,5);
f=sinc(x);
g=.5*sinc(x)+.25*sinc(x-1)+.25*sinc(x+1);
h=sinc(x)+sinc(x-1)+sinc(x+1);

plot(x,g,x,f,x,h);

j=cos(4*pi*x)
ws2=w2.*j
ws1=w1.*j
fws2=fft(ws2.*j)
fws1=fft(ws1)
fs=fft(j)
plot(x,real(fws),x,fs,x,real(fws2))
fws2=fft(ws2.*j)