I am trying to get an estimate for the induced 2-norm condition number $\kappa_{2}(M)$ [of this matrix $M$][1]:
$$M_{ij} = \frac{1}{(n-i)!(n-j)!(2n-i-j+1)} = \displaystyle\int_{0}^{1}\frac{x^{n-i}}{(n-i)!} \frac{x^{n-j}}{(n-j)!}\:{\rm{d}}x,$$
which is symmetric positive definite of size $n\times n$. The above link gives the elements of the inverse. We also know that its determinant is the reciprocal [of this OEIS sequence][2]. 

I want to obtain something like $\kappa_{2}(M) = \Omega(f(n))$ for some to-be-determined $f(n)$.

From numerics, it seems that $\kappa_2$ grows very fast. But I am not sure what might be an analytic way to obtain the asymptotic for $\lambda_{\max}/\lambda_{\min}$.


  [1]: https://mathoverflow.net/questions/370119/inverse-of-a-cauchy-like-matrix
  [2]: https://oeis.org/A107254