Partial Derivatives

⚡ Key Takeaways

One-line statement: A partial derivative measures the instantaneous rate of change of a multivariable function with respect to one variable while all other variables are held fixed. Key result in words: For f(x,y)f(x,y), the partial ∂f/∂x\partial f/\partial x is the ordinary derivative obtained by treating yy as a constant. Where it applies: Multivariable calculus, vector calculus, optimization, thermodynamics, PDEs, economics, machine learning (gradients). Who/when: Notation and systematic theory developed in the 18th–19th centuries; the curly ∂\partial symbol popularized by Legendre (1786) and Jacobi (1841).

💡 The Core Insight

Freeze everything but one variable, then differentiate normally. A partial derivative reduces a hard multivariable problem to a familiar single-variable one by slicing the domain along a coordinate axis. Geometrically, you intersect the surface z=f(x,y)z=f(x,y) with a plane parallel to a coordinate axis and measure the slope of the resulting curve. Every rule from single-variable calculus (product, quotient, chain) applies verbatim once the other variables are treated as constants.

Formal Statement

DefinitionPartial Derivative

Let f:U⊆Rn→Rf:U\subseteq\mathbb{R}^n\to\mathbb{R} be defined on an open set UU, and let a=(a1,…,an)∈U\mathbf{a}=(a_1,\dots,a_n)\in U. The partial derivative of ff with respect to xix_i at a\mathbf{a} is

1
∂f∂xi(a)=lim⁡h→0f(a1,…,ai+h,…,an)−f(a1,…,ai,…,an)h\quad \frac{\partial f}{\partial x_i}(\mathbf{a}) = \lim_{h\to 0}\frac{f(a_1,\dots,a_i+h,\dots,a_n)-f(a_1,\dots,a_i,\dots,a_n)}{h}

Derivation hint: This is the ordinary difference quotient in the single direction ei\mathbf{e}_i, holding all other coordinates fixed at their values aja_j (j≠ij\neq i). Here every symbol is: ff is the scalar-valued function; a\mathbf{a} is the base point; xix_i is the differentiation variable; hh is the increment; ei\mathbf{e}_i is the ii-th standard basis vector. Common notations include ∂f/∂xi\partial f/\partial x_i, fxif_{x_i}, ∂xif\partial_{x_i} f, and DifD_i f.

DefinitionEquivalent Directional Form

The partial derivative is the directional derivative along a coordinate axis:

2
∂f∂xi(a)=lim⁡h→0f(a+h ei)−f(a)h\quad \frac{\partial f}{\partial x_i}(\mathbf{a}) = \lim_{h\to 0}\frac{f(\mathbf{a}+h\,\mathbf{e}_i)-f(\mathbf{a})}{h}

Derivation hint: Replace the ii-th coordinate increment in [1] with the vector form a+hei\mathbf{a}+h\mathbf{e}_i, which changes only the ii-th component.

TheoremClairaut–Schwarz Theorem (Equality of Mixed Partials)

If f:U→Rf:U\to\mathbb{R} has continuous second-order partial derivatives ∂2f/∂x ∂y\partial^2 f/\partial x\,\partial y and ∂2f/∂y ∂x\partial^2 f/\partial y\,\partial x in a neighborhood of a\mathbf{a}, then they are equal:

3
∂2f∂x ∂y(a)=∂2f∂y ∂x(a)\quad \frac{\partial^2 f}{\partial x\,\partial y}(\mathbf{a}) = \frac{\partial^2 f}{\partial y\,\partial x}(\mathbf{a})

Derivation hint: Apply the Mean Value Theorem twice to the second difference Δ=f(a+h,b+k)−f(a+h,b)−f(a,b+k)+f(a,b)\Delta = f(a+h,b+k)-f(a+h,b)-f(a,b+k)+f(a,b), expressing it through both orders of differentiation, then take h,k→0h,k\to0 using continuity.

TheoremDifferentiability Implies Partials Exist

If ff is (totally) differentiable at a\mathbf{a} — meaning there exists a linear map LL with f(a+h)=f(a)+Lh+o(∥h∥)f(\mathbf{a}+\mathbf{h})=f(\mathbf{a})+L\mathbf{h}+o(\|\mathbf{h}\|) — then all partial derivatives exist at a\mathbf{a}, and the total differential is

