“How should I calculate and report uncertainty in my physics or chemistry lab report?” “Why do addition and subtraction combine absolute errors, while multiplication and division combine relative errors?” “How do I round the final value and uncertainty according to international standards without losing points?”

Experimental error analysis can feel confusing, but it resolves into three core principles:

  1. Additions and subtractions combine absolute uncertainties in quadrature; multiplications, divisions, and powers combine relative uncertainties in quadrature.
  2. Complex multi-variable formulas are resolved using the general partial derivative law: Δf=∑(∂f∂xiΔxi)2\Delta f = \sqrt{\sum \left(\frac{\partial f}{\partial x_i} \Delta x_i\right)^2}.
  3. In lab reports, round the uncertainty to 1 or 2 significant digits (keep 2 if the leading digit is 1), and align the main value to the exact same decimal place (e.g., 9.8009±0.0692→9.80±0.079.8009 \pm 0.0692 \to 9.80 \pm 0.07).

This guide walks through formula derivations, quadrature versus worst-case bounds, catastrophic cancellation warnings, and four practical worked examples (Ohm’s law, simple pendulum gravity, cylinder density, and temperature difference) verified with step-by-step arithmetic and a free browser tool.

Quick Reference: Error Propagation Formulas

When measurements x,yx, y carry independent uncertainties Δx,Δy\Delta x, \Delta y, the combined uncertainty Δf\Delta f of calculated result ff is:

Operation Formula Uncertainty Rule (Quadrature) Notes & Intuition
Addition f=x+yf = x + y Δf=(Δx)2+(Δy)2\Delta f = \sqrt{(\Delta x)^2 + (\Delta y)^2} Absolute uncertainties add in quadrature; errors never cancel out
Subtraction f=x−yf = x - y Δf=(Δx)2+(Δy)2\Delta f = \sqrt{(\Delta x)^2 + (\Delta y)^2} Identical to addition. Small differences amplify relative error
Constant Multiplier f=c⋅xf = c \cdot x Δf=∣c∣⋅Δx\Delta f = |c| \cdot \Delta x Error scales by ∣c∣|c|. Relative uncertainty is unchanged
Multiplication f=x⋅yf = x \cdot y Δf∣f∣=(Δx∣x∣)2+(Δy∣y∣)2\frac{\Delta f}{|f|} = \sqrt{\left(\frac{\Delta x}{|x|}\right)^2 + \left(\frac{\Delta y}{|y|}\right)^2} Expressed cleanly as relative uncertainties in quadrature
Division f=xyf = \frac{x}{y} Δf∣f∣=(Δx∣x∣)2+(Δy∣y∣)2\frac{\Delta f}{|f|} = \sqrt{\left(\frac{\Delta x}{|x|}\right)^2 + \left(\frac{\Delta y}{|y|}\right)^2} Same as multiplication. Denominator relative error adds in quadrature
Powers f=xnf = x^n Δf∣f∣=∣n∣⋅Δx∣x∣\frac{\Delta f}{|f|} = |n| \cdot \frac{\Delta x}{|x|} Relative error scales linearly by the exponent ∣n∣|n|
General Multivariable f(x1,x2,… )f(x_1, x_2, \dots) Δf=∑i(∂f∂xiΔxi)2\Delta f = \sqrt{\sum_{i} \left(\frac{\partial f}{\partial x_i} \Delta x_i\right)^2} Partial-derivative quadrature (GUM / NIST TN 1297 standard)

What is Relative Uncertainty? The absolute uncertainty $$\Delta x$$ divided by the magnitude of the measured value $$|x|$$ ($$\Delta x / |x|$$, often written as a percentage). For instance, $$5.00 \pm 0.05\text{ V}$$ has a relative uncertainty of $$0.05 / 5.00 = 0.010 = 1.0\%$$.

1. Mathematical Derivation via Partial Derivatives

Why can we compute error propagation using partial derivatives? The mathematical foundation lies in first-order Taylor series approximation (total differentials).

When measured quantities x,yx, y deviate from their true values by small perturbations Δx,Δy\Delta x, \Delta y, the deviation in f(x,y)f(x, y) is:

Δf≈∂f∂xΔx+∂f∂yΔy\Delta f \approx \frac{\partial f}{\partial x} \Delta x + \frac{\partial f}{\partial y} \Delta y

The partial derivative ∂f∂x\frac{\partial f}{\partial x} is called the sensitivity coefficient. It quantifies how steeply the output ff shifts per unit change in input xx.

Error Propagation Decision Flow Diagram

Why Root-Sum-of-Squares (Quadrature)?

