Put Anβ=βk=1nβakβ.
By hypothesis AnββA.
Set A0β:=0 so the formula below is also valid for m=1 (with Amβ1β=A0β=0).
Summation by parts (Abel's transform) gives, for mβ€n, βk=mnβakβbkβ=Anβbn+1ββAmβ1βbmβ+βk=mnβAkβ(bkββbk+1β), which is checked directly by substituting akβ=AkββAkβ1β and collecting the telescoping terms.
Write Akβ=A+Ξ΅kβ with Ξ΅kββ0. Then βk=mnβAkβ(bkββbk+1β) splits as A(bmββbn+1β)+βk=mnβΞ΅kβ(bkββbk+1β), the first part telescoping.
The key is the monotonicity of (bkβ): the differences have constant sign, so βk=mnββ£bkββbk+1ββ£=β£bmββbn+1ββ£.
Boundedness β£bkββ£β€M makes this at most 2M.
Hence β£βk=mnβΞ΅kβ(bkββbk+1β)β£β€maxmβ€kβ€nββ£Ξ΅kββ£β
2M.
For any Ξ·>0, taking m large makes this smaller than Ξ·β
2M.
Now the boundary term Anβbn+1ββAmβ1βbmβ.
Since bkββb (a monotone bounded sequence converges), this tends to AbβAb=0 in the Cauchy sense as m,nββ.
The leftover term A(bmββbn+1β) must also be estimated: β£A(bmββbn+1β)β£β€β£Aβ£β
β£bmββbn+1ββ£, and since bkββb, β£bmββbn+1ββ£β0 as m,nββ (crudely β£bmββbn+1ββ£β€2M).
Hence this term also tends to 0.
Every term in the decomposition is now estimated, so together with the previous estimate, the tails of βk=mnβakβbkβ become arbitrarily small, so the partial sums of βakβbkβ form a Cauchy sequence.
Hence the series converges.