sumpar.m

the Gibbs phenomenon graphic in matlab , you can change the number of armonics by changing the N val

function sumpar(n)

% f(x)= -1 si -pi<x<0
% 1 si 0<x<pi

x = -pi:0.001:pi;
f=-1*(x<0)+1*(x>0);

s = zeros(size(x));
for k=1:100
s=s+((1-(-1)^k)/k)*sin(k*x);
end
s = 2/pi*s;


plot(x, s, 'r', x, f, 'b'),grid;
title('Fenomeno de Gibbs');