evan's notes
cal 3discreteodeenv systems

On this page

  • Domains and graphs
  • Example: combine two domain restrictions
  • Level curves and contour maps
  • Example: the contours of xy\sqrt{xy}xy​
  • Functions of three variables
  • From one-variable limits to two variables
  • Formal definition
  • A warning from a surface
  • Example: prove a limit from the definition

Intermediate Calculus

§10 Functions and Multivariable Limits

Evan Luo · Oct 16, 2026

Intermediate Calculus

§10 Functions and Multivariable Limits

Evan Luo3 days ago

6 min read

A function of several variables takes a point as its input. The first job is to find which points are allowed. We can then study the function through its graph, level sets, and limits.

Domains and graphs

For a function f:R2→Rf:\mathbb R^2\to\mathbb Rf:R2→R, the domain is the set of pairs (x,y)(x,y)(x,y) for which the formula makes sense. Its graph is the surface

z=f(x,y)z=f(x,y) z=f(x,y)

in R3\mathbb R^3R3.

Example: combine two domain restrictions

Consider

f(x,y)=1−x2−y2+ln⁡(y+12).f(x,y)=\sqrt{1-x^2-y^2}+\ln\left(y+\frac12\right). f(x,y)=1−x2−y2​+ln(y+21​).

The square root requires

1−x2−y2≥0⟺x2+y2≤1.1-x^2-y^2\ge0 \quad\Longleftrightarrow\quad x^2+y^2\le1. 1−x2−y2≥0⟺x2+y2≤1.

The logarithm requires a strictly positive input:

y+12>0⟺y>−12.y+\frac12>0 \quad\Longleftrightarrow\quad y>-\frac12. y+21​>0⟺y>−21​.

Both conditions must hold at once, so

Dom⁡(f)={(x,y)∈R2:x2+y2≤1 and y>−12}.\boxed{ \operatorname{Dom}(f) =\left\{(x,y)\in\mathbb R^2: x^2+y^2\le1\text{ and }y>-\frac12 \right\}. } Dom(f)={(x,y)∈R2:x2+y2≤1 and y>−21​}.​

This is the part of the closed unit disk above y=−1/2y=-1/2y=−1/2. The circular boundary is included wherever it also satisfies the logarithm condition. The horizontal chord at y=−1/2y=-1/2y=−1/2 is excluded.

The domain is an intersection of two restrictions

The domain is an intersection of two restrictionsThe included part of the unit disk lies strictly above y equals negative one half. The circular boundary is solid and included, while the horizontal chord is dashed and excluded.
xxx
yyy
y=−12y=-\frac12y=−21​
x2+y2≤1x^2+y^2\le1x2+y2≤1

Level curves and contour maps

A level curve of f(x,y)f(x,y)f(x,y) is a set of inputs where the function has one fixed value:

f(x,y)=c.f(x,y)=c. f(x,y)=c.

On the graph z=f(x,y)z=f(x,y)z=f(x,y), this comes from the horizontal trace at height z=cz=cz=c. The trace is a curve in the plane z=cz=cz=c; the level curve records the same (x,y)(x,y)(x,y) coordinates in the input plane. For more background on traces, see Distance and Quadric Surfaces.

A collection of level curves is called a contour map. It shows how the output changes without requiring a three-dimensional surface drawing.

Example: the contours of xy\sqrt{xy}xy​

Let

f(x,y)=xy.f(x,y)=\sqrt{xy}. f(x,y)=xy​.

The square root is real when xy≥0xy\ge0xy≥0. Thus the domain is

(x≥0 and y≥0)or(x≤0 and y≤0).(x\ge0\text{ and }y\ge0) \quad\text{or}\quad (x\le0\text{ and }y\le0). (x≥0 and y≥0)or(x≤0 and y≤0).

It consists of quadrants I and III together with the coordinate axes.

Because a principal square root is nonnegative, a level value must satisfy c≥0c\ge0c≥0. The equation f(x,y)=cf(x,y)=cf(x,y)=c becomes

xy=c⟺xy=c2.\sqrt{xy}=c \quad\Longleftrightarrow\quad xy=c^2. xy​=c⟺xy=c2.

For the first few levels:

