Rsolve( { b[n+1]=b[n]+b[n]/(n+1) )