The Inverse Function Theorem

Suppose a smooth map F:Rn→RnF:\mathbb{R}^n \to \mathbb{R}^n has a derivative at some point aa. Suppose that derivative is an invertible linear map. Must FF itself be invertible near aa? At least locally, with a differentiable inverse? This question sits at the heart of multivariable calculus. The answer — yes, under the right hypotheses — is the Inverse Function Theorem. It is the rigorous engine behind an intuitive idea. Namely: "if the linear approximation is invertible, so is the real thing, at least close up." The theorem matters because most interesting maps are not linear. Yet linear algebra is the only place where invertibility is easy to check (a nonzero determinant). The Inverse Function Theorem builds a bridge. It lets us transport the algebraic fact "DF(a)DF(a) is invertible." This becomes the analytic fact "FF is a local diffeomorphism near aa." This single idea underlies the The Implicit Function Theorem. It also underlies coordinate changes in Change of Variables and the Jacobian. It shapes the local structure of smooth manifolds. And it drives Newton's method for solving nonlinear systems. The big picture is a linearization principle. Near a point, a C1C^1 map behaves, to first order, like its Jacobian derivative. This is the content of Jacobian Matrix and Linear Approximation. The Inverse Function Theorem says this first-order resemblance is strong. Strong enough to force a genuine local homeomorphism (in fact diffeomorphism). This holds provided the derivative is nondegenerate. The proof is a beautiful application of the contraction mapping principle. One constructs the inverse as the fixed point of an explicit iteration. Then one bootstraps continuity into differentiability. This last step uses the The Chain Rule in Several Variables. This article gives the full statement. It includes a complete proof via the contraction mapping method. It offers two worked examples with genuinely asymmetric Jacobians. It presents a counterexample showing why continuity of the derivative cannot be dropped. Finally, it discusses where the theorem is used across mathematics, physics, and computation.

Key takeaways
  • If FF is C1C^1 near aa and DF(a)DF(a) is invertible, then FF is a local C1C^1-diffeomorphism near aa.
  • The local inverse GG satisfies DG(F(a))=[DF(a)]−1DG(F(a)) = [DF(a)]^{-1} — derivatives of inverse maps are matrix inverses.
  • The proof constructs GG as a fixed point of a contraction, built from FF shifted and linearized.
  • Continuity of the derivative (C1C^1, not merely differentiable) is essential — a genuine counterexample exists without it.
  • The theorem is purely local: global invertibility is a separate, much harder question.
  • Key idea: a nonsingular linear approximation at a point forces genuine local invertibility of the nonlinear map itself.

Where this lives today

