support computation of Bernoulli polynomials