CALCULUS 3 • MULTIVARIABLE FUNCTIONS & PARTIAL DERIVATIVES

Differentials & Error Estimation — Differentials and Error Estimation (Multivariable)

Quantify how small measurement errors propagate through functions of several variables using total differentials.

Historical Context & Motivation

The need to understand how small changes in measured quantities affect computed results has driven mathematical inquiry since the birth of calculus itself. When Leibniz and Newton independently developed the differential calculus in the late seventeenth century, they recognized that the differential of a single-variable function provides a linear approximation to the actual change in the function's value. As scientific measurement grew more sophisticated and physical models began to depend on multiple independent quantities simultaneously—temperature and pressure in gas laws, latitude and longitude in geodesy—mathematicians needed a framework for tracking how errors in several inputs combine to produce error in a single output.

1684
Leibniz Publishes 'Nova Methodus'
Gottfried Wilhelm Leibniz introduces the notation dx and dy for infinitesimal changes, laying the symbolic groundwork for the differential. His framework makes it natural to write df = f′(x) dx for functions of one variable.
1797
Lagrange's Théorie des Fonctions
Joseph-Louis Lagrange formalizes derivatives without infinitesimals, but his work on functions of several variables inspires the rigorous treatment of partial derivatives that will underpin multivariable differentials.
1823
Cauchy Formalizes Limits and Differentials
Augustin-Louis Cauchy redefines the differential in terms of limits, giving the total differential of a multivariable function a precise analytical meaning that replaces the vague notion of infinitesimals.
1900s
Error Propagation in Experimental Science
Physicists and engineers adopt the total differential as the standard tool for error propagation, using it to estimate how uncertainties in measured inputs affect computed quantities such as resistance, volume, and energy.

Today the central question remains: if you measure several independent quantities and each measurement carries a small error, how large can the resulting error in the computed quantity be? The total differential answers this question by providing a first-order linear estimate of the change in a function due to small changes in all of its input variables. This section of the course connects single-variable differentials—which you already know well—to the richer multivariable setting where partial derivatives play the role of the ordinary derivative.

Core Principles & Definitions

Before diving into computations, it is essential to establish the conceptual pillars that support the entire framework of multivariable differentials and error estimation. Each of the following ideas generalizes a familiar single-variable concept: the ordinary derivative becomes a collection of partial derivatives, the single-variable differential becomes the total differential, and the linear approximation at a point becomes the tangent plane approximation. Understanding how these ideas relate is the key to applying differentials for error estimation.

1

Total Differential

For z = f(x, y), the total differential is dz = (∂f/∂x) dx + (∂f/∂y) dy. It captures the first-order sensitivity of f to changes in each independent variable.
2

Partial Derivatives as Sensitivity Coefficients

Each partial derivative ∂f/∂xᵢ measures the rate at which f changes per unit change in xᵢ while all other variables are held constant. These serve as weights in the differential.
3

Linear Approximation (Tangent Plane)

Near a point (a, b), f(x, y) ≈ f(a, b) + fₓ(a, b)(x − a) + f_y(a, b)(y − b). The total differential dz approximates the actual change Δz when Δx and Δy are small.
4

Error Propagation

Replace dx and dy with the maximum measurement errors |Δx| and |Δy|. The resulting |dz| provides an upper bound on the propagated error: |Δz| ≈ |∂f/∂x| |Δx| + |∂f/∂y| |Δy|.
5

Relative & Percentage Error

Dividing the absolute error estimate by the function value yields the relative error dz/z, and multiplying by 100 gives the percentage error. These dimensionless measures allow comparison across different quantities.
✦ KEY TAKEAWAY
Think of the total differential as a recipe for uncertainty. Imagine you are baking a cake whose taste depends on the amounts of sugar and butter. If your measuring cups are slightly off, the total differential tells you how those individual inaccuracies combine to alter the final flavor. Each partial derivative acts like a 'sensitivity weight'—the ingredient to which the recipe is most sensitive will dominate the error. In multivariable calculus, the partial derivatives play the same role: they weight each input's uncertainty according to its influence on the output.

Visualizing the Total Differential

