for kl=1:64
   u(kl)=(kl)*pi/64.0;
   x(kl)=sin((kl)*pi/64.0);
   x(kl)=x(kl)+1/2*sin(2*(kl)*pi/64.0)+1/3*sin(3*(kl)*pi/64.0)+1/4*sin(4*(kl)*pi/64.0);
   x(kl)=x(kl)+1/5*sin(5*(kl)*pi/64.0)+1/6*sin(6*(kl)*pi/64.0)+1/7*sin(7*(kl)*pi/64.0);
   x(kl)=x(kl)+1/8*sin(8*(kl)*pi/64.0)+1/9*sin(9*(kl)*pi/64.0)+1/10*sin(10*(kl)*pi/64.0);
 end

 plot (x)
 y=fft(x);
 z=fftshift(y);
 
 
 
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))