syms x; f=abs(x)/x;  % 定義方波信號
xx=[-pi:pi/200:pi]; xx=xx(xx~=0); xx=sort([xx,-eps,eps]); % 剔除零點
yy=subs(f,x,xx); plot(xx,yy), hold on  % 繪制出理論值并保持坐標系
for n=2:20
    [a,b,f1]=fseries(f,x,n), y1=subs(f1,x,xx); plot(xx,y1)
end
% 本例在 MATLAB 6.1, 6.5 和 7.0 中正常， 但6.5.1有問題

[a,b,f1]=fseries(f,x,14); latex(f1)