In single-variable calculus, the differential dy = f′(x) dx represents the change along the tangent line, which approximates the actual change Δy along the curve. In the multivariable setting, the surface z = f(x, y) has a tangent plane at each point, and the total differential dz represents the change along that plane. The diagram below illustrates this geometry: the gap between dz and Δz shrinks as the perturbations dx and dy decrease.

The curved purple line represents a cross-section of the surface z = f(x, y). The dashed green line is the tangent plane at the point (a, b). The pink segment Δz shows the actual change in f, while the solid green segment dz shows the differential approximation. As dx → 0, the gap between dz and Δz vanishes.

The geometric insight here is profound: the total differential replaces a potentially complicated nonlinear surface with a flat plane in a neighborhood of the point of interest. Every partial derivative contributes one 'tilt' direction to that plane—∂f/∂x tilts it in the x-direction and ∂f/∂y tilts it in the y-direction. For functions of three or more variables, the tangent plane becomes a tangent hyperplane, but the algebraic form of the total differential remains the same sum of partial derivative terms.

Mathematical Framework

We now formalize the total differential and its application to error estimation. Let f be a differentiable function of n independent variables x₁, x₂, …, xₙ. The total differential is the natural extension of the single-variable differential df = f′(x) dx to a setting where each variable contributes its own partial derivative and its own perturbation.

TOTAL DIFFERENTIAL (TWO VARIABLES)
dz = (∂f/∂x) dx + (∂f/∂y) dy
where z = f(x, y), and dx = Δx, dy = Δy are the actual (small) changes in x and y. The partial derivatives ∂f/∂x and ∂f/∂y are evaluated at the point of interest (a, b).
TOTAL DIFFERENTIAL (n VARIABLES)
df = (∂f/∂x₁) dx₁ + (∂f/∂x₂) dx₂ + ⋯ + (∂f/∂xₙ) dxₙ = Σᵢ₌₁ⁿ (∂f/∂xᵢ) dxᵢ
This extends naturally to any number of independent variables. Each term represents the first-order contribution of the i-th variable to the total change in f.
MAXIMUM ABSOLUTE ERROR
|Δf| ≈ |dz| ≤ |∂f/∂x| |Δx| + |∂f/∂y| |Δy|
By taking absolute values, we obtain a worst-case (maximum) estimate for the propagated error. This is conservative: it assumes all errors compound in the same direction.
RELATIVE AND PERCENTAGE ERROR
Relative error = |df|/|f|, Percentage error = (|df|/|f|) × 100%
The relative error normalizes the absolute error by the magnitude of the function value, making it a dimensionless quantity that facilitates comparison across different physical quantities or scales.
⚠️ When is the Approximation Valid?
The total differential provides a good approximation only when dx and dy are small relative to the scale on which the partial derivatives change appreciably. Formally, the error in the approximation is o(√(dx² + dy²)), meaning it shrinks faster than the perturbation itself. For practical purposes, if measurement errors are within a few percent of the measured values, the linear estimate is typically quite accurate.

Anatomy of Error Propagation

Understanding how each variable contributes to the total error is crucial for experimental design and engineering tolerance analysis. The diagram below dissects the error contributions for a function of three variables, showing how the magnitude of each partial derivative and each input error combine to produce the total propagated uncertainty.

This flow diagram shows how three input uncertainties combine into a total propagated error. Each input's contribution equals the product of its sensitivity (partial derivative magnitude) and its measurement error. Note that variable z dominates the total error even though its measurement uncertainty is the smallest—this is because the function is most sensitive to z.

This decomposition reveals a critical engineering insight: to reduce the total error most efficiently, you should invest effort in reducing the uncertainty of the variable whose error contribution is largest. That contribution depends on two factors—the size of the measurement error and the magnitude of the corresponding partial derivative. A large partial derivative amplifies even a tiny measurement error, while a small partial derivative can suppress a relatively large measurement error. In the example above, improving the precision of z-measurements would yield the greatest reduction in overall uncertainty.

💡 Logarithmic Differentiation for Products & Quotients
When f is a product or quotient of its variables (e.g., f = x²y/z), taking the total differential of ln|f| converts the problem into a sum of relative errors. The result is: |df/f| ≈ |2 dx/x| + |dy/y| + |dz/z|. The exponents become coefficients, making the relative error formula particularly elegant and easy to remember for power-law functions.

