function summe = partsum(x,vec) warning off; lgt = length(vec); i = sqrt(-1); summe = zeros(length(vec),length(x)); for k = 1:vec(1) summe(1,:) = summe(1,:) + exp((i*k).*x)./(i*k) + exp((i*(-k)).*x)./(i*(-k)); end for s = 2:lgt summe(s,:) = summe(s-1,:); for k = (vec(s-1)+1):vec(s) summe(s,:) = summe(s,:) + exp((i*k).*x)./(i*k) + exp((i*(-k)).*x)./(i*(-k)); end end %figure; %hold on; %plot(x,summe)