  A generalisation of the Faulhaber polynomials and Bernoulli numbers related
to elliptic curves is introduced and investigated. This is applied to compute
the density of states for the classical Lam\'e operators.