If Δx\Delta x and Δy\Delta y stem from independent instruments, random reading fluctuations, or independent calibration errors, the probability that both errors peak simultaneously in the same worst-case direction is statistically very low. Positive and negative fluctuations tend to partially offset each other.

Treating the uncertainties as independent zero-mean random variables with variances σx2,σy2\sigma_x^2, \sigma_y^2, the variance addition theorem yields:

σf2=(∂f∂x)2σx2+(∂f∂y)2σy2\sigma_f^2 = \left(\frac{\partial f}{\partial x}\right)^2 \sigma_x^2 + \left(\frac{\partial f}{\partial y}\right)^2 \sigma_y^2

Taking the square root establishes the combined standard uncertainty formula, mandated by BIPM’s Guide to the Expression of Uncertainty in Measurement (JCGM 100:2008 / GUM) and NIST Technical Note 1297 (Eq. A-3).

2. Quadrature vs. Worst-Case Maximum Error

Depending on your course syllabus or lab manual, you may be instructed to use either the root-sum-of-squares (quadrature) or linear worst-case error. Here is how they compare:

Feature Quadrature (Root-Sum-of-Squares / RMS) Worst-Case Maximum Error (Linear Bound)
Formula Δf=∑(∂f∂xiΔxi)2\Delta f = \sqrt{\sum \left(\frac{\partial f}{\partial x_i} \Delta x_i\right)^2} Δfmax=∑∣∂f∂xi∣Δxi\Delta f_{\text{max}} = \sum \left|\frac{\partial f}{\partial x_i}\right| \Delta x_i
Statistical Premise Errors are uncorrelated and independent All errors align in the same worst possible direction
Uncertainty Magnitude Realistic, expected 68% confidence interval Always larger; overestimates realistic uncertainty
Standard Usage Default in university physics, research & journals Critical structural safety engineering; some introductory labs
Example (Ohm’s Law) R=20.00±0.45 ΩR = 20.00 \pm 0.45\,\Omega R=20.00±0.60 ΩR = 20.00 \pm 0.60\,\Omega

3. Reporting Rules: Significant Figures and Rounding

The most common penalty in lab grading comes from copying all digits from a calculator screen.

International standards (GUM Section 7.2.6 and NIST TN 1297) dictate a strict two-step rounding procedure:

Step 1: Round the uncertainty to 1 or 2 significant digits

  • The standard rule is 1 significant digit (e.g., ±0.043→±0.04\pm 0.043 \to \pm 0.04).
  • If the leading digit of the uncertainty is 1 or 2, keep 2 significant digits to prevent excessive rounding error (e.g., ±0.142→±0.14\pm 0.142 \to \pm 0.14).

Step 2: Align the main value’s decimal place to the uncertainty

Round the calculated value to terminate at the same decimal place as the rounded uncertainty. Digits beyond the uncertainty have no physical significance.

[Incorrect / Penalized]
g = 9.800923 ± 0.06924 m/s²   <- Raw calculator output; too many digits
g = 9.8 ± 0.07 m/s²           <- Mismatched decimal precision (value is coarser)

[Correct / GUM & NIST Compliant]
g = (9.80 ± 0.07) m/s²         <- Uncertainty rounded to 0.07; value matches at 2nd decimal place

4. Four Practical Worked Examples for Lab Reports

Example 1: Ohm’s Law Resistance R=VIR = \frac{V}{I}

A voltmeter reads V=5.00±0.05 VV = 5.00 \pm 0.05\text{ V} and an ammeter reads I=0.250±0.005 AI = 0.250 \pm 0.005\text{ A}.

  1. Calculate Resistance: R=5.000.250=20.00 ΩR = \frac{5.00}{0.250} = 20.00\,\Omega
  2. Component Relative Uncertainties: ΔVV=0.055.00=0.010 (1.0%)\frac{\Delta V}{V} = \frac{0.05}{5.00} = 0.010\ (1.0\%) ΔII=0.0050.250=0.020 (2.0%)\frac{\Delta I}{I} = \frac{0.005}{0.250} = 0.020\ (2.0\%)
  3. Combined Relative Uncertainty (Quotient Rule): ΔRR=(0.010)2+(0.020)2=0.0005≈0.02236 (2.24%)\frac{\Delta R}{R} = \sqrt{(0.010)^2 + (0.020)^2} = \sqrt{0.0005} \approx 0.02236\ (2.24\%)
  4. Absolute Uncertainty & Report Format: ΔR=20.00×0.02236≈0.4472 Ω\Delta R = 20.00 \times 0.02236 \approx 0.4472\,\Omega Rounding to 1 digit yields (20.0±0.4) Ω(20.0 \pm 0.4)\,\Omega; keeping 2 digits yields (20.00±0.45) Ω(20.00 \pm 0.45)\,\Omega.

