PadeApproximant[r(1+r)^n/((1+r)^n-1),{r,0,{3,3}}]