Thursday, January 02, 2014
Euler-Maclaurin Summation
A classic formula for estimating sums via integrals. Suppose f is a function defined on [1,∞) whose derivatives exist and are absolutely integrable. Then
∑1≤k≤Nf(k)=∫N1f(x)dx+12f(N)+Cf+112f′(N)−1720f(3)(N)+⋯! |
HN=∫N11xdx+12⋅1N+γ−112⋅1N2+O(1N4)=lnN+12N+γ−112N2+O(1N4) |
lnN!=lnN+ln(N−1)+⋯=∫N1lnxdx+12lnN+12ln2π+112N+O(1N3)=NlnN−N+ln√2πN+112N+O(1N3)Stirling's approximation=(N+12)lnN−N+ln√2π+O(1N) |
Source: Analysis of Algorithm