echo on
fnct='normal';
a=-6;
b=6;
n=5;
tol=0.1;
xx=fseries(fnct,a,b,n,tol,0,1);
xx1=xx(n+1:-1:2);
xx1=[conj(xx1),xx];
absxx1=abs(xx1);
pause
n1=[-n:n];
stem(n1,absxx1)
title('the discrete magnitude spectrum')
phasexx1=angle(xx1);
pause
stem(n1,phasexx1);