Aby uzasadnić wzór (25) skorzystajmy ze wzoru Eulera-Maclaurina:
y*1 -j m-1
/ f(0^=2[f(°)+f(1)]+E>A2/c[^'1)(O)-f(2/f'1)(t)]-A2mf(2,n)(Co).
/c=1
gdzie £0 € (0.1). Związek współczynników A2r z liczbami Bernoulliego nie jest tu istotny. Przekształcając odpowiednio zmienną mamy
fcfi i
(26) Yt
*=1
gdzie £/ g (X/, x/+1 ). Sumując stronami (26) dla / = 0,1.....2n -1 i
zmieniając kolejność sumowania w podwójnej sumie (wszystkie składniki w niej zredukują się prócz skrajnych, w ostatniej sumie średni pochodnych zamieniamy pochodną w innym punkcie) otrzymamy