# Commutativity of iterated limits

The following is a weird result I've obtained with iterated limits. There must be a flaw somewhere in someone's reasoning but I can't discover what it is. The problem is that, in general, iterated limits are not supposed to be commutative. However, the purported proof below seems to indicate that they always are. Did I make a mistake somewhere?

Let $F$ be a real valued function of two real variables defined in some region around $(a,b)$. Then the standard limit [ correct? ] of $F$ as $(x,y)$ approaches $(a,b)$ equals $L\,$ if and only if for every $\epsilon > 0$ there exists $\delta > 0$ such that $F$ satisfies: $$| F(x,y) - L | < \epsilon$$ whenever the distance between $(x,y)$ and $(a,b)$ satisfies: $$0 < \sqrt{ (x-a)^2 + (y-b)^2 } < \delta$$ We will use the following notation for such limits of functions of two variables: $$\lim_{(x,y)\rightarrow (a,b)} F(x,y) = L$$ Note. At this moment, we do not wish to consider limits where (one of) the independent variable(s) approaches infinity.

Next we consider the following iterated limit: $$\lim_{y\rightarrow b} \left[ \lim_{x\rightarrow a} F(x,y) \right] = L$$ Theorem. Commutativity of iterated limits. $$\lim_{y\rightarrow b} \left[ \lim_{x\rightarrow a} F(x,y) \right] = \lim_{x\rightarrow a} \left[ \lim_{y\rightarrow b} F(x,y) \right] = \lim_{(x,y)\rightarrow (a,b)} F(x,y)$$ Proof. We split the first iterated limit in two pieces: $$\lim_{x\rightarrow a} F(x,y) = F_a(y)$$ And: $$\lim_{y\rightarrow b} F_a(y) = L$$ Thus it becomes evident that the (first) iterated limit is actually defined as follows.
For every $\epsilon_x > 0$ there is some $\delta_x > 0$ such that: $$| F(x,y) - F_a(y) | < \epsilon_x \quad \mbox{whenever} \quad 0 < | x - a | < \delta_x$$ For every $\epsilon_y > 0$ there is some $\delta_y > 0$ such that: $$| F_a(y) - L | < \epsilon_y \quad \mbox{whenever} \quad 0 < | y - b | < \delta_y$$ Applying the triangle inequality $|a| + |b| \ge |a + b|$ gives: $$| F(x,y) - F_a(y) | + | F_a(y) - L | \ge | F(x,y) - L |$$ Consequently: $$| F(x,y) - L | < \epsilon_x + \epsilon_y$$ On the other hand we have: $$0 < | x - a | < \delta_x \qquad \mbox{and} \qquad 0 < | y - b | < \delta_y$$ Hence: $$0 < \sqrt{ (x-a)^2 + (y-b)^2 } < \sqrt{ \delta_x^2 + \delta_y^2 }$$ This is exactly the definition of the above standard limit of a function of two variables if we put: $$\epsilon = \epsilon_y + \epsilon_y \qquad \mbox{and} \qquad \delta = \sqrt{\delta_x^2 + \delta_y^2}$$ Therefore: $$\lim_{y\rightarrow b} \left[ \lim_{x\rightarrow a} F(x,y) \right] = \lim_{(x,y)\rightarrow (a,b)} F(x,y)$$ In very much the same way we can prove that: $$\lim_{x\rightarrow a} \left[ \lim_{y\rightarrow b} F(x,y) \right] = \lim_{(x,y)\rightarrow (a,b)} F(x,y)$$ QED

• Before attempting to answer your question, it's worth noting that the current statement of your claim (stated as a "Theorem") could definitely be made clearer. Which of the three limits are you assuming converge? A far more precise formulation would be something like the following: "Claim: If both $\lim_{x\to a}\lim_{y\to b} F(x,y)$ and $\lim_{y\to b}\lim_{x\to a}F(x,y)$ exist and equal $L$, then $\lim_{(x,y)\to (a,b)}F(x,y)$ exists and equals $L$." I think that this is the result you're attempting to prove but can't be $100$% sure without you modifying your claim.
– Dan
Nov 1, 2013 at 21:22
• Somewhat weaker: if $\lim_{x\rightarrow a} \lim_{y \rightarrow b} F(x,y)$ and $\lim_{x\rightarrow b} \lim_{y \rightarrow a} F(x,y)$ and $\lim_{(x,y)\rightarrow (a,b)} F(x,y)$ do exist and that's all. Nov 2, 2013 at 11:08
– Dan
Nov 2, 2013 at 16:05

