Calculus on Euclidean space

In mathematics, calculus on Euclidean space is a generalization of calculus of functions in one or several variables to calculus of functions on Euclidean space $\mathbb {R} ^{n}$ as well as a finite-dimensional real vector space. This calculus is also known as advanced calculus, especially in the United States. It is similar to multivariable calculus but is somewhat more sophisticated in that it uses linear algebra (or some functional analysis) more extensively and covers some concepts from differential geometry such as differential forms and Stokes' formula in terms of differential forms. This extensive use of linear algebra also allows a natural generalization of multivariable calculus to calculus on Banach spaces or topological vector spaces.

Calculus on Euclidean space is also a local model of calculus on manifolds, a theory of functions on manifolds.

Basic notions

See also: Function of a real variable and Multivariable calculus

Functions in one real variable

This section is a brief review of function theory in one-variable calculus.

A real-valued function $f:\mathbb {R} \to \mathbb {R}$ is continuous at $a$ if it is approximately constant near $a$; i.e.,

$\lim _{h\to 0}(f(a+h)-f(a))=0.$

In contrast, the function $f$ is differentiable at $a$ if it is approximately linear near $a$; i.e., there is some real number $\lambda$ such that

$\lim _{h\to 0}{\frac {f(a+h)-f(a)-\lambda h}{h}}=0.$

(For simplicity, suppose $f(a)=0$. Then the above means that $f(a+h)=\lambda h+g(a,h)$ where $g(a,h)$ goes to 0 faster than h going to 0 and, in that sense, $f(a+h)$ behaves like $\lambda h$.)

The number $\lambda$ depends on $a$ and thus is denoted as $f'(a)$. If $f$ is differentiable on an open interval $U$ and if $f'$ is a continuous function on $U$, then $f$ is called a C¹ function. More generally, $f$ is called a C^k function if its derivative $f'$ is C^k-1 function. Taylor's theorem states that a C^k function is precisely a function that can be approximated by a polynomial of degree k.

If $f:\mathbb {R} \to \mathbb {R}$ is a C¹ function and $f'(a)\neq 0$ for some $a$, then either $f'(a)>0$ or $f'(a)<0$; i.e., either $f$ is strictly increasing or strictly decreasing in some open interval containing a. In particular, $f:f^{-1}(U)\to U$ is bijective for some open interval $U$ containing $f(a)$. The inverse function theorem then says that the inverse function $f^{-1}$ is differentiable on U with the derivatives: for $y\in U$

$(f^{-1})'(y)={1 \over f'(f^{-1}(y))}.$

Derivative of a map and chain rule

For functions $f$ defined in the plane or more generally on an Euclidean space $\mathbb {R} ^{n}$, it is necessary to consider functions that are vector-valued or matrix-valued. It is also conceptually helpful to do this in an invariant manner (i.e., a coordinate-free way). Derivatives of such maps at a point are then vectors or linear maps, not real numbers.

Let $f:X\to Y$ be a map from an open subset $X$ of $\mathbb {R} ^{n}$ to an open subset $Y$ of $\mathbb {R} ^{m}$. Then the map $f$ is said to be differentiable at a point $x$ in $X$ if there exists a (necessarily unique) linear transformation $f'(x):\mathbb {R} ^{n}\to \mathbb {R} ^{m}$, called the derivative of $f$ at $x$, such that

$\lim _{h\to 0}{\frac {1}{|h|}}|f(x+h)-f(x)-f'(x)h|=0$

where $f'(x)h$ is the application of the linear transformation $f'(x)$ to $h$. If $f$ is differentiable at $x$, then it is continuous at $x$ since

