sum[n^k/(r^n)] with n from 1 to infinity