[t,y,S,aquad,bquad]=fzzyquad(0,2,15);	%計算0到15次諧波		<1>
AA=abs(A_sym);AA(AA<1e-10)=NaN;	%在那些實際應為0的很小數處，考慮相對誤差
										%沒有意義，故予以排除。
BB=abs(B_sym);BB(BB<1e-10)=NaN;
A_quad=aquad(1:7)			%數值計算所得偶函數分量系數
B_quad=bquad(1:7)			%數值計算所得奇函數分量系數
asq=abs(A_quad-A_sym)./AA	%偶函數系數與符號計算值的相對誤差
bsq=abs(B_quad-B_sym)./BB	%奇函數系數與符號計算值的相對誤差  

SS=[S;y]';ribbon(t,SS);%把原始波形列在圖形最前方
xlabel('諧波次數\rightarrow'),ylabel('t\rightarrow')
view([46,38]),colormap(jet),shading flat,light,lighting gouraud