def u(n): U = [None] * (n+1) U[0] = (1,0) for k in range(1,n+1): S = U[k-1][1] + k * U[k-1][0] U[k] = (S/k, S) return U[n][0] for n in range(20): print (u(n))