Navigation

Concept Review

Differential Calculus

Linear Approximation and Differentials

Linear approximation is the practical payoff of the derivative: at a point where ff is differentiable, the graph of ff is indistinguishable from its tangent line, so a hard function can be replaced locally by an easy linear one. This topic turns that idea into a formula (the linearisation), a notation for small changes (differentials), and — the part students skip — a way to say how wrong the answer is and in which direction. It is the first-order Taylor polynomial you will meet again in integral calculus, and it is the engine behind every "propagation of measurement error" calculation in physics and engineering.

The linearisation and the tangent-line approximation

Definition. If ff is differentiable at aa, the linearisation of ff at aa is L(x)=f(a)+f(a)(xa).L(x) = f(a) + f'(a)(x-a). The graph y=L(x)y = L(x) is exactly the tangent line to y=f(x)y=f(x) at the point (a,f(a))\left(a, f(a)\right). The linear approximation (tangent-line approximation) is the statement f(x)L(x)for x near a.f(x) \approx L(x) \qquad \text{for } x \text{ near } a.

Two hypotheses carry all the weight: f(a)f'(a) must exist, and xx must be near aa. In fact differentiability at aa is precisely the statement that f(x)=f(a)+f(a)(xa)+E(x)(xa),E(x)0 as xa,f(x) = f(a) + f'(a)(x-a) + E(x)(x-a), \qquad E(x) \to 0 \text{ as } x \to a, so the error is not merely small — it is small compared with xax-a. Nothing here promises anything far from aa: sinxx\sin x \approx x is excellent at x=0.05x = 0.05 and useless at x=3x = 3.

Choosing aa. Pick the nearest point at which f(a)f(a) and f(a)f'(a) are computable exactly: a perfect square or cube, x=1x=1, x=0x=0, π/6\pi/6, π/4\pi/4. To estimate 26\sqrt{26} use a=25a=25, not a=16a=16.

Standard linearisations at a=0a = 0

f(x)f(x)L(x)L(x) near 00Sign of ff'' near 00LL over- or under-estimates
(1+x)k(1+x)^{k} (needs x>1x>-1)1+kx1+kxsign of k(k1)k(k-1)over if 0<k<10<k<1; under if k<0k<0 or k>1k>1; exact if k=0k=0 or k=1k=1
1+x\sqrt{1+x}1+x21+\dfrac{x}{2}negativeover
11+x\dfrac{1}{1+x}1x1-xpositiveunder
exe^{x}1+x1+xpositiveunder
ln(1+x)\ln(1+x)xxnegativeover
sinx\sin xxxopposite sign to xxover for x>0x>0, under for x<0x<0
tanx\tan xxxsame sign as xxunder for x>0x>0, over for x<0x<0
arctanx\arctan xxxopposite sign to xxover for x>0x>0, under for x<0x<0
cosx\cos x11negativeover

Note the last row: f(0)=0f'(0)=0, so the linearisation of cos\cos at 00 is the constant 11. That is not a failure — it is why cosx1\cos x \approx 1 is unusually accurate for small xx (the error is of order x2x^{2}, not xx).

Differentials

Definition. Let y=f(x)y = f(x) with ff differentiable. The symbol dxdx is an independent variable: it may be given any real value. The differential dydy is then defined by dy=f(x)dx.dy = f'(x) dx.

Compare it with the true change. Fix xx and a step Δx\Delta x, and set dx=Δxdx = \Delta x. Then

  • Δy=f(x+Δx)f(x)\Delta y = f(x + \Delta x) - f(x) is the actual rise along the curve;
  • dy=f(x)Δxdy = f'(x)\Delta x is the rise along the tangent line.

So Δydy\Delta y \approx dy, and Δydy0\Delta y - dy \to 0 faster than Δx\Delta x does. Evaluating the differential at the base point — the derivative taken at aa, with dx=xadx = x-a — recovers the linearisation: f(x)f(a)+dy=L(x)f(x) \approx f(a) + dy = L(x). Once dydy and dxdx are separate objects, the Leibniz symbol dydx=f(x)\dfrac{dy}{dx} = f'(x) really is a quotient — which is why chain-rule and separable-equation notation behaves so well.

Estimating a small change: the recipe

  1. Identify the quantity y=f(x)y = f(x) and the base point aa where f(a)f(a), f(a)f'(a) are exact.
  2. Write dx=xadx = x - a (the small change or the measurement error).
  3. Compute dy=f(a)dxdy = f'(a) dx.
  4. Report f(x)f(a)+dyf(x) \approx f(a) + dy, or "the change in yy is about dydy".
  5. Say whether it is an over- or under-estimate (concavity) and, if asked, bound the error.

Absolute, relative and percentage error

For a quantity qq with estimated error dqdq:

  • absolute error =dq= \left| dq \right| (carries units);
  • relative error =dqq= \dfrac{\left| dq \right|}{\left| q \right|} (dimensionless);
  • percentage error =100= 100 \cdot relative error.

Power rule for error propagation. If y=Cxny = Cx^{n} with C0C \ne 0 constant, and x0x \ne 0 (so neither denominator below vanishes), then dy=Cnxn1dxdy = Cnx^{n-1} dx, so dyy=Cnxn1dxCxn=ndxx.\frac{dy}{y} = \frac{Cnx^{n-1} dx}{Cx^{n}} = n \cdot \frac{dx}{x}. The percentage error in yy is about n\left| n \right| times the percentage error in xx. Hence a 1%1\% error in a radius gives about 2%2\% in an area and 3%3\% in a volume, while T=2πL/gT = 2\pi\sqrt{L/g} turns a 1%1\% error in LL into only 0.5%0.5\% in TT.

Over- or under-estimate: read it off the concavity

