polygonal(3, 1:10) # 1 3 6 10 15 21 28 36 45 55
polygonal(4, 1:10) # 1 4 9 16 25 36 49 64 81 100
polygonal(5, 1:10) # 1 5 12 22 35 51 70 92 117 145
polygonal(6, 1:10) # 1 6 15 28 45 66 91 120 153 190
# Sums of reciprocals of polygonal numbers:
n = 1000
sum(1/polygonal(3, 1:n)) # 1.998002 -> 2.0
sum(1/polygonal(4, 1:n)) # 1.643935 -> pi^2/6
sum(1/polygonal(5, 1:n)) # 1.481371 -> 3*log(3)-pi/sqrt(3)
sum(1/polygonal(6, 1:n)) # 1.385794 -> 2*log(2)
Run the code above in your browser using DataLab