4
df=∑i=1n∂f∂xi dxi\quad df = \sum_{i=1}^{n}\frac{\partial f}{\partial x_i}\,dx_i

Derivation hint: Evaluate the differentiability condition along h=hei\mathbf{h}=h\mathbf{e}_i; the linear map's action on ei\mathbf{e}_i recovers exactly ∂f/∂xi\partial f/\partial x_i. The converse is FALSE — see the counterexample below.

Why It Holds / Derivation

ProofClairaut's Theorem via the Mean Value Theorem

We prove [3] for f(x,y)f(x,y) near (a,b)(a,b).

Step 1. Define the second difference

Δ(h,k)=f(a+h,b+k)−f(a+h,b)−f(a,b+k)+f(a,b).\Delta(h,k)=f(a+h,b+k)-f(a+h,b)-f(a,b+k)+f(a,b).

Step 2. Set φ(x)=f(x,b+k)−f(x,b)\varphi(x)=f(x,b+k)-f(x,b), so Δ=φ(a+h)−φ(a)\Delta=\varphi(a+h)-\varphi(a). By the Mean Value Theorem there is ξ\xi between aa and a+ha+h with

Δ=h φ′(ξ)=h[fx(ξ,b+k)−fx(ξ,b)].\Delta=h\,\varphi'(\xi)=h\big[f_x(\xi,b+k)-f_x(\xi,b)\big].

Step 3. Apply the MVT again to y↦fx(ξ,y)y\mapsto f_x(\xi,y) on [b,b+k][b,b+k]: there is η\eta between bb and b+kb+k with

Δ=h k fyx(ξ,η).\Delta = h\,k\,f_{yx}(\xi,\eta).

Step 4. By the symmetric argument (subtracting in the other order, defining ψ(y)=f(a+h,y)−f(a,y)\psi(y)=f(a+h,y)-f(a,y)), we obtain Δ=h k fxy(ξ′,η′)\Delta = h\,k\,f_{xy}(\xi',\eta') for some (ξ′,η′)(\xi',\eta').

Step 5. Equate: fyx(ξ,η)=fxy(ξ′,η′)f_{yx}(\xi,\eta)=f_{xy}(\xi',\eta'). As h,k→0h,k\to0, all four sample points tend to (a,b)(a,b). Continuity of both mixed partials forces fyx(a,b)=fxy(a,b)f_{yx}(a,b)=f_{xy}(a,b). ■\blacksquare

∎
mathematician · Carl Gustav Jacob Jacobi | 1804–1851

Jacobi systematized the use of partial derivatives in his work on determinants of derivative arrays — what we now call the Jacobian matrix and Jacobian determinant. His 1841 treatise De determinantibus functionalibus fixed the modern ∂\partial notation in the mathematical mainstream and demonstrated the power of organizing partial derivatives into matrices for studying coordinate transformations, implicit functions, and systems of equations. The chain rule for compositions of multivariable maps is most naturally stated through Jacobian matrix multiplication.

Building Intuition

🧠 Intuition

Imagine a mountainous terrain whose height is z=f(x,y)z=f(x,y), with xx pointing east and yy pointing north. Standing at a point, ∂f/∂x\partial f/\partial x is the slope you feel if you walk due east — you ignore any north–south change because yy is held fixed. Likewise ∂f/∂y\partial f/\partial y is the slope walking due north. Each partial is a one-dimensional slice: cut the surface with a vertical plane aligned to an axis, and you get an ordinary curve whose slope is the partial derivative. The full gradient ∇f=(∂f/∂x,∂f/∂y)\nabla f=(\partial f/\partial x,\partial f/\partial y) assembles these slices into the single direction of steepest ascent.

Figure
Slicing a surface: the partial derivative ∂f/∂x is the slope of the curve cut by the plane y = b
x z y plane y = b z = f(x, b) slope = ∂f/∂x a

Worked Examples

ExamplePolynomial Partials with Numbers