c=0:x=0 or y=0,c=1:xy=1,c=2:xy=4.\begin{aligned} c=0&:\quad x=0\text{ or }y=0,\\ c=1&:\quad xy=1,\\ c=2&:\quad xy=4. \end{aligned} c=0c=1c=2​:x=0 or y=0,:xy=1,:xy=4.​

For c>0c>0c>0, the level curve is a hyperbola with one branch in quadrant I and one in quadrant III. Larger values of ccc lie farther from the axes.

Contour map for the square root of xy

Contour map for the square root of xyThe level zero is the pair of coordinate axes. Levels one and two are hyperbolas in quadrants one and three. The level two branches lie farther from the axes, and arrows point toward increasing function values.
xxx
yyy
000
111
222
higher z\text{higher }zhigher z

Functions of three variables

For a function f:R3→Rf:\mathbb R^3\to\mathbb Rf:R3→R, the full graph

w=f(x,y,z)w=f(x,y,z) w=f(x,y,z)

lives in R4\mathbb R^4R4, so we cannot draw it directly in ordinary space. We can still study the three-dimensional level surfaces

f(x,y,z)=c.f(x,y,z)=c. f(x,y,z)=c.

For example, take

f(x,y,z)=1000−x2−y2−z2.f(x,y,z)=1000-x^2-y^2-z^2. f(x,y,z)=1000−x2−y2−z2.

The level equation is

x2+y2+z2=1000−c.x^2+y^2+z^2=1000-c. x2+y2+z2=1000−c.

This gives

c=0:x2+y2+z2=1000,a sphere of radius 1000,c=500:x2+y2+z2=500,a sphere of radius 500,c=1000:x2+y2+z2=0,the single point (0,0,0).\begin{aligned} c=0&:\quad x^2+y^2+z^2=1000, &&\text{a sphere of radius }\sqrt{1000},\\ c=500&:\quad x^2+y^2+z^2=500, &&\text{a sphere of radius }\sqrt{500},\\ c=1000&:\quad x^2+y^2+z^2=0, &&\text{the single point }(0,0,0). \end{aligned} c=0c=500c=1000​:x2+y2+z2=1000,:x2+y2+z2=500,:x2+y2+z2=0,​​a sphere of radius 1000​,a sphere of radius 500​,the single point (0,0,0).​

The last level set is a degenerate sphere of radius 000. A level set need not be an ordinary surface; it can collapse to a point or even be empty.

From one-variable limits to two variables

For a one-variable function, the statement

lim⁡x→af(x)=L\lim_{x\to a}f(x)=L x→alim​f(x)=L

means that every requested output tolerance ε>0\varepsilon>0ε>0 can be met by choosing an input tolerance δ>0\delta>0δ>0 such that

0<∣x−a∣<δ⟹∣f(x)−L∣<ε.0<|x-a|<\delta \quad\Longrightarrow\quad |f(x)-L|<\varepsilon. 0<∣x−a∣<δ⟹∣f(x)−L∣<ε.

In two variables, replace the one-dimensional distance ∣x−a∣|x-a|∣x−a∣ with the Euclidean distance between (x,y)(x,y)(x,y) and (a,b)(a,b)(a,b):

∥(x,y)−(a,b)∥=(x−a)2+(y−b)2.\lVert(x,y)-(a,b)\rVert =\sqrt{(x-a)^2+(y-b)^2}. ∥(x,y)−(a,b)∥=(x−a)2+(y−b)2​.

Formal definition

Let f:D→Rf:D\to\mathbb Rf:D→R, where D⊆R2D\subseteq\mathbb R^2D⊆R2, and let (a,b)(a,b)(a,b) be an accumulation point of DDD. We write

lim⁡(x,y)→(a,b)f(x,y)=L\lim_{(x,y)\to(a,b)}f(x,y)=L (x,y)→(a,b)lim​f(x,y)=L

if, for every ε>0\varepsilon>0ε>0, there is a δ>0\delta>0δ>0 such that every (x,y)∈D(x,y)\in D(x,y)∈D satisfies

0<(x−a)2+(y−b)2<δ⟹∣f(x,y)−L∣<ε.\boxed{ 0<\sqrt{(x-a)^2+(y-b)^2}<\delta \quad\Longrightarrow\quad |f(x,y)-L|<\varepsilon. } 0<(x−a)2+(y−b)2​<δ⟹∣f(x,y)−L∣<ε.​

