Chain rule

This article is about the chain rule in calculus. For the chain rule in probability theory, see Chain rule (probability). For other uses, see Chain rule (disambiguation).

Part of a series of articles about
Calculus
Fundamental theorem Limits of functions Continuity Mean value theorem Rolle's theorem
Differential Definitions Derivative (generalizations) Differential infinitesimal of a function total Concepts Differentiation notation Second derivative Third derivative Change of variables Implicit differentiation Related rates Taylor's theorem Rules and identities Sum Product Chain Power Quotient Inverse General Leibniz Faà di Bruno's formula
Integral Lists of integrals Definitions Antiderivative Integral (improper) Riemann integral Lebesgue integration Contour integration Integration by Parts Discs Cylindrical shells Substitution (trigonometric) Partial fractions Order Reduction formulae
Series Geometric (arithmetico-geometric) Harmonic Alternating Power Binomial Taylor Convergence tests Summand limit (term test) Ratio Root Integral Direct comparison Limit comparison Alternating series Cauchy condensation Dirichlet Abel
Vector Gradient Divergence Curl Laplacian Directional derivative Identities Theorems Divergence Gradient Green's Kelvin–Stokes Stokes
Multivariable Formalisms Matrix Tensor Exterior Geometric Definitions Partial derivative Multiple integral Line integral Surface integral Volume integral Jacobian Hessian
Specialized Fractional Malliavin Stochastic Variations
Glossary of calculus Glossary of calculus

In calculus, the chain rule is a formula for computing the derivative of the composition of two or more functions. That is, if f and g are functions, then the chain rule expresses the derivative of their composition f ∘ g (the function which maps x to f(g(x)) ) in terms of the derivatives of f and g and the product of functions as follows:

$(f\circ g)'=(f'\circ g)\cdot g'.$

This may equivalently be expressed in terms of the variable. Let F = f ∘ g, or equivalently, F(x) = f(g(x)) for all x. Then one can also write

$F'(x)=f'(g(x))g'(x).$

The chain rule may be written in Leibniz's notation in the following way. If a variable z depends on the variable y, which itself depends on the variable x, so that y and z are therefore dependent variables, then z, via the intermediate variable of y, depends on x as well. The chain rule then states,

${\frac {dz}{dx}}={\frac {dz}{dy}}\cdot {\frac {dy}{dx}}.$

The two versions of the chain rule are related; if $z=f(y)$ and $y=g(x)$, then

${\frac {dz}{dx}}={\frac {dz}{dy}}\cdot {\frac {dy}{dx}}=f'(y)g'(x)=f'(g(x))g'(x).$

In integration, the counterpart to the chain rule is the substitution rule.

History

The chain rule seems to have first been used by Gottfried Wilhelm Leibniz. He used it to calculate the derivative of ${\sqrt {a+bz+cz^{2}}}$ as the composite of the square root function and the function $a+bz+cz^{2}$. He first mentioned it in a 1676 memoir (with a sign error in the calculation). The common notation of chain rule is due to Leibniz. Guillaume de l'Hôpital used the chain rule implicitly in his Analyse des infiniment petits. The chain rule does not appear in any of Leonhard Euler's analysis books, even though they were written over a hundred years after Leibniz's discovery.

One dimension

First example

Suppose that a skydiver jumps from an aircraft. Assume that t seconds after his jump, his height above sea level in meters is given by g(t) = 4000 − 4.9t². One model for the atmospheric pressure at a height h is f(h) = 101325 e^−0.0001h. These two equations can be differentiated and combined in various ways to produce the following data:

• g′(t) = −9.8t is the velocity of the skydiver at time t.
• f′(h) = −10.1325e^−0.0001h is the rate of change in atmospheric pressure with respect to height at the height h and is proportional to the buoyant force on the skydiver at h meters above sea level. (The true buoyant force depends on the volume of the skydiver.)
• (f ∘ g)(t) is the atmospheric pressure the skydiver experiences t seconds after his jump.
• (f ∘ g)′(t) is the rate of change in atmospheric pressure with respect to time at t seconds after the skydiver's jump and is proportional to the buoyant force on the skydiver at t seconds after his jump.

