function gp = gibbs_perc(mvec)
if max(mvec)<=100
    h = .01;
else
    h = max(mvec)^(-1)/2;
end

xn = -pi:.01:0;
xp = 0:.01:pi;
zn = -pi-xn;
zp = pi-xp;

x1 = -pi:h:0;
x2 = 0:h:pi;
x=[x1 x2];
s = partsum(x,mvec);

sprung = max(x)-min(x);
[nz,ns] = size(s);

for k = 1:nz
    gp(k) = max(s(k,:)) - min(s(k,:));
end
gp = gp/sprung*100-100;

plot(x,s,xn,zn,'k',xp,zp,'k')
title('m-te Fourier-Partialsumme(n) zum ersten Bernoulli-Polynom')
% Beispielinput
% gibbs_perc([10,20,200,400])