Theorem. Suppose ff'' exists on an open interval II containing aa.

  • If f(x)0f''(x) \ge 0 for all xx in II (ff concave up), then f(x)L(x)f(x) \ge L(x) on II: the tangent line lies below the graph, so LL under-estimates.
  • If f(x)0f''(x) \le 0 for all xx in II (ff concave down), then f(x)L(x)f(x) \le L(x) on II: the tangent lies above, so LL over-estimates.

It makes no difference whether x>ax > a or x<ax < a — only the sign of ff'' matters. If ff'' changes sign between aa and xx (an inflection point in between), no conclusion is available.

Error bounds

Theorem (Taylor's inequality, first order). Suppose ff' is continuous on the closed interval with endpoints aa and xx, and ff'' exists on the open interval between them. Then there is a number cc strictly between aa and xx with f(x)L(x)=f(c)2(xa)2.f(x) - L(x) = \frac{f''(c)}{2}(x-a)^{2}. Consequently, if f(t)M\left| f''(t) \right| \le M for every tt between aa and xx, then f(x)L(x)M2(xa)2.\left| f(x) - L(x) \right| \le \frac{M}{2}(x-a)^{2}.

Three consequences worth memorising: (i) MM must bound f\left| f'' \right| across the whole interval, not just at aa; (ii) any valid MM gives a valid bound, so a clean over-estimate of MM is perfectly acceptable; (iii) the bound is quadratic in the step, so halving xa\left| x-a \right| quarters it.

Worked example 1 — a cube root, with direction and error bound

Estimate 8.063\sqrt[3]{8.06}, decide whether the estimate is too big or too small, and bound the error.

Take f(x)=x1/3f(x) = x^{1/3} and a=8a = 8 (the nearest perfect cube).

f(8)=2,f(x)=13x2/3,f(8)=134=112,L(x)=2+112(x8).\begin{aligned} f(8) &= 2, \\ f'(x) &= \frac{1}{3}x^{-2/3}, \qquad f'(8) = \frac{1}{3 \cdot 4} = \frac{1}{12}, \\ L(x) &= 2 + \frac{1}{12}(x-8). \end{aligned}

With x=8.06x = 8.06, so dx=0.06dx = 0.06: L(8.06)=2+0.0612=2+0.005=2.005.L(8.06) = 2 + \frac{0.06}{12} = 2 + 0.005 = 2.005.

Direction. f(x)=29x5/3<0f''(x) = -\dfrac{2}{9}x^{-5/3} < 0 for x>0x>0, so ff is concave down and 2.0052.005 is an over-estimate.

Bound. For tt in [8,8.06][8, 8.06] we have t8t \ge 8, hence t5/385/3=32t^{5/3} \ge 8^{5/3} = 32 and f(t)=29t5/329132=1144=M.\left| f''(t) \right| = \frac{2}{9}t^{-5/3} \le \frac{2}{9} \cdot \frac{1}{32} = \frac{1}{144} = M. Therefore 8.0632.0051288(0.06)2=0.0036288=1.25×105.\left| \sqrt[3]{8.06} - 2.005 \right| \le \frac{1}{288}(0.06)^{2} = \frac{0.0036}{288} = 1.25 \times 10^{-5}. Combining with the direction: 2.00498758.0632.0052.0049875 \le \sqrt[3]{8.06} \le 2.005. (The true value is 2.00498762.0049876, so both statements are correct.)

Worked example 2 — measurement error

The radius of a circular disk is measured as 2424 cm with a maximum error of 0.020.02 cm. Estimate the resulting maximum error, relative error and percentage error in the computed area.

A=πr2A = \pi r^{2}, so dA=2πrdrdA = 2\pi r dr. With r=24r = 24 and dr0.02\left| dr \right| \le 0.02: dA2π(24)(0.02)=0.96π3.02 cm2.\left| dA \right| \le 2\pi (24)(0.02) = 0.96\pi \approx 3.02 \text{ cm}^{2}. Relative error: dAA0.96ππ(24)2=0.96576=16000.00167,\frac{\left| dA \right|}{A} \le \frac{0.96\pi}{\pi (24)^{2}} = \frac{0.96}{576} = \frac{1}{600} \approx 0.00167, i.e. about 0.167%0.167\%. The shortcut dAA=2drr=2(0.02)24=1600\dfrac{dA}{A} = 2\dfrac{dr}{r} = \dfrac{2(0.02)}{24} = \dfrac{1}{600} gives the same thing with no arithmetic. (Exactly: A(24.02)A(24)=0.9604π3.017A(24.02) - A(24) = 0.9604\pi \approx 3.017, against the estimate 3.0163.016.)

Common mistakes

WrongRight
L(x)=f(a)+f(x)(xa)L(x) = f(a) + f'(x)(x-a)L(x)=f(a)+f(a)(xa)L(x) = f(a) + f'(a)(x-a) — the derivative is frozen at aa, so LL is genuinely linear in xx
Estimating sin31\sin 31^{\circ} with dx=1dx = 1Convert to radians first: dx=π/1800.01745dx = \pi/180 \approx 0.01745
Using a=8.06a = 8.06 to estimate 8.063\sqrt[3]{8.06}Then xa=0x-a=0 and you learn nothing; aa must be a nearby point you can evaluate exactly
"x>ax > a, so the estimate is too small"The side of aa is irrelevant; the sign of ff'' decides over- vs under-estimate
Treating Δy\Delta y and dydy as equaldydy is the tangent-line rise; Δydy\Delta y - dy is nonzero. Differentiability alone gives only Δydy0\Delta y - dy \to 0 faster than Δx\Delta x; the familiar order (Δx)2(\Delta x)^{2} needs ff'' to exist and be bounded near xx (for f(x)=x4/3f(x)=x^{4/3} at x=0x=0 the gap is (Δx)4/3(\Delta x)^{4/3}, far bigger than (Δx)2(\Delta x)^{2})
Writing dy=f(x)dy = f'(x)dy=f(x)dxdy = f'(x) dx — without the dxdx the equation is dimensionally meaningless
Taking $M = \leftf''(a) \right
Reporting 0.001670.00167 as "the percentage error"Percentage error is 100×100 \times the relative error: 0.167%0.167\%
Linearising $\leftx \right

Key terms

  • linearisation L(x)
  • tangent-line approximation
  • base point a
  • differential dy = f'(x) dx
  • dx as an independent variable
  • actual change Delta y vs. differential dy
  • first-order Taylor polynomial
  • standard small-x approximations
  • absolute error
  • relative error
  • percentage error
  • error propagation
  • power rule for relative error
  • concavity and over-estimation
  • concavity and under-estimation
  • tangent line below/above the graph
  • Taylor's inequality (first order)
  • Lagrange remainder f''(c)(x-a)^2/2
  • bound M on the second derivative
  • quadratic decay of the error bound

Practice Problems

Find the linearisation of f(x)=xf(x)=\sqrt{x} at a=9a=9 and use it to approximate 9.1\sqrt{9.1}. Is your estimate too large or too small? Justify the direction.

Show hint

You need f(9)f(9) and f(9)f'(9) exactly — 99 is a perfect square, which is why it was chosen. For the direction, look at the sign of ff''.

Show answer

Set up. f(x)=x=x1/2f(x)=\sqrt{x}=x^{1/2}, base point a=9a=9.

f(9)=9=3,f(x)=12x,f(9)=123=16.\begin{aligned} f(9) &= \sqrt{9} = 3, \\ f'(x) &= \frac{1}{2\sqrt{x}}, \qquad f'(9) = \frac{1}{2 \cdot 3} = \frac{1}{6}. \end{aligned}

Linearisation. L(x)=f(9)+f(9)(x9)=3+16(x9).L(x) = f(9) + f'(9)(x-9) = 3 + \frac{1}{6}(x-9).

Approximate. Put x=9.1x = 9.1, so x9=0.1x - 9 = 0.1: 9.1L(9.1)=3+0.16=3+160=3.0166667 (7 d.p.).\sqrt{9.1} \approx L(9.1) = 3 + \frac{0.1}{6} = 3 + \frac{1}{60} = 3.0166667 \text{ (7 d.p.)}.

Direction. f(x)=14x3/2=14x3/2<0f''(x) = -\dfrac{1}{4}x^{-3/2} = -\dfrac{1}{4x^{3/2}} < 0 for x>0x>0, so ff is concave down on (0,)(0,\infty) and the tangent line lies above the graph. The estimate is therefore too large.

Independent check. Square the estimate: (3+160)2=9+23160+13600=9+0.1+0.0002778=9.1002778>9.1,\left(3 + \frac{1}{60}\right)^{2} = 9 + 2 \cdot 3 \cdot \frac{1}{60} + \frac{1}{3600} = 9 + 0.1 + 0.0002778 = 9.1002778 > 9.1, confirming that 3.0166667>9.13.0166667 > \sqrt{9.1}.

Answer: L(x)=3+x96L(x) = 3 + \dfrac{x-9}{6} and 9.1181603.01667\sqrt{9.1} \approx \dfrac{181}{60} \approx 3.01667, an over-estimate. (True value 3.01662063.0166206; the error is about 4.6×1054.6\times 10^{-5}.)

Let y=x32xy = x^{3} - 2x.

(a) Find the differential dydy.

(b) Evaluate dydy and the actual change Δy\Delta y when xx changes from 22 to 2.12.1.

(c) By how much do they differ, and what happens to that difference if the step is 0.010.01 instead?

Show hint

dy=f(x)dxdy = f'(x)dx with dxdx equal to the step; Δy=f(x+Δx)f(x)\Delta y = f(x+\Delta x) - f(x) is computed from the function itself, not from the tangent line.

Show answer

(a) With f(x)=x32xf(x) = x^{3}-2x we have f(x)=3x22f'(x) = 3x^{2}-2, so dy=(3x22)dx.dy = \left(3x^{2}-2\right) dx.

(b) Here x=2x = 2 and dx=Δx=0.1dx = \Delta x = 0.1: dy=(3(2)22)(0.1)=(122)(0.1)=10(0.1)=1.dy = \left(3(2)^{2}-2\right)(0.1) = (12-2)(0.1) = 10(0.1) = 1.

For the actual change,

f(2.1)=(2.1)32(2.1)=9.2614.2=5.061,f(2)=84=4,Δy=5.0614=1.061.\begin{aligned} f(2.1) &= (2.1)^{3} - 2(2.1) = 9.261 - 4.2 = 5.061, \\ f(2) &= 8 - 4 = 4, \\ \Delta y &= 5.061 - 4 = 1.061 . \end{aligned}

(c) Δydy=1.0611=0.061\Delta y - dy = 1.061 - 1 = 0.061.

To see the pattern, expand with a general step hh:

Δy=[(2+h)32(2+h)][84]=[8+12h+6h2+h342h]4=10h+6h2+h3.\begin{aligned} \Delta y &= \left[(2+h)^{3} - 2(2+h)\right] - \left[8-4\right] \\ &= \left[8 + 12h + 6h^{2} + h^{3} - 4 - 2h\right] - 4 \\ &= 10h + 6h^{2} + h^{3}. \end{aligned}

The differential is exactly the linear part, dy=10hdy = 10h, so Δydy=6h2+h3.\Delta y - dy = 6h^{2} + h^{3}. With h=0.1h = 0.1 this is 0.06+0.001=0.0610.06 + 0.001 = 0.061 (matching the direct computation). With h=0.01h = 0.01 it is 6(0.0001)+0.000001=0.0006016(0.0001) + 0.000001 = 0.000601 — roughly 100100 times smaller for a step 1010 times smaller, exactly the quadratic behaviour predicted by the error bound.

Answer: dy=(3x22)dxdy = (3x^{2}-2)dx; at x=2x=2, dx=0.1dx = 0.1: dy=1dy = 1, Δy=1.061\Delta y = 1.061, difference 0.0610.061. For dx=0.01dx = 0.01 the difference drops to 0.0006010.000601.

Use a linear approximation at a suitable point to estimate 14.002\dfrac{1}{4.002}. State whether the estimate is above or below the true value.

Show hint

Take f(x)=1/xf(x)=1/x and choose the nearby point where ff and ff' are exact. Concavity of 1/x1/x on x>0x>0 settles the direction.

Show answer

Set up. f(x)=1x=x1f(x) = \dfrac{1}{x} = x^{-1}, base point a=4a = 4 (nearby and easy).

f(4)=14=0.25,f(x)=x2=1x2,f(4)=116=0.0625,L(x)=14116(x4).\begin{aligned} f(4) &= \frac{1}{4} = 0.25, \\ f'(x) &= -x^{-2} = -\frac{1}{x^{2}}, \qquad f'(4) = -\frac{1}{16} = -0.0625, \\ L(x) &= \frac{1}{4} - \frac{1}{16}(x-4). \end{aligned}

Approximate. With x=4.002x = 4.002, dx=0.002dx = 0.002: 14.002L(4.002)=0.250.00216=0.250.000125=0.249875.\frac{1}{4.002} \approx L(4.002) = 0.25 - \frac{0.002}{16} = 0.25 - 0.000125 = 0.249875 .

Direction. f(x)=2x3=2x3>0f''(x) = 2x^{-3} = \dfrac{2}{x^{3}} > 0 for x>0x>0, so ff is concave up there, the tangent line lies below the graph, and 0.2498750.249875 is an under-estimate.

Independent check. Multiply back: 0.249875×4.002=0.9995+0.00049975=0.99999975<1,0.249875 \times 4.002 = 0.9995 + 0.00049975 = 0.99999975 < 1, so 0.249875<14.0020.249875 < \dfrac{1}{4.002}, confirming the under-estimate.

Answer: 14.0020.249875\dfrac{1}{4.002} \approx 0.249875, slightly below the true value 0.249875060.24987506 (error about 6.2×1086.2 \times 10^{-8}).

Find the linearisation of f(x)=(1+x)1/3f(x) = (1+x)^{1/3} at a=0a=0 and use it to approximate 1.033\sqrt[3]{1.03} and 0.973\sqrt[3]{0.97}. Are these over- or under-estimates?

Show hint

This is the binomial approximation (1+x)k1+kx(1+x)^{k} \approx 1 + kx. For the direction, note that the sign of ff'' does not depend on which side of 00 you are on.

Show answer

Linearisation. f(x)=(1+x)1/3f(x) = (1+x)^{1/3}, so

f(0)=1,f(x)=13(1+x)2/3,f(0)=13,L(x)=1+x3.\begin{aligned} f(0) &= 1, \\ f'(x) &= \frac{1}{3}(1+x)^{-2/3}, \qquad f'(0) = \frac{1}{3}, \\ L(x) &= 1 + \frac{x}{3}. \end{aligned}

(This is the case k=13k = \frac{1}{3} of (1+x)k1+kx(1+x)^{k} \approx 1+kx.)

First estimate. 1.033\sqrt[3]{1.03} corresponds to 1+x=1.031 + x = 1.03, i.e. x=0.03x = 0.03: 1.0331+0.033=1.01.\sqrt[3]{1.03} \approx 1 + \frac{0.03}{3} = 1.01 .

Second estimate. 0.973\sqrt[3]{0.97} corresponds to x=0.03x = -0.03: 0.9731+0.033=0.99.\sqrt[3]{0.97} \approx 1 + \frac{-0.03}{3} = 0.99 .

Direction. f(x)=13(23)(1+x)5/3=29(1+x)5/3<0for x>1,f''(x) = \frac{1}{3}\left(-\frac{2}{3}\right)(1+x)^{-5/3} = -\frac{2}{9}(1+x)^{-5/3} < 0 \quad \text{for } x > -1, so ff is concave down on (1,)(-1,\infty) and the tangent line lies above the graph on both sides of 00. Both estimates are over-estimates.

Independent check. Cube each estimate:

(1.01)3=1.030301>1.031.01>1.033,(0.99)3=0.970299>0.970.99>0.973.\begin{aligned} (1.01)^{3} &= 1.030301 > 1.03 \quad \Rightarrow \quad 1.01 > \sqrt[3]{1.03}, \\ (0.99)^{3} &= 0.970299 > 0.97 \quad \Rightarrow \quad 0.99 > \sqrt[3]{0.97}. \end{aligned}

Answer: L(x)=1+x3L(x) = 1 + \dfrac{x}{3}; 1.0331.01\sqrt[3]{1.03} \approx 1.01 and 0.9730.99\sqrt[3]{0.97} \approx 0.99, both over-estimates (true values 1.0099021.009902 and 0.9898980.989898).

Use a linear approximation to estimate cos29\cos 29^{\circ}. Give the exact expression and a decimal to 66 decimal places, say whether it over- or under-estimates, and bound the error.

Show hint

Work in radians throughout and base the approximation at 3030^{\circ}, where the cosine and sine are known exactly. The step dxdx is one degree expressed in radians.

Show answer

Convert to radians. 29=29π18029^{\circ} = \dfrac{29\pi}{180} and the base point is a=30=π6=30π180a = 30^{\circ} = \dfrac{\pi}{6} = \dfrac{30\pi}{180}. Hence dx=xa=29π18030π180=π1800.0174533.dx = x - a = \frac{29\pi}{180} - \frac{30\pi}{180} = -\frac{\pi}{180} \approx -0.0174533 .

Set up. f(x)=cosxf(x) = \cos x, so f(x)=sinxf'(x) = -\sin x and f(π6)=32,f(π6)=sinπ6=12.f\left(\frac{\pi}{6}\right) = \frac{\sqrt{3}}{2}, \qquad f'\left(\frac{\pi}{6}\right) = -\sin\frac{\pi}{6} = -\frac{1}{2}.

Linearisation. L(x)=3212(xπ6).L(x) = \frac{\sqrt{3}}{2} - \frac{1}{2}\left(x - \frac{\pi}{6}\right).

Evaluate.

cos293212(π180)=32+π3600.8660254+0.0087266=0.8747520.\begin{aligned} \cos 29^{\circ} &\approx \frac{\sqrt{3}}{2} - \frac{1}{2}\left(-\frac{\pi}{180}\right) \\ &= \frac{\sqrt{3}}{2} + \frac{\pi}{360} \\ &\approx 0.8660254 + 0.0087266 \\ &= 0.8747520 . \end{aligned}

To six decimal places, cos290.874752\cos 29^{\circ} \approx 0.874752.

Direction. f(x)=cosxf''(x) = -\cos x, and cosx>0\cos x > 0 throughout [29,30]\left[29^{\circ}, 30^{\circ}\right], so f<0f'' < 0 there: ff is concave down, the tangent lies above the graph, and this is an over-estimate.

Error bound. For tt in that interval, f(t)=costcos29<0.875\left| f''(t) \right| = \cos t \le \cos 29^{\circ} < 0.875, so take M=0.875M = 0.875: cos290.8747520.8752(π180)2=0.4375×0.000304620.000133.\left| \cos 29^{\circ} - 0.874752 \right| \le \frac{0.875}{2}\left(\frac{\pi}{180}\right)^{2} = 0.4375 \times 0.00030462 \approx 0.000133 .

Check. The true value is cos29=0.8746197\cos 29^{\circ} = 0.8746197, so the actual error is 0.8747520.874620=0.0001320.874752 - 0.874620 = 0.000132, comfortably inside the bound, and indeed the estimate is too large.

Answer: cos2932+π3600.874752\cos 29^{\circ} \approx \dfrac{\sqrt{3}}{2} + \dfrac{\pi}{360} \approx 0.874752; an over-estimate, with error at most 1.34×1041.34 \times 10^{-4}.

The radius of a sphere is measured as 2121 cm with a possible error of at most 0.050.05 cm. Use differentials to estimate the maximum error in the calculated volume, and find the relative and percentage error.

Show hint

Differentiate V=43πr3V = \frac{4}{3}\pi r^{3} to get dVdV in terms of drdr. For the relative error, divide dVdV by VV before substituting numbers — most of the mess cancels.

Show answer

Model. The volume of a sphere of radius rr is V=43πr3,dVdr=4πr2,dV=4πr2dr.V = \frac{4}{3}\pi r^{3}, \qquad \frac{dV}{dr} = 4\pi r^{2}, \qquad dV = 4\pi r^{2} dr .

Maximum error in VV. With r=21r = 21 and dr0.05\left| dr \right| \le 0.05:

dV4π(21)2(0.05)=4π(441)(0.05)=88.2π277.09 cm3.\begin{aligned} \left| dV \right| &\le 4\pi (21)^{2}(0.05) \\ &= 4\pi (441)(0.05) \\ &= 88.2\pi \\ &\approx 277.09 \text{ cm}^{3}. \end{aligned}

Relative error. Divide the differential by the volume symbolically first: dVV=4πr2dr43πr3=3drr.\frac{dV}{V} = \frac{4\pi r^{2} dr}{\frac{4}{3}\pi r^{3}} = 3 \cdot \frac{dr}{r}. Hence dVV30.0521=0.1521=11400.0071429.\frac{\left| dV \right|}{V} \le 3 \cdot \frac{0.05}{21} = \frac{0.15}{21} = \frac{1}{140} \approx 0.0071429 .

Percentage error. 100×0.00714290.714%100 \times 0.0071429 \approx 0.714\%.

Independent check. V=43π(21)3=43π(9261)=12348π38792.4V = \frac{4}{3}\pi(21)^{3} = \frac{4}{3}\pi(9261) = 12348\pi \approx 38792.4 cm3^{3}, and 88.2π12348π=88.212348=0.0071429\frac{88.2\pi}{12348\pi} = \frac{88.2}{12348} = 0.0071429 which agrees with the shortcut.

Answer: maximum error in the volume 88.2π277\approx 88.2\pi \approx 277 cm3^{3}; relative error 11400.00714\dfrac{1}{140} \approx 0.00714; percentage error 0.71%\approx 0.71\%. (Note the general rule: a 0.238%0.238\% error in rr produces three times that, 0.714%0.714\%, in VV.)

(a) Show that if y=Cxny = Cx^{n} (CC and nn constants, C0C \ne 0, x0x \ne 0), the relative error in yy satisfies dyy=ndxx\dfrac{dy}{y} = n\,\dfrac{dx}{x}.

(b) The period of a simple pendulum is T=2πL/gT = 2\pi\sqrt{L/g}. If LL is measured with a percentage error of at most 1%1\%, estimate the resulting percentage error in TT. Illustrate with L=1L = 1 m, g=9.8g = 9.8 m/s2^{2}.

Show hint

For (a) compute dydy and divide by yy — every constant cancels. For (b), write TT as a constant times a power of LL and read off the exponent.

Show answer

(a) Differentiate y=Cxny = Cx^{n}: dydx=Cnxn1dy=Cnxn1dx.\frac{dy}{dx} = Cnx^{n-1} \qquad \Rightarrow \qquad dy = Cnx^{n-1} dx . Divide by y=Cxny = Cx^{n} (nonzero because C0C \ne 0 and x0x \ne 0): dyy=Cnxn1dxCxn=ndxx.\frac{dy}{y} = \frac{Cnx^{n-1} dx}{Cx^{n}} = n \cdot \frac{dx}{x}. So the relative error is multiplied by the exponent nn, and the percentage error in yy is about n\left| n \right| times that in xx. (Equivalently: lny=lnC+nlnx\ln \left| y \right| = \ln \left| C \right| + n\ln \left| x \right|, and differentiating gives dyy=ndxx\dfrac{dy}{y} = n\dfrac{dx}{x}.)

(b) Write the period as a constant times a power of LL: T=2πLg=2πgL1/2,T = 2\pi \sqrt{\frac{L}{g}} = \frac{2\pi}{\sqrt{g}} L^{1/2}, so part (a) applies with n=12n = \dfrac{1}{2}: dTT=12dLL.\frac{dT}{T} = \frac{1}{2} \cdot \frac{dL}{L}. A 1%1\% error in LL therefore gives dTT12(0.01)=0.005about 0.5% in T.\left| \frac{dT}{T} \right| \le \frac{1}{2}(0.01) = 0.005 \quad \Rightarrow \quad \text{about } 0.5\% \text{ in } T .

Illustration. With L=1L = 1 and g=9.8g = 9.8: T=2π19.8=2π(0.3194383)2.007090 s.T = 2\pi\sqrt{\frac{1}{9.8}} = 2\pi(0.3194383) \approx 2.007090 \text{ s}. Taking dL=0.01dL = 0.01 m (that is 1%1\%), dT=12dLLT=12(0.01)(2.007090)0.010035 s.dT = \frac{1}{2} \cdot \frac{dL}{L} \cdot T = \frac{1}{2}(0.01)(2.007090) \approx 0.010035 \text{ s}.

Independent check. Exactly, T(1.01)=2π1.019.8=2.0070901.01=2.007090(1.0049876)2.017100,T(1.01) = 2\pi\sqrt{\frac{1.01}{9.8}} = 2.007090\sqrt{1.01} = 2.007090(1.0049876) \approx 2.017100, so ΔT0.010010\Delta T \approx 0.010010 s, against the differential estimate 0.0100350.010035 s — agreement to three significant figures.

Answer: dyy=ndxx\dfrac{dy}{y} = n\dfrac{dx}{x}; a 1%1\% error in LL produces roughly a 0.5%0.5\% error in TT (about 0.0100.010 s when T2.007T \approx 2.007 s).

Use the linearisation of f(x)=exf(x) = e^{x} at a=0a=0 to estimate e0.1e^{0.1}. Then use Taylor's inequality to produce a rigorous bound on the error, and state whether the estimate is too big or too small.

Show hint

The linearisation is the familiar 1+x1+x. For the bound you need a number MM with f(t)M\left| f''(t) \right| \le M for all tt in [0,0.1][0, 0.1] — find one you can justify without already knowing e0.1e^{0.1}.

Show answer

Linearisation. f(x)=exf(x) = e^{x}, f(0)=1f(0) = 1, f(x)=exf'(x) = e^{x}, f(0)=1f'(0) = 1, so L(x)=1+xe0.1L(0.1)=1.1.L(x) = 1 + x \qquad \Rightarrow \qquad e^{0.1} \approx L(0.1) = 1.1 .

Direction. f(x)=ex>0f''(x) = e^{x} > 0 everywhere, so ff is concave up and the tangent line lies below the graph: 1.11.1 is an under-estimate, i.e. e0.1>1.1e^{0.1} > 1.1.

A usable bound MM. On [0,0.1][0, 0.1], f(t)=ete0.1\left| f''(t) \right| = e^{t} \le e^{0.1} (the exponential is increasing). To avoid circularity, bound e0.1e^{0.1} crudely: since e<3e < 3, e0.1<30.1,e^{0.1} < 3^{0.1}, and 30.1<1.23^{0.1} < 1.2 because 1.210=6.1917>31.2^{10} = 6.1917 > 3. So M=1.2M = 1.2 works.

Taylor's inequality. ff' is continuous on [0,0.1][0,0.1] and ff'' exists on (0,0.1)(0,0.1), so

e0.11.1M2(0.10)2=1.22(0.01)=0.006.\begin{aligned} \left| e^{0.1} - 1.1 \right| &\le \frac{M}{2}(0.1 - 0)^{2} \\ &= \frac{1.2}{2}(0.01) \\ &= 0.006 . \end{aligned}

Conclusion. Combining the two facts, 1.1<e0.11.1+0.006=1.106.1.1 < e^{0.1} \le 1.1 + 0.006 = 1.106 .

Independent check. The true value is e0.1=1.1051709e^{0.1} = 1.1051709, so the actual error is 0.00517090.0051709, which is indeed positive (under-estimate) and less than the bound 0.0060.006.

Answer: e0.11.1e^{0.1} \approx 1.1, an under-estimate, with error0.006\left| \text{error} \right| \le 0.006; in fact 1.1<e0.1<1.1061.1 < e^{0.1} < 1.106.

Estimate 263\sqrt[3]{26} using a linear approximation. State the estimate as an exact fraction, decide whether it is an over- or under-estimate, and give a rigorous error bound. Compare with the actual error.

Show hint

Base the approximation at the nearest perfect cube. For the bound you must dominate f\left| f'' \right| on the whole interval between the base point and 2626; a crude lower bound on t1/3t^{1/3} there is enough.

Show answer

Set up. f(x)=x1/3f(x) = x^{1/3}, base point a=27a = 27 (the nearest perfect cube).

f(27)=3,f(x)=13x2/3,f(27)=131272/3=139=127,L(x)=3+127(x27).\begin{aligned} f(27) &= 3, \\ f'(x) &= \frac{1}{3}x^{-2/3}, \qquad f'(27) = \frac{1}{3} \cdot \frac{1}{27^{2/3}} = \frac{1}{3 \cdot 9} = \frac{1}{27}, \\ L(x) &= 3 + \frac{1}{27}(x - 27). \end{aligned}

Estimate. With x=26x = 26, so xa=1x - a = -1: 263L(26)=3127=81127=80272.9629630.\sqrt[3]{26} \approx L(26) = 3 - \frac{1}{27} = \frac{81 - 1}{27} = \frac{80}{27} \approx 2.9629630 .

Direction. f(x)=29x5/3<0f''(x) = -\dfrac{2}{9}x^{-5/3} < 0 for x>0x > 0, so ff is concave down: the tangent lies above the graph and 8027\dfrac{80}{27} is an over-estimate.

Error bound. Let tt lie in [26,27][26, 27]. Then t1/3261/3>2.9t^{1/3} \ge 26^{1/3} > 2.9, because 2.93=24.389<262.9^{3} = 24.389 < 26. Hence t5/3=(t1/3)5>2.95=205.11149,t^{5/3} = \left(t^{1/3}\right)^{5} > 2.9^{5} = 205.11149, and therefore f(t)=29t5/3<291205.11149=0.0010834.\left| f''(t) \right| = \frac{2}{9}t^{-5/3} < \frac{2}{9} \cdot \frac{1}{205.11149} = 0.0010834 . Take M=0.0011M = 0.0011. Taylor's inequality gives 26380270.00112(2627)2=0.00055.\left| \sqrt[3]{26} - \frac{80}{27} \right| \le \frac{0.0011}{2}(26-27)^{2} = 0.00055 .

Conclusion. 2.96241302632.96296302.9624130 \le \sqrt[3]{26} \le 2.9629630.

Independent check. Cube the estimate: (8027)3=51200019683=26.0122949>26,\left(\frac{80}{27}\right)^{3} = \frac{512000}{19683} = 26.0122949 > 26, which re-confirms the over-estimate. The true value is 263=2.9624961\sqrt[3]{26} = 2.9624961, so the actual error is 2.96296302.9624961=0.0004669<0.00055.2.9629630 - 2.9624961 = 0.0004669 < 0.00055 .

Answer: 26380272.96296\sqrt[3]{26} \approx \dfrac{80}{27} \approx 2.96296, an over-estimate, with error at most 0.000550.00055 (actual error 0.000470.00047).

Suppose ff is twice differentiable on the whole real line with f(2)=5f(2) = 5, f(2)=3f'(2) = -3, and f(x)>0f''(x) > 0 for every xx.

(a) Use the linearisation at a=2a=2 to estimate f(1.9)f(1.9) and f(2.2)f(2.2).

(b) Is each estimate too large or too small? Explain why the answer is the same for both.

(c) If in addition f(x)4\left| f''(x) \right| \le 4 on [1.9,2.2][1.9, 2.2], give an interval that certainly contains f(2.2)f(2.2).

Show hint

Write down L(x)L(x) from the given data, then remember that the over/under question is decided by the sign of ff'' alone, not by which side of aa the point lies on.

Show answer

(a) Linearisation at a=2a = 2. L(x)=f(2)+f(2)(x2)=53(x2).L(x) = f(2) + f'(2)(x-2) = 5 - 3(x-2).

f(1.9)L(1.9)=53(1.92)=53(0.1)=5+0.3=5.3,f(2.2)L(2.2)=53(2.22)=53(0.2)=50.6=4.4.\begin{aligned} f(1.9) &\approx L(1.9) = 5 - 3(1.9 - 2) = 5 - 3(-0.1) = 5 + 0.3 = 5.3, \\ f(2.2) &\approx L(2.2) = 5 - 3(2.2 - 2) = 5 - 3(0.2) = 5 - 0.6 = 4.4 . \end{aligned}

(b) Since f(x)>0f''(x) > 0 for all xx, ff is concave up everywhere, so its graph lies above every tangent line (strictly, except at the point of tangency): f(x)>L(x)for all x2.f(x) > L(x) \qquad \text{for all } x \ne 2 . Therefore f(1.9)>5.3f(1.9) > 5.3 and f(2.2)>4.4f(2.2) > 4.4: both estimates are too small.

They agree in direction because the tangent line is below the curve on both sides of x=2x=2 — the relevant fact is the sign of ff'', not whether x>2x > 2 or x<2x < 2. (The sign of ff' is irrelevant here too: ff is decreasing at 22, yet both estimates are still under-estimates.)

(c) Taylor's inequality with a=2a = 2, x=2.2x = 2.2 and M=4M = 4: f(2.2)4.4M2(2.22)2=42(0.2)2=2(0.04)=0.08.\left| f(2.2) - 4.4 \right| \le \frac{M}{2}(2.2-2)^{2} = \frac{4}{2}(0.2)^{2} = 2(0.04) = 0.08 . So 4.32f(2.2)4.484.32 \le f(2.2) \le 4.48. Combining with part (b), which says f(2.2)>4.4f(2.2) > 4.4, the sharper conclusion is 4.4<f(2.2)4.48.4.4 < f(2.2) \le 4.48 .

Remark on accuracy. The same bound at x=1.9x = 1.9 gives f(1.9)5.32(0.1)2=0.02\left| f(1.9) - 5.3 \right| \le 2(0.1)^{2} = 0.02 — four times smaller, because the error bound is quadratic in the distance from aa and 1.91.9 is half as far from 22 as 2.22.2 is.

Answer: (a) f(1.9)5.3f(1.9) \approx 5.3, f(2.2)4.4f(2.2) \approx 4.4; (b) both are under-estimates, because f>0f'' > 0 puts the graph above the tangent line on both sides of 22; (c) 4.4<f(2.2)4.484.4 < f(2.2) \le 4.48.

The area of a circular disk is to be computed from a measurement of its radius, which is nominally 2424 cm. The area must be correct to within a percentage error of 0.5%0.5\%. Using differentials, how accurately must the radius be measured — as a percentage, and in centimetres?

Show hint

This is the error-propagation relation run backwards: express the relative error in AA in terms of the relative error in rr, then solve the inequality for dr\left| dr \right|.

Show answer

Relation between the relative errors. For A=πr2A = \pi r^{2}, dA=2πrdrdAA=2πrdrπr2=2drr.dA = 2\pi r dr \qquad \Rightarrow \qquad \frac{dA}{A} = \frac{2\pi r dr}{\pi r^{2}} = 2 \cdot \frac{dr}{r}. (The general power rule with n=2n=2.)

Impose the tolerance. We require dAA0.005(that is, 0.5%).\left| \frac{dA}{A} \right| \le 0.005 \quad \text{(that is, } 0.5\% \text{)}. Substituting, 2drr0.005drr0.0025.2\left| \frac{dr}{r} \right| \le 0.005 \qquad \Rightarrow \qquad \left| \frac{dr}{r} \right| \le 0.0025 .

So the radius must be measured to within 0.25%0.25\%half the allowed percentage error in the area, because the exponent is 22.

Convert to centimetres. With r=24r = 24 cm, dr0.0025×24=0.06 cm.\left| dr \right| \le 0.0025 \times 24 = 0.06 \text{ cm}.

Independent check (and a caution). An error of exactly +0.06+0.06 cm gives A(24.06)A(24)=(24.0624)2=(1.0025)2=1.00500625,\frac{A(24.06)}{A(24)} = \left(\frac{24.06}{24}\right)^{2} = (1.0025)^{2} = 1.00500625, a true percentage error of 0.5006%0.5006\% — a hair over 0.5%0.5\%, because the differential drops the (dr)2\left(dr\right)^{2} term. An error of 0.06-0.06 cm gives (0.9975)2=0.99500625(0.9975)^{2} = 0.99500625, i.e. 0.4994%0.4994\%. The linear answer is thus off by only about 0.00060.0006 percentage points (0.5006%0.5%=0.0006250.5006\% - 0.5\% = 0.000625 points); if the specification is a hard limit, tighten slightly (for example to dr0.0599\left| dr \right| \le 0.0599 cm).

Answer: the radius must be measured with a percentage error of at most about 0.25%0.25\%, i.e. dr0.06\left| dr \right| \le 0.06 cm.

Use differentials to estimate the amount of paint needed to apply a coat of paint 0.050.05 cm thick to a hemispherical dome of diameter 5050 m. Give the answer in cubic metres and in litres.

Show hint

The paint forms a thin shell on the curved surface. Model its volume as the change dVdV in the volume of a solid hemisphere when the radius increases by the paint thickness — and be ruthless about units.

Show answer

Model. A solid hemisphere of radius rr has volume V(r)=23πr3,V(r)=2πr2,dV=2πr2dr.V(r) = \frac{2}{3}\pi r^{3}, \qquad V'(r) = 2\pi r^{2}, \qquad dV = 2\pi r^{2} dr . The coat of paint is the thin shell added when rr increases by the paint thickness, so its volume is approximately dVdV. (Sanity check on the formula: 2πr22\pi r^{2} is exactly the curved surface area of the hemisphere, so dV=area×thicknessdV = \text{area} \times \text{thickness}, as it should be.)

Units. Work entirely in centimetres: r=502=25 m=2500 cm,dr=0.05 cm.r = \frac{50}{2} = 25 \text{ m} = 2500 \text{ cm}, \qquad dr = 0.05 \text{ cm}.

Compute.

dV=2π(2500)2(0.05)=2π(6,250,000)(0.05)=2π(312,500)=625,000π cm31,963,495 cm3.\begin{aligned} dV &= 2\pi (2500)^{2}(0.05) \\ &= 2\pi (6{,}250{,}000)(0.05) \\ &= 2\pi (312{,}500) \\ &= 625{,}000\pi \text{ cm}^{3} \\ &\approx 1{,}963{,}495 \text{ cm}^{3}. \end{aligned}

Convert. Since 1 m3=106 cm31 \text{ m}^{3} = 10^{6} \text{ cm}^{3} and 11 litre =1000= 1000 cm3^{3}: dV1.963 m31963 litres.dV \approx 1.963 \text{ m}^{3} \approx 1963 \text{ litres}.

Independent check. The exact shell volume is

V(2500.05)V(2500)=23π[(2500.05)325003]=23π[3(2500)2(0.05)+3(2500)(0.05)2+(0.05)3]=23π[937,500+18.75+0.000125]=23π(937,518.750125)625,012.5π1,963,535 cm3,\begin{aligned} V(2500.05) - V(2500) &= \frac{2}{3}\pi\left[(2500.05)^{3} - 2500^{3}\right] \\ &= \frac{2}{3}\pi\left[3(2500)^{2}(0.05) + 3(2500)(0.05)^{2} + (0.05)^{3}\right] \\ &= \frac{2}{3}\pi\left[937{,}500 + 18.75 + 0.000125\right] \\ &= \frac{2}{3}\pi \left(937{,}518.750125\right) \\ &\approx 625{,}012.5\pi \approx 1{,}963{,}535 \text{ cm}^{3}, \end{aligned}

so the differential estimate is low by about 3939 cm3^{3} — a relative error of roughly 2×1052 \times 10^{-5}, utterly negligible next to the uncertainty in how evenly anyone paints.

Answer: about 625,000π1.96×106625{,}000\pi \approx 1.96 \times 10^{6} cm3^{3}, i.e. roughly 1.961.96 m3^{3} (about 19631963 litres) of paint.