Some new GP features
harmonic
harmonic(n,k=1) computes the sum
Pn
i=1 1/ik as a rational
number.
? harmonic(10)
%62 = 7381/2520
? sum(i=1,10,1/i)
%63 = 7381/2520
? harmonic(10,3)
%64 = 19164113947/16003008000
? sum(i=1,10,1/i^3)
%65 = 19164113947/16003008000