Worked Example

Consider a right circular cylinder whose radius r is measured as 5.00 cm with a maximum error of ±0.03 cm, and whose height h is measured as 12.00 cm with a maximum error of ±0.05 cm. We wish to estimate the maximum error and the percentage error in the computed volume V = πr²h.

Maximum Error in the Volume of a Cylinder
1
Step 1 — Identify the Function and Given ValuesThe volume is V = πr²h. The measured values are r = 5.00 cm and h = 12.00 cm. The maximum measurement errors are |Δr| = 0.03 cm and |Δh| = 0.05 cm.
2
Step 2 — Compute the Partial DerivativesDifferentiate V with respect to each variable: ∂V/∂r = 2πrh and ∂V/∂h = πr². Evaluate at r = 5, h = 12: ∂V/∂r = 2π(5)(12) = 120π ≈ 376.99 cm² ∂V/∂h = π(5)² = 25π ≈ 78.54 cm².
∂V/∂r = 120π, ∂V/∂h = 25π
3
Step 3 — Form the Total DifferentialThe total differential is dV = (∂V/∂r) dr + (∂V/∂h) dh = 120π dr + 25π dh.
4
Step 4 — Estimate the Maximum Absolute ErrorTake absolute values and substitute the maximum errors: |dV| ≤ |120π|(0.03) + |25π|(0.05) = 3.6π + 1.25π = 4.85π ≈ 15.24 cm³.
Maximum absolute error ≈ 15.24 cm³
5
Step 5 — Compute the Percentage ErrorThe exact volume is V = π(5)²(12) = 300π ≈ 942.48 cm³. The relative error is |dV|/V = 4.85π / 300π = 4.85/300 ≈ 0.01617. Multiplying by 100 gives a percentage error of approximately 1.62%. Notice that the radius contributes (3.6π)/(4.85π) ≈ 74.2% of the total error, even though its measurement error (0.03 cm) is smaller than that of the height (0.05 cm). This is because V depends on r² — the sensitivity to r is amplified by the exponent.
Percentage error ≈ 1.62%
🔍 OBSERVATION
Using the logarithmic differentiation shortcut for V = πr²h: dV/V = 2(dr/r) + (dh/h) = 2(0.03/5) + (0.05/12) = 0.012 + 0.00417 ≈ 0.01617, which gives the same 1.62% relative error with less computation. This shortcut works because V is a product of powers of its variables.

Strengths, Limitations & Comparisons

The total differential method for error estimation is remarkably versatile, but it is important to understand both its power and its boundaries. The table below compares the differential approach with the two most common alternative methods for error propagation.

Comparison of three common error propagation methods
FeatureTotal Differential (Linear)Root-Sum-of-Squares (RSS)Monte Carlo Simulation
AssumptionsErrors are small; function is differentiable; worst-case (all errors compound)Errors are small, independent, and normally distributedRequires known probability distributions for inputs; no linearity needed
Formula|df| = Σ |∂f/∂xᵢ| |Δxᵢ|σ_f = √(Σ (∂f/∂xᵢ)² σᵢ²)Numerical sampling — no closed-form formula
Type of estimateWorst-case upper boundProbable (1σ) estimateFull probability distribution
Computational costVery low — pencil-and-paperLow — pencil-and-paperHigh — requires software and many samples
Accuracy for large errorsPoor — linear approximation breaks downModerate — still first-orderExcellent — handles nonlinearity and correlations
✦ KEY TAKEAWAY
The total differential method is the conservative choice: it gives a worst-case error bound by assuming all individual errors conspire to push the output in the same direction. In practice, random errors partially cancel, which is why the root-sum-of-squares method often gives a more realistic estimate. Think of it this way: if you are navigating by adding several compass bearings, the total differential assumes every bearing error steers you in the worst direction, while RSS assumes the errors are random and partially cancel. For safety-critical applications (bridge design, pharmaceutical dosing), the conservative bound from the total differential is preferred.

Connection to Advanced Theory

