Unit 3, lesson 19

Implicit and Inverse Function Theorems

When an equation can be solved for some variables near a point, when a map has a local inverse, how to differentiate the solutions without finding them, and why the answers are only local.

  • Theory part 37
  • Practice part 38
  • Derivatives and Approximation

The unit circle is not the graph of one function, yet near the point (0.6,0.8)(0.6, 0.8) a small arc of it is the graph of y=1−x2y = \sqrt{1 - x^2}. Near (1,0)(1, 0) no such description works, because the circle turns vertical there. The implicit function theorem gives a simple test that guarantees the first kind of behavior, using one partial derivative at that point. The inverse function theorem answers the same kind of question for maps that send points to points.

Key terms

C1 regularity on a neighborhood

Definition Let WW be an open subset of Rk\mathbb{R}^k and F ⁣:W→RmF \colon W \to \mathbb{R}^m. The map FF is of class C1C^1 on WW when every first partial derivative of every component of FF exists at every point of WW and is continuous on WW. A map is C1C^1 near a point pp when it is C1C^1 on some open set containing pp. Such a map is differentiable at every point of WW, and its derivative is the Jacobian matrix DFDF.

In plain words The function has partial derivatives everywhere in a small open region around the point, and they change continuously. Both theorems in this lesson need this on a whole neighborhood, not only at one point.

Example Every polynomial, such as F(x,y)=x2+y2−1F(x, y) = x^2 + y^2 - 1, is C1C^1 on all of R2\mathbb{R}^2. The function H(x,y)=x+∣y∣H(x, y) = x + \lvert y \rvert is continuous everywhere, but ∂H/∂y\partial H / \partial y does not exist at any point with y=0y = 0, so HH is not C1C^1 on any open set that meets the xx-axis.

Implicit equations

Definition An implicit equation is an equation F(x,y)=0F(x, y) = 0, where FF is defined on an open set of Rn+m\mathbb{R}^{n + m} with values in Rm\mathbb{R}^m, x=(x1,…,xn)x = (x_1, \dots, x_n) and y=(y1,…,ym)y = (y_1, \dots, y_m). It is a system of mm scalar equations in n+mn + m unknowns. Its solution set is the set of all points (x,y)(x, y) with F(x,y)=0F(x, y) = 0, and the question is whether the mm variables yy can be regarded as functions of the nn variables xx on that set.

In plain words You have as many equations as unknowns you want to solve for, plus some leftover variables that are free to move. The equations tie the unknowns to the free variables without handing you a formula.

Example The single equation xz+y2+z3=3xz + y^2 + z^3 = 3 has n=2n = 2 free variables x,yx, y and m=1m = 1 unknown zz. The point (1,1,1)(1, 1, 1) lies on it, since 1+1+1=31 + 1 + 1 = 3. The pair of equations

xu+yv=2,u2−v+x−y=0\begin{gathered} xu + yv = 2, \\ u^2 - v + x - y = 0 \end{gathered}

has n=2n = 2 free variables x,yx, y and m=2m = 2 unknowns u,vu, v, and (x,y,u,v)=(1,1,1,1)(x, y, u, v) = (1, 1, 1, 1) satisfies both.

Local solvability

Definition Let F(a,b)=0F(a, b) = 0. The equation F(x,y)=0F(x, y) = 0 is locally solvable for yy near (a,b)(a, b) when there are open sets U∋aU \ni a in Rn\mathbb{R}^n and V∋bV \ni b in Rm\mathbb{R}^m and a function g ⁣:U→Vg \colon U \to V with g(a)=bg(a) = b such that, for (x,y)(x, y) in U×VU \times V,

F(x,y)=0  ⟺  y=g(x).F(x, y) = 0 \iff y = g(x).

The function gg is called a local branch of the solution set through (a,b)(a, b).

In plain words Inside a small box around the point, the solution set is exactly the graph of one function. Points of the solution set outside the box do not matter.

Example For x2+y2−1=0x^2 + y^2 - 1 = 0 near (0.6,0.8)(0.6, 0.8), take U=(0,1)U = (0, 1), V=(0,2)V = (0, 2) and g(x)=1−x2g(x) = \sqrt{1 - x^2}. Near (1,0)(1, 0) there is no local branch y=g(x)y = g(x), because for xx slightly less than 11 every box around the point contains the two solutions y=±1−x2y = \pm\sqrt{1 - x^2}, and for xx slightly more than 11 it contains none.