Great question! The error in your proof is somewhat subtle.

Let's assume that $F_a(y) = \lim_{x\to a}F(x,y)$ exists for all $y$ sufficiently close to $b$. Let $\epsilon > 0$. Then by definition there exists $\delta(y) > 0$ possibly dependent on $y$ such that

$$|x - a| < \delta(y) \implies |F(x,y) - F_a(y)| < \epsilon/2.$$

Assume further that $L = \lim_{y\to b}F_a(y)$ exists. Then there exists $\delta > 0$ such that $$|y - b| < \delta \implies |F_a(y) - L| < \epsilon/2.$$

So if $|y-b| < \delta$ and $|x-a| < \delta(y)$, then we do in fact have $$|F(x,y) - L| \leq \epsilon.$$ However, this is only for $|y-b| < \delta$ and $|x-a| < \delta(y)$, which (depending on the function $\delta(y)$) may not be a neighborhood of $(x,y)$ in $\mathbb{R}^2$. Accordingly, this is not the same as saying $\lim_{(x,y)\to(a,b)}F(x,y) = L$. Similar statements can be made about the second half of the argument.

Hope this makes sense!

Edit: Thanks for clarifying your assumptions in your comment. Good news: you're right! In fact, all you need to assume is that $\lim_{(x,y)\to(a,b)}F(x,y)$ exists. Intuitively, this is because $\lim_{(x,y)\to(a,b)}F(x,y) = L$ means that $F(x,y) \to L$ as $(x,y)$ approaches $(a,b)$ along any path, and you can think of the iterated limits as limits of $F(x,y)$ as $(x,y)$ approaches $(a,b)$ along particular paths. So the assumption $F(x,y) \to L$ as $(x,y) \to (a,b)$ includes the statements $\lim_{x\to a}\lim_{y\to b}F(x,y) = L$ and $\lim_{y\to b}\lim_{x\to a}F(x,y)$ as special cases.

A precise proof of this claim is probably easiest to write by starting with the definition of the limit $L = \lim_{(x,y) \to (a,b)}F(x,y)$ and showing (in a similar fashion to what you were doing) that the iterated limits both equal $L$.

Edit 2: This is in response to your last comment:

Assume for example that a>0 and b>0 and construct two paths e.g. from the origin $(0,0)$ to the limit point $(a,b)$. Let $(x_1(t),y_1(t))=(t,0)$ for $0\leq t\leq a$ and $=(a,t−a)$ for $a≤t≤a+b$. Define $G_1(t)=F(x_1(t),y_1(t))$ and the "path limit" $\lim_{t\to a+b}G_1(t)$. Let $(x_2(t),y_2(t))=(0,t)$ for $0\leq t\leq b$ and $=(t−b,b)$ for $b\leq t\leq a+b$. Define $G_2(t)=F(x_2(t),y_2(t))$ and the "path limit" $\lim_{t\to a+b}G_2(t)$. Then I agree with you that, in general, $\lim_{t\to a+b}G_1(t)\neq \lim_{t\to a+b}G_2(t)$. But these are not iterated limits.

It doesn't matter that they're not iterated limits. If you reread my last comment, my point is that even if the iterated limits converge to the same value, i.e., $\lim_{x\to a}\lim_{y\to b}F(x,y) = L = \lim_{y\to b}\lim_{x\to a}F(x,y)$, it may very well be the case that the limit $\lim_{(x,y)\to(a,b)}F(x,y)$ doesn't exist, since it might not converge to $L$ along other paths to $(a,b)$ (which by the definition of $\lim_{(x,y)\to(a,b)}F(x,y)$ must all converge to the same value). A standard example of this is $F(x,y) = \frac{xy}{x^2+y^2}$ as $(x,y) \to (0,0)$. It's easy to see that $$\lim_{x\to 0}\lim_{y\to0}\frac{xy}{x^2+y^2} = \lim_{y\to 0}\lim_{x\to 0}\frac{xy}{x^2+y^2} = 0.$$ In fact $F(x,0) = 0 = F(0,y)$ for all $x,y \neq 0$, and so every neighborhood of $(0,0)$ contains a point $(x,y)$ at which $F(x,y) = 0$.

