N=10;n=0:1:N-1;x1=(0.8).^n;k=0:1:N-1;WN=exp(-j*2*pi/N);nk=n'*k;WNnk=WN.^nk;Xk=x1*WNnk;subplot(2,1,1)stem(n,x1);subplot(2,1,2)stem(k,abs(Xk));