I want to understand the asymptotics of a function which depends on a parameter $N$ which is a large integer.
The function is of the form (variable $x$ is of order 1) $$ F_N(x)=\sum_{n=0}^\infty f_{N,n}x^n.$$
I suspect that the coefficient $f_{N,N}$ dominates the asymptotics, so that $F_N(x)\approx f_{N,N}x^N$, in some sense. But I would like to actually prove this, and estimate the order of the error.
The coefficients are given by $$ f_{N,n}=\sum_{k=0}^{{\rm min}(n,N)-1}\frac{k!(n-1-k)!(N-1-k)!}{(N+n-1-k)!}.$$
What I have done: I can see that if $n$ is small, then $f_{N,n}\sim N^{-n+1}$. On the other hand, if $n=N$ then the contribution to the sum from $k=N-1-m$ is of order $N^{-2m}$, so the sum is dominated by $m=0$ and we get a value of order 1.
Using the Stirling approximation blindly, without considering whether each factorial is large or not, I get an approximation for the summand which is $$ \frac{k^k(n-k-1)^{n-k-1/2}(1-\frac{k+1}{N})^{N-k-1/2}}{N^n(1+\frac{n-k-1}{N})^{N+n-k-1/2}}$$ which could be useful. But this approximation is not valid in particular for the case $k=n-1=N-1$ which I think is the most important one.
How do I carry out a rigorous asymptotic analysis here?