The chain rule gives a method for computing (f ∘ g)′(t) in terms of f′ and g′. While it is always possible to directly apply the definition of the derivative to compute the derivative of a composite function, this is usually very difficult. The utility of the chain rule is that it turns a complicated derivative into several easy derivatives.

The chain rule states that, under appropriate conditions,

$(f\circ g)'(t)=f'(g(t))\cdot g'(t).$

In this example, this equals

$(f\circ g)'(t)={\big (}{\mathord {-}}10.1325e^{-0.0001(4000-4.9t^{2})}{\big )}\cdot {\big (}{\mathord {-}}9.8t{\big )}.$

In the statement of the chain rule, f and g play slightly different roles because f′ is evaluated at g(t) whereas g′ is evaluated at t. This is necessary to make the units work out correctly. For example, suppose that we want to compute the rate of change in atmospheric pressure ten seconds after the skydiver jumps. This is (f ∘ g)′(10) and has units of pascals per second. The factor g′(10) in the chain rule is the velocity of the skydiver ten seconds after his jump, and it is expressed in meters per second. f′(g(10)) is the change in pressure with respect to height at the height g(10) and is expressed in pascals per meter. The product of f′(g(10)) and g′(10) therefore has the correct units of pascals per second. It is not possible to evaluate f anywhere else. For instance, because the 10 in the problem represents ten seconds, the expression f′(10) represents the change in pressure at a height of ten seconds, which is nonsense. Similarly, because g′(10) = −98 meters per second, the expression f′(g′(10)) represents the change in pressure at a height of −98 meters per second, which is also nonsense. However, g(10) is 3020 meters above sea level, the height of the skydiver ten seconds after his jump. This has the correct units for an input to f.

Statement

The simplest form of the chain rule is for real-valued functions of one real variable. It says that if g is a function that is differentiable at a point c (i.e. the derivative g′(c) exists) and f is a function that is differentiable at g(c), then the composite function f ∘ g is differentiable at c, and the derivative is

$(f\circ g)'(c)=f'(g(c))\cdot g'(c).$

The rule is sometimes abbreviated as

$(f\circ g)'=(f'\circ g)\cdot g'.$

If y = f(u) and u = g(x), then this abbreviated form is written in Leibniz notation as:

${\frac {dy}{dx}}={\frac {dy}{du}}\cdot {\frac {du}{dx}}.$

The points where the derivatives are evaluated may also be stated explicitly:

$\left.{\frac {dy}{dx}}\right|_{x=c}=\left.{\frac {dy}{du}}\right|_{u=g(c)}\cdot \left.{\frac {du}{dx}}\right|_{x=c}.$

Further examples

Absence of formulas

It may be possible to apply the chain rule even when there are no formulas for the functions which are being differentiated. This can happen when the derivatives are measured directly. Suppose that a car is driving up a tall mountain. The car's speedometer measures its speed directly. If the grade is known, then the rate of ascent can be calculated using trigonometry. Suppose that the car is ascending at 6999694444444444444♠2.5 km/h. Standard models for the Earth's atmosphere imply that the temperature drops about 7002279649999999999♠6.5 °C per kilometer ascended (called the lapse rate). To find the temperature drop per hour, we apply the chain rule. Let the function g(t) be the altitude of the car at time t, and let the function f(h) be the temperature h kilometers above sea level. f and g are not known exactly: For example, the altitude where the car starts is not known and the temperature on the mountain is not known. However, their derivatives are known: f′ is 3002349999999999999♠−6.5 °C/km, and g′ is 6999694444444444444♠2.5 km/h. The chain rule says that the derivative of the composite function is the product of the derivative of f and the derivative of g. This is 3002349999999999999♠−6.5 °C/km ⋅ 6999694444444444444♠2.5 km/h = 3002548611111111111♠−16.25 °C/h.

One of the reasons why this computation is possible is because f′ is a constant function. This is because the above model is very simple. A more accurate description of how the temperature near the car varies over time would require an accurate model of how the temperature varies at different altitudes. This model may not have a constant derivative. To compute the temperature change in such a model, it would be necessary to know g and not just g′, because without knowing g it is not possible to know where to evaluate f′.

