Higher-Order Derivatives

Suppose a car's position at time tt is s(t)s(t). Its velocity is s′(t)s'(t). But what tells you whether the car is speeding up or slowing down? You need the derivative of the velocity itself — a "derivative of a derivative." This simple idea, iterating the operation of differentiation, is the subject of this article. It turns out to carry enormous structural weight in analysis. Higher-order derivatives let us classify functions by how much smoothness they possess. Some functions are differentiable only once. Others are differentiable any finite number of times. Some are differentiable infinitely often. This hierarchy is organized into the classes CnC^n and C∞C^\infty. It is not a pedantic bookkeeping device. It is exactly the hypothesis needed to control the error in polynomial approximations (via Taylor's Theorem with Remainder). It also helps detect concavity and convexity of a graph. It lets us distinguish local maxima from minima beyond the first-derivative test. And it helps make sense of physical quantities like acceleration and jerk. The construction is inductive. Define f′′f'' as the derivative of f′f' wherever that derivative exists. Then define f′′′f''' as the derivative of f′′f'', and so on. Each step requires the previous derivative to exist on an open set (or at least a punctured neighborhood). This ensures differentiating it again makes sense. This recursive dependency is the source of every subtlety in the theory. A function can be differentiable without its derivative being differentiable. The resulting gaps in smoothness are exactly what the classes CnC^n are designed to measure. This article gives the precise inductive definition of f(n)f^{(n)}. It proves the basic algebraic facts (linearity, the product rule in its Leibniz form). It defines and studies the smoothness classes Cn(I)C^n(I) and C∞(I)C^\infty(I). It also exhibits sharp counterexamples showing why each hypothesis in these definitions is load-bearing.

Key takeaways
  • The nn-th derivative f(n)f^{(n)} is defined inductively: f(0)=ff^{(0)} = f, and f(n)=(f(n−1))′f^{(n)} = (f^{(n-1)})' wherever the right side exists.
  • Existence of f(n)(a)f^{(n)}(a) requires f(n−1)f^{(n-1)} to exist and be differentiable at aa, hence f(n−1)f^{(n-1)} must exist on a neighborhood of aa.
  • f∈Cn(I)f \in C^n(I) means f(n)f^{(n)} exists on II and is continuous there; f∈C∞(I)f \in C^\infty(I) means derivatives of all orders exist.
  • The Leibniz rule generalizes the product rule: (fg)(n)=∑k=0n(nk)f(k)g(n−k)(fg)^{(n)} = \sum_{k=0}^n \binom{n}{k} f^{(k)} g^{(n-k)}.
  • Smoothness classes are strictly nested: Cn+1(I)⊊Cn(I)C^{n+1}(I) \subsetneq C^n(I), witnessed by explicit functions like xn∣x∣x^n|x|.
  • Key idea: differentiation is an operator you can iterate, and the n-th time you can iterate it and stay continuous measures a function's smoothness.

Where this lives today

Higher-order derivatives are the backbone of Taylor approximation: the coefficients of the Taylor polynomial of ff at a point are literally f(k)(a)/k!f^{(k)}(a)/k!, and the size of the remainder term is controlled by a bound on f(n+1)f^{(n+1)} — the entire subject of Taylor's Theorem with Remainder is a direct continuation of this article. In geometry and optimization, the sign of f′′f'' characterizes Convexity and the Second Derivative of a function, which underlies convex optimization algorithms throughout machine learning and operations research. In physics, f′′,f′′′f'', f''', and beyond appear as acceleration, jerk, snap, and higher kinematic quantities used in robotics motion planning and aerospace trajectory design, where smoothness of a planned path (often demanded to be C2C^2 or higher) prevents mechanically damaging jerks in actuators. In numerical analysis, finite-difference schemes for solving differential equations require explicit error bounds built from higher derivatives (via Taylor expansions), and the stability of these schemes depends on which CnC^n class the true solution belongs to. In differential geometry and PDE theory, the classes CkC^k and C∞C^\infty (and their infinite-dimensional generalizations, Sobolev spaces) are the standard language for stating exactly how much regularity a solution to a differential equation is assumed or proved to have.

Ask the tutor: Why does f∈C∞f \in C^\infty at a point not guarantee that ff equals its Taylor series near that point?

DefinitionHigher-order derivatives and the classes CnC^n, C∞C^\infty

Let I⊆RI \subseteq \mathbb{R} be an open interval and let f:I→Rf: I \to \mathbb{R}.

Inductive definition of f(n)f^{(n)}. Define f(0):=ff^{(0)} := f. Suppose n≥1n \geq 1 and that f(n−1)f^{(n-1)} has been defined on some open set Un−1⊆IU_{n-1} \subseteq I. We say ff is nn times differentiable at a point a∈Un−1a \in U_{n-1} if f(n−1)f^{(n-1)} is differentiable at aa, i.e. if

f(n)(a):=lim⁡h→0f(n−1)(a+h)−f(n−1)(a)hf^{(n)}(a) := \lim_{h \to 0} \frac{f^{(n-1)}(a+h) - f^{(n-1)}(a)}{h}

exists (which in particular requires f(n−1)f^{(n-1)} to be defined on an open interval around aa). If this limit exists for every aa in an open set Un⊆Un−1U_n \subseteq U_{n-1}, we call the resulting function f(n):Un→Rf^{(n)}: U_n \to \mathbb{R} the nn-th derivative of ff. We write f′′:=f(2)f'' := f^{(2)}, f′′′:=f(3)f''' := f^{(3)}, and f(n)f^{(n)} for n≥4n \geq 4.

Smoothness classes. For n≥0n \geq 0, we say f∈Cn(I)f \in C^n(I) (f is of class CnC^n on II) if f(n)f^{(n)} exists at every point of II and f(n):I→Rf^{(n)}: I \to \mathbb{R} is continuous on II. We say f∈C∞(I)f \in C^\infty(I) if f∈Cn(I)f \in C^n(I) for every n≥0n \geq 0; such ff is called smooth. By convention C0(I)C^0(I) is the set of continuous functions on II, consistent with f(0)=ff^{(0)} = f.

Intuition

Think of differentiation as a machine: feed it a function, it returns another function (the derivative), provided the limit defining that derivative exists everywhere you need it. Iterating the machine is exactly what f′′f'', f′′′f''', etc. record. The subtlety is that the machine can jam: if f′f' itself is not differentiable at some point, the process simply stops there, and f′′f'' fails to exist at that point — even though ff and f′f' might both be perfectly well-behaved (continuous, even differentiable) nearby. A good mental picture is a tower of functions f,f′,f′′,f′′′,…f, f', f'', f''', \dots, each one built from the last by the same limiting operation, with the tower's height (how many rungs you can climb before the process breaks) measuring the smoothness of ff. The classes CnC^n slice this tower at height nn and additionally demand that the top function you reach is continuous — not just that the limit exists pointwise, but that it exists and varies continuously. This extra continuity requirement is what makes CnC^n well-behaved enough to support approximation theorems. Geometrically, f′′f'' measures the curvature or bending of the graph of ff: where f′′>0f'' > 0 the graph curves upward (convex), where f′′<0f'' < 0 it curves downward (concave), foreshadowing the link between Second Derivative Test and convexity. Physically, if ff is position, f′f' is velocity, f′′f'' is acceleration, and f′′′f''' is "jerk" — the rate of change of acceleration, which is what makes a smooth elevator ride feel smooth rather than jarring.

Figure
The tower of derivatives: each rung requires the one below to be differentiable
smoothness f (given) f′ f″ f‴, f⁽⁴⁾, … exists if f differentiable exists if f′ differentiable exists if f″ differentiable exists if f‴ differentiable, … C⁰ if f continuous C¹ if f′ continuous C² if f″ continuous C³, …, C∞
TheoremAlgebra of higher derivatives, and the Leibniz formula

Let I⊆RI \subseteq \mathbb{R} be an open interval, let n≥1n \geq 1 be an integer, and suppose f,g:I→Rf, g: I \to \mathbb{R} are both nn times differentiable at every point of II (equivalently, f(n)f^{(n)} and g(n)g^{(n)} exist on all of II).

(a) Linearity. For all α,β∈R\alpha, \beta \in \mathbb{R}, the function αf+βg\alpha f + \beta g is nn times differentiable on II, and

(αf+βg)(n)=αf(n)+βg(n).(\alpha f + \beta g)^{(n)} = \alpha f^{(n)} + \beta g^{(n)}.

(b) Leibniz's rule. The product fgfg is nn times differentiable on II, and

(fg)(n)=∑k=0n(nk)f(k)g(n−k).(fg)^{(n)} = \sum_{k=0}^{n} \binom{n}{k} f^{(k)} g^{(n-k)}.

(c) Class closure. If f,g∈Cn(I)f, g \in C^n(I), then αf+βg∈Cn(I)\alpha f + \beta g \in C^n(I) and fg∈Cn(I)fg \in C^n(I).

Proof

We use induction on nn, building on the Algebraic Rules of Differentiation (linearity and product rule for a single derivative), which we take as known base machinery.

Part (a): Linearity, by induction on nn. Base case n=1n=1. This is exactly the linearity rule for a single derivative: if f,gf, g are differentiable at a∈Ia \in I, then αf+βg\alpha f + \beta g is differentiable at aa with (αf+βg)′(a)=αf′(a)+βg′(a)(\alpha f + \beta g)'(a) = \alpha f'(a) + \beta g'(a). This is a standard consequence of the linearity of limits applied to the difference quotient of αf+βg\alpha f + \beta g. Inductive step. Suppose the claim holds for n−1≥1n - 1 \geq 1: whenever f,gf, g are (n−1)(n-1) times differentiable on II, (αf+βg)(n−1)=αf(n−1)+βg(n−1)(\alpha f + \beta g)^{(n-1)} = \alpha f^{(n-1)} + \beta g^{(n-1)} on II. Now assume f,gf, g are nn times differentiable on II; in particular f(n−1),g(n−1)f^{(n-1)}, g^{(n-1)} exist on II and are themselves differentiable there (since f(n),g(n)f^{(n)}, g^{(n)} exist). By the inductive hypothesis, αf(n−1)+βg(n−1)\alpha f^{(n-1)} + \beta g^{(n-1)} equals (αf+βg)(n−1)(\alpha f + \beta g)^{(n-1)} as a function on II. Since f(n−1)f^{(n-1)} and g(n−1)g^{(n-1)} are each differentiable on II, the base case (n=1n=1 case of linearity, now applied to f(n−1),g(n−1)f^{(n-1)}, g^{(n-1)} in place of f,gf, g) gives that αf(n−1)+βg(n−1)\alpha f^{(n-1)} + \beta g^{(n-1)} is differentiable on II with derivative αf(n)+βg(n)\alpha f^{(n)} + \beta g^{(n)}. Combining,

(αf+βg)(n)=(αf(n−1)+βg(n−1))′=αf(n)+βg(n).(\alpha f + \beta g)^{(n)} = \big(\alpha f^{(n-1)} + \beta g^{(n-1)}\big)' = \alpha f^{(n)} + \beta g^{(n)}.

This completes the induction.

Part (b): Leibniz's rule, by induction on nn. Base case n=1n = 1. This is the ordinary product rule: (fg)′=f′g+fg′(fg)' = f'g + fg', which is exactly ∑k=01(1k)f(k)g(1−k)=(10)fg′+(11)f′g=fg′+f′g\sum_{k=0}^1 \binom{1}{k} f^{(k)} g^{(1-k)} = \binom{1}{0}f g' + \binom{1}{1} f' g = fg' + f'g. This matches the standard product rule from Algebraic Rules of Differentiation. Inductive step. Suppose the formula holds for n−1≥1n-1 \geq 1 whenever both functions are (n−1)(n-1) times differentiable:

(fg)(n−1)=∑k=0n−1(n−1k)f(k)g(n−1−k).(fg)^{(n-1)} = \sum_{k=0}^{n-1} \binom{n-1}{k} f^{(k)} g^{(n-1-k)}.

Now let f,gf, g be nn times differentiable on II. Then each f(k)f^{(k)} for 0≤k≤n−10 \leq k \leq n-1 and each g(n−1−k)g^{(n-1-k)} for the corresponding index is differentiable on II (since f(n),g(n)f^{(n)}, g^{(n)} exist, all lower derivatives exist and are differentiable). Differentiate the right-hand side of the (n−1)(n-1)-formula term by term using linearity (Part (a), one-step case) and the ordinary product rule:

ddx[f(k)g(n−1−k)]=f(k+1)g(n−1−k)+f(k)g(n−k).\frac{d}{dx}\Big[f^{(k)} g^{(n-1-k)}\Big] = f^{(k+1)} g^{(n-1-k)} + f^{(k)} g^{(n-k)}.

Hence

(fg)(n)=∑k=0n−1(n−1k)[f(k+1)g(n−1−k)+f(k)g(n−k)].(fg)^{(n)} = \sum_{k=0}^{n-1} \binom{n-1}{k}\Big[f^{(k+1)} g^{(n-1-k)} + f^{(k)} g^{(n-k)}\Big].

Split this into two sums, S1=∑k=0n−1(n−1k)f(k+1)g(n−1−k)S_1 = \sum_{k=0}^{n-1}\binom{n-1}{k} f^{(k+1)} g^{(n-1-k)} and S2=∑k=0n−1(n−1k)f(k)g(n−k)S_2 = \sum_{k=0}^{n-1}\binom{n-1}{k} f^{(k)} g^{(n-k)}. Re-index S1S_1 by j=k+1j = k+1 (so jj runs from 11 to nn):

S1=∑j=1n(n−1j−1)f(j)g(n−j).S_1 = \sum_{j=1}^{n} \binom{n-1}{j-1} f^{(j)} g^{(n-j)}.

Renaming j→kj \to k in S1S_1 and adding to S2S_2:

(fg)(n)=∑k=1n(n−1k−1)f(k)g(n−k)+∑k=0n−1(n−1k)f(k)g(n−k).(fg)^{(n)} = \sum_{k=1}^n \binom{n-1}{k-1} f^{(k)}g^{(n-k)} + \sum_{k=0}^{n-1}\binom{n-1}{k} f^{(k)}g^{(n-k)}.

The k=0k=0 term appears only in the second sum (coefficient (n−10)=1=(n0)\binom{n-1}{0}=1=\binom{n}{0}), and the k=nk=n term appears only in the first sum (coefficient (n−1n−1)=1=(nn)\binom{n-1}{n-1} = 1 = \binom{n}{n}). For 1≤k≤n−11 \leq k \leq n-1, both sums contribute, and Pascal's identity gives

(n−1k−1)+(n−1k)=(nk).\binom{n-1}{k-1} + \binom{n-1}{k} = \binom{n}{k}.

Therefore

(fg)(n)=∑k=0n(nk)f(k)g(n−k),(fg)^{(n)} = \sum_{k=0}^n \binom{n}{k} f^{(k)} g^{(n-k)},

which is exactly the claimed formula. This completes the induction, using Pascal's identity (a purely combinatorial fact, proved directly from the definition of binomial coefficients, hence not circular).

Part (c): Class closure. Suppose f,g∈Cn(I)f, g \in C^n(I): so f(n),g(n)f^{(n)}, g^{(n)} exist on II and are continuous. By Part (a), (αf+βg)(n)=αf(n)+βg(n)(\alpha f + \beta g)^{(n)} = \alpha f^{(n)} + \beta g^{(n)} exists on II, and it is continuous since it is a linear combination of continuous functions (sum and scalar multiple of continuous functions are continuous, a basic fact about limits). Hence αf+βg∈Cn(I)\alpha f + \beta g \in C^n(I). For the product: by Part (b), (fg)(n)=∑k=0n(nk)f(k)g(n−k)(fg)^{(n)} = \sum_{k=0}^n \binom{n}{k} f^{(k)} g^{(n-k)}. Since f,g∈Cn(I)f, g \in C^n(I), in particular f(k)f^{(k)} exists on II for all 0≤k≤n0 \le k \le n (differentiability of order nn implies existence of all lower-order derivatives, immediate from the inductive definition — each f(k)f^{(k)} for k<nk<n must exist and be differentiable for f(n)f^{(n)} to exist at all). Moreover each f(k)f^{(k)} for k<nk < n is differentiable, hence continuous (differentiability implies continuity, a standard one-variable fact), and f(n)f^{(n)} is continuous by hypothesis. The same holds for gg. So each product f(k)g(n−k)f^{(k)}g^{(n-k)} is a product of continuous functions, hence continuous, and (fg)(n)(fg)^{(n)} is a finite sum of continuous functions, hence continuous. Thus fg∈Cn(I)fg \in C^n(I). ■\blacksquare

∎
CorollaryStrict nesting of smoothness classes

For every n≥0n \geq 0, Cn+1(I)⊆Cn(I)C^{n+1}(I) \subseteq C^n(I), and this inclusion is strict: there exists f∈Cn(I)∖Cn+1(I)f \in C^n(I) \setminus C^{n+1}(I).

We prove the inclusion here; the strictness is exhibited explicitly in the Worked Examples and the Counterexample section below.

Proof

Let f∈Cn+1(I)f \in C^{n+1}(I), so f(n+1)f^{(n+1)} exists on II and is continuous. In particular f(n+1)=(f(n))′f^{(n+1)} = (f^{(n)})' exists on II, so f(n)f^{(n)} is differentiable on II, and differentiability implies continuity (standard one-variable fact: if gg is differentiable at aa, then g(a+h)−g(a)=g(a+h)−g(a)h⋅h→g′(a)⋅0=0g(a+h) - g(a) = \frac{g(a+h)-g(a)}{h}\cdot h \to g'(a)\cdot 0 = 0 as h→0h \to 0, so gg is continuous at aa). Hence f(n)f^{(n)} is continuous on II, which is exactly the statement f∈Cn(I)f \in C^n(I). ■\blacksquare

∎

Worked examples

Example 1: A function that is exactly C2C^2 but not C3C^3, computed explicitly. Let f:R→Rf: \mathbb{R} \to \mathbb{R} be defined by

f(x)=x3∣x∣={x4x≥0−x4x<0.f(x) = x^3 |x| = \begin{cases} x^4 & x \geq 0 \\ -x^4 & x < 0. \end{cases}

We compute the derivatives directly and locate exactly where smoothness fails. Computing f′f'. For x>0x > 0, f(x)=x4f(x) = x^4, so f′(x)=4x3f'(x) = 4x^3. For x<0x < 0, f(x)=−x4f(x) = -x^4, so f′(x)=−4x3f'(x) = -4x^3. At x=0x=0, using the definition of the derivative,

f′(0)=lim⁡h→0f(h)−f(0)h=lim⁡h→0h3∣h∣h=lim⁡h→0h2∣h∣=0,f'(0) = \lim_{h \to 0} \frac{f(h) - f(0)}{h} = \lim_{h\to 0} \frac{h^3|h|}{h} = \lim_{h \to 0} h^2 |h| = 0,

since ∣h2∣h∣∣≤h2→0|h^2|h|| \le h^2 \to 0 (squeeze). So f′(x)=4∣x∣x2=4x2∣x∣f'(x) = 4|x|x^2 = 4x^2|x| for all xx (check: at x>0x>0 this is 4x34x^3 ✓, at x<0x<0 this is −4x3-4x^3 ✓, at 00 it's 00 ✓). Note f′f' is continuous everywhere (it is 4x2∣x∣4x^2|x|, a product/composition of continuous functions), so f∈C1(R)f \in C^1(\mathbb{R}). Computing f′′f''. For x>0x>0: f′(x)=4x3⇒f′′(x)=12x2f'(x) = 4x^3 \Rightarrow f''(x) = 12x^2. For x<0x<0: f′(x)=−4x3⇒f′′(x)=−12x2f'(x) = -4x^3 \Rightarrow f''(x) = -12x^2. At x=0x=0:

f′′(0)=lim⁡h→0f′(h)−f′(0)h=lim⁡h→04h2∣h∣h=lim⁡h→04h∣h∣=0.f''(0) = \lim_{h\to 0} \frac{f'(h) - f'(0)}{h} = \lim_{h \to 0} \frac{4h^2|h|}{h} = \lim_{h\to 0} 4h|h| = 0.

So f′′(x)=12x∣x∣f''(x) = 12x|x| for all xx (check at x>0x>0: 12x212x^2 ✓; at x<0x<0: −12x2-12x^2 ✓). This is continuous everywhere (again a product of continuous functions), so f∈C2(R)f \in C^2(\mathbb{R}). Computing f′′′f'''. For x>0x > 0: f′′(x)=12x2⇒f′′′(x)=24xf''(x) = 12x^2 \Rightarrow f'''(x) = 24x. For x<0x<0: f′′(x)=−12x2⇒f′′′(x)=−24xf''(x) = -12x^2 \Rightarrow f'''(x) = -24x. At x=0x = 0:

f′′′(0)=lim⁡h→0f′′(h)−f′′(0)h=lim⁡h→012h∣h∣h=lim⁡h→012∣h∣=0.f'''(0) = \lim_{h \to 0} \frac{f''(h) - f''(0)}{h} = \lim_{h\to0} \frac{12h|h|}{h} = \lim_{h \to 0} 12|h| = 0.

So f′′′f''' does exist everywhere, and equals f′′′(x)=24∣x∣f'''(x) = 24|x| (check: x>0x>0 gives 24x24x ✓, x<0x<0 gives −24x=24∣x∣-24x = 24|x| ✓, and at 00 gives 00 ✓). This function 24∣x∣24|x| is continuous, so in fact f∈C3(R)f \in C^3(\mathbb{R}) too! Computing f(4)f^{(4)}. Now differentiate f′′′(x)=24∣x∣f'''(x) = 24|x| at x=0x=0:

lim⁡h→0+24∣h∣−0h=lim⁡h→0+24hh=24,lim⁡h→0−24∣h∣−0h=lim⁡h→0−−24hh=−24.\lim_{h \to 0^+} \frac{24|h| - 0}{h} = \lim_{h\to 0^+} \frac{24h}{h} = 24, \qquad \lim_{h \to 0^-} \frac{24|h|-0}{h} = \lim_{h \to 0^-}\frac{-24h}{h} = -24.

The one-sided limits disagree (24≠−2424 \neq -24), so f(4)(0)f^{(4)}(0) does not exist. Hence f∉C4(R)f \notin C^4(\mathbb{R}) — in fact f(4)f^{(4)} fails to exist at even a single point, so ff is not 44 times differentiable on any open interval containing 00.

Conclusion: f(x)=x3∣x∣f(x) = x^3|x| satisfies f∈C3(R)f \in C^3(\mathbb{R}) but f∉C4(R)f \notin C^4(\mathbb{R}) (indeed f(4)(0)f^{(4)}(0) does not exist at all). This shows the chain C0⊋C1⊋C2⊋C3⊋C4C^0 \supsetneq C^1 \supsetneq C^2 \supsetneq C^3 \supsetneq C^4 is strict at this level, generalizing the classic x∣x∣x|x| example (which is the n=1n=1 analogue) to a higher order.

Example 2: Leibniz's rule applied to an asymmetric product, n=4n=4. Let f(x)=e2xf(x) = e^{2x} and g(x)=x3−xg(x) = x^3 - x on I=RI = \mathbb{R}; we compute (fg)(4)(fg)^{(4)} two ways to cross-check Leibniz's formula. Direct route. h(x)=e2x(x3−x)h(x) = e^{2x}(x^3 - x). We could expand and differentiate four times directly, but let's instead verify via Leibniz and then spot-check numerically at x=1x=1. Derivatives of ff. f(x)=e2xf(x) = e^{2x}, so f(k)(x)=2ke2xf^{(k)}(x) = 2^k e^{2x} for all k≥0k \geq 0. Thus f(1)=e2f(1) = e^2, f′(1)=2e2f'(1) = 2e^2, f′′(1)=4e2f''(1) = 4e^2, f′′′(1)=8e2f'''(1) = 8e^2, f(4)(1)=16e2f^{(4)}(1) = 16e^2. Derivatives of gg. g(x)=x3−xg(x) = x^3 - x: g′(x)=3x2−1g'(x) = 3x^2 - 1, g′′(x)=6xg''(x) = 6x, g′′′(x)=6g'''(x) = 6, g(4)(x)=0g^{(4)}(x) = 0. At x=1x=1: g(1)=0g(1) = 0, g′(1)=2g'(1) = 2, g′′(1)=6g''(1) = 6, g′′′(1)=6g'''(1) = 6, g(4)(1)=0g^{(4)}(1) = 0. Applying Leibniz's rule with n=4n=4:

(fg)(4)(1)=∑k=04(4k)f(k)(1) g(4−k)(1).(fg)^{(4)}(1) = \sum_{k=0}^{4} \binom{4}{k} f^{(k)}(1)\, g^{(4-k)}(1).

Compute each term:

  • k=0k=0: (40)f(1)g(4)(1)=1⋅e2⋅0=0\binom{4}{0} f(1) g^{(4)}(1) = 1 \cdot e^2 \cdot 0 = 0.
  • k=1k=1: (41)f′(1)g(3)(1)=4⋅2e2⋅6=48e2\binom{4}{1} f'(1) g^{(3)}(1) = 4 \cdot 2e^2 \cdot 6 = 48e^2.
  • k=2k=2: (42)f′′(1)g′′(1)=6⋅4e2⋅6=144e2\binom{4}{2} f''(1) g''(1) = 6 \cdot 4e^2 \cdot 6 = 144 e^2.
  • k=3k=3: (43)f′′′(1)g′(1)=4⋅8e2⋅2=64e2\binom{4}{3} f'''(1) g'(1) = 4 \cdot 8 e^2 \cdot 2 = 64e^2.
  • k=4k=4: (44)f(4)(1)g(1)=1⋅16e2⋅0=0\binom{4}{4} f^{(4)}(1) g(1) = 1 \cdot 16 e^2 \cdot 0 = 0.

Summing: (fg)(4)(1)=0+48e2+144e2+64e2+0=256e2(fg)^{(4)}(1) = 0 + 48e^2 + 144e^2 + 64e^2 + 0 = 256 e^2. Cross-check by direct expansion. Write h(x)=e2xg(x)h(x) = e^{2x} g(x) with g=x3−xg = x^3 - x. One can show by induction (or repeated product rule) that

h(4)(x)=e2x∑k=04(4k)2kg(4−k)(x)=e2x[g(4)+8g(3)+24g′′+32g′+16g](x),h^{(4)}(x) = e^{2x}\sum_{k=0}^4 \binom{4}{k} 2^k g^{(4-k)}(x) = e^{2x}\Big[g^{(4)} + 8g^{(3)} + 24 g'' + 32 g' + 16 g\Big](x),

matching the Leibniz coefficients (4k)2k\binom{4}{k}2^k for k=0,…,4k=0,\dots,4: namely 16,32,24,8,116,32,24,8,1 paired with g,g′,g′′,g(3),g(4)g, g', g'', g^{(3)}, g^{(4)} respectively (same computation, reorganized). At x=1x=1: g(1)=0,g′(1)=2,g′′(1)=6,g′′′(1)=6,g(4)(1)=0g(1){=}0, g'(1){=}2, g''(1){=}6, g'''(1){=}6, g^{(4)}(1){=}0, giving

h(4)(1)=e2[0+8(6)+24(6)+32(2)+16(0)]=e2[48+144+64]=256e2,h^{(4)}(1) = e^2[0 + 8(6) + 24(6) + 32(2) + 16(0)] = e^2[48 + 144 + 64] = 256e^2,

matching exactly. This numerically confirms Leibniz's rule in a genuinely asymmetric case (different function types, non-palindromic coefficients all engaged).

CounterexampleWhy "f′f' exists on II" is not enough — you need f′f' differentiable, not just ff twice "formally" iterable

The definition of f(n)(a)f^{(n)}(a) requires f(n−1)f^{(n-1)} to exist not just at aa but on a full neighborhood of aa, and then to be differentiable at aa. One might be tempted to weaken this to: "f(n−1)(a)f^{(n-1)}(a) exists, and some difference quotient involving only values of f(n−2)f^{(n-2)} near aa converges" — hoping to skip verifying that f(n−1)f^{(n-1)} is genuinely differentiable. Here is why that shortcut is unsound, via a classical example due to the structure of oscillating functions. Define f:R→Rf: \mathbb{R} \to \mathbb{R} by

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

Step 1: ff is differentiable everywhere, but f′f' is discontinuous at 00. For x≠0x \neq 0, the product and chain rules give

f′(x)=2xsin⁡(1/x)−cos⁡(1/x).f'(x) = 2x\sin(1/x) - \cos(1/x).

At x=0x=0: f′(0)=lim⁡h→0h2sin⁡(1/h)−0h=lim⁡h→0hsin⁡(1/h)=0f'(0) = \lim_{h\to0} \frac{h^2\sin(1/h) - 0}{h} = \lim_{h\to 0} h\sin(1/h) = 0 (squeeze, since ∣hsin⁡(1/h)∣≤∣h∣|h\sin(1/h)| \le |h|). So f′(0)=0f'(0) = 0 exists, and ff is differentiable on all of R\mathbb{R} — this is a genuine, everywhere-defined derivative f′f'. But as x→0x \to 0, f′(x)=2xsin⁡(1/x)−cos⁡(1/x)f'(x) = 2x\sin(1/x) - \cos(1/x) does not converge, because cos⁡(1/x)\cos(1/x) oscillates between −1-1 and 11 without settling (e.g. along xn=1/(2πn)→0x_n = 1/(2\pi n) \to 0, cos⁡(1/xn)=1\cos(1/x_n) = 1, while along xn=1/(π+2πn)→0x_n = 1/(\pi + 2\pi n) \to 0, cos⁡(1/xn)=−1\cos(1/x_n) = -1). So lim⁡x→0f′(x)\lim_{x\to 0} f'(x) does not exist, hence f′f' is discontinuous at 00, even though f′(0)f'(0) itself is a perfectly well-defined real number.

Step 2: Consequently f′′(0)f''(0) does not exist, because differentiability at 00 of f′f' would force f′f' to be continuous at 00 (differentiability implies continuity, as used in the Corollary's proof) — but we just showed f′f' is discontinuous at 00. So the inductive chain to define f′′(0)f''(0) breaks down at the very first step: f′f' exists everywhere (so "iterating the operation" seems to make sense symbolically) but is not differentiable at 00, hence f′′(0)f''(0) simply does not exist.

What this shows. If we dropped the requirement in the definition that f(n−1)f^{(n-1)} be genuinely differentiable at aa (as opposed to merely defined near aa), and instead tried some naive symbol-pushing — e.g., differentiating the formula "2xsin⁡(1/x)−cos⁡(1/x)2x\sin(1/x) - \cos(1/x)" term by term while ignorning that this formula is not valid at x=0x=0 and that the resulting expression's limit at 00 must be checked — we would get nonsensical or ill-defined results. The hypothesis "f(n−1)f^{(n-1)} differentiable at aa" in the definition of f(n)(a)f^{(n)}(a) is not a formality: dropping it (settling for "the formula for f(n−1)f^{(n-1)} makes sense off a point") produces exactly this kind of failure, where f∈C0(R)f \in C^0(\mathbb{R}), ff is differentiable everywhere (so we might loosely say f∈D1(R)f \in \mathcal{D}^1(\mathbb{R})), but f∉C1(R)f \notin C^1(\mathbb{R}) and f′′(0)f''(0) does not exist at all. This is precisely why the classes CnC^n (continuity of the derivative), not mere "nn-fold differentiability," are the robust objects used in analysis, e.g. as hypotheses for Taylor's Theorem with Remainder.

Scope, edge cases & common traps

SituationWhat holdsWhat can fail
f(n)(a)f^{(n)}(a) exists at a single point aaf(n−1)f^{(n-1)} must exist on a neighborhood of aa and be differentiable at aaf(n)f^{(n)} need not exist at nearby points, and need not be continuous even at aa
f∈Cn(I)f \in C^n(I)f(k)f^{(k)} exists and is continuous on II for every 0≤k≤n0 \le k \le nf(n+1)f^{(n+1)} may fail to exist anywhere
f,g∈Cn(I)f, g \in C^n(I)αf+βg, fg∈Cn(I)\alpha f + \beta g,\ fg \in C^n(I) (Theorem, parts a,c)f/g∈Cnf/g \in C^n only where g≠0g \neq 0; must exclude zeros of gg
f∈C∞(I)f \in C^\infty(I)All derivatives exist and are continuousDoes not imply ff is analytic (equal to its Taylor series) — see FAQ
Composition f∘gf\circ gIf f,g∈Cnf,g \in C^n, then f∘g∈Cnf \circ g \in C^n (provable by induction using the chain rule, stated here without full proof)The formula for (f∘g)(n)(f\circ g)^{(n)} (Faà di Bruno's formula) is far more complex than Leibniz's rule and is easy to misstate

The mathematician behind it

Origins · 1 minGottfried Wilhelm Leibniz and how this was discovered
GW
Gottfried Wilhelm Leibniz
1646–1716 · Co-inventor of calculus, Hanover and Paris

Leibniz developed the differential calculus independently of Newton in the 1670s–80s, and it is his notation — dy/dxdy/dx, ∫\int, and the systematic use of dny/dxnd^n y / dx^n for higher derivatives — that survives in nearly all of modern analysis. His 1695 correspondence with Johann Bernoulli contains the first explicit statement of what we now call Leibniz's rule for the nn-th derivative of a product, discovered while investigating fractional (non-integer order) derivatives, an idea decades ahead of its formal development. Leibniz's insistence on symbolic manipulation and rigorous notation, rather than purely geometric argument, is precisely what makes the inductive machinery of f(n)f^{(n)} tractable and teachable today.

History
- 1675–1676: Leibniz introduces differential notation dxdx, dydy and begins systematic symbolic calculus.
- 1695: Leibniz states the general product-rule formula for (fg)(n)(fg)^{(n)} in a letter to Johann Bernoulli, prompted by a question about fractional derivatives.
- 1797: Lagrange, in Théorie des fonctions analytiques, popularizes the prime notation f′,f′′,f′′′f', f'', f''' for successive derivatives.
- 1821–1823: Cauchy, in his Cours d'Analyse and related lectures, gives the first fully rigorous ε\varepsilon–δ\delta treatment of limits underlying the modern definition of f(n)f^{(n)}.
- 1872: Weierstrass exhibits a function continuous everywhere but differentiable nowhere, sharpening the community's understanding of exactly how fragile smoothness properties are — motivating the careful class hierarchy C0⊋C1⊋C2⊋⋯C^0 \supsetneq C^1 \supsetneq C^2 \supsetneq \cdots used today.
What is the difference between "ff is nn times differentiable" and "f∈Cnf \in C^n"?

"ff is nn times differentiable at a point (or on a set)" only asserts that the limit defining f(n)f^{(n)} exists there — nothing is said about whether f(n)f^{(n)} varies continuously. "f∈Cn(I)f \in C^n(I)" additionally requires f(n)f^{(n)} to be continuous on II. These are genuinely different: the counterexample f(x)=x2sin⁡(1/x)f(x)=x^2\sin(1/x) (extended by f(0)=0f(0)=0) is differentiable everywhere (so "1 time differentiable" everywhere), yet f′f' is discontinuous at 00, so f∉C1(R)f \notin C^1(\mathbb{R}). Every CnC^n function is nn times differentiable, but not conversely. This distinction matters because theorems like Taylor's theorem with the Lagrange remainder typically require the stronger CnC^n (or Cn+1C^{n+1}) hypothesis, not merely "n times differentiable," to control error terms rigorously.

Does f∈C∞f \in C^\infty mean ff equals its own Taylor series?

No — this is one of the most important subtleties in analysis. A function all of whose derivatives exist and are continuous (smooth) need not be analytic (equal to the sum of its Taylor series in a neighborhood of each point). The standard counterexample is f(x)=e−1/x2f(x) = e^{-1/x^2} for x≠0x \ne 0, f(0)=0f(0)=0: one can show by induction that every derivative f(n)(0)=0f^{(n)}(0) = 0, yet ff is not identically zero near 00. So its Taylor series at 00 is the zero function, which does not equal ff anywhere except at x=0x=0 itself. Thus C∞C^\infty is a strictly larger class than the analytic functions CωC^\omega; this gap is exactly what Taylor's Theorem with Remainder is designed to quantify and control.

Why does the definition require f(n−1)f^{(n-1)} to exist on a whole neighborhood, not just at the point aa?

Differentiability at a point is inherently a two-sided limiting statement: f(n−1)f^{(n-1)} must be defined for a+ha+h for all sufficiently small hh (both positive and negative, unless working on a one-sided endpoint) in order for the difference quotient f(n−1)(a+h)−f(n−1)(a)h\frac{f^{(n-1)}(a+h)-f^{(n-1)}(a)}{h} to even make sense as h→0h \to 0. If f(n−1)f^{(n-1)} were defined only at aa itself, the difference quotient could not be formed for h≠0h \ne 0, and no derivative could ever be extracted. This is not a technicality specific to higher derivatives — it is inherited directly from the ordinary definition of a single derivative, which already requires the function to be defined near the point in question.

Is every CnC^n function's nn-th derivative automatically differentiable?

No, and this is exactly the gap between CnC^n and Cn+1C^{n+1}. Being in Cn(I)C^n(I) only guarantees f(n)f^{(n)} is continuous, not differentiable. The function f(x)=xn+1∣x∣f(x) = x^{n+1}|x| (generalizing Worked Example 1's x3∣x∣x^3|x|, which is the n=2n=2 case) is a clean example: it lies in Cn(R)C^n(\mathbb{R}) but its nn-th derivative fails to be differentiable at 00, so f∉Cn+1(R)f \notin C^{n+1}(\mathbb{R}). This is precisely how one proves the class inclusions Cn+1(I)⊊Cn(I)C^{n+1}(I) \subsetneq C^n(I) are strict for every nn, not just an isolated low-order phenomenon.

How does Leibniz's rule relate to the binomial theorem?

The formula (fg)(n)=∑k=0n(nk)f(k)g(n−k)(fg)^{(n)} = \sum_{k=0}^n \binom{n}{k} f^{(k)} g^{(n-k)} has exactly the same combinatorial structure as the binomial theorem (x+y)n=∑k=0n(nk)xkyn−k(x+y)^n = \sum_{k=0}^n \binom{n}{k}x^k y^{n-k}, and the proofs are structurally parallel: both proceed by induction using Pascal's identity (n−1k−1)+(n−1k)=(nk)\binom{n-1}{k-1}+\binom{n-1}{k}=\binom{n}{k} to merge two sums into one with binomial coefficients. Formally, one can think of the differentiation operator DD as behaving like "D(fg)=(Df)g+f(Dg)D(fg) = (Df)g + f(Dg)," which is the same recursive structure that generates binomial coefficients when iterated nn times — this is why the same numbers (nk)\binom{n}{k} appear in both settings.

Do higher-order derivatives make sense for functions of several variables?

Yes, but the theory becomes richer: for f:Rm→Rf: \mathbb{R}^m \to \mathbb{R}, one first defines partial derivatives ∂f/∂xi\partial f/\partial x_i, then second-order partials ∂2f/∂xi∂xj\partial^2 f/\partial x_i \partial x_j, and the analogue of CnC^n smoothness additionally requires all mixed partials up to order nn to exist and be continuous. A celebrated theorem (Clairaut's/Schwarz's theorem) states that under continuity hypotheses, mixed partials commute: ∂2f/∂xi∂xj=∂2f/∂xj∂xi\partial^2 f/\partial x_i\partial x_j = \partial^2 f/\partial x_j \partial x_i. This multivariable extension is essential for defining the Hessian matrix, which generalizes the single-variable second derivative test used to classify critical points.

Why does the problem statement distinguish CnC^n from merely writing f(n)f^{(n)} "exists everywhere"?

Because "exists everywhere" is a strictly weaker, less useful hypothesis for proving theorems. Many key results in analysis — including error bounds in Taylor's theorem, and criteria for convexity — require not just that f(n)f^{(n)} exists pointwise but that it does not jump around unpredictably, which is exactly what continuity guarantees. Working with CnC^n instead of the weaker pointwise-existence notion lets you invoke the Extreme Value Theorem and Intermediate Value Theorem on f(n)f^{(n)} itself (since those theorems require continuity), which is essential machinery for bounding remainder terms and locating critical points of higher order.

Before this

Where this leads