xx= linspace(0,2*pi,33);
for i=1:32
     x(i) = xx(i)
  end
  y = (cos(14..*x) + 2*sin(15.*x) + 3*cos(16.*x))/16.;
  plot(x,y);
  z = fft(y)
  zz=fftshift(z)