The implicit function theorem

Definition Let W⊆Rn+mW \subseteq \mathbb{R}^{n + m} be open, let F ⁣:W→RmF \colon W \to \mathbb{R}^m be C1C^1 on WW, and let (a,b)(a, b) be a point of WW with F(a,b)=0F(a, b) = 0. Write DyFD_y F for the m×mm \times m matrix of partial derivatives of the components of FF with respect to y1,…,ymy_1, \dots, y_m. If

det⁡DyF(a,b)≠0,\det D_y F(a, b) \ne 0,

then there are open sets U∋aU \ni a and V∋bV \ni b with U×V⊆WU \times V \subseteq W and a unique function g ⁣:U→Vg \colon U \to V with g(a)=bg(a) = b such that F(x,y)=0F(x, y) = 0 with (x,y)(x, y) in U×VU \times V holds exactly when y=g(x)y = g(x). Moreover gg is C1C^1 on UU, and UU can be chosen so small that DyF(x,g(x))D_y F(x, g(x)) is invertible for every xx in UU. For one equation (m=1m = 1) the condition is simply that the relevant partial derivative is nonzero, for example Fz(a,b)≠0F_z(a, b) \ne 0 when solving for zz.

In plain words Three things are needed. The function is C1C^1 near the point, the point solves the equation, and the derivative with respect to the unknowns you want is invertible there. Then the solution set near the point is the graph of a smooth function, even if you can never write it down.

Example Take F(x,y,z)=xz+y2+z3−3F(x, y, z) = xz + y^2 + z^3 - 3 at (1,1,1)(1, 1, 1). It is a polynomial, so it is C1C^1 on R3\mathbb{R}^3, and F(1,1,1)=0F(1, 1, 1) = 0. Since

Fz=x+3z2,Fz(1,1,1)=4≠0,\begin{gathered} F_z = x + 3z^2, \\ F_z(1, 1, 1) = 4 \ne 0, \end{gathered}

the theorem gives a C1C^1 function z=g(x,y)z = g(x, y) near (1,1)(1, 1) with g(1,1)=1g(1, 1) = 1. For the system above in the unknowns u,vu, v, the relevant Jacobian at (1,1,1,1)(1, 1, 1, 1) is

D(u,v)F=[xy2u−1]=[112−1],D_{(u,v)} F = \begin{bmatrix} x & y \\ 2u & -1 \end{bmatrix} = \begin{bmatrix} 1 & 1 \\ 2 & -1 \end{bmatrix},

with determinant −3≠0-3 \ne 0. So uu and vv are C1C^1 functions of (x,y)(x, y) near (1,1)(1, 1).

Derivatives of local branches

Definition Under the hypotheses of the implicit function theorem, differentiating the identity F(x,g(x))=0F(x, g(x)) = 0 with the chain rule gives DxF+DyF Dg=0D_x F + D_y F \, Dg = 0. Since DyFD_y F is invertible along the branch, for every xx in UU

Dg(x)=−[DyF]−1DxF,Dg(x) = -\big[ D_y F \big]^{-1} D_x F,

with both matrices evaluated at the point (x,g(x))(x, g(x)). Here DgDg is the m×nm \times n Jacobian of gg and DxFD_x F is the m×nm \times n matrix of partial derivatives of FF in x1,…,xnx_1, \dots, x_n.

For one equation F(x,y,z)=0F(x, y, z) = 0 solved for zz this reads

∂z∂x=−FxFz,∂z∂y=−FyFz,Fz≠0.\begin{gathered} \frac{\partial z}{\partial x} = -\frac{F_x}{F_z}, \qquad \frac{\partial z}{\partial y} = -\frac{F_y}{F_z}, \\ F_z \ne 0. \end{gathered}

In plain words You get the derivatives of the hidden function without solving for it. Differentiate the equations, then solve a small linear system for the rates you want.

Example On xz+y2+z3=3xz + y^2 + z^3 = 3 near (1,1,1)(1, 1, 1) we have Fx=zF_x = z, Fy=2yF_y = 2y and Fz=x+3z2F_z = x + 3z^2, so at that point

∂z∂x=−14,∂z∂y=−24=−12.\frac{\partial z}{\partial x} = -\frac{1}{4}, \qquad \frac{\partial z}{\partial y} = -\frac{2}{4} = -\frac{1}{2}.

