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.