Find a formula for $b_m$ by evaluating both sides for $f(x)=e^{\lambda x}$ where $\lambda$ is a parameter. The formula is $\int_0^1 f(x)dx=1/2(f(0)+f(1))+\sum_{m=1}^\inf b_m(f^{(2m-1)}(1)-f^{(2m-1)}(0))$ for some unknown constants $b_m$ independent of f.
I plug in $f(x)=e^{\lambda x}$ and get that $\sum_{m=1}^\inf b_m=(2e^\lambda -2-\lambda e^\lambda-\lambda)/2\lambda e^\lambda$. I don't know how to get rid of the sum. Any thought?