Given: f(x,y)=3x2y+4xy3−5y+7f(x,y)=3x^2y + 4xy^3 - 5y + 7. Find ∂f/∂x\partial f/\partial x and ∂f/∂y\partial f/\partial y at the point (2,1)(2,1).

Step 1 — Differentiate w.r.t. xx (treat yy as constant):

∂f∂x=6xy+4y3.\frac{\partial f}{\partial x}=6xy + 4y^3.

Step 2 — Evaluate at (2,1)(2,1):

∂f∂x(2,1)=6(2)(1)+4(1)3=12+4=16.\frac{\partial f}{\partial x}(2,1)=6(2)(1)+4(1)^3 = 12+4 = 16.

Step 3 — Differentiate w.r.t. yy (treat xx as constant):

∂f∂y=3x2+12xy2−5.\frac{\partial f}{\partial y}=3x^2 + 12xy^2 - 5.

Step 4 — Evaluate at (2,1)(2,1):

∂f∂y(2,1)=3(4)+12(2)(1)−5=12+24−5=31.\frac{\partial f}{\partial y}(2,1)=3(4)+12(2)(1)-5 = 12+24-5 = 31.

Boxed answer: fx(2,1)=16,fy(2,1)=31\boxed{f_x(2,1)=16,\quad f_y(2,1)=31} (dimensionless slopes; if x,yx,y carry units, slopes carry units of ff per unit of the respective variable).

ExampleIdeal Gas — Partial Derivatives in Thermodynamics

Given: The ideal gas law P=nRTVP=\dfrac{nRT}{V}, with n=2 moln=2\ \text{mol}, R=8.314 J ⋅ mol−1K−1R=8.314\ \text{J}\,\cdot\,\text{mol}^{-1}\text{K}^{-1}, T=300 KT=300\ \text{K}, V=0.05 m3V=0.05\ \text{m}^3. Find how pressure changes with temperature (at fixed VV) and with volume (at fixed TT).

Step 1 — Partial w.r.t. TT at constant VV:

(∂P∂T)V=nRV=(2)(8.314)0.05=332.6 Pa ⋅ K−1.\left(\frac{\partial P}{\partial T}\right)_V=\frac{nR}{V}=\frac{(2)(8.314)}{0.05}=332.6\ \text{Pa}\,\cdot\,\text{K}^{-1}.

Step 2 — Partial w.r.t. VV at constant TT:

(∂P∂V)T=−nRTV2=−(2)(8.314)(300)(0.05)2=−1.996×106 Pa ⋅ m−3.\left(\frac{\partial P}{\partial V}\right)_T=-\frac{nRT}{V^2}=-\frac{(2)(8.314)(300)}{(0.05)^2}=-1.996\times10^{6}\ \text{Pa}\,\cdot\,\text{m}^{-3}.

Boxed answer: (∂P/∂T)V=+332.6 Pa/K,(∂P/∂V)T=−2.00×106 Pa/m3\boxed{(\partial P/\partial T)_V=+332.6\ \text{Pa/K},\quad (\partial P/\partial V)_T=-2.00\times10^{6}\ \text{Pa/m}^3}. The subscripts VV and TT explicitly state what is held constant — essential in thermodynamics.

Common Misconceptions

⚠️ Common mistake

The subscript "held constant" is not optional in physics. Writing ∂P/∂T\partial P/\partial T is ambiguous in thermodynamics because PP can be expressed in terms of different variable sets. The notation (∂P/∂T)V(\partial P/\partial T)_V versus (∂P/∂T)S(\partial P/\partial T)_S (constant volume vs. constant entropy) gives genuinely different numbers. Always declare which variables are frozen.

⚠️ Common mistake