Composites of more than two functions

The chain rule can be applied to composites of more than two functions. To take the derivative of a composite of more than two functions, notice that the composite of f, g, and h (in that order) is the composite of f with g ∘ h. The chain rule says that to compute the derivative of f ∘ g ∘ h, it is sufficient to compute the derivative of f and the derivative of g ∘ h. The derivative of f can be calculated directly, and the derivative of g ∘ h can be calculated by applying the chain rule again.

For concreteness, consider the function

$y=e^{\sin(x^{2})}.$

This can be decomposed as the composite of three functions:

${\begin{aligned}y&=f(u)=e^{u},\\[6pt]u&=g(v)=\sin v=\sin(x^{2}),\\[6pt]v&=h(x)=x^{2}.\end{aligned}}$

Their derivatives are:

${\begin{aligned}{\frac {dy}{du}}&=f'(u)=e^{u}=e^{\sin(x^{2})},\\[6pt]{\frac {du}{dv}}&=g'(v)=\cos v=\cos(x^{2}),\\[6pt]{\frac {dv}{dx}}&=h'(x)=2x.\end{aligned}}$

The chain rule says that the derivative of their composite at the point x = a is:

${\begin{aligned}(f\circ g\circ h)'(a)&=f'((g\circ h)(a))\cdot (g\circ h)'(a)\\[10pt]&=f'((g\circ h)(a))\cdot g'(h(a))\cdot h'(a)=(f'\circ g\circ h)(a)\cdot (g'\circ h)(a)\cdot h'(a).\end{aligned}}$

In Leibniz notation, this is:

${\frac {dy}{dx}}=\left.{\frac {dy}{du}}\right|_{u=g(h(a))}\cdot \left.{\frac {du}{dv}}\right|_{v=h(a)}\cdot \left.{\frac {dv}{dx}}\right|_{x=a},$

or for short,

${\frac {dy}{dx}}={\frac {dy}{du}}\cdot {\frac {du}{dv}}\cdot {\frac {dv}{dx}}.$

The derivative function is therefore:

${\frac {dy}{dx}}=e^{\sin(x^{2})}\cdot \cos(x^{2})\cdot 2x.$

Another way of computing this derivative is to view the composite function f ∘ g ∘ h as the composite of f ∘ g and h. Applying the chain rule to this situation gives:

$(f\circ g\circ h)'(a)=(f\circ g)'(h(a))\cdot h'(a)=f'(g(h(a)))\cdot g'(h(a))\cdot h'(a).$

This is the same as what was computed above. This should be expected because (f ∘ g) ∘ h = f ∘ (g ∘ h).

Sometimes it is necessary to differentiate an arbitrarily long composition of the form $f_{1}\circ f_{2}\circ \cdots \circ f_{n-1}\circ f_{n}$. In this case, define

$f_{a\,.\,.\,b}=f_{a}\circ f_{a+1}\circ \cdots \circ f_{b-1}\circ f_{b}$

where $f_{a\,.\,.\,a}=f_{a}$ and $f_{a\,.\,.\,b}(x)=x$ when $b<a$. Then the chain rule takes the form

$Df_{1\,.\,.\,n}=(Df_{1}\circ f_{2\,.\,.\,n})(Df_{2}\circ f_{3\,.\,.\,n})\cdots (Df_{n-1}\circ f_{n\,.\,.\,n})Df_{n}=\prod _{k=1}^{n}\left[Df_{k}\circ f_{(k+1)\,.\,.\,n}\right]$

or, in the Lagrange notation,

$f_{1\,.\,.\,n}'(x)=f_{1}'\left(f_{2\,.\,.\,n}(x)\right)\;f_{2}'\left(f_{3\,.\,.\,n}(x)\right)\cdots f_{n-1}'\left(f_{n\,.\,.\,n}(x)\right)\;f_{n}'(x)=\prod _{k=1}^{n}f_{k}'\left(f_{(k+1\,.\,.\,n)}(x)\right)$

Quotient rule

See also: Quotient rule

The chain rule can be used to derive some well-known differentiation rules. For example, the quotient rule is a consequence of the chain rule and the product rule. To see this, write the function f(x)/g(x) as the product f(x) · 1/g(x). First apply the product rule:

${\begin{aligned}{\frac {d}{dx}}\left({\frac {f(x)}{g(x)}}\right)&={\frac {d}{dx}}\left(f(x)\cdot {\frac {1}{g(x)}}\right)\\&=f'(x)\cdot {\frac {1}{g(x)}}+f(x)\cdot {\frac {d}{dx}}\left({\frac {1}{g(x)}}\right).\end{aligned}}$

To compute the derivative of 1/g(x), notice that it is the composite of g with the reciprocal function, that is, the function that sends x to 1/x. The derivative of the reciprocal function is −1/x². By applying the chain rule, the last expression becomes:

$f'(x)\cdot {\frac {1}{g(x)}}+f(x)\cdot \left(-{\frac {1}{g(x)^{2}}}\cdot g'(x)\right)={\frac {f'(x)g(x)-f(x)g'(x)}{g(x)^{2}}},$

which is the usual formula for the quotient rule.

Derivatives of inverse functions

Main article: Inverse functions and differentiation

Suppose that y = g(x) has an inverse function. Call its inverse function f so that we have x = f(y). There is a formula for the derivative of f in terms of the derivative of g. To see this, note that f and g satisfy the formula

$f(g(x))=x.$

Because the functions f(g(x)) and x are equal, their derivatives must be equal. The derivative of x is the constant function with value 1, and the derivative of f(g(x)) is determined by the chain rule. Therefore, we have:

$f'(g(x))g'(x)=1.$

To express f′ as a function of an independent variable y, we substitute f(y) for x wherever it appears. Then we can solve for f′.

${\begin{aligned}f'(g(f(y)))g'(f(y))&=1\\[5pt]f'(y)g'(f(y))&=1\\[5pt]f'(y)={\frac {1}{g'(f(y))}}.\end{aligned}}$