For the system, differentiate both equations with respect to xx, holding yy fixed. At (1,1,1,1)(1, 1, 1, 1) this gives

1+ux+vx=0,2ux−vx+1=0.\begin{gathered} 1 + u_x + v_x = 0, \\ 2u_x - v_x + 1 = 0. \end{gathered}

Adding the two equations gives 3ux+2=03u_x + 2 = 0, so ux=−23u_x = -\frac{2}{3} and then vx=−13v_x = -\frac{1}{3}. The matrix formula gives the same values together with uy=0u_y = 0 and vy=−1v_y = -1.

The inverse function theorem

Definition Let U⊆RnU \subseteq \mathbb{R}^n be open, let f ⁣:U→Rnf \colon U \to \mathbb{R}^n be C1C^1 on UU, and let aa be a point of UU with

det⁡Df(a)≠0.\det Df(a) \ne 0.

Then there are open sets V∋aV \ni a and W∋f(a)W \ni f(a) with V⊆UV \subseteq U such that ff maps VV one-to-one onto WW, the inverse f−1 ⁣:W→Vf^{-1} \colon W \to V is C1C^1 on WW, and

D(f−1)(f(a))=[Df(a)]−1.D(f^{-1})(f(a)) = \big[ Df(a) \big]^{-1}.

The same formula holds at every point of WW, with aa replaced by the point of VV that ff sends there.

In plain words If the derivative at a point is an invertible matrix, then near that point the map itself can be undone smoothly, and the derivative of the undoing map is the inverse matrix.

Example Let f(u,v)=⟨u2−v2, 2uv⟩f(u, v) = \langle u^2 - v^2, \, 2uv \rangle, which is C1C^1 on R2\mathbb{R}^2. Its Jacobian has determinant

det⁡[2u−2v2v2u]=4(u2+v2),\det \begin{bmatrix} 2u & -2v \\ 2v & 2u \end{bmatrix} = 4(u^2 + v^2),

which is nonzero away from the origin. At (2,1)(2, 1) we get f(2,1)=(3,4)f(2, 1) = (3, 4) and

D(f−1)(3,4)=[4−224]−1=120[42−24].\begin{aligned} D(f^{-1})(3, 4) &= \begin{bmatrix} 4 & -2 \\ 2 & 4 \end{bmatrix}^{-1} \\[0.4em] &= \frac{1}{20} \begin{bmatrix} 4 & 2 \\ -2 & 4 \end{bmatrix}. \end{aligned}

So if x=u2−v2x = u^2 - v^2 and y=2uvy = 2uv, the local inverse near (3,4)(3, 4) has ∂u/∂x=15\partial u / \partial x = \frac{1}{5} and ∂v/∂x=−110\partial v / \partial x = -\frac{1}{10} there.

Local versus global inverses

Definition A map f ⁣:U→Rnf \colon U \to \mathbb{R}^n has a global inverse on UU when it is one-to-one on all of UU. The inverse function theorem gives only local inverses. A C1C^1 map with det⁡Df≠0\det Df \ne 0 at every point of an open set need not be one-to-one on that set. In the other direction, det⁡Df(a)=0\det Df(a) = 0 does not rule out an inverse, but if an inverse exists it cannot be differentiable at f(a)f(a).

In plain words An invertible derivative guarantees that the map can be undone near each point separately. Gluing those local inverses into one inverse on the whole domain can fail, because the map may wrap around and hit the same point twice.

Example The polar map f(r,θ)=⟨rcos⁡θ, rsin⁡θ⟩f(r, \theta) = \langle r \cos \theta, \, r \sin \theta \rangle on the open set r>0r > 0 has det⁡Df=r≠0\det Df = r \ne 0 everywhere, yet f(1,0)=f(1,2π)=(1,0)f(1, 0) = f(1, 2\pi) = (1, 0), so it has no global inverse. The squaring map above behaves the same way, since f(−u,−v)=f(u,v)f(-u, -v) = f(u, v). In one variable, f(x)=x3f(x) = x^3 is one-to-one with f′(0)=0f'(0) = 0, and its inverse x3\sqrt[3]{x} is not differentiable at 00.

Tangent planes to level surfaces