The total differential is the entry point to several powerful generalizations that appear later in multivariable calculus and beyond. Understanding where the differential fits within this broader landscape helps solidify the current material and prepares you for what follows.

From differentials to advanced theory
Concept in This LessonAdvanced GeneralizationWhere You'll See It
Total differential dz = fₓ dx + f_y dyGradient vector ∇f = ⟨fₓ, f_y⟩ and the dot product df = ∇f · d𝐫Directional derivatives, gradient fields
Tangent plane approximationTaylor series in several variables (second-order and higher terms)Optimization, Hessian matrix, second derivative test
Linear error propagation |df| = Σ |∂f/∂xᵢ| |Δxᵢ|Jacobian matrix for vector-valued functions; covariance propagationCoordinate transformations, statistical modeling, Kalman filtering
Differentiability as existence of total differentialFréchet derivative on Banach spacesFunctional analysis, differential geometry

The most immediate connection is with the gradient vector. Writing dz = fₓ dx + f_y dy is equivalent to dz = ∇f · ⟨dx, dy⟩, so the total differential is simply the dot product of the gradient with the displacement vector d𝐫. This interpretation gives the differential a geometric meaning: it measures the component of the function's steepest ascent direction along the direction of the perturbation. When you later encounter the Jacobian matrix for vector-valued functions 𝐅: ℝⁿ → ℝᵐ, each row of the Jacobian is precisely the gradient of one component function, so the Jacobian generalizes the total differential to the vector-valued setting.

Practice Problems

PROBLEM 1 — CONCEPTUAL
Explain in your own words why the total differential dz is generally not equal to the actual change Δz. Under what condition does dz become an increasingly accurate approximation to Δz?
PROBLEM 2 — BASIC CALCULATION
Let f(x, y) = x³y − 2xy². Compute the total differential df and evaluate it at (1, 2) with dx = 0.01 and dy = −0.02.
PROBLEM 3 — INTERMEDIATE
The resistance R of a wire is given by R = ρL/A, where the resistivity ρ = 1.7 × 10⁻⁸ Ω·m (exact), the length L is measured as 50.0 m with |ΔL| = 0.2 m, and the cross-sectional area A is measured as 2.00 × 10⁻⁶ m² with |ΔA| = 0.05 × 10⁻⁶ m². Use the total differential to estimate the maximum absolute error and the percentage error in R.
PROBLEM 4 — APPLIED
An engineer calculates the power dissipated in a resistor using P = V²/R. The voltage is measured as V = 120.0 V ± 0.5 V and the resistance as R = 60.0 Ω ± 0.4 Ω. Estimate the maximum absolute error in P, the percentage error, and determine which measurement contributes more to the uncertainty in P.
PROBLEM 5 — CRITICAL THINKING
Let f(x, y, z) = ln(x² + y² + z²). Show that the relative error |df/f| depends not only on the measurement errors but also on the point (x, y, z) at which the function is evaluated. Then discuss: for a fixed set of measurement errors, at which types of points (x, y, z) will the relative error be largest, and at which will it be smallest? Assume x² + y² + z² > 1 and f > 0.

Lesson Summary

The total differential of a function z = f(x₁, …, xₙ) is dz = Σ (∂f/∂xᵢ) dxᵢ, and it provides a first-order linear approximation to the actual change Δz. Each partial derivative acts as a sensitivity coefficient, weighting the corresponding input perturbation. The tangent plane (or tangent hyperplane in higher dimensions) is the geometric realization of this linear approximation, and the quality of the approximation improves as the perturbations shrink.

For error estimation, replace dx and dy with maximum measurement errors and take absolute values: |df| ≤ Σ |∂f/∂xᵢ| |Δxᵢ|. This yields a conservative, worst-case maximum absolute error bound. Dividing by |f| gives the relative error. For functions that are products or quotients of their variables, logarithmic differentiation provides an efficient shortcut. To reduce total propagated error most effectively, focus on the input variable whose error contribution — the product of its partial derivative magnitude and its measurement uncertainty — is largest.

Varsity Tutors • Calculus 3 • Differentials & Error Estimation — Differentials and Error Estimation (Multivariable)