function cesaromittel = cesaro(x,cvec)

clgt = length(cvec);
cesaromittel = zeros(clgt,length(x));

for j = 1:clgt
    y = partsum(x,1);
    for k = 2:cvec(j)
        y(k,:) = y(k-1,:) + partsum(x,k);
    end
    cesaromittel(j,:) = y(end,:)/cvec(j);
end