First, the integrals exist: a monotone function on [a,b] is Riemann integrable (standard theorem), so the decreasing f is integrable on every [n,n+1] and every β«nn+1βf exists.
Since f is decreasing, f(n+1)β€f(x)β€f(n) for xβ[n,n+1].
Integrating over [n,n+1] (length 1) gives f(n+1)β€β«nn+1βfβ€f(n) by monotonicity of the integral.
Sum over n=1,β¦,N.
The left inequalities sum to βn=1Nβf(n+1)=βk=2N+1βf(k)β€β«1N+1βf; the right ones to β«1N+1βfβ€βn=1Nβf(n).
With SNβ=βn=1Nβf(n), the correct identity is βk=2N+1βf(k)=SN+1ββf(1)=SNββf(1)+f(N+1).
Hence the left side gives SN+1ββf(1)β€β«1N+1βf, i.e. SN+1ββ€f(1)+β«1N+1βf, and since f>0, SNββ€SN+1β, so SNββ€f(1)+β«1N+1βf as well.
Chaining with the right side gives β«1N+1βfβ€βn=1Nβf(n)β€f(1)+β«1N+1βf.
(β) If β«1ββf converges (i.e. β«1Rβf has a finite limit as Rββ), then the partial sums βn=1Nβf(n) are bounded above by f(1)+β«1ββf and, the terms being nonnegative, converge by the monotone convergence theorem.
(β) By contraposition: if β«1ββf diverges (to +β, the integrand being positive), then βn=1Nβf(n)β₯β«1N+1βfββ, so the series diverges too.
Hence βf(n) converges βΊ β«1ββf converges.
QED.