However, the limit $\lim_{(x,y)\to(0,0)}\frac{xy}{x^2+y^2}$ does not exist. For example, along the line $y=x$, the function equals $$F(x,x) = \frac{x^2}{2x^2} = \frac{1}{2}.$$ So every neighborhood of $(0,0)$ contains a point $(x,y)\neq (0,0)$ at which $F(x,y) = 1/2$. (For any $\delta > 0$, choose $(x,y) = (x,x)$ where $|x| < \delta/\sqrt2$. Then $\|(x,x)-(0,0)\| = \sqrt{2x^2} < \delta$, and we've shown $F(x,x) = 1/2$.) We've shown that $F(x,y)$ takes on the values $0$ and $1/2$ in every neighborhood of $(0,0)$, and so $\lim_{(x,y)\to(0,0)}F(x,y)$ does not exist.

• Thanks! But what if we only assume that the iterated limits $\lim_{x\to a}\lim_{y\to b}F(x,y)$ and $\lim_{y\to b}\lim_{x\to a}F(x,y)$ both exist? That's the other way around. Nov 8, 2013 at 18:13
• @HandeBruijn: Then the statement is false, even if we make the additional assumption that these two limits are equal. Intuitively: Just because a function approaches $L$ when moving on these two rectilinear paths to $(a,b)$ doesn't mean it will approach $L$ along all paths to $(a,b)$.
– Dan
Nov 8, 2013 at 20:19
• Assume for example that $a>0$ and $b>0$ and construct two paths e.g. from the origin $(0,0)$ to the limit point $(a,b)$. Let $(x_1(t),y_1(t)) = (t,0)$ for $0 \le t \le a$ and $= (a,t-a)$ for $a \le t \le a+b$. Define $G_1(t) = F(x_1(t),y_1(t))$ and the "path limit" $\lim_{t \to a+b} G_1(t)$. Let $(x_2(t),y_2(t)) = (0,t)$ for $0 \le t \le b$ and $= (t-b,b)$ for $b \le t \le a+b$. Define $G_2(t) = F(x_2(t),y_2(t))$ and the "path limit" $\lim_{t \to a+b} G_2(t)$. Then I agree with you that, in general, $\lim_{t \to a+b} G_1(t) \ne \lim_{t \to a+b} G_2(t)$. But these are not iterated limits. Nov 10, 2013 at 10:37
– Dan
Nov 10, 2013 at 17:59
• Yes, and another (quite similar and quite standard) example is $f(x,y) = (x^2 - y^2)/(x^2 + y^2)$ where the iterated limits are supposed to be $+1$ and $-1$ and the standard limit for $(x,y) \to (0,0)$ doesn't exist either. We all know that. Nov 10, 2013 at 19:23

There is a property called "uniform continuity" and is stronger than "side-way" continuity. Consider the sequence of functions $$f_n(x)=n^2xe^{-nx}, \; x>0,\; n=1,2,\dots$$ It's easy to show that for each $n$, the maximum of $f_n(x)$ over all $x\ge0$ is $\;n/e$ and is achieved at $x_n=1/n$. Also for any fixed $x\ge0$, $f_n(x) \to 0$ as $n \to \infty$. The notation $\lim_{n \to \infty}f_n(x)=0$ is somewhat misleading in this context. How can we have that and still have a sequence $x_n$ such that $f_n(x_n)=n/e$ ? The definition goes: $\forall \epsilon \gt 0, \exists \; N \gt 0$ such that $f_n(x) \lt \epsilon \; \forall \;n \gt N$. But this $N$ may be (and it is in this example) a (unbounded) function of $x$. The definition for uniform convergence at say $\bar{x}$ goes: $\forall \epsilon \gt 0, \exists \; N \gt 0$ and an $\alpha \gt 0$ such that $f_n(x) \lt \epsilon \; \forall \;n \gt N$ and $\forall x \in(\bar{x}- \alpha,\bar{x}+\alpha)$. Here we have continuity and uniform continuity everywhere but we don't have uniform continuity at $x=0$. To cut the long story short, if you have the existence of a "side-ways" limit over e.g. $x$ and uniform continuity in $y$ then you have it all. Similarly starting with $y$.

• Your answer (+1) is very much appreciated. Apr 16, 2014 at 13:20

Rather a comment than an answer. The only thing missing in Georgy's answer is an integral: $$\int_0^\infty n^2xe^{-nx} dx = -\frac{1}{n} \int_0^\infty n^2x\, d e^{-nx} = \left[ - \frac{n^2x e^{-nx}}{n} \right]_0^\infty + \int_0^\infty n e^{-nx} dx = 0-\left[e^{-nx}\right]_0^\infty = 1$$ And a picture, where the grey area is independent of $\,n$ , i.e. the one just calculated: It is seen that Georgy's function converges to a Dirac delta at $x=0$.
This is similar to the behaviour of a simpler function at : Iterated Limits Schizophrenia .