$$ \sum_{i=1}^{\infty}\left[\sum_{j=i}^{\infty}f(i,j)\right]=\sum_{j=1}^{\infty}\left[\sum_{i=1}^{j}f(i,j)\right] $$
I have trouble understanding the change of order of this summation. I know that $$ 1≤i≤j≤\infty $$ however, I can't seem to make sense of this even with manually comparing each sum of both double summations.
In the left side, we always have the sums: $$ S_1=f(1,1)+\dots f(1,\infty) $$ $$ S_2=f(2,2)+\dots f(2,\infty) $$ $$ \vdots $$ so the total sum
$$ S=\sum S_n $$ looks like summing each element of an $N\text{x}N$ matrix as $N\to\infty$
In the right side, we have the sums:
$$ S_1=f(1,1) $$ $$ S_2=f(1,2) + f(2,2) $$ $$ S_3=f(1,3)+f(2,3) + f(3,3) $$ $$ \vdots $$ The total sum $S$ looks like summing the top-right to bottom-left entries of a "pyramid" structure, which looks like the earlier matrix but deformed and slanted. However, I do not want to just assume that this is the case since I am working with infinite sums.
Is there a better intuition on how changing order of summations work so that I may extend it to more difficult scenarios particularly in the proof of the convolution property of Z transforms.