$|f(x+h)-f(x)|\leq (|h|^{-1}|f(x+h)-f(x)-f'(x)h|)|h|+|f'(x)h|\to 0$ as $h\to 0$.

As in the one-variable case, there is

Chain rule —  Let $f$ be as above and $g:Y\to Z$ a map for some open subset $Z$ of $\mathbb {R} ^{l}$. If $f$ is differentiable at $x$ and $g$ differentiable at $y=f(x)$, then the composition $g\circ f$ is differentiable at $x$ with the derivative

$(g\circ f)'(x)=g'(y)\circ f'(x).$

This is proved exactly as for functions in one variable. Indeed, with the notation ${\widetilde {h}}=f(x+h)-f(x)$, we have:

${\begin{aligned}&{\frac {1}{|h|}}|g(f(x+h))-g(y)-g'(y)f'(x)h|\\&\leq {\frac {1}{|h|}}|g(y+{\widetilde {h}})-g(y)-g'(y){\widetilde {h}}|+{\frac {1}{|h|}}|g'(y)(f(x+h)-f(x)-f'(x)h)|.\end{aligned}}$

Here, since $f$ is differentiable at $x$, the second term on the right goes to zero as $h\to 0$. As for the first term, it can be written as:

${\begin{cases}{\frac {|{\widetilde {h}}|}{|h|}}|g(y+{\widetilde {h}})-g(y)-g'(y){\widetilde {h}}|/|{\widetilde {h}}|,&{\widetilde {h}}\neq 0,\\0,&{\widetilde {h}}=0.\end{cases}}$

Now, by the argument showing the continuity of $f$ at $x$, we see ${\frac {|{\widetilde {h}}|}{|h|}}$ is bounded. Also, ${\widetilde {h}}\to 0$ as $h\to 0$ since $f$ is continuous at $x$. Hence, the first term also goes to zero as $h\to 0$ by the differentiability of $g$ at $y$. $\square$

The map $f$ as above is called continuously differentiable or $C^{1}$ if it is differentiable on the domain and also the derivatives vary continuously; i.e., $x\mapsto f'(x)$ is continuous.

Corollary — If $f,g$ are continuously differentiable, then $g\circ f$ is continuously differentiable.

As a linear transformation, $f'(x)$ is represented by an $m\times n$-matrix, called the Jacobian matrix $Jf(x)$ of $f$ at $x$ and we write it as:

$(Jf)(x)={\begin{bmatrix}{\frac {\partial f_{1}}{\partial x_{1}}}(x)&\cdots &{\frac {\partial f_{1}}{\partial x_{n}}}(x)\\\vdots &\ddots &\vdots \\{\frac {\partial f_{m}}{\partial x_{1}}}(x)&\cdots &{\frac {\partial f_{m}}{\partial x_{n}}}(x)\end{bmatrix}}.$

Taking $h$ to be $he_{j}$, $h$ a real number and $e_{j}=(0,\cdots ,1,\cdots ,0)$ the j-th standard basis element, we see that the differentiability of $f$ at $x$ implies:

$\lim _{h\to 0}{\frac {f_{i}(x+he_{j})-f_{i}(x)}{h}}={\frac {\partial f_{i}}{\partial x_{j}}}(x)$

where $f_{i}$ denotes the i-th component of $f$. That is, each component of $f$ is differentiable at $x$ in each variable with the derivative ${\frac {\partial f_{i}}{\partial x_{j}}}(x)$. In terms of Jacobian matrices, the chain rule says $J(g\circ f)(x)=Jg(y)Jf(x)$; i.e., as $(g\circ f)_{i}=g_{i}\circ f$,

${\frac {\partial (g_{i}\circ f)}{\partial x_{j}}}(x)={\frac {\partial g_{i}}{\partial y_{1}}}(y){\frac {\partial f_{1}}{\partial x_{j}}}(x)+\cdots +{\frac {\partial g_{i}}{\partial y_{m}}}(y){\frac {\partial f_{m}}{\partial x_{j}}}(x),$

which is the form of the chain rule that is often stated.

A partial converse to the above holds. Namely, if the partial derivatives ${\partial f_{i}}/{\partial x_{j}}$ are all defined and continuous, then $f$ is continuously differentiable. This is a consequence of the mean value inequality:

Mean value inequality —  Given the map $f:X\to Y$ as above and points $x,y$ in $X$ such that the line segment between $x,y$ lies in $X$, if $t\mapsto f(x+ty)$ is continuous on $[0,1]$ and is differentiable on the interior, then, for any vector $v\in \mathbb {R} ^{m}$,

$|\Delta _{y}f(x)-v|\leq \sup _{0<t<1}\left|{\frac {d}{dt}}f(x+ty)-v\right|$

where $\Delta _{y}f(x)=f(x+y)-f(x).$

(This version of mean value inequality follows from mean value inequality in Mean value theorem § Mean value theorem for vector-valued functions applied to the function $[0,1]\to \mathbb {R} ^{m},\,t\mapsto f(x+ty)-tv$, where the proof on mean value inequality is given.)

Indeed, let $g(x)=(Jf)(x)$. We note that, if $y=y_{i}e_{i}$, then

${\frac {d}{dt}}f(x+ty)={\frac {\partial f}{\partial x_{i}}}(x+ty)y=g(x+ty)(y_{i}e_{i}).$

For simplicity, assume $n=2$ (the argument for the general case is similar). Then, by mean value inequality, with the operator norm $\|\cdot \|$,

${\begin{aligned}&|\Delta _{y}f(x)-g(x)y|\\&\leq |\Delta _{y_{1}e_{1}}f(x_{1},x_{2}+y_{2})-g(x)(y_{1}e_{1})|+|\Delta _{y_{2}e_{2}}f(x_{1},x_{2})-g(x)(y_{2}e_{2})|\\&\leq |y_{1}|\sup _{0<t<1}\|g(x_{1}+ty_{1},x_{2}+y_{2})-g(x)\|+|y_{2}|\sup _{0<t<1}\|g(x_{1},x_{2}+ty_{2})-g(x)\|,\end{aligned}}$

which implies $|\Delta _{y}f(x)-g(x)y|/|y|\to 0$ as required. $\square$

Example: Let $U$ be the set of all invertible real square matrices of size n. Note $U$ can be identified as an open subset of $\mathbb {R} ^{n^{2}}$ with coordinates $x_{ij},0\leq i,j\neq n$. Consider the function $f(g)=g^{-1}$ = the inverse matrix of $g$ defined on $U$. To guess its derivatives, assume $f$ is differentiable and consider the curve $c(t)=ge^{tg^{-1}h}$ where $e^{A}$ means the matrix exponential of $A$. By the chain rule applied to $f(c(t))=e^{-tg^{-1}h}g^{-1}$, we have:

$f'(c(t))\circ c'(t)=-g^{-1}he^{-tg^{-1}h}g^{-1}$.

Taking $t=0$, we get:

$f'(g)h=-g^{-1}hg^{-1}$.

Now, we then have:

$\|(g+h)^{-1}-g^{-1}+g^{-1}hg^{-1}\|\leq \|(g+h)^{-1}\|\|h\|\|g^{-1}hg^{-1}\|.$

Since the operator norm is equivalent to the Euclidean norm on $\mathbb {R} ^{n^{2}}$ (any norms are equivalent to each other), this implies $f$ is differentiable. Finally, from the formula for $f'$, we see the partial derivatives of $f$ are smooth (infinitely differentiable); whence, $f$ is smooth too.

Higher derivatives and Taylor formula

If $f:X\to \mathbb {R} ^{m}$ is differentiable where $X\subset \mathbb {R} ^{n}$ is an open subset, then the derivatives determine the map $f':X\to \operatorname {Hom} (\mathbb {R} ^{n},\mathbb {R} ^{m})$, where $\operatorname {Hom}$ stands for homomorphisms between vector spaces; i.e., linear maps. If $f'$ is differentiable, then $f'':X\to \operatorname {Hom} (\mathbb {R} ^{n},\operatorname {Hom} (\mathbb {R} ^{n},\mathbb {R} ^{m}))$. Here, the codomain of $f''$ can be identified with the space of bilinear maps by:

$\operatorname {Hom} (\mathbb {R} ^{n},\operatorname {Hom} (\mathbb {R} ^{n},\mathbb {R} ^{m})){\overset {\varphi }{\underset {\sim }{\to }}}\{(\mathbb {R} ^{n})^{2}\to \mathbb {R} ^{m}{\text{ bilinear}}\}$

where $\varphi (g)(x,y)=g(x)y$ and $\varphi$ is bijective with the inverse $\psi$ given by $(\psi (g)x)y=g(x,y)$. In general, $f^{(k)}=(f^{(k-1)})'$ is a map from $X$ to the space of $k$-multilinear maps $(\mathbb {R} ^{n})^{k}\to \mathbb {R} ^{m}$.

Just as $f'(x)$ is represented by a matrix (Jacobian matrix), when $m=1$ (a bilinear map is a bilinear form), the bilinear form $f''(x)$ is represented by a matrix called the Hessian matrix of $f$ at $x$; namely, the square matrix $H$ of size $n$ such that $f''(x)(y,z)=(Hy,z)$, where the paring refers to an inner product of $\mathbb {R} ^{n}$, and $H$ is none other than the Jacobian matrix of $f':X\to (\mathbb {R} ^{n})^{*}\simeq \mathbb {R} ^{n}$. The $(i,j)$-th entry of $H$ is thus given explicitly as $H_{ij}={\frac {\partial ^{2}f}{\partial x_{i}\partial x_{j}}}(x)$.

Moreover, if $f''$ exists and is continuous, then the matrix $H$ is symmetric, the fact known as the symmetry of second derivatives. This is seen using the mean value inequality. For vectors $u,v$ in $\mathbb {R} ^{n}$, using mean value inequality twice, we have:

$|\Delta _{v}\Delta _{u}f(x)-f''(x)(u,v)|\leq \sup _{0<t_{1},t_{2}<1}|f''(x+t_{1}u+t_{2}v)(u,v)-f''(x)(u,v)|,$

which says

$f''(x)(u,v)=\lim _{s,t\to 0}(\Delta _{tv}\Delta _{su}f(x)-f(x))/(st).$

Since the right-hand side is symmetric in $u,v$, so is the left-hand side: $f''(x)(u,v)=f''(x)(v,u)$. By induction, if $f$ is $C^{k}$, then the k-multilinear map $f^{(k)}(x)$ is symmetric; i.e., the order of taking partial derivatives does not matter.

As in the case of one variable, the Taylor series expansion can then be proved by integration by parts:

$f(z+(h,k))=\sum _{a+b<n}\partial _{x}^{a}\partial _{y}^{b}f(z){h^{a}k^{b} \over a!b!}+n\int _{0}^{1}(1-t)^{n-1}\sum _{a+b=n}\partial _{x}^{a}\partial _{y}^{b}f(z+t(h,k)){h^{a}k^{b} \over a!b!}\,dt.$

Taylor's formula has an effect of dividing a function by variables, which can be illustrated by the next typical theoretical use of the formula.

Example: Let $T:{\mathcal {S}}\to {\mathcal {S}}$ be a linear map between the vector space ${\mathcal {S}}$ of smooth functions on $\mathbb {R} ^{n}$ with rapidly decreasing derivatives; i.e., $\sup |x^{\beta }\partial ^{\alpha }\varphi |<\infty$ for any multi-index $\alpha ,\beta$. (The space ${\mathcal {S}}$ is called a Schwartz space.) For each $\varphi$ in ${\mathcal {S}}$, Taylor's formula implies we can write:

$\varphi -\psi \varphi (y)=\sum _{j=1}^{n}(x_{j}-y_{j})\varphi _{j}$

with $\varphi _{j}\in {\mathcal {S}}$, where $\psi$ is a smooth function with compact support and $\psi (y)=1$. Now, assume $T$ commutes with coordinates; i.e., $T(x_{j}\varphi )=x_{j}T\varphi$. Then

$T\varphi -\varphi (y)T\psi =\sum _{j=1}^{n}(x_{j}-y_{j})T\varphi _{j}$.

Evaluating the above at $y$, we get $T\varphi (y)=\varphi (y)T\psi (y).$ In other words, $T$ is a multiplication by some function $m$; i.e., $T\varphi =m\varphi$. Now, assume further that $T$ commutes with partial differentiations. We then easily see that $m$ is a constant; $T$ is a multiplication by a constant.

(Aside: the above discussion almost proves the Fourier inversion formula. Indeed, let $F,R:{\mathcal {S}}\to {\mathcal {S}}$ be the Fourier transform and the reflection; i.e., $(R\varphi )(x)=\varphi (-x)$. Then, dealing directly with the integral that is involved, one can see $T=RF^{2}$ commutes with coordinates and partial differentiations; hence, $T$ is a multiplication by a constant. This is almost a proof since one still has to compute this constant.)

A partial converse to the Taylor formula also holds; see Borel's lemma and Whitney extension theorem.

Inverse function theorem and submersion theorem

Inverse function theorem — Let $f:X\to Y$ be a map between open subsets $X,Y$ in $\mathbb {R} ^{n},\mathbb {R} ^{m}$. If $f$ is continuously differentiable (or more generally $C^{k}$) and $f'(x)$ is bijective, there exists neighborhoods $U,V$ of $x,f(x)$ and the inverse $f^{-1}:V\to U$ that is continuously differentiable (or respectively $C^{k}$).

A $C^{k}$-map with the $C^{k}$-inverse is called a $C^{k}$-diffeomorphism. Thus, the theorem says that, for a map $f$ satisfying the hypothesis at a point $x$, $f$ is a diffeomorphism near $x,f(x).$ For a proof, see Inverse function theorem § A proof using successive approximation.

The implicit function theorem says: given a map $f:\mathbb {R} ^{n}\times \mathbb {R} ^{m}\to \mathbb {R} ^{m}$, if $f(a,b)=0$, $f$ is $C^{k}$ in a neighborhood of $(a,b)$ and the derivative of $y\mapsto f(a,y)$ at $b$ is invertible, then there exists a differentiable map $g:U\to V$ for some neighborhoods $U,V$ of $a,b$ such that $f(x,g(x))=0$. The theorem follows from the inverse function theorem; see Inverse function theorem § Implicit function theorem.

Another consequence is the submersion theorem.

Integrable functions on Euclidean spaces

A partition of an interval $[a,b]$ is a finite sequence $a=t_{0}\leq t_{1}\leq \cdots \leq t_{k}=b$. A partition $P$ of a rectangle $D$ (product of intervals) in $\mathbb {R} ^{n}$ then consists of partitions of the sides of $D$; i.e., if $D=\prod _{1}^{n}[a_{i},b_{i}]$, then $P$ consists of $P_{1},\dots ,P_{n}$ such that $P_{i}$ is a partition of $[a_{i},b_{i}]$.

Given a function $f$ on $D$, we then define the upper Riemann sum of it as:

$U(f,P)=\sum _{Q\in P}(\sup _{Q}f)\operatorname {vol} (Q)$

where

• $Q$ is a partition element of $P$; i.e., $Q=\prod _{i=1}^{n}[t_{i,j_{i}},t_{i,j_{i}+1}]$ when $P_{i}:a_{i}=t_{i,0}\leq \dots \cdots \leq t_{i,k_{i}}=b_{i}$ is a partition of $[a_{i},b_{i}]$.
• The volume $\operatorname {vol} (Q)$ of $Q$ is the usual Euclidean volume; i.e., $\operatorname {vol} (Q)=\prod _{1}^{n}(t_{i,j_{i}+1}-t_{i,j_{i}})$.

The lower Riemann sum $L(f,P)$ of $f$ is then defined by replacing $\sup$ by $\inf$. Finally, the function $f$ is called integrable if it is bounded and $\sup\{L(f,P)\mid P\}=\inf\{U(f,P)\mid P\}$. In that case, the common value is denoted as $\int _{D}f\,dx$.

A subset of $\mathbb {R} ^{n}$ is said to have measure zero if for each $\epsilon >0$, there are some possibly infinitely many rectangles $D_{1},D_{2},\dots ,$ whose union contains the set and $\sum _{i}\operatorname {vol} (D_{i})<\epsilon .$

A key theorem is

Theorem —  A bounded function $f$ on a closed rectangle is integrable if and only if the set $\{x|f{\text{ is not continuous at }}x\}$ has measure zero.

The next theorem allows us to compute the integral of a function as the iteration of the integrals of the function in one-variables:

Fubini's theorem — If $f$ is a continuous function on a closed rectangle $D=\prod [a_{i},b_{i}]$ (in fact, this assumption is too strong), then

$\int _{D}f\,dx=\int _{a_{n}}^{b_{n}}\cdots \left(\int _{a_{1}}^{b_{1}}f(x_{1},\dots ,x_{n})dx_{1}\right)dx_{2}\cdots dx_{n}.$

In particular, the order of integrations can be changed.

Finally, if $M\subset \mathbb {R} ^{n}$ is a bounded open subset and $f$ a function on $M$, then we define $\int _{M}f\,dx:=\int _{D}\chi _{M}f\,dx$ where $D$ is a closed rectangle containing $M$ and $\chi _{M}$ is the characteristic function on $M$; i.e., $\chi _{M}(x)=1$ if $x\in M$ and $=0$ if $x\not \in M,$ provided $\chi _{M}f$ is integrable.

Surface integral

If a bounded surface $M$ in $\mathbb {R} ^{3}$ is parametrized by ${\textbf {r}}={\textbf {r}}(u,v)$ with domain $D$, then the surface integral of a measurable function $F$ on $M$ is defined and denoted as:

$\int _{M}F\,dS:=\int \int _{D}(F\circ {\textbf {r}})|{\textbf {r}}_{u}\times {\textbf {r}}_{v}|\,dudv$

If $F:M\to \mathbb {R} ^{3}$ is vector-valued, then we define

$\int _{M}F\cdot dS:=\int _{M}(F\cdot {\textbf {n}})\,dS$

where ${\textbf {n}}$ is an outward unit normal vector to $M$. Since ${\textbf {n}}={\frac {{\textbf {r}}_{u}\times {\textbf {r}}_{v}}{|{\textbf {r}}_{u}\times {\textbf {r}}_{v}|}}$, we have:

$\int _{M}F\cdot dS=\int \int _{D}(F\circ {\textbf {r}})\cdot ({\textbf {r}}_{u}\times {\textbf {r}}_{v})\,dudv=\int \int _{D}\det(F\circ {\textbf {r}},{\textbf {r}}_{u},{\textbf {r}}_{v})\,dudv.$

Vector analysis

Tangent vectors and vector fields

Let $c:[0,1]\to \mathbb {R} ^{n}$ be a differentiable curve. Then the tangent vector to the curve $c$ at $t$ is a vector $v$ at the point $c(t)$ whose components are given as:

$v=(c_{1}'(t),\dots ,c_{n}'(t))$.

For example, if $c(t)=(a\cos(t),a\sin(t),bt),a>0,b>0$ is a helix, then the tangent vector at t is:

$c'(t)=(-a\sin(t),a\cos(t),b).$

It corresponds to the intuition that the a point on the helix moves up in a constant speed.

If $M\subset \mathbb {R} ^{n}$ is a differentiable curve or surface, then the tangent space to $M$ at a point p is the set of all tangent vectors to the differentiable curves $c:[0,1]\to M$ with $c(0)=p$.

A vector field X is an assignment to each point p in M a tangent vector $X_{p}$ to M at p such that the assignment varies smoothly.

Differential forms

The dual notion of a vector field is a differential form. Given an open subset $M$ in $\mathbb {R} ^{n}$, by definition, a differential 1-form (often just 1-form) $\omega$ is an assignment to a point $p$ in $M$ a linear functional $\omega _{p}$ on the tangent space $T_{p}M$ to $M$ at $p$ such that the assignment varies smoothly. For a (real or complex-valued) smooth function $f$, define the 1-form $df$ by: for a tangent vector $v$ at $p$,

$df_{p}(v)=v(f)$

where $v(f)$ denotes the directional derivative of $f$ in the direction $v$ at $p$. For example, if $x_{i}$ is the $i$-th coordinate function, then $dx_{i,p}(v)=v_{i}$; i.e., $dx_{i,p}$ are the dual basis to the standard basis on $T_{p}M$. Then every differential 1-form $\omega$ can be written uniquely as

$\omega =f_{1}\,dx_{1}+\cdots +f_{n}\,dx_{n}$

for some smooth functions $f_{1},\dots ,f_{n}$ on $M$ (since, for every point $p$, the linear functional $\omega _{p}$ is a unique linear combination of $dx_{i}$ over real numbers). More generally, a differential k-form is an assignment to a point $p$ in $M$ a vector $\omega _{p}$ in the $k$-th exterior power $\bigwedge ^{k}T_{p}^{*}M$ of the dual space $T_{p}^{*}M$ of $T_{p}M$ such that the assignment varies smoothly. In particular, a 0-form is the same as a smooth function. Also, any $k$-form $\omega$ can be written uniquely as:

$\omega =\sum _{i_{1}<\cdots <i_{k}}f_{i_{1}\dots i_{k}}\,dx_{i_{1}}\wedge \cdots \wedge dx_{i_{k}}$

for some smooth functions $f_{i_{1}\dots i_{k}}$.

Like a smooth function, we can differentiate and integrate differential forms. If $f$ is a smooth function, then $df$ can be written as:

$df=\sum _{i=1}^{n}{\frac {\partial f}{\partial x_{i}}}\,dx_{i}$

since, for $v=\partial /\partial x_{j}|_{p}$, we have: $df_{p}(v)={\frac {\partial f}{\partial x_{j}}}(p)=\sum _{i=1}^{n}{\frac {\partial f}{\partial x_{i}}}(p)\,dx_{i}(v)$. Note that, in the above expression, the left-hand side (whence the right-hand side) is independent of coordinates $x_{1},\dots ,x_{n}$; this property is called the invariance of differential.

The operation $d$ is called the exterior derivative and it extends to any differential forms inductively by the requirement (Leibniz rule)

$d(\alpha \wedge \beta )=d\alpha \wedge \beta +(-1)^{p}\alpha \wedge d\beta .$

where $\alpha ,\beta$ are a p-form and a q-form.

The exterior derivative has the important property that $d\circ d=0$; that is, the exterior derivative $d$ of a differential form $d\omega$ is zero. This property is a consequence of the symmetry of second derivatives (mixed partials are equal).

Boundary and orientation

A circle can be oriented clockwise or counterclockwise. Mathematically, we say that a subset $M$ of $\mathbb {R} ^{n}$ is oriented if there is a consistent choice of normal vectors to $M$ that varies continuously. For example, a circle or, more generally, an n-sphere can be oriented; i.e., orientable. On the other hand, a Möbius strip (a surface obtained by identified by two opposite sides of the rectangle in a twisted way) cannot oriented: if we start with a normal vector and travel around the strip, the normal vector at end will point to the opposite direction.

Proposition — A bounded differentiable region $M$ in $\mathbb {R} ^{n}$ of dimension $k$ is oriented if and only if there exists a nowhere-vanishing $k$-form on $M$ (called a volume form).

The proposition is useful because it allows us to give an orientation by giving a volume form.

Integration of differential forms

If $\omega =f\,dx_{1}\wedge \cdots \wedge dx_{n}$ is a differential n-form on an open subset M in $\mathbb {R} ^{n}$ (any n-form is that form), then the integration of it over $M$ with the standard orientation is defined as:

$\int _{M}\omega =\int _{M}f\,dx_{1}\cdots dx_{n}.$

If M is given the orientation opposite to the standard one, then $\int _{M}\omega$ is defined as the negative of the right-hand side.

Then we have the fundamental formula relating exterior derivative and integration:

Stokes' formula — For a bounded region $M$ in $\mathbb {R} ^{n}$ of dimension $k$ whose boundary is a union of finitely many $C^{1}$-subsets, if $M$ is oriented, then

$\int _{\partial M}\omega =\int _{M}d\omega$

for any differential $(k-1)$-form $\omega$ on the boundary $\partial M$ of $M$.

Here is a sketch of proof of the formula. If $f$ is a smooth function on $\mathbb {R} ^{n}$ with compact support, then we have:

$\int d(f\omega )=0$

(since, by the fundamental theorem of calculus, the above can be evaluated on boundaries of the set containing the support.) On the other hand,

$\int d(f\omega )=\int df\wedge \omega +\int f\,d\omega .$

Let $f$ approach the characteristic function on $M$. Then the second term on the right goes to $\int _{M}d\omega$ while the first goes to $-\int _{\partial M}\omega$, by the argument similar to proving the fundamental theorem of calculus. $\square$

The formula generalizes the fundamental theorem of calculus as well as Stokes' theorem in multivariable calculus. Indeed, if $M=[a,b]$ is an interval and $\omega =f$, then $d\omega =f'\,dx$ and the formula says:

$\int _{M}f'\,dx=f(b)-f(a)$.

Similarly, if $M$ is an oriented bounded surface in $\mathbb {R} ^{3}$ and $\omega =f\,dx+g\,dy+h\,dz$, then $d(f\,dx)=df\wedge dx={\frac {\partial f}{\partial y}}\,dy\wedge dx+{\frac {\partial f}{\partial z}}\,dz\wedge dx$ and similarly for $d(g\,dy)$ and $d(g\,dy)$. Collecting the terms, we thus get:

$d\omega =\left({\frac {\partial h}{\partial y}}-{\frac {\partial g}{\partial z}}\right)dy\wedge dz+\left({\frac {\partial f}{\partial z}}-{\frac {\partial h}{\partial x}}\right)dz\wedge dx+\left({\frac {\partial g}{\partial x}}-{\frac {\partial f}{\partial y}}\right)dx\wedge dy.$

Then, from the definition of the integration of $\omega$, we have $\int _{M}d\omega =\int _{M}(\nabla \times F)\cdot dS$ where $F=(f,g,h)$ is the vector-valued function and $\nabla =\left({\frac {\partial }{\partial x}},{\frac {\partial }{\partial y}},{\frac {\partial }{\partial z}}\right)$. Hence, Stokes’ formula becomes

$\int _{M}(\nabla \times F)\cdot dS=\int _{\partial M}(f\,dx+g\,dy+h\,dz),$

which is the usual form of the Stokes' theorem on surfaces. Green’s theorem is also a special case of Stokes’ formula.

Stokes' formula also yields a general version of Cauchy's integral formula. To state and prove it, for the complex variable $z=x+iy$ and the conjugate ${\bar {z}}$, let us introduce the operators

${\frac {\partial }{\partial z}}={\frac {1}{2}}\left({\frac {\partial }{\partial x}}-i{\frac {\partial }{\partial y}}\right),\,{\frac {\partial }{\partial {\bar {z}}}}={\frac {1}{2}}\left({\frac {\partial }{\partial x}}+i{\frac {\partial }{\partial y}}\right).$

In these notations, a function $f$ is holomorphic (complex-analytic) if and only if ${\frac {\partial f}{\partial {\bar {z}}}}=0$ (the Cauchy–Riemann equations). Also, we have:

$df={\frac {\partial f}{\partial z}}dz+{\frac {\partial f}{\partial {\bar {z}}}}d{\bar {z}}.$

Let $D_{\epsilon }=\{z\in \mathbb {C} \mid \epsilon <|z-z_{0}|<r\}$ be a punctured disk with center $z_{0}$. Since $1/(z-z_{0})$ is holomorphic on $D_{\epsilon }$, We have:

$d\left({\frac {f}{z-z_{0}}}dz\right)={\frac {\partial f}{\partial {\bar {z}}}}{\frac {d{\bar {z}}\wedge dz}{z-z_{0}}}$.

By Stokes’ formula,

$\int _{D_{\epsilon }}{\frac {\partial f}{\partial {\bar {z}}}}{\frac {d{\bar {z}}\wedge dz}{z-z_{0}}}=\left(\int _{|z-z_{0}|=r}-\int _{|z-z_{0}|=\epsilon }\right){\frac {f}{z-z_{0}}}dz.$

Letting $\epsilon \to 0$ we then get:

$2\pi i\,f(z_{0})=\int _{|z-z_{0}|=r}{\frac {f}{z-z_{0}}}dz+\int _{|z-z_{0}|\leq r}{\frac {\partial f}{\partial {\bar {z}}}}{\frac {dz\wedge d{\bar {z}}}{z-z_{0}}}.$

Winding numbers and Poincaré lemma

A differential form $\omega$ is called closed if $d\omega =0$ and is called exact if $\omega =d\eta$ for some differential form $\eta$ (often called a potential). Since $d\circ d=0$, an exact form is closed. But the converse does not hold in general; there might be a non-exact closed form. A classic example of such a form is:

$\omega ={\frac {-y}{x^{2}+y^{2}}}+{\frac {x}{x^{2}+y^{2}}}$,

which is a differential form on $\mathbb {R} ^{2}-0$. Suppose we switch to polar coordinates: $x=r\cos \theta ,y=r\sin \theta$ where $r={\sqrt {x^{2}+y^{2}}}$. Then

$\omega =r^{-2}(-r\sin \theta \,dx+r\cos \theta \,dy)=d\theta .$

This does not show that $\omega$ is exact: the trouble is that $\theta$ is not a well-defined continuous function on $\mathbb {R} ^{2}-0$. Since any function $f$ on $\mathbb {R} ^{2}-0$ with $df=\omega$ differ from $\theta$ by constant, this means that $\omega$ is not exact. The calculation, however, shows that $\omega$ is exact, for example, on $\mathbb {R} ^{2}-\{x=0\}$ since we can take $\theta =\arctan(y/x)$ there.

There is a result (Poincaré lemma) that gives a condition that guarantees closed forms are exact. To state it, we need some notions from topology. Given two continuous maps $f,g:X\to Y$ between subsets of $\mathbb {R} ^{m},\mathbb {R} ^{n}$ (or more generally topological spaces), a homotopy from $f$ to $g$ is a continuous function $H:X\times [0,1]\to Y$ such that $f(x)=H(x,0)$ and $g(x)=H(x,1)$. Intuitively, a homotopy is a continuous variation of one function to another. A loop in a set $X$ is a curve whose starting point coincides with the end point; i.e., $c:[0,1]\to X$ such that $c(0)=c(1)$. Then a subset of $\mathbb {R} ^{n}$ is called simply connected if every loop is homotopic to a constant function. A typical example of a simply connected set is a disk $D=\{(x,y)\mid {\sqrt {x^{2}+y^{2}}}\leq r\}\subset \mathbb {R} ^{2}$. Indeed, given a loop $c:[0,1]\to D$, we have the homotopy $H:[0,1]^{2}\to D,\,H(x,t)=(1-t)c(x)+tc(0)$ from $c$ to the constant function $c(0)$. A punctured disk, on the other hand, is not simply connected.

Poincaré lemma — If $M$ is a simply connected open subset of $\mathbb {R} ^{n}$, then each closed 1-form on $M$ is exact.

Geometry of curves and surfaces

Moving frame

Vector fields $E_{1},\dots ,E_{3}$ on $\mathbb {R} ^{3}$ are called a frame field if they are orthogonal to each other at each point; i.e., $E_{i}\cdot E_{j}=\delta _{ij}$ at each point. The basic example is the standard frame $U_{i}$; i.e., $U_{i}(x)$ is a standard basis for each point $x$ in $\mathbb {R} ^{3}$. Another example is the cylindrical frame

$E_{1}=\cos \theta U_{1}+\sin \theta U_{2},\,E_{2}=-\sin \theta U_{1}+\cos \theta U_{2},\,E_{3}=U_{3}.$

For the study of the geometry of a curve, the important frame to use is a Frenet frame $T,N,B$ on a unit-speed curve $\beta :I\to \mathbb {R} ^{3}$ given as:

The Gauss–Bonnet theorem

The Gauss–Bonnet theorem relates the topology of a surface and its geometry.

The Gauss–Bonnet theorem —  For each bounded surface $M$ in $\mathbb {R} ^{3}$, we have:

$2\pi \,\chi (M)=\int _{M}K\,dS$

where $\chi (M)$ is the Euler characteristic of $M$ and $K$ the curvature.

Calculus of variations

Method of Lagrange multiplier

Lagrange multiplier —  Let $g:U\to \mathbb {R} ^{r}$ be a differentiable function from an open subset of $\mathbb {R} ^{n}$ such that $g'$ has rank $r$ at every point in $g^{-1}(0)$. For a differentiable function $f:\mathbb {R} ^{n}\to \mathbb {R}$, if $f$ attains either a maximum or minimum at a point $p$ in $g^{-1}(0)$, then there exists real numbers $\lambda _{1},\dots ,\lambda _{r}$ such that

$\nabla f(p)=\lambda _{i}\sum _{i=1}^{r}\nabla g_{i}(p)$.

In other words, $p$ is a stationary point of $f-\sum _{1}^{r}\lambda _{i}g_{i}$.

The set $g^{-1}(0)$ is usually called a constraint.

Example: Suppose we want to find the minimum distance between the circle $x^{2}+y^{2}=1$ and the line $x+y=4$. That means that we want to minimize the function $f(x,y,u,v)=(x-u)^{2}+(y-v)^{2}$, the square distance between a point $(x,y)$ on the circle and a point $(u,v)$ on the line, under the constraint $g=(x^{2}+y^{2}-1,u+v-4)$. We have:

$\nabla f=(2(x-u),2(y-v),-2(x-u),-2(y-v)).$

$\nabla g_{1}=(2x,2y,0,0),\nabla g_{2}=(0,0,1,1).$

Since the Jacobian matrix of $g$ has rank 2 everywhere on $g^{-1}(0)$, the Lagrange multiplier gives:

$x-u=\lambda _{1}x,\,y-v=\lambda _{1}y,\,2(x-u)=-\lambda _{2},\,2(y-v)=-\lambda _{2}.$

If $\lambda _{1}=0$, then $x=u,y=v$, not possible. Thus, $\lambda _{1}\neq 0$ and

$x={\frac {x-u}{\lambda _{1}}},\,y={\frac {y-v}{\lambda _{1}}}.$

From this, it easily follows that $x=y=1/{\sqrt {2}}$ and $u=v=2$. Hence, the minimum distance is $2{\sqrt {2}}-1$ (as a minimum distance clearly exists).

Here is an application to linear algebra. Let $V$ be a finite-dimensional real vector space and $T:V\to V$ a self-adjoint operator. We shall show $V$ has a basis consisting of eigenvectors of $T$ (i.e., $T$ is diagonalizable) by induction on the dimension of $V$. Choosing a basis on $V$ we can identify $V=\mathbb {R} ^{n}$ and $T$ is represented by the matrix $[a_{ij}]$. Consider the function $f(x)=(Tx,x)$, where the bracket means the inner product. Then $\nabla f=2(\sum a_{1i}x_{i},\dots ,\sum a_{ni}x_{i})$. On the other hand, for $g=\sum x_{i}^{2}-1$, since $g^{-1}(0)$ is compact, $f$ attains a maximum or minimum at a point $u$ in $g^{-1}(0)$. Since $\nabla g=2(x_{1},\dots ,x_{n})$, by Lagrange multiplier, we find a real number $\lambda$ such that $2\sum _{i}a_{ji}u_{i}=2\lambda u_{j},1\leq j\leq n.$ But that means $Tu=\lambda u$. By inductive hypothesis, the self-adjoint operator $T:W\to W$, $W$ the orthogonal complement to $u$, has a basis consisting of eigenvectors. Hence, we are done. $\square$.

Weak derivatives

Up to measure-zero sets, two functions can be determined to be equal or not by means of integration against other functions (called test functions). Namely, the following sometimes called the fundamental lemma of calculus of variations:

Lemma — If $f,g$ are locally integrable functions on an open subset $M\subset \mathbb {R} ^{n}$ such that

$\int (f-g)\varphi \,dx=0$

for every $\varphi \in C_{c}^{\infty }(M)$ (called a test function). Then $f=g$ almost everywhere. If, in addition, $f,g$ are continuous, then $f=g$.

Given a continuous function $f$, by the lemma, a continuously differentiable function $u$ is such that ${\frac {\partial u}{\partial x_{i}}}=f$ if and only if

$\int {\frac {\partial u}{\partial x_{i}}}\varphi \,dx=\int f\varphi \,dx$

for every $\varphi \in C_{c}^{\infty }(M)$. But, by integration by parts, the partial derivative on the left-hand side of $u$ can be moved to that of $\varphi$; i.e.,

$-\int u{\frac {\partial \varphi }{\partial x_{i}}}\,dx=\int f\varphi \,dx$

where there is no boundary term since $\varphi$ has compact support. Now the key point is that this expression makes sense even if $u$ is not necessarily differentiable and thus can be used to give sense to a derivative of such a function.

Note each locally integrable function $u$ defines the linear functional $\varphi \mapsto \int u\varphi \,dx$ on $C_{c}^{\infty }(M)$ and, moreover, each locally integrable function can be identified with such linear functional, because of the early lemma. Hence, quite generally, if $u$ is a linear functional on $C_{c}^{\infty }(M)$, then we define ${\frac {\partial u}{\partial x_{i}}}$ to be the linear functional $\varphi \mapsto -\left\langle u,{\frac {\partial \varphi }{\partial x_{i}}}\right\rangle$ where the bracket means $\langle \alpha ,\varphi \rangle =\alpha (\varphi )$. It is then called the weak derivative of $u$ with respect to $x_{i}$. If $u$ is continuously differentiable, then the weak derivate of it coincides with the usual one; i.e., the linear functional ${\frac {\partial u}{\partial x_{i}}}$ is the same as the linear functional determined by the usual partial derivative of $u$ with respect to $x_{i}$. A usual derivative is often then called a classical derivative. When a linear functional on $C_{c}^{\infty }(M)$ is continuous with respect to a certain topology on $C_{c}^{\infty }(M)$, such a linear functional is called a distribution, an example of a generalized function.

A classic example of a weak derivative is that of the Heaviside function $H$, the characteristic function on the interval $(0,\infty )$. For every test function $\varphi$, we have:

$\langle H',\varphi \rangle =-\int _{0}^{\infty }\varphi '\,dx=\varphi (0).$

Let $\delta _{a}$ denote the linear functional $\varphi \mapsto \varphi (a)$, called the Dirac delta function (although not exactly a function). Then the above can be written as:

$H'=\delta _{0}.$

Cauchy's integral formula has a similar interpretation in terms of weak derivatives. For the complex variable $z=x+iy$, let $E_{z_{0}}(z)={\frac {1}{\pi (z-z_{0})}}$. For a test function $\varphi$, if the disk $|z-z_{0}|\leq r$ contains the support of $\varphi$, by Cauchy's integral formula, we have:

$\varphi (z_{0})={1 \over 2\pi i}\int {\frac {\partial \varphi }{\partial {\bar {z}}}}{\frac {dz\wedge d{\bar {z}}}{z-z_{0}}}.$

Since $dz\wedge d{\bar {z}}=-2idx\wedge dy$, this means:

$\varphi (z_{0})=-\int E_{z_{0}}{\frac {\partial \varphi }{\partial {\bar {z}}}}dxdy=\left\langle {\frac {\partial E_{z_{0}}}{\partial {\bar {z}}}},\varphi \right\rangle ,$

or

${\frac {\partial E_{z_{0}}}{\partial {\bar {z}}}}=\delta _{z_{0}}.$

In general, a generalized function is called a fundamental solution for a linear partial differential operator if the application of the operator to it is the Dirac delta. Hence, the above says $E_{z_{0}}$ is the fundamental solution for the differential operator $\partial /\partial {\bar {z}}$.

See also: Limit of distributions

Hamilton–Jacobi theory

Main article: Hamilton–Jacobi equation

Calculus on manifolds

Definition of a manifold

This section requires some background in general topology.

A manifold is a Hausdorff topological space that is locally modeled by an Euclidean space. By definition, an atlas of a topological space $M$ is a set of maps $\varphi _{i}:U_{i}\to \mathbb {R} ^{n}$, called charts, such that

• $U_{i}$ are an open cover of $M$; i.e., each $U_{i}$ is open and $M=\cup _{i}U_{i}$,
• $\varphi _{i}:U_{i}\to \varphi _{i}(U_{i})$ is a homeomorphism and
• $\varphi _{j}\circ \varphi _{i}^{-1}:\varphi _{i}(U_{i}\cap U_{j})\to \varphi _{j}(U_{i}\cap U_{j})$ is smooth; thus a diffeomorphism.

By definition, a manifold is a second-countable Hausdorff topological space with a maximal atlas (called a differentiable structure); "maximal" means that it is not contained in strictly larger atlas. The dimension of the manifold $M$ is the dimension of the model Euclidean space $\mathbb {R} ^{n}$; namely, $n$ and a manifold is called an n-manifold when it has dimension n. A function on a manifold $M$ is said to be smooth if $f|_{U}\circ \varphi ^{-1}$ is smooth on $\varphi (U)$ for each chart $\varphi :U\to \mathbb {R} ^{n}$ in the differentiable structure.

A manifold is paracompact; this has an implication that it admits a partition of unity subordinate to a given open cover.

If $\mathbb {R} ^{n}$ is replaced by an upper half-space $\mathbb {H} ^{n}$, then we get the notion of a manifold-with-boundary. The set of points that map to the boundary of $\mathbb {H} ^{n}$ under charts is denoted by $\partial M$ and is called the boundary of $M$. This boundary may not be the topological boundary of $M$. Since the interior of $\mathbb {H} ^{n}$ is diffeomorphic to $\mathbb {R} ^{n}$, a manifold is a manifold-with-boundary with empty boundary.

The next theorem furnishes many examples of manifolds.

Theorem —  Let $g:U\to \mathbb {R} ^{r}$ be a differentiable map from an open subset $U\subset \mathbb {R} ^{n}$ such that $g'(p)$ has rank $r$ for every point $p$ in $g^{-1}(0)$. Then the zero set $g^{-1}(0)$ is an $(n-r)$-manifold.

For example, for $g(x)=x_{1}^{2}+\cdots +x_{n+1}^{2}-1$, the derivative $g'(x)={\begin{bmatrix}2x_{1}&2x_{2}&\cdots &2x_{n+1}\end{bmatrix}}$ has rank one at every point $p$ in $g^{-1}(0)$. Hence, the n-sphere $g^{-1}(0)$ is an n-manifold.

The theorem is proved as a corollary of the inverse function theorem.

Many familiar manifolds are subsets of $\mathbb {R} ^{n}$. The next theoretically important result says that there is no other kind of manifolds. An immersion is a smooth map whose differential is injective. An embedding is an immersion that is homeomorphic (thus diffeomorphic) to the image.

Whitney's embedding theorem — Each $k$-manifold can be embedded into $\mathbb {R} ^{2k}$.

The proof that a manifold can be embedded into $\mathbb {R} ^{N}$ for some N is considerably easier and can be readily given here. It is known that a manifold has a finite atlas $\{\varphi _{i}:U_{i}\to \mathbb {R} ^{n}\mid 1\leq i\leq r\}$. Let $\lambda _{i}$ be smooth functions such that $\operatorname {Supp} (\lambda _{i})\subset U_{i}$ and $\{\lambda _{i}=1\}$ cover $M$ (e.g., a partition of unity). Consider the map

$f=(\lambda _{1}\varphi _{1},\dots ,\lambda _{r}\varphi _{r},\lambda _{1},\dots ,\lambda _{r}):M\to \mathbb {R} ^{(k+1)r}$

It is easy to see that $f$ is an injective immersion. It may not be an embedding. To fix that, we shall use:

$(f,g):M\to \mathbb {R} ^{(k+1)r+1}$

where $g$ is a smooth proper map. The existence of a smooth proper map is a consequence of a partition of unity. See [1] for the rest of the proof in the case of an immersion. $\square$

Nash's embedding theorem says that, if $M$ is equipped with a Riemannian metric, then the embedding can be taken to be isometric with an expense of increasing $2k$; for this, see this T. Tao's blog.

Tubular neighborhood and transversality

A technically important result is:

Tubular neighborhood theorem — Let M be a manifold and $N\subset M$ a compact closed submanifold. Then there exists a neighborhood $U$ of $N$ such that $U$ is diffeomorphic to the normal bundle $\nu _{N}=TM|_{N}/TN$ to $i:N\hookrightarrow M$ and $N$ corresponds to the zero section of $\nu _{i}$ under the diffeomorphism.

This can be proved by putting a Riemannian metric on the manifold $M$. Indeed, the choice of metric makes the normal bundle $\nu _{i}$ a complementary bundle to $TN$; i.e., $TM|_{N}$ is the direct sum of $TN$ and $\nu _{N}$. Then, using the metric, we have the exponential map $\exp :U\to V$ for some neighborhood $U$ of $N$ in the normal bundle $\nu _{N}$ to some neighborhood $V$ of $N$ in $M$. The exponential map here may not be injective but it is possible to make it injective (thus diffeomorphic) by shrinking $U$ (for now, see [2]).

Integration on manifolds and distribution densities

The starting point for the topic of integration on manifolds is that there is no invariant way to integrate functions on manifolds. This may be obvious if we asked: what is an integration of functions on a finite-dimensional real vector space? (In contrast, there is an invariant way to do differentiation since, by definition, a manifold comes with a differentiable structure). There are several ways to introduce integration theory to manifolds:

• Integrate differential forms.
• Do integration against some measure.
• Equip a manifold with a Riemannian metric and do integration against such a metric.

For example, if a manifold is embedded into an Euclidean space $\mathbb {R} ^{n}$, then it acquires the Lebesgue measure restricting from the ambient Euclidean space and then the second approach works. The first approach is fine in many situations but it requires the manifold to be oriented (and there is a non-orientable manifold that is not pathological). The third approach generalizes and that gives rise to the notion of a density.

Generalizations

Extensions to infinite-dimensional normed spaces

The notions like differentiability extend to normed spaces.

See also

• Differential geometry of surfaces
• Integration along fibers
• Lusin's theorem
• Density on a manifold

References

• do Carmo, Manfredo P. (1976), Differential Geometry of Curves and Surfaces, Prentice-Hall, ISBN 978-0-13-212589-5
• Edwards, Charles Henry (1994) [1973], Advanced Calculus of Several Variables, Mineola, New York: Dover Publications, ISBN 0-486-68336-2
• Folland, Gerald, Real Analysis: Modern Techniques and Their Applications (2nd ed.)
• Cartan, Henri (1971), Calcul Differentiel (in French), Hermann, ISBN 9780395120330
• Hirsch, Morris (1994), Differential Topology (2nd ed.), Springer-Verlag
• Hörmander, Lars (2015), The Analysis of Linear Partial Differential Operators I: Distribution Theory and Fourier Analysis, Classics in Mathematics (2nd ed.), Springer, ISBN 9783642614972
• Loomis, Lynn Harold; Sternberg, Shlomo (1968), Advanced Calculus, Addison-Wesley (revised 1990, Jones and Bartlett; reprinted 2014, World Scientific) [this text in particular discusses density]
• O'Neill, Barrett (2006), Elementary Differential Geometry (revised 2nd ed.), Amsterdam: Elsevier/Academic Press, ISBN 0-12-088735-5
• Rudin, Walter (1976) [1953], Principles of Mathematical Analysis (3rd ed.), New York: McGraw Hill, pp. 204–299, ISBN 978-0-07-054235-8
• Spivak, Michael (1965). Calculus on Manifolds: A Modern Approach to Classical Theorems of Advanced Calculus. San Francisco: Benjamin Cummings. ISBN 0-8053-9021-9.