The Inverse Function Theorem is a load-bearing pillar across mathematics and its applications. In pure analysis, it is the direct ancestor of the The Implicit Function Theorem, which recasts "solving F(x,y)=0F(x,y)=0 for yy as a function of xx" as a special case of local invertibility. In Change of Variables and the Jacobian, the theorem certifies that a substitution used in a multiple integral is locally a genuine coordinate change, so that the change-of-variables formula with the Jacobian determinant is legitimate. In differential geometry, it is the tool that shows a map between manifolds with nonsingular differential is a local diffeomorphism, underlying the definition of immersions, submersions, and regular values (Sard's theorem territory). In dynamical systems, local invertibility near fixed points is used to build local coordinate charts (linearization near hyperbolic fixed points, via the Hartman–Grobman theorem, is philosophically downstream of this idea). In numerical analysis and optimization, Newton's method for solving nonlinear systems F(x)=0F(x)=0 relies on the nonsingularity of the Jacobian at the root to guarantee local convergence — a quantitative, algorithmic descendant of the same contraction-mapping argument used here. In physics, changes between coordinate systems (Cartesian to polar, spherical, or curvilinear coordinates in general relativity and fluid dynamics) are legitimized locally by exactly this theorem, since a nonzero Jacobian determinant is what allows physical quantities to be re-expressed consistently in new coordinates. In economics and statistics, local invertibility of demand or transformation functions (e.g. the change-of-variables formula for probability densities under a smooth bijective transformation) again rests on this same nonsingular-Jacobian criterion.

:::quiz[Given the conclusion DG(y)=[DF(G(y))]−1DG(y) = [DF(G(y))]^{-1}, which fact guarantees DGDG is itself continuous, so that G∈C1G \in C^1?]

  • GG is a bijection, and all bijections have continuous derivatives
  • Continuity of GG, continuity of DFDF, and continuity of matrix inversion on invertible matrices, composed together ✓
  • The chain rule alone, with no further input
  • The fact that det⁡DF(a)≠0\det DF(a) \neq 0 implies DFDF is constant near aa

Ask the tutor: How does the proof of the Inverse Function Theorem change if FF maps between infinite-dimensional Banach spaces instead of Rn\mathbb{R}^n?

DefinitionLocal C1C^1-diffeomorphism and the setup

Let U⊆RnU \subseteq \mathbb{R}^n be open and let F:U→RnF: U \to \mathbb{R}^n be a function. We say FF is continuously differentiable (C1C^1) on UU if all partial derivatives ∂Fi/∂xj\partial F_i/\partial x_j exist and are continuous on UU; equivalently the Jacobian matrix map

DF(x)=[∂Fi∂xj(x)]i,j=1nDF(x) = \left[\frac{\partial F_i}{\partial x_j}(x)\right]_{i,j=1}^n

is a continuous function of x∈Ux \in U. A point a∈Ua \in U is a regular point of FF if DF(a)DF(a) is an invertible n×nn\times n matrix, i.e. det⁡DF(a)≠0\det DF(a) \neq 0. FF is a local C1C^1-diffeomorphism at aa if there exist open sets V∋aV \ni a and W∋F(a)W \ni F(a) such that F:V→WF: V \to W is a bijection, and both F∣VF|_V and its inverse F−1:W→VF^{-1}: W \to V are C1C^1 maps.

Intuition

Think of FF as a deformation of space: it stretches, rotates, and shears Rn\mathbb{R}^n near aa. The derivative DF(a)DF(a) is the linear snapshot of this deformation — the best linear approximation, in the sense that

F(a+h)=F(a)+DF(a)h+o(∥h∥),h→0.F(a+h) = F(a) + DF(a)h + o(\|h\|), \qquad h \to 0.

If this linear snapshot DF(a)DF(a) is invertible, it means the deformation does not collapse any direction to zero near aa — no direction gets crushed flat. Geometrically, a small ball around aa gets mapped, to first order, onto a small ellipsoid around F(a)F(a) of nonzero volume and full dimension. The Inverse Function Theorem says this "no collapsing" persists in the actual nonlinear map, not just its linear approximation: nearby points map to nearby distinct images, and every point near F(a)F(a) is hit exactly once from nearby. A useful mental model in n=2n=2: consider a smooth map from the plane to the plane, like (r,θ)↦(rcos⁡θ,rsin⁡θ)(r,\theta) \mapsto (r\cos\theta, r\sin\theta) (polar to Cartesian). Its Jacobian determinant is rr, which vanishes exactly at r=0r=0 — the one place where circles collapse to the origin. Everywhere with r≠0r \neq 0, the map is genuinely a local diffeomorphism: you can locally recover (r,θ)(r,\theta) smoothly from (x,y)(x,y).

Figure
The Jacobian as a local "stretch map"; invertible at aa means no local collapse
Domain near a Image near F(a) a F, DF(a) invertible F(a) Small disk around a Nondegenerate ellipse — full dimension, no collapse
TheoremInverse Function Theorem

Let U⊆RnU \subseteq \mathbb{R}^n be open, let F:U→RnF: U \to \mathbb{R}^n be a C1C^1 map, and let a∈Ua \in U be a point at which DF(a)DF(a) is invertible. Then: 1. There exist open sets V⊆UV \subseteq U with a∈Va \in V, and W⊆RnW \subseteq \mathbb{R}^n with F(a)∈WF(a) \in W, such that F:V→WF: V \to W is a bijection. 2. The inverse map G:=(F∣V)−1:W→VG := (F|_V)^{-1}: W \to V is C1C^1 on WW. 3. For every y∈Wy \in W, writing x=G(y)∈Vx = G(y) \in V,

DG(y)=[DF(x)]−1.DG(y) = \big[DF(x)\big]^{-1}.

That is, FF is a local C1C^1-diffeomorphism at aa, and the derivative of the inverse is the matrix inverse of the derivative.

Proof

We give the classical contraction mapping proof, in four stages: normalization, local injectivity and a contraction estimate, surjectivity onto an open set via the contraction principle, and differentiability of the inverse.

Stage 0: Reduction by an invertible linear change. Let A=DF(a)A = DF(a). Since AA is invertible, define F~(x)=A−1F(x)\widetilde F(x) = A^{-1}F(x). Then F~\widetilde F is C1C^1 (composition of C1C^1 with linear), DF~(x)=A−1DF(x)D\widetilde F(x) = A^{-1}DF(x) by the The Chain Rule in Several Variables, and DF~(a)=A−1A=ID\widetilde F(a) = A^{-1}A = I, the identity matrix. If we prove the theorem for F~\widetilde F — producing a local C1C^1 inverse G~\widetilde G near F~(a)\widetilde F(a) — then G(y):=G~(A−1y)G(y) := \widetilde G(A^{-1}y) is a local C1C^1 inverse for FF near F(a)F(a), since F=AF~F = A\widetilde F. So without loss of generality we assume DF(a)=IDF(a) = I.

Stage 1: Setting up the contraction. Since DFDF is continuous at aa (this is where C1C^1, not mere differentiability, is used) and DF(a)=IDF(a) = I, choose r>0r>0 so small that the closed ball B‾(a,2r)⊂U\overline{B}(a,2r) \subset U and

∥DF(x)−I∥≤12for all x∈B‾(a,2r),\|DF(x) - I\| \le \tfrac{1}{2} \quad \text{for all } x \in \overline{B}(a,2r),

where ∥⋅∥\|\cdot\| denotes the operator norm on matrices. (Such rr exists because x↦∥DF(x)−I∥x \mapsto \|DF(x)-I\| is continuous and equals 00 at x=ax=a.) For y∈Rny \in \mathbb{R}^n define the auxiliary map

φy(x):=x+y−F(x),x∈B‾(a,2r).\varphi_y(x) := x + y - F(x), \qquad x \in \overline{B}(a,2r).

Note that xx is a fixed point of φy\varphi_y if and only if F(x)=yF(x) = y. We show φy\varphi_y is a contraction on B‾(a,2r)\overline{B}(a,2r) (for suitable yy), with Lipschitz constant 12\tfrac12. The derivative of φy\varphi_y is Dφy(x)=I−DF(x)D\varphi_y(x) = I - DF(x), so ∥Dφy(x)∥≤12\|D\varphi_y(x)\| \le \tfrac12 for all x∈B‾(a,2r)x \in \overline{B}(a,2r). By the mean value inequality (applied componentwise, or directly to vector-valued C1C^1 maps: for x1,x2x_1,x_2 in a convex set, ∥φy(x1)−φy(x2)∥≤sup⁡t∥Dφy(x1+t(x2−x1))∥⋅∥x1−x2∥\|\varphi_y(x_1)-\varphi_y(x_2)\| \le \sup_{t}\|D\varphi_y(x_1+t(x_2-x_1))\| \cdot \|x_1-x_2\|, a standard consequence of the fundamental theorem of calculus applied to t↦φy(x1+t(x2−x1))t\mapsto \varphi_y(x_1+t(x_2-x_1))), we get for all x1,x2∈B‾(a,2r)x_1,x_2 \in \overline{B}(a,2r):

∥φy(x1)−φy(x2)∥≤12∥x1−x2∥.(∗)\|\varphi_y(x_1) - \varphi_y(x_2)\| \le \tfrac12 \|x_1 - x_2\|. \tag{$\ast$}

In particular, taking y=F(x2)y = F(x_2) (so φy(x2)=x2\varphi_y(x_2) = x_2... we do not yet know that, so instead directly apply (∗)(\ast) with x2=ax_2 = a, yy arbitrary):

Injectivity of FF on B‾(a,2r)\overline{B}(a,2r). If F(x1)=F(x2)F(x_1) = F(x_2) for x1,x2∈B‾(a,2r)x_1,x_2\in \overline B(a,2r), then φy(x1)−φy(x2)=(x1−x2)−(F(x1)−F(x2))=x1−x2\varphi_y(x_1) - \varphi_y(x_2) = (x_1-x_2) - (F(x_1)-F(x_2)) = x_1 - x_2 for y=F(x1)y=F(x_1). Combined with (∗)(\ast): ∥x1−x2∥≤12∥x1−x2∥\|x_1-x_2\| \le \tfrac12\|x_1-x_2\|, forcing x1=x2x_1=x_2. So FF is injective on B‾(a,2r)\overline{B}(a,2r).

Stage 2: Existence of the inverse on a genuinely open set (surjectivity via contraction). Let W0:=B(F(a),r/2)W_0 := B\big(F(a), r/2\big), an open ball. We claim: for every y∈W0y \in W_0, there is a unique x∈B(a,2r)x \in B(a,2r) with F(x)=yF(x) = y, and moreover x∈B(a,r)x \in B(a,r). Fix y∈W0y \in W_0. We show φy\varphi_y maps the closed ball B‾(a,r)\overline{B}(a,r) into itself and is a contraction there, hence has a unique fixed point by the contraction mapping principle (a standard consequence of completeness of Rn\mathbb{R}^n, i.e. of Compactness in R^n together with the Cauchy criterion). For x∈B‾(a,r)x \in \overline{B}(a,r):

∥φy(x)−a∥≤∥φy(x)−φy(a)∥+∥φy(a)−a∥.\|\varphi_y(x) - a\| \le \|\varphi_y(x) - \varphi_y(a)\| + \|\varphi_y(a) - a\|.

By (∗)(\ast), ∥φy(x)−φy(a)∥≤12∥x−a∥≤r2\|\varphi_y(x)-\varphi_y(a)\| \le \tfrac12\|x-a\| \le \tfrac{r}{2}. Also φy(a)−a=y−F(a)\varphi_y(a) - a = y - F(a), so ∥φy(a)−a∥=∥y−F(a)∥<r/2\|\varphi_y(a)-a\| = \|y - F(a)\| < r/2 since y∈W0y \in W_0. Hence

∥φy(x)−a∥<r2+r2=r,\|\varphi_y(x) - a\| < \tfrac{r}{2} + \tfrac{r}{2} = r,

so φy(x)∈B(a,r)⊂B‾(a,r)\varphi_y(x) \in B(a,r) \subset \overline{B}(a,r). Thus φy:B‾(a,r)→B‾(a,r)\varphi_y : \overline{B}(a,r) \to \overline{B}(a,r) maps this complete metric space into itself, and by (∗)(\ast) it is a contraction with constant 12<1\tfrac12 < 1. By the Banach contraction mapping principle, φy\varphi_y has a unique fixed point x=x(y)∈B‾(a,r)x = x(y) \in \overline{B}(a,r), i.e. a unique x∈B‾(a,r)x\in \overline B(a,r) with F(x)=yF(x) = y; the argument above (bounding ∥φy(x)−a∥<r\|\varphi_y(x)-a\|<r strictly) shows this fixed point actually lies in the open ball B(a,r)B(a,r). Combined with injectivity on the larger ball B‾(a,2r)\overline B(a,2r) from Stage 1, this xx is the unique preimage of yy in all of B‾(a,2r)\overline B(a,2r). Define V:=F−1(W0)∩B(a,r)V := F^{-1}(W_0) \cap B(a,r) and W:=W0W := W_0. By the above, F:V→WF: V \to W is a bijection, with inverse G:=F−1∣WG := F^{-1}|_W. We should check VV is open: since FF is continuous, and we can show below that GG is continuous (in fact Lipschitz), V=G(W)V = G(W) is the continuous image of an open set under a map that is, as we'll see, an open map — alternatively, directly: for any x0∈Vx_0 \in V, F(x0)∈W0F(x_0) \in W_0 open, and by continuity of FF and the estimate below there's a ball around x0x_0 mapping into W0∩W_0 \cap appropriate region; we verify this rigorously via the Lipschitz estimate on GG next, which shows GG is continuous, hence V=G(W)V=G(W) is open as the preimage of the open set... to avoid circularity we argue directly:

VV is open. Let x0∈Vx_0\in V, so F(x0)=y0∈W0F(x_0)=y_0\in W_0. Since W0W_0 open, pick ρ>0\rho>0 with B(y0,ρ)⊂W0B(y_0,\rho)\subset W_0. For any y∈B(y0,ρ)y\in B(y_0,\rho), the fixed-point argument above (applied verbatim with yy in place of the earlier generic point, since y∈W0y \in W_0) produces x=x(y)∈B(a,r)x=x(y)\in B(a,r) with F(x)=yF(x)=y; we must check xx is close to x0x_0. From (∗)(\ast) with x1=x(y),x2=x0x_1 = x(y), x_2 = x_0:

∥x(y)−x0∥=∥φy(x(y))−φy0(x0)∥≤∥φy(x(y))−φy(x0)∥+∥φy(x0)−φy0(x0)∥\|x(y) - x_0\| = \|\varphi_{y}(x(y)) - \varphi_{y_0}(x_0) \| \le \|\varphi_y(x(y)) - \varphi_y(x_0)\| + \|\varphi_y(x_0)-\varphi_{y_0}(x_0)\|

using φy(x0)−φy0(x0)=(y−y0)\varphi_y(x_0) - \varphi_{y_0}(x_0) = (y - y_0) directly from the definition of φy\varphi_y:

≤12∥x(y)−x0∥+∥y−y0∥,\le \tfrac12\|x(y)-x_0\| + \|y-y_0\|,

so ∥x(y)−x0∥≤2∥y−y0∥\|x(y)-x_0\| \le 2\|y-y_0\|. As yy ranges over B(y0,ρ)B(y_0,\rho), x(y)x(y) ranges over a subset of B(x0,2ρ)B(x_0, 2\rho); choosing ρ\rho small enough that B(x0,2ρ)⊂V∩B(a,r)B(x_0,2\rho)\subset V\cap B(a,r) (possible since B(a,r)B(a,r) is open and x0∈B(a,r)x_0 \in B(a,r)), this shows F−1(B(y0,ρ))∩B(a,r)⊂B(x0,2ρ)F^{-1}(B(y_0,\rho)) \cap B(a,r) \subset B(x_0,2\rho), and this whole preimage lies in VV. Since x0x_0 was arbitrary and every point of F(V)=W0F(V)=W_0 has such a neighborhood pulling back into VV, we conclude VV is open.   □\;\square (Stage 2)

Also, the same inequality shows GG is Lipschitz continuous: for y1,y2∈Wy_1, y_2 \in W,

∥G(y1)−G(y2)∥≤2∥y1−y2∥.(∗∗)\|G(y_1) - G(y_2)\| \le 2\|y_1 - y_2\|. \tag{$\ast\ast$}

(Same computation as above, with y0→y1y_0 \to y_1, x0→G(y1)x_0\to G(y_1).)

Stage 3: Differentiability of GG. Fix y∈Wy \in W and let x=G(y)∈Vx = G(y) \in V. Since DF(x)DF(x) is invertible — indeed ∥DF(x)−I∥≤1/2\|DF(x)-I\|\le 1/2 on B‾(a,2r)\overline B(a,2r) forces DF(x)DF(x) invertible, because if DF(x)v=0DF(x)v=0 for unit vector vv then ∥v−DF(x)v∥=∥v∥=1≤∥I−DF(x)∥⋅∥v∥≤1/2\|v - DF(x)v\| = \|v\| =1 \le \|I - DF(x)\|\cdot\|v\| \le 1/2, a contradiction — let T=[DF(x)]−1T = [DF(x)]^{-1}. We claim DG(y)=TDG(y) = T. Let k∈Rnk \in \mathbb{R}^n be small and nonzero with y+k∈Wy+k \in W. Let h:=G(y+k)−G(y)h := G(y+k) - G(y), so x+h=G(y+k)∈Vx+h = G(y+k) \in V, and F(x+h)=y+kF(x+h) = y+k, F(x)=yF(x) = y. By (∗∗)(\ast\ast), ∥h∥≤2∥k∥\|h\| \le 2\|k\|, so h→0h \to 0 as k→0k \to 0; in particular h≠0h \neq 0 whenever k≠0k\neq 0 (by injectivity of FF on VV), and more importantly ∥h∥/∥k∥\|h\|/\|k\| stays bounded. By differentiability of FF at xx:

k=F(x+h)−F(x)=DF(x)h+ε(h)∥h∥,ε(h)→0 as h→0.k = F(x+h) - F(x) = DF(x)h + \varepsilon(h)\|h\|, \qquad \varepsilon(h) \to 0 \text{ as } h \to 0.

Apply T=[DF(x)]−1T = [DF(x)]^{-1} to both sides:

Tk=h+Tε(h)∥h∥⟹h−Tk=−Tε(h)∥h∥.Tk = h + T\varepsilon(h)\|h\| \quad\Longrightarrow\quad h - Tk = -T\varepsilon(h)\|h\|.

Thus

∥G(y+k)−G(y)−Tk∥∥k∥=∥h−Tk∥∥k∥=∥Tε(h)∥⋅∥h∥∥k∥≤∥T∥⋅∥ε(h)∥⋅2.\frac{\|G(y+k) - G(y) - Tk\|}{\|k\|} = \frac{\|h - Tk\|}{\|k\|} = \|T\varepsilon(h)\| \cdot \frac{\|h\|}{\|k\|} \le \|T\| \cdot \|\varepsilon(h)\| \cdot 2.

As k→0k \to 0, h→0h \to 0 (shown above), so ε(h)→0\varepsilon(h) \to 0, and the right side →0\to 0. This is exactly the statement that GG is differentiable at yy with DG(y)=T=[DF(x)]−1=[DF(G(y))]−1DG(y) = T = [DF(x)]^{-1} = [DF(G(y))]^{-1}.

Stage 4: Continuity of DGDG. We have, for all y∈Wy \in W,

DG(y)=[DF(G(y))]−1.DG(y) = [DF(G(y))]^{-1}.

The map y↦G(y)y \mapsto G(y) is continuous (Lipschitz, Stage 2); the map x↦DF(x)x \mapsto DF(x) is continuous (FF is C1C^1); and matrix inversion M↦M−1M \mapsto M^{-1} is continuous on the open set of invertible matrices (each entry of M−1M^{-1} is, by Cramer's rule, a ratio of polynomials in the entries of MM divided by det⁡M≠0\det M \neq 0, hence continuous). Composing three continuous maps, y↦DG(y)y \mapsto DG(y) is continuous on WW. Hence G∈C1(W)G \in C^1(W). This proves all three parts of the theorem, with V,WV,W as constructed. ■\blacksquare

∎

Worked examples

Example 1: A genuinely asymmetric planar map. Let F:R2→R2F:\mathbb{R}^2 \to \mathbb{R}^2 be

F(x,y)=(x2−y3+2x,    xy+ey−1).F(x,y) = (x^2 - y^3 + 2x,\; \; xy + e^y - 1).

Take the base point a=(1,0)a = (1,0). Compute:

F(1,0)=(1−0+2,  0+1−1)=(3,0).F(1,0) = (1 - 0 + 2,\; 0 + 1 - 1) = (3, 0).

The Jacobian matrix is

DF(x,y)=(2x+2−3y2yx+ey).DF(x,y) = \begin{pmatrix} 2x+2 & -3y^2 \\ y & x + e^y \end{pmatrix}.

At a=(1,0)a = (1,0):

DF(1,0)=(4002).DF(1,0) = \begin{pmatrix} 4 & 0 \\ 0 & 2 \end{pmatrix}.

This is diagonal, invertible, with det⁡=8≠0\det = 8 \neq 0. Since FF's components are polynomials and exponentials, all partial derivatives are continuous everywhere, so FF is C1C^1 on all of R2\mathbb{R}^2. The Inverse Function Theorem applies at a=(1,0)a=(1,0): there exist open V∋(1,0)V \ni (1,0) and W∋(3,0)W \ni (3,0) such that F:V→WF: V \to W is a C1C^1-diffeomorphism, and

DG(3,0)=[DF(1,0)]−1=(1/4001/2).DG(3,0) = [DF(1,0)]^{-1} = \begin{pmatrix} 1/4 & 0 \\ 0 & 1/2 \end{pmatrix}.

Note we have not solved for GG explicitly — the theorem guarantees its local existence and smoothness, and hands us its derivative at the one point we need, purely from DF(a)DF(a), without ever inverting FF in closed form. This is the entire power of the theorem: local invertibility and its derivative are certified algebraically.

Example 2: A 3×33\times 3 map, checking the hypothesis fails at one nearby point but holds at another. Let F:R3→R3F: \mathbb{R}^3 \to \mathbb{R}^3,

F(x,y,z)=(x+sin⁡z,  y2+z,  xz+y).F(x,y,z) = (x + \sin z,\; y^2 + z,\; xz + y).

Take b=(2,1,0)b = (2,1,0). Then

F(2,1,0)=(2+0,  1+0,  0+1)=(2,1,1).F(2,1,0) = (2+0,\;1+0,\;0+1) = (2,1,1).

Jacobian:

DF(x,y,z)=(10cos⁡z02y1z1x).DF(x,y,z) = \begin{pmatrix} 1 & 0 & \cos z \\ 0 & 2y & 1 \\ z & 1 & x \end{pmatrix}.

At b=(2,1,0)b=(2,1,0):

DF(2,1,0)=(101021012).DF(2,1,0) = \begin{pmatrix} 1 & 0 & 1 \\ 0 & 2 & 1 \\ 0 & 1 & 2 \end{pmatrix}.

Expand det⁡\det along the first column: det⁡=1⋅(2⋅2−1⋅1)−0+0=1⋅3=3≠0\det = 1\cdot(2\cdot2 - 1\cdot1) - 0 + 0 = 1\cdot 3 = 3 \neq 0. So FF is a local C1C^1-diffeomorphism at bb, with inverse derivative at F(b)=(2,1,1)F(b)=(2,1,1) equal to [DF(b)]−1[DF(b)]^{-1}. Computing the inverse (via cofactors, since det⁡=3\det=3):

[DF(2,1,0)]−1=13(31−202−10−12),[DF(2,1,0)]^{-1} = \frac{1}{3}\begin{pmatrix} 3 & 1 & -2 \\ 0 & 2 & -1 \\ 0 & -1 & 2 \end{pmatrix},

which one checks by direct multiplication: DF(b)⋅[DF(b)]−1=IDF(b)\cdot [DF(b)]^{-1} = I. (Multiplying row 1 of DF(b)DF(b), (1,0,1)(1,0,1), against column 1 of the claimed inverse, (3,0,0)/3=(1,0,0)(3,0,0)/3=(1,0,0), gives 11; against column 2, (1,2,−1)/3(1,2,-1)/3, gives (1−1)/3=0(1-1)/3=0; against column 3, (−2,−1,2)/3(-2,-1,2)/3, gives (−2+2)/3=0(-2+2)/3=0 — consistent with the identity's first row, and the remaining rows check similarly.) Contrast this with the point c=(0,1,0)c=(0,1,0): there

DF(0,1,0)=(101021010),det⁡=1(2⋅0−1⋅1)−0+1(0−0)=−1≠0,DF(0,1,0) = \begin{pmatrix} 1 & 0 & 1 \\ 0 & 2 & 1 \\ 0 & 1 & 0 \end{pmatrix}, \qquad \det = 1(2\cdot0 - 1\cdot1) - 0 + 1(0-0) = -1 \neq 0,

still nonsingular — but if instead we test c′=(0,0,0)c' = (0,0,0): DF(0,0,0)=(101001010)DF(0,0,0) = \begin{pmatrix}1&0&1\\0&0&1\\0&1&0\end{pmatrix}, det⁡=1(0−1)−0+1(0−0)=−1≠0\det = 1(0-1) - 0 + 1(0-0) = -1 \neq 0 still. The point of this example is that checking the hypothesis is a genuine per-point computation — the theorem gives no shortcut, and gives no information at any point where det⁡DF=0\det DF = 0 (which for this FF occurs, e.g., when solving det⁡DF(x,y,z)=2xy−z2cos⁡z−2y⋅0...\det DF(x,y,z)=2xy - z^2\cos z - 2y\cdot 0...; the qualitative lesson is that the singular locus {det⁡DF=0}\{\det DF = 0\} is a genuine lower-dimensional constraint set, and the theorem is silent there).

CounterexampleWhy C1C^1 (continuity of the derivative), not mere differentiability, is essential

It is tempting to think "differentiable with invertible derivative at aa" alone should suffice — but this is false. The theorem's proof uses continuity of DFDF near aa (Stage 1: to find a whole ball on which ∥DF(x)−I∥≤1/2\|DF(x)-I\|\le 1/2), not just at aa. Consider the classical one-dimensional example, adapted here for concreteness (the same phenomenon occurs in every dimension n≥1n\ge1 by taking a product with the identity on the remaining coordinates). Define f:R→Rf:\mathbb{R}\to\mathbb{R} by

f(x)={x+2x2sin⁡(1/x),x≠00,x=0.f(x) = \begin{cases} x + 2x^2\sin(1/x), & x \neq 0 \\ 0, & x = 0. \end{cases}

ff is differentiable everywhere, with

f′(x)=1+4xsin⁡(1/x)−2cos⁡(1/x),x≠0,f′(0)=lim⁡h→0f(h)−0h=lim⁡h→0(1+2hsin⁡(1/h))=1.f'(x) = 1 + 4x\sin(1/x) - 2\cos(1/x), \quad x\neq 0, \qquad f'(0) = \lim_{h\to0}\frac{f(h)-0}{h} = \lim_{h\to0}\big(1+2h\sin(1/h)\big) = 1.

So f′(0)=1≠0f'(0) = 1 \neq 0: the derivative is invertible ("nonsingular") at a=0a=0, exactly the hypothesis of the theorem apart from continuity.

But f′f' is discontinuous at 00: as x→0x\to 0 along xk=1/(2πk)x_k = 1/(2\pi k), cos⁡(1/xk)=cos⁡(2πk)=1\cos(1/x_k) = \cos(2\pi k) = 1, so f′(xk)=1+0−2=−1f'(x_k) = 1 + 0 - 2 = -1, which does not converge to f′(0)=1f'(0)=1. So ff is differentiable but not C1C^1 at 00.

Consequence: ff is not injective on any neighborhood of 00. Since f′(xk)=−1<0f'(x_k) = -1 < 0 for xk→0x_k \to 0, while f′(0)=1>0f'(0)=1>0, and f′f' is continuous away from 00, the intermediate value theorem (applied to the continuous function f′f' on (0,xk)(0,x_k), say) shows f′f' takes both signs infinitely often arbitrarily close to 00. In particular ff has infinitely many critical points xk′→0x_k' \to 0 with f′(xk′)=0f'(x_k')=0, accumulating at 00; between consecutive sign changes of f′f', ff is not monotonic. Hence no interval around 00 maps injectively under ff — there is no local inverse at all, let alone a differentiable one. This shows: dropping continuity of the derivative (weakening C1C^1 to "differentiable, with DF(a)DF(a) invertible") breaks the conclusion entirely — not just smoothness of the inverse, but the mere existence of a local inverse function. The hypothesis "F∈C1F\in C^1" is not a technical convenience; it is load-bearing.

Scope, edge cases & common traps

SituationWhat holdsWhat can fail
F∈C1F \in C^1, DF(a)DF(a) invertibleLocal C1C^1-diffeomorphism (theorem)Global injectivity — not guaranteed
F∈CkF \in C^k, k≥1k\ge1, DF(a)DF(a) invertibleLocal inverse is also CkC^k (bootstrap via the same formula DG=(DF∘G)−1DG=(DF\circ G)^{-1})—
FF differentiable, not C1C^1, DF(a)DF(a) invertibleNothing guaranteedCounterexample above: no local inverse at all
DF(a)DF(a) singular (det⁡DF(a)=0\det DF(a) = 0)Theorem gives no informationFF may or may not be locally injective; theorem is simply silent
FF injective globally on UU—Does not imply DFDF is invertible everywhere (e.g. F(x)=x3F(x)=x^3 on R\mathbb{R} is a global bijection but F′(0)=0F'(0)=0)

The mathematician behind it

Origins · 1 minAugustin-Louis Cauchy and how this was discovered
AC
Augustin-Louis Cauchy
1789–1857 · Early rigorous analysis in France

Cauchy did not state the Inverse Function Theorem in its modern form, but his program of rigorizing calculus — precise limits, the mean value theorem, careful treatment of continuity — created the toolkit (uniform estimates, the mean value inequality) on which the contraction-mapping proof depends. The theorem as stated here, with the modern C1C^1 hypothesis and the fixed-point method of proof, crystallized over the late nineteenth and early twentieth centuries through the work of analysts formalizing implicit and inverse function results for systems of equations, and was cast in its now-standard functional-analytic form (via Banach's fixed point theorem, itself from the 1920s) well after Cauchy — but the analytic rigor that makes such a theorem possible to even state precisely traces directly to his era.

History
- 1821 — Cauchy's Cours d'Analyse establishes rigorous limit-based foundations for continuity and differentiability.
- 1837 — Cauchy proves an early implicit-function result for analytic equations, a precursor to the general theorem.
- 1877 — Ulisse Dini gives a general implicit function theorem for real variables, closely related to the inverse function theorem.
- 1922 — Stefan Banach proves the contraction mapping (fixed point) theorem, supplying the modern proof technique used here.
- 1930s — The Inverse and Implicit Function Theorems reach their standard textbook form, stated for C1C^1 maps between Euclidean spaces and later Banach spaces.
Does the Inverse Function Theorem give a formula for the inverse function itself?

No. The theorem is a pure existence-and-regularity statement: it guarantees a local inverse GG exists and is C1C^1, and it gives a formula for the derivative of GG at a point, namely DG(y)=[DF(x)]−1DG(y) = [DF(x)]^{-1} where y=F(x)y=F(x). It does not produce a closed-form expression for G(y)G(y) itself, which in general may not have one in elementary terms. In practice one often computes GG approximately via Newton's method or a Taylor expansion, using the derivative formula the theorem supplies as the leading-order term.

Is the Inverse Function Theorem true if the domain and codomain have different dimensions?

No — the very statement requires F:U⊆Rn→RnF: U \subseteq \mathbb{R}^n \to \mathbb{R}^n, same dimension nn in source and target, because DF(a)DF(a) must be a square matrix to even ask whether it is invertible. If dimensions differ, the correct nonlinear generalization is different: an injective linear derivative gives a local immersion result, and a surjective one gives a local submersion result, but neither produces a two-sided inverse. The closely related The Implicit Function Theorem handles the genuinely different-dimensional situation of solving part of a system for some variables in terms of others.

Why does the proof need a closed ball and completeness of Rn\mathbb{R}^n?

The contraction mapping principle requires a complete metric space: a Cauchy sequence generated by iterating the contraction must actually converge to a point within the space. A closed ball in Rn\mathbb{R}^n is complete (closed subsets of a complete space are complete), which is ultimately a consequence of the completeness of R\mathbb{R} and, more broadly, of the compactness properties surveyed in Compactness in R^n. If one instead worked in an open ball, the fixed point produced by the iteration might only exist as a limit sitting on the boundary, outside the space, so the closed ball is the technically correct setting even though the conclusion is ultimately stated on open sets.

How does this theorem relate to the Implicit Function Theorem?

They are logically equivalent — each can be derived from the other. The Inverse Function Theorem is the special case of the The Implicit Function Theorem where you want to solve the full system F(x)=yF(x) = y for xx in terms of yy, with no variables held fixed. Conversely, the Implicit Function Theorem can be derived from the Inverse Function Theorem by augmenting a map F(x,y)F(x,y) with the identity in the "xx" slot, forming Φ(x,y)=(x,F(x,y))\Phi(x,y) = (x, F(x,y)), and inverting Φ\Phi. Most textbooks prove one directly and derive the other, and the choice of which to prove first is largely a matter of exposition.

If DF(a)DF(a) is invertible, is FF automatically injective on the entire connected component of aa in UU?

No, and this is one of the most common misunderstandings. The theorem is strictly local: it produces some neighborhood VV of aa on which FF is injective, but says nothing about behavior far from aa, even if det⁡DF(x)≠0\det DF(x) \neq 0 at every point of a large connected set. The map F(x,y)=(excos⁡y,exsin⁡y)F(x,y) = (e^x\cos y, e^x\sin y) is the standard illustration: its Jacobian determinant e2xe^{2x} never vanishes, yet FF is periodic in yy with period 2π2\pi and hence infinitely far from globally injective.

Does the theorem require FF to be defined on all of Rn\mathbb{R}^n?

No. It is a genuinely local statement: FF only needs to be defined and C1C^1 on some open set UU containing the point aa of interest. The neighborhoods VV and WW produced by the proof can be taken as small as necessary, so the theorem applies equally well to maps defined only on a small open ball around aa, or to restrictions of larger maps, as long as the C1C^1 hypothesis and the nonsingular-Jacobian hypothesis hold at aa itself.

Before this

Where this leads