Example 2: Pendulum Acceleration Due to Gravity g=4π2LT2g = \frac{4\pi^2 L}{T^2}

A pendulum of string length L=1.000±0.005 mL = 1.000 \pm 0.005\text{ m} swings with period T=2.007±0.005 sT = 2.007 \pm 0.005\text{ s}.

  1. Calculate Gravity: g=4π2×1.000(2.007)2≈9.8009 m/s2g = \frac{4 \pi^2 \times 1.000}{(2.007)^2} \approx 9.8009\text{ m/s}^2
  2. Relative Uncertainties:
    • Length LL (power 1): ΔLL=0.0051.000=0.50%\frac{\Delta L}{L} = \frac{0.005}{1.000} = 0.50\%
    • Period TT (power 2 in denominator): By the power rule, its contribution scales by 2: 2×ΔTT=2×0.0052.007≈2×0.2491%=0.4982%2 \times \frac{\Delta T}{T} = 2 \times \frac{0.005}{2.007} \approx 2 \times 0.2491\% = 0.4982\%
  3. Combined Uncertainty in Quadrature: Δgg=(0.005)2+(0.004982)2≈0.007058 (0.71%)\frac{\Delta g}{g} = \sqrt{(0.005)^2 + (0.004982)^2} \approx 0.007058\ (0.71\%) Δg=9.8009×0.007058≈0.0692 m/s2\Delta g = 9.8009 \times 0.007058 \approx 0.0692\text{ m/s}^2
  4. Final Reported Value: Rounding to 1 significant digit yields (9.80±0.07) m/s2(9.80 \pm 0.07)\text{ m/s}^2. Notice how squaring TT doubled its impact on uncertainty.

Example 3: Cylinder Density ρ=mπd2h4\rho = \frac{m}{\frac{\pi d^2 h}{4}}

A solid metal cylinder has mass m=52.3±0.1 gm = 52.3 \pm 0.1\text{ g}, diameter d=2.50±0.02 cmd = 2.50 \pm 0.02\text{ cm}, and height h=10.2±0.1 cmh = 10.2 \pm 0.1\text{ cm}.

  1. Volume and Density: V=π×(2.50)2×10.24≈50.069 cm3V = \frac{\pi \times (2.50)^2 \times 10.2}{4} \approx 50.069\text{ cm}^3 ρ=52.350.069≈1.0446 g/cm3\rho = \frac{52.3}{50.069} \approx 1.0446\text{ g/cm}^3

  2. Relative Uncertainty Breakdown:

    • Mass: 0.152.3≈0.19%\frac{0.1}{52.3} \approx 0.19\%
    • Diameter (squared): 2×0.022.50=1.60%2 \times \frac{0.02}{2.50} = 1.60\%
    • Height: 0.110.2≈0.98%\frac{0.1}{10.2} \approx 0.98\%

  3. Variance Contribution Analysis: Because the diameter is squared, it contributes over 70% of the total variance. In lab discussions, note that improving caliper precision on diameter will yield the greatest accuracy gain.

  4. Combined Result:

    Δρρ=(0.00191)2+(0.01600)2+(0.00980)2≈1.89%\frac{\Delta \rho}{\rho} = \sqrt{(0.00191)^2 + (0.01600)^2 + (0.00980)^2} \approx 1.89\%

    Δρ=1.0446×0.0189≈0.0197 g/cm3\Delta \rho = 1.0446 \times 0.0189 \approx 0.0197\text{ g/cm}^3 With a leading digit of 1, keep two digits: (1.045±0.020) g/cm3(1.045 \pm 0.020)\text{ g/cm}^3.

Example 4: Catastrophic Cancellation in Temperature Difference ΔT=T2−T1\Delta T = T_2 - T_1

In a calorimetry run, initial temperature is T1=34.9±0.1∘CT_1 = 34.9 \pm 0.1^\circ\text{C} and final temperature is T2=35.2±0.1∘CT_2 = 35.2 \pm 0.1^\circ\text{C}.

  1. Temperature Rise: ΔT=35.2−34.9=0.30 K\Delta T = 35.2 - 34.9 = 0.30\text{ K}
  2. Combined Absolute Uncertainty: Δ(ΔT)=(0.1)2+(0.1)2=0.02≈0.1414 K\Delta(\Delta T) = \sqrt{(0.1)^2 + (0.1)^2} = \sqrt{0.02} \approx 0.1414\text{ K}
  3. Catastrophic Relative Error Explosion: While each temperature carried a tiny relative error of 0.1/35≈0.3%0.1 / 35 \approx 0.3\%, the difference carries: 0.14140.30≈47.1%\frac{0.1414}{0.30} \approx 47.1\% A nearly 50% relative uncertainty! This is catastrophic cancellation. When measuring subtle shifts, direct differential sensors (such as differential thermocouples) are required.

