%tex1

L=16;N=64;
n=1:N;
xn=[ones(1,L),zeros(1,N-L)];
   Xk=dfs(xn,N);
   magXk=abs([Xk(N/2+1:N) Xk(1:N/2+1)]);
   k=[-N/2:N/2];
   figure(1)
   subplot(2,1,1);
      stem(n,xn);xlabel('n');ylabel('xtide(n)');
   title('DFS of SQ.wave:L=16,N=64');
   subplot(2,1,2);
   stem(k,magXk);
   axis([-N/2,N/2,0,16]);
   xlabel('k');ylabel('Xtide(k)');