Throughout, Ξ£ denotes the set of strictly increasing functions NβN, a subset of the set of all functions NβN; it contains the identity id and is closed under composition by Monotone Sequences and Subsequences: Comparison of All Terms, Growth of the Indices, and Subsequences of Subsequences (iii), and for ΟβΞ£ we use kβ€Ο(k) and the equivalence of Ο(j)<Ο(k) with j<k, from Monotone Sequences and Subsequences: Comparison of All Terms, Growth of the Indices, and Subsequences of Subsequences (ii).
Step 1: a canonical convergent subsequence. Let b=(bnβ) be a bounded sequence in R and let Lbβ=limsupnβββbnβ be its limit superior. We first show: for every real Ξ΅>0 and every uβN there is nβN with u<n and β£bnββLbββ£<Ξ΅. By Properties of the Limit Superior and the Limit Inferior: Tail Formula, Order, Comparison, Bounds, Eventual and Frequent Inequalities, Subsequences, Negation and Convergence (v) there is NβN with bnβ<Lbβ+Ξ΅ for every nβ₯N. By Properties of the Limit Superior and the Limit Inferior: Tail Formula, Order, Comparison, Bounds, Eventual and Frequent Inequalities, Subsequences, Negation and Convergence (vi), applied with the index N+u, there is nβ₯N+u with LbββΞ΅<bnβ. Then nβ₯N and n>u, so LbββΞ΅<bnβ<Lbβ+Ξ΅, as claimed.
Hence, for j,uβN, the set
Ebβ(j,u)={nβN:u<nΒ andΒ β£bnββLbββ£<1/j}
is nonempty and so has a least element by Arithmetic and Order of the Natural Numbers (xv); call it ΞΌbβ(j,u). By Recursion on the Natural Numbers Starting at One (i) there is exactly one function Οbβ:NβN with
Οbβ(1)=ΞΌbβ(1,1),Οbβ(j+1)=ΞΌbβ(j+1,Οbβ(j))(jβN).
By the definition of Ebβ, Οbβ(j)<Οbβ(j+1), so ΟbββΞ£, and β£bΟbβ(j)ββLbββ£<1/j for every jβN (for j=1 directly, for a successor j from the recursion). Since 1/jβ0 by Completeness of the Real Numbers for Sequences: Monotone Convergence, the Bolzano-Weierstrass Theorem and Cauchy Sequences (i), the subsequence bβΟbβ converges to Lbβ by Limits of Sequences of Real Numbers: Uniqueness, Boundedness, Constants, Tails, Arithmetic, Quotients, Absolute Values, Finite Sums, Order, Squeezing, Domination and Subsequences (xii). The function Οbβ is determined by b alone.
Step 2: nested index sequences. For kβN and ΟβΞ£ let ck,Ο be the sequence iβ¦aΟ(i),kβ. It is bounded: if β£am,kββ£β€M for every m, then β£aΟ(i),kββ£β€M for every i. So, by Step 1, (k,Ο)β¦ΟβΟck+1,Οβ is a function NΓΞ£βΞ£, and Recursion on the Natural Numbers Starting at One (i) gives a sequence (Οkβ)kβNβ in Ξ£ with
Ο1β=Οc1,idβ,Οk+1β=ΟkββΟkβ,whereΒ Οkβ=Οck+1,ΟkβββΞ£.
For every kβN the sequence (aΟkβ(i),kβ)iβNβ converges: for k=1 it is c1,idβΟc1,idβ, and for k=l+1 it is cl+1,ΟlββΟcl+1,Οlββ, both convergent by Step 1. Let Akβ be its limit.
Step 3: nesting. We show by induction on jβN: for every kβ€j there is ΟβΞ£ with Οjβ=ΟkββΟ. For j=1 the only k is 1, and Ο=id works. Suppose the assertion holds for j, and let kβ€j+1. If k=j+1, take Ο=id. Otherwise kβ€j, and with Ο from the hypothesis, Οj+1β=ΟjββΟjβ=Οkββ(ΟβΟjβ) with ΟβΟjββΞ£.
Step 4: the diagonal. Put mjβ=Οjβ(j) for jβN. Since j<j+1β€Οjβ(j+1), we get
mjβ=Οjβ(j)<Οjβ(Οjβ(j+1))=Οj+1β(j+1)=mj+1β,
so (mjβ) is strictly increasing. Fix kβN; we show that (amjβ,kβ)jβ converges to Akβ. Let Ξ΅>0. First choose, by Convergent Sequences of Real Numbers (i) and Step 2, IβN with β£aΟkβ(i),kββAkββ£<Ξ΅ for every iβ₯I; then put J=k+I. Let jβ₯J. Then kβ€j, so Step 3 gives ΟβΞ£ with Οjβ=ΟkββΟ, whence mjβ=Οkβ(Ο(j)); and Ο(j)β₯jβ₯I. Therefore
β£amjβ,kββAkββ£=β£aΟkβ(Ο(j)),kββAkββ£<Ξ΅(jβ₯J).
So (amjβ,kβ)jβ converges to Akβ by Convergent Sequences of Real Numbers (i). As k was arbitrary, (mjβ) is as required by (i).