For example, consider the function g(x) = e^x. It has an inverse f(y) = ln y. Because g′(x) = e^x, the above formula says that

${\frac {d}{dy}}\ln y={\frac {1}{e^{\ln y}}}={\frac {1}{y}}.$

This formula is true whenever g is differentiable and its inverse f is also differentiable. This formula can fail when one of these conditions is not true. For example, consider g(x) = x³. Its inverse is f(y) = y^1/3, which is not differentiable at zero. If we attempt to use the above formula to compute the derivative of f at zero, then we must evaluate 1/g′(f(0)). f(0) = 0 and g′(0) = 0, so we must evaluate 1/0, which is undefined. Therefore, the formula fails in this case. This is not surprising because f is not differentiable at zero.

Higher derivatives

Faà di Bruno's formula generalizes the chain rule to higher derivatives. Assuming that y = f(u) and u = g(x), then the first few derivatives are:

${\begin{aligned}{\frac {dy}{dx}}&={\frac {dy}{du}}{\frac {du}{dx}}\\[4pt]{\frac {d^{2}y}{dx^{2}}}&={\frac {d^{2}y}{du^{2}}}\left({\frac {du}{dx}}\right)^{2}+{\frac {dy}{du}}{\frac {d^{2}u}{dx^{2}}}\\[4pt]{\frac {d^{3}y}{dx^{3}}}&={\frac {d^{3}y}{du^{3}}}\left({\frac {du}{dx}}\right)^{3}+3\,{\frac {d^{2}y}{du^{2}}}{\frac {du}{dx}}{\frac {d^{2}u}{dx^{2}}}+{\frac {dy}{du}}{\frac {d^{3}u}{dx^{3}}}\\[4pt]{\frac {d^{4}y}{dx^{4}}}&={\frac {d^{4}y}{du^{4}}}\left({\frac {du}{dx}}\right)^{4}+6\,{\frac {d^{3}y}{du^{3}}}\left({\frac {du}{dx}}\right)^{2}{\frac {d^{2}u}{dx^{2}}}+{\frac {d^{2}y}{du^{2}}}\left(4\,{\frac {du}{dx}}{\frac {d^{3}u}{dx^{3}}}+3\,\left({\frac {d^{2}u}{dx^{2}}}\right)^{2}\right)+{\frac {dy}{du}}{\frac {d^{4}u}{dx^{4}}}.\end{aligned}}$