5. Common Pitfalls Checklist

  • Are angular uncertainties in radians? Trigonometric derivatives (ddθsin⁡θ=cos⁡θ\frac{d}{d\theta}\sin\theta = \cos\theta) require radians. If your instrument measures Δθ=1∘\Delta\theta = 1^\circ, convert it: Δθ=1×π180≈0.0175 rad\Delta\theta = 1 \times \frac{\pi}{180} \approx 0.0175\text{ rad}.
  • Are variables repeated across terms? In f=xx+yf = \frac{x}{x + y}, treating xx in the numerator and denominator as independent overestimates error. Differentiate the consolidated formula: ∂f∂x=y(x+y)2\frac{\partial f}{\partial x} = \frac{y}{(x+y)^2}.
  • Are you rounding before the final step? Intermediate rounding compounds error. Keep all floating-point digits until the final reporting line.

6. Verifying Results with the Error Propagation Calculator

You can verify any arbitrary formula with automatic symbolic/algorithmic differentiation using the Error Propagation Calculator.

[Error Propagation Calculator Layout]
┌────────────────────────────────────────────────────────┐
│ Formula: g = 4*pi^2*L / T^2                            │
│ [Presets: Pendulum / Ohm's Law / Cylinder / Delta T]   │
├────────────────────────────────────────────────────────┤
│ Quantity L: Value [ 1.000 ]  Uncertainty ± [ 0.005 ] m │
│ Quantity T: Value [ 2.007 ]  Uncertainty ± [ 0.005 ] s │
├────────────────────────────────────────────────────────┤
│ Rounding: [ Auto (2 digits if leading 1) ]  [ RMS ]    │
├────────────────────────────────────────────────────────┤
│ Report Format: (9.80 ± 0.07) m/s² (Relative 0.71%)     │
│ • Dominant contributor: T (49.8% variance), L (50.2%)  │
│ • Step-by-step partial derivative expansion shown      │
│ [ Copy Report Text ]  [ Reset ]                        │
└────────────────────────────────────────────────────────┘
  1. Input Your Formula: Type g = 4*pi^2*L/T^2 or click a quick preset.
  2. Enter Values and Uncertainties: Provide measured values and absolute or percentage uncertainties.
  3. Inspect Sensitivity & Working: Instantly view the variance contribution breakdown to identify which measurement limits experimental accuracy.
  4. Copy Formatted Report: Export formatted Value ± Uncertainty text directly into your lab notebook.

7. Frequently Asked Questions FAQ

Q1. What should I do if no uncertainty is stated for an instrument?

A. If unspecified, assume half of the smallest scale division (0.05 mm0.05\text{ mm} for vernier calipers, 0.5 mm0.5\text{ mm} for millimeter rulers, or ±0.5\pm 0.5 on the least significant digital digit). If the manufacturer provides an accuracy tolerance (e.g., ±0.5%+2 digits\pm 0.5\% + 2\text{ digits}), use that specification.

Q2. Can systematic tolerances and random errors be combined in quadrature?

A. Yes. Under the GUM framework, Type A uncertainties (evaluated statistically) and Type B uncertainties (evaluated from specifications or calibration certificates) are converted into equivalent standard deviations and combined via the same root-sum-of-squares formula.

Q3. Should I report quadrature or worst-case maximum error?

A. Quadrature is the recognized international scientific standard. If your instructor explicitly asks for maximum error, state both: e.g., “g=(9.80±0.07) m/s2g = (9.80 \pm 0.07)\text{ m/s}^2 (independent quadrature) vs. ±0.10 m/s2\pm 0.10\text{ m/s}^2 (worst-case maximum error).”

Standards & Literature

  1. NIST Technical Note 1297: Guidelines for Evaluating and Expressing the Uncertainty of NIST Measurement Results (B. N. Taylor and C. E. Kuyatt)
  2. JCGM 100:2008 (GUM): Guide to the expression of uncertainty in measurement (BIPM, IEC, IFCC, ILAC, ISO, IUPAC, IUPAP, OIML)
  3. OpenStax University Physics Volume 1: Chapter 1 Units and Measurement