Bernoulli numbers and solitons
  We present a new formula for the Bernoulli numbers as the following integral $$B_{2m} =\frac{(-1)^{m-1}}{2^{2m+1}} \int_{-\infty}^{+\infty} (\frac{d^{m-1}}{dx^{m-1}} {sech}^2 x)^2dx. $$ This formula is motivated by the results of Fairlie and Veselov, who discovered the relation of Bernoulli polynomials with soliton theory. 
M-P. Grosset and A.P. Veselov