Is there a good way to compute the ratio ( B[n] / n! ) that occurs so often in power series coefficients? Good in the sense that you get an answer that does not overflow a double; the largest n such that B[n] fits in a double is relatively small, but the ratio in question should never overflow, although I guess it will eventually underflow.


The Wikipedia article on Bernoulli numbers has a simple algorithm for the Bernoulli numbers, more than adequate for this range, but it's not obvious how to modify it to answer this question. To clarify, I'm not asking for an exact computation, but a decimal computation.



computing Bernoulli numbers
these are socalled "scaled" Bernoulli numbers; an efficient algorithm is discussed by Brent & Harvey, [arXiv:1108.0286](http://arxiv.org/abs/1108.0286), Eq. 8.


alternatively, you can use any method that computes the Riemann zeta function (say via the Euler product), because of the identity


$$\frac{B\_{2n}}{(2n)!}=(-1)^{n+1}2\zeta(2n)(2\pi)^{-2n}$$


(note that $B\_{2n+1}=0$ for $n\geq 1$)