Proofs

First proof

One proof of the chain rule begins with the definition of the derivative:

$(f\circ g)'(a)=\lim _{x\to a}{\frac {f(g(x))-f(g(a))}{x-a}}.$

Assume for the moment that g(x) does not equal g(a) for any x near a. Then the previous expression is equal to the product of two factors:

$\lim _{x\to a}{\frac {f(g(x))-f(g(a))}{g(x)-g(a)}}\cdot {\frac {g(x)-g(a)}{x-a}}.$

If g oscillates near a, then it might happen that no matter how close one gets to a, there is always an even closer x such that g(x) equals g(a). For example, this happens for g(x) = x²sin(1 / x) near the point a = 0. Whenever this happens, the above expression is undefined because it involves division by zero. To work around this, introduce a function Q as follows:

$Q(y)={\begin{cases}{\frac {f(y)-f(g(a))}{y-g(a)}},&y\neq g(a),\\f'(g(a)),&y=g(a).\end{cases}}$

We will show that the difference quotient for f ∘ g is always equal to:

$Q(g(x))\cdot {\frac {g(x)-g(a)}{x-a}}.$

Whenever g(x) is not equal to g(a), this is clear because the factors of g(x) − g(a) cancel. When g(x) equals g(a), then the difference quotient for f ∘ g is zero because f(g(x)) equals f(g(a)), and the above product is zero because it equals f′(g(a)) times zero. So the above product is always equal to the difference quotient, and to show that the derivative of f ∘ g at a exists and to determine its value, we need only show that the limit as x goes to a of the above product exists and determine its value.

To do this, recall that the limit of a product exists if the limits of its factors exist. When this happens, the limit of the product of these two factors will equal the product of the limits of the factors. The two factors are Q(g(x)) and (g(x) − g(a)) / (x − a). The latter is the difference quotient for g at a, and because g is differentiable at a by assumption, its limit as x tends to a exists and equals g′(a).

It remains to study Q(g(x)). Q is defined wherever f is. Furthermore, because f is differentiable at g(a) by assumption, Q is continuous at g(a). g is continuous at a because it is differentiable at a, and therefore Q ∘ g is continuous at a. So its limit as x goes to a exists and equals Q(g(a)), which is f′(g(a)).

This shows that the limits of both factors exist and that they equal f′(g(a)) and g′(a), respectively. Therefore, the derivative of f ∘ g at a exists and equals f′(g(a))g′(a).

Second proof

Another way of proving the chain rule is to measure the error in the linear approximation determined by the derivative. This proof has the advantage that it generalizes to several variables. It relies on the following equivalent definition of differentiability at a point: A function g is differentiable at a if there exists a real number g′(a) and a function ε(h) that tends to zero as h tends to zero, and furthermore

$g(a+h)-g(a)=g'(a)h+\varepsilon (h)h.$

Here the left-hand side represents the true difference between the value of g at a and at a + h, whereas the right-hand side represents the approximation determined by the derivative plus an error term.

In the situation of the chain rule, such a function ε exists because g is assumed to be differentiable at a. Again by assumption, a similar function also exists for f at g(a). Calling this function η, we have

$f(g(a)+k)-f(g(a))=f'(g(a))k+\eta (k)k.$

The above definition imposes no constraints on η(0), even though it is assumed that η(k) tends to zero as k tends to zero. If we set η(0) = 0, then η is continuous at 0.

Proving the theorem requires studying the difference f(g(a + h)) − f(g(a)) as h tends to zero. The first step is to substitute for g(a + h) using the definition of differentiability of g at a:

$f(g(a+h))-f(g(a))=f(g(a)+g'(a)h+\varepsilon (h)h)-f(g(a)).$

The next step is to use the definition of differentiability of f at g(a). This requires a term of the form f(g(a) + k) for some k. In the above equation, the correct k varies with h. Set k_h = g′(a) h + ε(h) h and the right hand side becomes f(g(a) + k_h) − f(g(a)). Applying the definition of the derivative gives:

$f(g(a)+k_{h})-f(g(a))=f'(g(a))k_{h}+\eta (k_{h})k_{h}.$

To study the behavior of this expression as h tends to zero, expand k_h. After regrouping the terms, the right-hand side becomes:

$f'(g(a))g'(a)h+[f'(g(a))\varepsilon (h)+\eta (k_{h})g'(a)+\eta (k_{h})\varepsilon (h)]h.$

Because ε(h) and η(k_h) tend to zero as h tends to zero, the first two bracketed terms tend to zero as h tends to zero. Applying the same theorem on products of limits as in the first proof, the third bracketed term also tends zero. Because the above expression is equal to the difference f(g(a + h)) − f(g(a)), by the definition of the derivative f ∘ g is differentiable at a and its derivative is f′(g(a)) g′(a).

The role of Q in the first proof is played by η in this proof. They are related by the equation:

$Q(y)=f'(g(a))+\eta (y-g(a)).$

The need to define Q at g(a) is analogous to the need to define η at zero.

Third proof

Constantin Carathéodory's alternative definition of the differentiability of a function can be used to give an elegant proof of the chain rule.

Under this definition, a function f is differentiable at a point a if and only if there is a function q, continuous at a and such that f(x) − f(a) = q(x)(x − a). There is at most one such function, and if f is differentiable at a then f '(a) = q(a).

Given the assumptions of the chain rule and the fact that differentiable functions and compositions of continuous functions are continuous, we have that there exist functions q, continuous at g(a) and r, continuous at a and such that,

$f(g(x))-f(g(a))=q(g(x))(g(x)-g(a))$

and

$g(x)-g(a)=r(x)(x-a).$

Therefore,

$f(g(x))-f(g(a))=q(g(x))r(x)(x-a),$

but the function given by h(x) = q(g(x))r(x) is continuous at a, and we get, for this a

$(f(g(a)))'=q(g(a))r(a)=f'(g(a))g'(a).$

A similar approach works for continuously differentiable (vector-)functions of many variables. This method of factoring also allows a unified approach to stronger forms of differentiability, when the derivative is required to be Lipschitz continuous, Hölder continuous, etc. Differentiation itself can be viewed as the polynomial remainder theorem (the little Bézout theorem, or factor theorem), generalized to an appropriate class of functions.

Proof via infinitesimals

If $y=f(x)$ and $x=g(t)$ then choosing infinitesimal $\Delta t\not =0$ we compute the corresponding $\Delta x=g(t+\Delta t)-g(t)$ and then the corresponding $\Delta y=f(x+\Delta x)-f(x)$, so that

${\frac {\Delta y}{\Delta t}}={\frac {\Delta y}{\Delta x}}{\frac {\Delta x}{\Delta t}}$

and applying the standard part we obtain

${\frac {dy}{dt}}={\frac {dy}{dx}}{\frac {dx}{dt}}$

which is the chain rule.

Higher dimensions

The simplest generalization of the chain rule to higher dimensions uses the total derivative. The total derivative is a linear transformation that captures how the function changes in all directions. Fix differentiable functions f : R^m → R^k and g : Rⁿ → R^m and a point a in Rⁿ. Let D_ag denote the total derivative of g at a and D_g(a)f denote the total derivative of f at g(a). These two derivatives are linear transformations Rⁿ → R^m and R^m → R^k, respectively, so they can be composed. The chain rule for total derivatives says that their composite is the total derivative of f ∘ g at a:

$D_{\mathbf {a} }(f\circ g)=D_{g(\mathbf {a} )}f\circ D_{\mathbf {a} }g,$

or for short,

$D(f\circ g)=Df\circ Dg.$

The higher-dimensional chain rule can be proved using a technique similar to the second proof given above.

Because the total derivative is a linear transformation, the functions appearing in the formula can be rewritten as matrices. The matrix corresponding to a total derivative is called a Jacobian matrix, and the composite of two derivatives corresponds to the product of their Jacobian matrices. From this perspective the chain rule therefore says:

$J_{f\circ g}(\mathbf {a} )=J_{f}(g(\mathbf {a} ))J_{g}(\mathbf {a} ),$

or for short,

$J_{f\circ g}=(J_{f}\circ g)J_{g}.$

That is, the Jacobian of a composite function is the product of the Jacobians of the composed functions (evaluated at the appropriate points).

The higher-dimensional chain rule is a generalization of the one-dimensional chain rule. If k, m, and n are 1, so that f : R → R and g : R → R, then the Jacobian matrices of f and g are 1 × 1. Specifically, they are:

${\begin{aligned}J_{g}(a)&={\begin{pmatrix}g'(a)\end{pmatrix}},\\J_{f}(g(a))&={\begin{pmatrix}f'(g(a))\end{pmatrix}}.\end{aligned}}$

The Jacobian of f ∘ g is the product of these 1 × 1 matrices, so it is f′(g(a))⋅g′(a), as expected from the one-dimensional chain rule. In the language of linear transformations, D_a(g) is the function which scales a vector by a factor of g′(a) and D_g(a)(f) is the function which scales a vector by a factor of f′(g(a)). The chain rule says that the composite of these two linear transformations is the linear transformation D_a(f ∘ g), and therefore it is the function that scales a vector by f′(g(a))⋅g′(a).

Another way of writing the chain rule is used when f and g are expressed in terms of their components as y = f(u) = (f₁(u), …, f_k(u)) and u = g(x) = (g₁(x), …, g_m(x)). In this case, the above rule for Jacobian matrices is usually written as:

${\frac {\partial (y_{1},\ldots ,y_{k})}{\partial (x_{1},\ldots ,x_{n})}}={\frac {\partial (y_{1},\ldots ,y_{k})}{\partial (u_{1},\ldots ,u_{m})}}{\frac {\partial (u_{1},\ldots ,u_{m})}{\partial (x_{1},\ldots ,x_{n})}}.$

The chain rule for total derivatives implies a chain rule for partial derivatives. Recall that when the total derivative exists, the partial derivative in the ith coordinate direction is found by multiplying the Jacobian matrix by the ith basis vector. By doing this to the formula above, we find:

${\frac {\partial (y_{1},\ldots ,y_{k})}{\partial x_{i}}}={\frac {\partial (y_{1},\ldots ,y_{k})}{\partial (u_{1},\ldots ,u_{m})}}{\frac {\partial (u_{1},\ldots ,u_{m})}{\partial x_{i}}}.$

Since the entries of the Jacobian matrix are partial derivatives, we may simplify the above formula to get:

${\frac {\partial (y_{1},\ldots ,y_{k})}{\partial x_{i}}}=\sum _{\ell =1}^{m}{\frac {\partial (y_{1},\ldots ,y_{k})}{\partial u_{\ell }}}{\frac {\partial u_{\ell }}{\partial x_{i}}}.$

More conceptually, this rule expresses the fact that a change in the x_i direction may change all of g₁ through g_m, and any of these changes may affect f.

In the special case where k = 1, so that f is a real-valued function, then this formula simplifies even further:

${\frac {\partial y}{\partial x_{i}}}=\sum _{\ell =1}^{m}{\frac {\partial y}{\partial u_{\ell }}}{\frac {\partial u_{\ell }}{\partial x_{i}}}.$

This can be rewritten as a dot product. Recalling that u = (g₁, …, g_m), the partial derivative ∂u / ∂x_i is also a vector, and the chain rule says that:

${\frac {\partial y}{\partial x_{i}}}=\nabla y\cdot {\frac {\partial \mathbf {u} }{\partial x_{i}}}.$

Example

Given u(x, y) = x² + 2y where x(r, t) = r sin(t) and y(r,t) = sin²(t), determine the value of ∂u / ∂r and ∂u / ∂t using the chain rule.

${\frac {\partial u}{\partial r}}={\frac {\partial u}{\partial x}}{\frac {\partial x}{\partial r}}+{\frac {\partial u}{\partial y}}{\frac {\partial y}{\partial r}}=(2x)(\sin(t))+(2)(0)=2r\sin ^{2}(t),$

and

${\begin{aligned}{\frac {\partial u}{\partial t}}&={\frac {\partial u}{\partial x}}{\frac {\partial x}{\partial t}}+{\frac {\partial u}{\partial y}}{\frac {\partial y}{\partial t}}\\&=(2x)(r\cos(t))+(2)(2\sin(t)\cos(t))\\&=(2r\sin(t))(r\cos(t))+4\sin(t)\cos(t)\\&=2(r^{2}+2)\sin(t)\cos(t)\\&=(r^{2}+2)\sin(2t).\end{aligned}}$

Higher derivatives of multivariable functions

Main article: Faà di Bruno's formula § Multivariate version

Faà di Bruno's formula for higher-order derivatives of single-variable functions generalizes to the multivariable case. If y = f(u) is a function of u = g(x) as above, then the second derivative of f ∘ g is:

${\frac {\partial ^{2}y}{\partial x_{i}\partial x_{j}}}=\sum _{k}\left({\frac {\partial y}{\partial u_{k}}}{\frac {\partial ^{2}u_{k}}{\partial x_{i}\partial x_{j}}}\right)+\sum _{k,\ell }\left({\frac {\partial ^{2}y}{\partial u_{k}\partial u_{\ell }}}{\frac {\partial u_{k}}{\partial x_{i}}}{\frac {\partial u_{\ell }}{\partial x_{j}}}\right).$

Further generalizations

All extensions of calculus have a chain rule. In most of these, the formula remains the same, though the meaning of that formula may be vastly different.

One generalization is to manifolds. In this situation, the chain rule represents the fact that the derivative of f ∘ g is the composite of the derivative of f and the derivative of g. This theorem is an immediate consequence of the higher dimensional chain rule given above, and it has exactly the same formula.

The chain rule is also valid for Fréchet derivatives in Banach spaces. The same formula holds as before. This case and the previous one admit a simultaneous generalization to Banach manifolds.

In abstract algebra, the derivative is interpreted as a morphism of modules of Kähler differentials. A ring homomorphism of commutative rings f : R → S determines a morphism of Kähler differentials Df : Ω_R → Ω_S which sends an element dr to d(f(r)), the exterior differential of f(r). The formula D(f ∘ g) = Df ∘ Dg holds in this context as well.

The common feature of these examples is that they are expressions of the idea that the derivative is part of a functor. A functor is an operation on spaces and functions between them. It associates to each space a new space and to each function between two spaces a new function between the corresponding new spaces. In each of the above cases, the functor sends each space to its tangent bundle and it sends each function to its derivative. For example, in the manifold case, the derivative sends a C^r-manifold to a C^r−1-manifold (its tangent bundle) and a C^r-function to its total derivative. There is one requirement for this to be a functor, namely that the derivative of a composite must be the composite of the derivatives. This is exactly the formula D(f ∘ g) = Df ∘ Dg.

There are also chain rules in stochastic calculus. One of these, Itō's lemma, expresses the composite of an Itō process (or more generally a semimartingale) dX_t with a twice-differentiable function f. In Itō's lemma, the derivative of the composite function depends not only on dX_t and the derivative of f but also on the second derivative of f. The dependence on the second derivative is a consequence of the non-zero quadratic variation of the stochastic process, which broadly speaking means that the process can move up and down in a very rough way. This variant of the chain rule is not an example of a functor because the two functions being composed are of different types.

See also

• Integration by substitution
• Leibniz integral rule
• Quotient rule
• Triple product rule
• Product rule
• Automatic differentiation, a computational method that makes heavy use of the chain rule to compute exact numerical derivatives.

External links

• Hazewinkel, Michiel, ed. (2001) [1994], "Leibniz rule", Encyclopedia of Mathematics, Springer Science+Business Media B.V. / Kluwer Academic Publishers, ISBN 978-1-55608-010-4
• Weisstein, Eric W. "Chain Rule". MathWorld.