The input condition describes a punctured disk: points within radius δ\deltaδ of (a,b)(a,b)(a,b), with the center removed. The output condition requires the corresponding heights to lie in the band

L−ε<f(x,y)<L+ε.L-\varepsilon<f(x,y)<L+\varepsilon. L−ε<f(x,y)<L+ε.

A punctured input disk controls an output band

A punctured input disk controls an output bandA punctured disk of radius delta around the input point maps into the horizontal output band from L minus epsilon to L plus epsilon. The center is excluded because a limit does not use the function value at the center.Input plane
(a,b)(a,b)(a,b)
δ\deltaδ
fff
Output values
L+εL+\varepsilonL+ε
LLL
L−εL-\varepsilonL−ε
∣f(x,y)−L∣<ε|f(x,y)-L|<\varepsilon∣f(x,y)−L∣<ε

The puncture matters. A limit depends on nearby values, not on the value at the center. The function may have a different value at (a,b)(a,b)(a,b) or be undefined there.

A warning from a surface

The function

g(x,y)=x2−y2x2+y2g(x,y)=\frac{x^2-y^2}{x^2+y^2} g(x,y)=x2+y2x2−y2​

is undefined at (0,0)(0,0)(0,0). Its graph near the origin suggests that approaching from different directions may produce incompatible values, so a single limit should not be assumed from a picture. A proof of nonexistence must exhibit the contradiction formally; that argument comes after the definition-based example below.

Example: prove a limit from the definition

We will prove

lim⁡(x,y)→(0,0)3x2yx2+y2=0.\lim_{(x,y)\to(0,0)} \frac{3x^2y}{x^2+y^2}=0. (x,y)→(0,0)lim​x2+y23x2y​=0.

The expression is undefined at the origin, but the punctured-neighborhood condition excludes that point.

To determine how small δ\deltaδ must be for a given ε\varepsilonε, first bound the error in terms of the distance x2+y2\sqrt{x^2+y^2}x2+y2​. For (x,y)≠(0,0)(x,y)\ne(0,0)(x,y)=(0,0),

∣3x2yx2+y2−0∣=3∣y∣x2x2+y2≤3∣y∣≤3x2+y2.\begin{aligned} \left|\frac{3x^2y}{x^2+y^2}-0\right| &=3|y|\frac{x^2}{x^2+y^2}\\ &\le3|y|\\ &\le3\sqrt{x^2+y^2}. \end{aligned} ​x2+y23x2y​−0​​=3∣y∣x2+y2x2​≤3∣y∣≤3x2+y2​.​

The first inequality uses 0≤x2/(x2+y2)≤10\le x^2/(x^2+y^2)\le10≤x2/(x2+y2)≤1. The second uses y2≤x2+y2y^2\le x^2+y^2y2≤x2+y2.

Now let ε>0\varepsilon>0ε>0 and choose

δ=ε3.\delta=\frac{\varepsilon}{3}. δ=3ε​.

If

0<x2+y2<δ,0<\sqrt{x^2+y^2}<\delta, 0<x2+y2​<δ,

then

∣3x2yx2+y2−0∣≤3x2+y2<3δ=ε.\begin{aligned} \left|\frac{3x^2y}{x^2+y^2}-0\right| &\le3\sqrt{x^2+y^2}\\ &<3\delta\\ &=\varepsilon. \end{aligned} ​x2+y23x2y​−0​​≤3x2+y2​<3δ=ε.​

This is exactly the implication required by the definition. Therefore,

lim⁡(x,y)→(0,0)3x2yx2+y2=0.\boxed{ \lim_{(x,y)\to(0,0)} \frac{3x^2y}{x^2+y^2}=0. } (x,y)→(0,0)lim​x2+y23x2y​=0.​

Source: https://notes.ohevan.com/notes/intermediate-calculus/10-functions-and-multivariable-limits

© 2026 Evan Luo. All rights reserved.

Back to Intermediate Calculus

0 reads

·Last edited Today
  • newsletter

  • about me

  • sponsor

© 2026 Evan Luo. All rights reserved.