% p.85 ex.2 x1 = (-1:0.01:0)'; x2 = (0:0.01:1)'; x = [x1; x2]; y1 = 0*x1; y2 = 0*x2+1; y = [y1; y2]; m = size(y, 1); f = 1/2 * ones(m,1); for n=1:30 bn = (1 - (-1)^n)/(n*pi); f = f + bn * sin(n*pi*x); plot(x, y, 'b', x, f, 'r'); n pause; end set(findobj('type','line'),'linewidth',4); print -depsc ex2.ps;