Maxima Function
burn (n)
Returns the n'th Bernoulli number for integer n.
burn may be more efficient than bern for large, isolated n
(perhaps n greater than 105 or so),
as bern computes all the Bernoulli numbers up to index n before returning.
burn exploits the observation that (rational) Bernoulli numbers can be
approximated by (transcendental) zetas with tolerable efficiency.
load ("bffac") loads this function.