Definition Let FF be C1C^1 on an open set of R3\mathbb{R}^3, let F(p)=cF(p) = c and suppose ∇F(p)≠0\nabla F(p) \ne 0. Then some partial derivative of FF is nonzero at pp, so by the implicit function theorem the level set F=cF = c is, near pp, the graph of a C1C^1 function of two of the variables. Its tangent plane at p=(x0,y0,z0)p = (x_0, y_0, z_0) is

Fx(p)(x−x0)+Fy(p)(y−y0)+Fz(p)(z−z0)=0,\begin{aligned} F_x(p)(x - x_0) &+ F_y(p)(y - y_0) \\ &+ F_z(p)(z - z_0) = 0, \end{aligned}

and ∇F(p)\nabla F(p) is a normal vector. If ∇F(p)=0\nabla F(p) = 0 the theorem says nothing, and the level set may have no tangent plane there.

In plain words A nonzero gradient both proves that the level set is a smooth surface near the point and hands you its normal vector.

Example For F(x,y,z)=xz+y2+z3F(x, y, z) = xz + y^2 + z^3 at (1,1,1)(1, 1, 1), the gradient is ⟨z,2y,x+3z2⟩=⟨1,2,4⟩\langle z, 2y, x + 3z^2 \rangle = \langle 1, 2, 4 \rangle. The tangent plane is

(x−1)+2(y−1)+4(z−1)=0,(x - 1) + 2(y - 1) + 4(z - 1) = 0,

which simplifies to x+2y+4z=7x + 2y + 4z = 7. Solving for zz gives slopes −14-\frac{1}{4} and −12-\frac{1}{2}, matching the branch derivatives found earlier. By contrast the cone x2+y2−z2=0x^2 + y^2 - z^2 = 0 has zero gradient at the origin, and it has no tangent plane at its tip.

Check yourself

  1. Show that xey+y=1x e^y + y = 1 defines yy as a C1C^1 function of xx near (1,0)(1, 0), and find dydx\frac{dy}{dx} there.
Show the answer

The function F(x,y)=xey+y−1F(x, y) = x e^y + y - 1 is C1C^1 on R2\mathbb{R}^2 and F(1,0)=0F(1, 0) = 0. Its partial Fy=xey+1F_y = x e^y + 1 equals 2≠02 \ne 0 at (1,0)(1, 0), so the theorem applies. Then

dydx=−FxFy=−eyxey+1=−12.\frac{dy}{dx} = -\frac{F_x}{F_y} = -\frac{e^y}{x e^y + 1} = -\frac{1}{2}.
  1. For f(u,v)=⟨u2−v2, 2uv⟩f(u, v) = \langle u^2 - v^2, \, 2uv \rangle near (u,v)=(2,1)(u, v) = (2, 1), find ∂u/∂y\partial u / \partial y and ∂v/∂y\partial v / \partial y of the local inverse at (x,y)=(3,4)(x, y) = (3, 4). Why is there no local inverse near (0,0)(0, 0)?
Show the answer

The second column of the inverse matrix found above gives ∂u/∂y=220=110\partial u / \partial y = \frac{2}{20} = \frac{1}{10} and ∂v/∂y=420=15\partial v / \partial y = \frac{4}{20} = \frac{1}{5}. At the origin the Jacobian determinant 4(u2+v2)4(u^2 + v^2) is 00, so the theorem does not apply. In fact no local inverse exists, because every disk around (0,0)(0, 0) contains a point (u,v)≠(0,0)(u, v) \ne (0, 0) together with (−u,−v)(-u, -v), and these two different points have the same image.

  1. Find the tangent plane to x2+2y2+3z2=6x^2 + 2y^2 + 3z^2 = 6 at (1,1,1)(1, 1, 1).
Show the answer

The gradient of F=x2+2y2+3z2F = x^2 + 2y^2 + 3z^2 is ⟨2x,4y,6z⟩=⟨2,4,6⟩≠0\langle 2x, 4y, 6z \rangle = \langle 2, 4, 6 \rangle \ne 0. The tangent plane is 2(x−1)+4(y−1)+6(z−1)=02(x - 1) + 4(y - 1) + 6(z - 1) = 0, which simplifies to x+2y+3z=6x + 2y + 3z = 6.

How this lesson connects

It also uses Partial derivatives and C1 regularity, Chain rule and Jacobian matrices and matrix inverses.