for i=1:1000
   x(i)=0;
   end
>> for i=480:520
   x(i)=sin((i-480)*pi/40);
  end
>> plot (x)
 y=fft(x);
 z=fftshift(y);
 bar(z(400:600))
New
 for i=1:1000
   x(i)=0;
   end
 for i=480:520
   x(i)=sin((i-480)*pi/8);
  end
  print -deps Plot2.eps
0
plot(x(400:600))