% mm2203.mN = 25;	% number of harmonicsTo = 0.2;	% choose periodn = 2*N;t = linspace(0,To,n+1); % (n+1)th point is one period awayt(end) = []; % throw away undesired last pointf = sawtooth(t,To); % compute sawtoothFn = fft(f);% compute FFTFn = [conj(Fn(N+1)) Fn(N+2:end) Fn(1:N+1)]; % rearrange valuesFn = Fn/n; % scale results
Bn=-2*imag(Fn(N+2:end));
idx=-N:N;
Fna = 5j./(idx*pi);
Fna(N+1)=5;
Bna=-2*imag(Fna(N+2:end));

Bn_error=(Bn-Bna)./Bna;

stem(idx,abs(Fn))xlabel('Harmonic Index')title('Figure 22.3: Sawtooth Harmonic Content')
axis tight