Mixed partials are NOT always equal. The interchange fxy=fyxf_{xy}=f_{yx} requires continuity of the mixed second derivatives (Clairaut's hypothesis). Without it, the order of differentiation can matter — see the counterexample.

CounterexampleExistence of Partials Does Not Imply Continuity

Consider

f(x,y)={xyx2+y2,(x,y)≠(0,0)0,(x,y)=(0,0).f(x,y)=\begin{cases}\dfrac{xy}{x^2+y^2}, & (x,y)\neq(0,0)\\[4pt] 0, & (x,y)=(0,0).\end{cases}

At the origin both partials exist: along the xx-axis f(x,0)=0f(x,0)=0, so fx(0,0)=0f_x(0,0)=0; similarly fy(0,0)=0f_y(0,0)=0. Yet ff is not continuous at (0,0)(0,0) — along the line y=xy=x, f(x,x)=12f(x,x)=\tfrac12 for all x≠0x\neq0, so the limit is 12≠0=f(0,0)\tfrac12\neq0=f(0,0). Thus a function can have all first-order partials at a point and still be discontinuous there. Partial differentiability is strictly weaker than total differentiability.

CounterexampleMixed Partials That Disagree

For

g(x,y)={xy(x2−y2)x2+y2,(x,y)≠(0,0)0,(x,y)=(0,0),g(x,y)=\begin{cases}\dfrac{xy(x^2-y^2)}{x^2+y^2}, & (x,y)\neq(0,0)\\[4pt] 0,& (x,y)=(0,0),\end{cases}

one computes gxy(0,0)=−1g_{xy}(0,0)=-1 but gyx(0,0)=+1g_{yx}(0,0)=+1. The mixed second partials exist but are unequal because they fail to be continuous at the origin — exactly the situation Clairaut's theorem excludes.

Connections

connection

Partial derivatives are the components of the Gradient, whose direction gives steepest ascent and whose magnitude gives the maximal rate of change. Assembling partials of a vector-valued map produces the Jacobian Matrix, which governs the multivariable Chain Rule and the Inverse Function Theorem. When set to zero, partials locate stationary points in Multivariable Optimization, with the Hessian Matrix of second partials classifying them. Partial derivatives are the alphabet of Partial Differential Equations, and the Total Differential in [4] links them to Differentiability in Several Variables. They also underlie the Directional Derivative, of which they are the axis-aligned special case.

Real-World Applications

🔬 Where this appears

Historical Development

⏳ History
- 1690s — Leibniz and the Bernoullis differentiate expressions in several variables informally, treating extra variables as parameters.
- 1734 — Leonhard Euler develops systematic methods for functions of several variables and conditions for exactness of differentials.
- 1740s — Alexis Clairaut and Euler independently establish the symmetry of mixed second derivatives for smooth functions.
- 1770s — Lagrange advances the theory in the calculus of variations, treating partial derivatives of action-like functionals.
- 1786 — Adrien-Marie Legendre introduces a rounded "∂\partial" symbol to distinguish partial from total differentiation.
- 1841 — Jacobi's De determinantibus functionalibus firmly standardizes the ∂\partial notation and the use of derivative matrices (Jacobians).
- 1860s–1870s — Weierstrass and others place limits and differentiability on rigorous ε\varepsilon–δ\delta foundations, clarifying the gap between partial and total differentiability.
- 1880s — Hermann Schwarz gives the careful modern statement and proof of the equality-of-mixed-partials theorem (hence "Clairaut–Schwarz").
- 20th century — The notion generalizes to the Fréchet and Gâteaux derivatives in infinite-dimensional spaces, and partial derivatives become central to the theory of distributions and PDEs.

Subtleties & Limits

open-problem

The frontier is not whether partials exist, but how weakly differentiable a function can be while still admitting meaningful "partial derivatives." Sobolev spaces replace classical partials with weak derivatives defined through integration by parts, allowing differential equations to have solutions that are not differentiable in the classical sense. Determining the optimal regularity of solutions to nonlinear PDEs — for instance, whether solutions to the 3D incompressible Navier–Stokes equations stay smooth (a Clay Millennium Problem) — hinges on controlling partial derivatives. Even at the undergraduate boundary, the precise gap between "all partials exist," "all directional derivatives exist," and "fully differentiable" reveals how subtle multivariable smoothness truly is.

See Also

Gradient · Jacobian Matrix · Hessian Matrix · Directional Derivative · Multivariable Chain Rule · Total Differential · Partial Differential Equations · Clairaut's Theorem · Implicit Function Theorem · Maxwell Relations