sump(n, p)

    function s = sump(n, p)

    s = 0;
    for i = 1:n
      s = s + i^p;
    end;