pdl> $n = sequence(15)+1

pdl> p $a = 1 + 1/$n
[2 1.5 1.33333333333333 1.25 1.2 1.16666666666667 1.14285714285714 1.125 1.11111111111111 1.1 1.09090909090909 1.08333333333333 1.07692307692308 1.07142857142857 1.06666666666667]
pdl> p $s = $a->sumover
18.318228993229