Put T(M,N)=βmβ€Mββnβ€Nβam,nβ (a finite sum, β₯0) and S=supM,NβT(M,N)β[0,β].
T is monotone in each of M,N.
For each fixed m, the inner series βn=1ββam,nββ[0,β] converges (finitely or to +β) since the terms are nonnegative.
Put R=βm=1βββn=1ββam,nβ.
For each M,N: T(M,N)=βmβ€Mββnβ€Nβam,nββ€βmβ€Mββn=1ββam,nββ€R (inner finite sum β€ inner infinite sum, outer finite β€ outer infinite β all terms nonnegative).
Taking the supremum over M,N gives Sβ€R.
The other direction: for each M, βmβ€Mββn=1ββam,nβ=supNβT(M,N)β€S (the inner infinite sum is the supremum of the inner finite sums).
This equality needs justification: for fixed M, each inner sum βnβ€Nβam,nβ is increasing in N (since am,nββ₯0) and converges to βn=1ββam,nβ.
For a finite sum (M terms), supremum and summation interchange β a finite sum of increasing sequences is increasing, and its supremum equals the sum of the suprema (checked termwise).
Hence βmβ€MβsupNββnβ€Nβam,nβ=supNββmβ€Mββnβ€Nβam,nβ=supNβT(M,N).
Taking the supremum over M gives R=supMββmβ€Mββn=1ββam,nββ€S.
Hence R=S.
Swapping the roles of rows and columns, the identical argument gives C=βn=1βββm=1ββam,nβ=S (the same supremum of the same T(M,N)).
Hence R=C=S β the two iterated sums agree.
If S<β the common value is finite; if S=β both sides diverge to +β.
QED.