MATH1231 4,946 words·25 min read

Averages, Arc Length, Speed and Surface Area

In this chapter we point calculus at four geometric problems:

Every one of them is attacked in exactly the same way: chop the object into nn small pieces, approximate each piece by something we already know how to measure (a rectangle, a straight line segment, the frustum of a cone), add the approximations up to get a Riemann sum, and let n→∞n\to\infty. Basically, the whole chapter is the fundamental theorem of calculus being aimed at geometry; the only new content is deciding what the little pieces should be in each case.

5.1 The average value of a function#

Suppose a cable is suspended between two poles, so that each end of the cable is 33 metres above the ground and the lowest point of the cable is 22 metres above the ground. What is the average height of the cable above the ground? Clearly the answer is somewhere between 22 and 33 metres, but 'somewhere between' is not an answer. Recall from MATH1131 that a suspended cable is always the graph of a function of the form

f(x)=1ccosh⁡(cx)f(x)=\frac{1}{c}\cosh(cx)

over some interval [a,b][a,b], where cc is a constant depending on the tension in and the mass of the cable, and where the coordinate system has been chosen conveniently. So the question becomes: what is the average value of ff on [a,b][a,b]?

Averaging finitely many numbers is easy, but a function takes uncountably many values on [a,b][a,b], so we have to say what we even mean. Divide [a,b][a,b] into nn subintervals of equal length, sample the height of the graph in the kkth subinterval by choosing a point ckc_k in it, and take the ordinary arithmetic mean of the sampled heights:

an=1n(f(c1)+f(c2)+⋯+f(cn))=1n∑k=1nf(ck).a_n=\frac{1}{n}\big(f(c_1)+f(c_2)+\dots+f(c_n)\big)=\frac{1}{n}\sum_{k=1}^{n}f(c_k).

As nn increases, ana_n should get closer and closer to what we intuitively call the average height of the graph. Now comes the trick that turns this into an integral: multiply and divide by b−ab-a,

an=1n∑k=1nf(ck)=1n∑k=1nf(ck)⋅b−ab−a=1b−a∑k=1nf(ck) b−an.\begin{align*} a_n &= \frac{1}{n}\sum_{k=1}^{n}f(c_k) \\ &= \frac{1}{n}\sum_{k=1}^{n}f(c_k)\cdot\frac{b-a}{b-a} \\ &= \frac{1}{b-a}\sum_{k=1}^{n}f(c_k)\,\frac{b-a}{n}. \end{align*}

Notice what the last sum is: f(ck)f(c_k) is the height and b−an\frac{b-a}{n} is the width of the kkth rectangle in the usual picture, so it is a Riemann sum for ff on [a,b][a,b]. Hence, provided ff is Riemann integrable,

lim⁡n→∞an=1b−alim⁡n→∞∑k=1nf(ck) b−an=1b−a∫abf(x) dx.\lim_{n\to\infty}a_n=\frac{1}{b-a}\lim_{n\to\infty}\sum_{k=1}^{n}f(c_k)\,\frac{b-a}{n}=\frac{1}{b-a}\int_a^bf(x)\,dx.

This limit is what we take as the definition.

Note

Definition
Suppose that ff is integrable on a closed interval [a,b][a,b]. Then the average value fˉ\bar{f} of ff on [a,b][a,b] is defined by the formula

fˉ=1b−a∫abf(x) dx.\bar{f}=\frac{1}{b-a}\int_a^bf(x)\,dx.

Rearranging this formula gives a second way to read it: fˉ\bar{f} is the unique constant with

(b−a)fˉ=∫abf(x) dx.(b-a)\bar{f}=\int_a^bf(x)\,dx.

Geometrically, fˉ\bar{f} is the unique yy-value for which the rectangle of width b−ab-a and height fˉ\bar{f} has the same area as the region under the graph of ff. Basically, the average value is the height you would flatten the graph down to if you were allowed to redistribute the area but not create or destroy any of it.

Because the integral measures signed area, so does the average value. For example the average value of sin⁡\sin over [0,2π][0,2\pi] is 00, not because sin⁡\sin is small but because the hump below the axis exactly cancels the hump above it. Over half a period the answer is much more interesting.

Example. Find the average value of sin⁡\sin on [0,π][0,\pi].
By the definition,

fˉ=1π−0∫0πsin⁡x dx=1π[−cos⁡x]0π=1π(1−(−1))=2π.\begin{align*} \bar{f} &= \frac{1}{\pi-0}\int_0^\pi\sin x\,dx \\ &= \frac{1}{\pi}\Big[-\cos x\Big]_0^\pi \\ &= \frac{1}{\pi}\big(1-(-1)\big) \\ &= \frac{2}{\pi}. \end{align*}

Therefore the average value of sin⁡\sin on [0,π][0,\pi] is 2π≈0.6366\frac{2}{\pi}\approx0.6366. As a sanity check, the graph sits between 00 and 11 and spends most of its time up near the top, so an average a bit under 23\frac{2}{3} is exactly what you would expect; a common wrong guess is 12\frac{1}{2}, which is the average of the endpoint and peak values, not of the function.

Example. The cable illustrated at the start of this section is the curve

y=2cosh⁡(x/2),−a≤x≤a,y=2\cosh(x/2), \qquad -a\leq x\leq a,

where the xx-axis runs along the ground, the yy-axis passes through the vertex of the curve, and a=2cosh⁡−1(3/2)a=2\cosh^{-1}(3/2). Find, to the nearest centimetre, the average height of the cable above the ground.
Set f(x)=2cosh⁡(x/2)f(x)=2\cosh(x/2). The interval has length a−(−a)=2aa-(-a)=2a, so

fˉ=12a∫−aa2cosh⁡(x/2) dx=22a∫0a2cosh⁡(x/2) dx(since cosh⁡ is even)=1a[4sinh⁡(x/2)]0a=4asinh⁡(a/2)=2cosh⁡−1(3/2)sinh⁡(cosh⁡−1(3/2)).\begin{align*} \bar{f} &= \frac{1}{2a}\int_{-a}^{a}2\cosh(x/2)\,dx \\ &= \frac{2}{2a}\int_{0}^{a}2\cosh(x/2)\,dx \quad (\text{since } \cosh \text{ is even}) \\ &= \frac{1}{a}\Big[4\sinh(x/2)\Big]_0^a \\ &= \frac{4}{a}\sinh(a/2) \\ &= \frac{2}{\cosh^{-1}(3/2)}\sinh\big(\cosh^{-1}(3/2)\big). \end{align*}

The remaining sinh⁡\sinh of an inverse cosh⁡\cosh is unwound with the identity cosh⁡2t−sinh⁡2t=1\cosh^2t-\sinh^2t=1:

sinh⁡(cosh⁡−1(3/2))=cosh⁡2(cosh⁡−1(3/2))−1=94−1=52,\sinh\big(\cosh^{-1}(3/2)\big)=\sqrt{\cosh^2\big(\cosh^{-1}(3/2)\big)-1}=\sqrt{\frac{9}{4}-1}=\frac{\sqrt{5}}{2},

where the positive root is correct because cosh⁡−1(3/2)>0\cosh^{-1}(3/2)>0 and sinh⁡\sinh is positive there. Hence

fˉ=2cosh⁡−1(3/2)×52=5cosh⁡−1(3/2)≈2.32.\bar{f}=\frac{2}{\cosh^{-1}(3/2)}\times\frac{\sqrt5}{2}=\frac{\sqrt5}{\cosh^{-1}(3/2)}\approx2.32.

Therefore the average height of the cable above the ground is about 2.322.32 metres, which does indeed sit between 22 and 33 as predicted. Note that the fact that the xx-axis runs along the ground is special to this example; for a general hanging cable you have to work out where the ground is before you start.

Example. The air temperature T(t)T(t), measured in degrees Celsius tt hours after noon, is given by

T(t)=25+2t−t23.T(t)=25+2t-\frac{t^2}{3}.

Find the average temperature between noon and 5 p.m.
Noon is t=0t=0 and 5 p.m. is t=5t=5, so

Tˉ=15−0∫05(25+2t−t23)dt=15[25t+t2−t39]05=15(125+25−1259)=15⋅1350−1259=122545=2459.\begin{align*} \bar{T} &= \frac{1}{5-0}\int_0^5\left(25+2t-\frac{t^2}{3}\right)dt \\ &= \frac{1}{5}\left[25t+t^2-\frac{t^3}{9}\right]_0^5 \\ &= \frac{1}{5}\left(125+25-\frac{125}{9}\right) \\ &= \frac{1}{5}\cdot\frac{1350-125}{9} \\ &= \frac{1225}{45} \\ &= \frac{245}{9}. \end{align*}

Therefore the average temperature between noon and 5 p.m. is 2459=2729≈27.2∘\frac{245}{9}=27\frac{2}{9}\approx27.2^\circC. Do not average the endpoint values T(0)=25T(0)=25 and T(5)=2623T(5)=26\frac{2}{3} and call it a day; that gives 25.825.8, which is wrong, because the temperature does not change linearly. The endpoint average is only ever correct for a straight line.

We saw in MATH1131 that a continuous function on a closed interval attains its maximum and minimum values (the maximum-minimum theorem). The next result says such a function also attains its average value.

Note

Theorem (The mean value theorem for integrals)
Suppose that ff is continuous on [a,b][a,b]. Then there is a number cc in (a,b)(a,b) such that

∫abf(t) dt=f(c)(b−a).\int_a^bf(t)\,dt=f(c)(b-a).

Proof. Define F:[a,b]→RF:[a,b]\to\mathbb{R} by

F(x)=∫axf(t) dt.F(x)=\int_a^xf(t)\,dt.

By the fundamental theorem of calculus, FF is continuous on [a,b][a,b], differentiable on (a,b)(a,b), and F′(x)=f(x)F'(x)=f(x). So FF satisfies the hypotheses of the ordinary mean value theorem, which supplies a c∈(a,b)c\in(a,b) with

F(b)−F(a)b−a=F′(c).\frac{F(b)-F(a)}{b-a}=F'(c).

But F(a)=0F(a)=0, F(b)=∫abf(t) dtF(b)=\int_a^bf(t)\,dt and F′(c)=f(c)F'(c)=f(c), so this reads

1b−a∫abf(t) dt=f(c),\frac{1}{b-a}\int_a^bf(t)\,dt=f(c),

which is the claim (in the form fˉ=f(c)\bar{f}=f(c)). ■\blacksquare

Basically, the mean value theorem for integrals says the horizontal line y=fˉy=\bar{f} from the rectangle picture must actually cut the graph somewhere strictly inside the interval. Continuity is doing all the work here; drop it and the theorem is false. For instance the step function equal to 00 on [0,1)[0,1) and 11 at x=1x=1 has average value 00... but a function equal to 00 on [0,1)[0,1) and 22 on [1,2][1,2] has average 11 on [0,2][0,2] and never takes the value 11 anywhere.

Example. Verify the mean value theorem for integrals for f(x)=x2f(x)=x^2 on [0,3][0,3] by finding the point cc explicitly.
The function is continuous, and its average value is

fˉ=13∫03x2 dx=13[x33]03=13×9=3.\bar{f}=\frac{1}{3}\int_0^3x^2\,dx=\frac{1}{3}\left[\frac{x^3}{3}\right]_0^3=\frac{1}{3}\times9=3.

So we need c∈(0,3)c\in(0,3) with c2=3c^2=3, i.e. c=3≈1.732c=\sqrt3\approx1.732 (rejecting c=−3c=-\sqrt3, which is not in the interval). Therefore the theorem holds with c=3c=\sqrt3. Notice that cc is not the midpoint 1.51.5 of the interval; the parabola is steeper on the right, so the average value is achieved a little past halfway.

A more general version of this theorem appears in the tutorial problems for Chapter 4, and it is what is used to prove the Lagrange form of the remainder in Taylor's theorem — so this innocent-looking result is doing real work later in the course.

5.2 The arc length of a curve#

Suppose P0(x0,y0)P_0(x_0,y_0) and P1(x1,y1)P_1(x_1,y_1) are two points in R2\mathbb{R}^2. The distance between them is

dist⁡(P0,P1)=(x1−x0)2+(y1−y0)2,\operatorname{dist}(P_0,P_1)=\sqrt{(x_1-x_0)^2+(y_1-y_0)^2},

and if P0P_0 and P1P_1 are the endpoints of a straight line segment then we define the length of that segment to be this distance. A line segment is a very special curve; the point of this section is to measure the length of curves that bend. The idea is the obvious one — approximate the curve by a chain of little straight segments, whose lengths we can already compute, and take a limit as the segments get short.

5.2.1 An intuitive derivation of the arc length formula#

Suppose CC is a curve in R2\mathbb{R}^2 which can be expressed in parametric form as

C={(x(t),y(t))∈R2:a≤t≤b},C=\big\{(x(t),y(t))\in\mathbb{R}^2:a\leq t\leq b\big\},

where xx and yy are differentiable functions of tt. We also assume the parametrisation is chosen so that the moving point (x(t),y(t))(x(t),y(t)) does not retrace its steps, either forwards or backwards. For a≤s≤ba\leq s\leq b, let ℓ(s)\ell(s) denote the arc length of the partial curve

Cs={(x(t),y(t))∈R2:a≤t≤s},C_s=\big\{(x(t),y(t))\in\mathbb{R}^2:a\leq t\leq s\big\},

so that ℓ(a)=0\ell(a)=0 and ℓ(b)\ell(b) is the number we want.

Take a small piece of the curve and approximate it by a secant. Fix tt with a<t<ba<t<b, let hh be a small nonzero real number, and consider the two points P(x(t),y(t))P(x(t),y(t)) and Q(x(t+h),y(t+h))Q(x(t+h),y(t+h)). The length of the arc from PP to QQ is ℓ(t+h)−ℓ(t)\ell(t+h)-\ell(t), and it is approximately the length of the chord PQPQ, which the distance formula gives as

ℓ(t+h)−ℓ(t)≈[x(t+h)−x(t)]2+[y(t+h)−y(t)]2.\ell(t+h)-\ell(t)\approx\sqrt{[x(t+h)-x(t)]^2+[y(t+h)-y(t)]^2}.

Dividing both sides by hh (and pushing the hh inside the square root, where it becomes h2h^2),

ℓ(t+h)−ℓ(t)h≈(x(t+h)−x(t)h)2+(y(t+h)−y(t)h)2.\frac{\ell(t+h)-\ell(t)}{h}\approx\sqrt{\left(\frac{x(t+h)-x(t)}{h}\right)^2+\left(\frac{y(t+h)-y(t)}{h}\right)^2}.

The approximation improves as hh shrinks, so assuming that ℓ\ell is a differentiable function, letting h→0h\to0 gives

ℓ′(t)=[x′(t)]2+[y′(t)]2.\ell'(t)=\sqrt{[x'(t)]^2+[y'(t)]^2}.

By the fundamental theorem of calculus,

ℓ(s)=∫as[x′(t)]2+[y′(t)]2 dt+K\ell(s)=\int_a^s\sqrt{[x'(t)]^2+[y'(t)]^2}\,dt+K

for some constant KK, and putting s=as=a gives 0=ℓ(a)=0+K0=\ell(a)=0+K, so K=0K=0. Taking s=bs=b recovers the length of the whole curve.

This is a heuristic derivation, not a proof. A rigorous treatment would first define arc length (as the supremum of the lengths of inscribed polygons) and then evaluate some genuinely awkward Riemann sums; we do not do that here. Three hypotheses were quietly used and are worth naming, because they are exactly the ones that get violated in exam questions: xx and yy must be differentiable, the parametrisation must trace the curve once and only once, and ℓ\ell itself is assumed differentiable.

5.2.2 Arc length for a parametrised curve#

Collecting the derivation above: if a curve CC is given parametrically by

C={(x(t),y(t))∈R2:a≤t≤b},C=\big\{(x(t),y(t))\in\mathbb{R}^2:a\leq t\leq b\big\},

where xx and yy are differentiable functions of tt, then its arc length ℓ\ell is

ℓ=∫ab[x′(t)]2+[y′(t)]2 dt.\boxed{\ell=\int_a^b\sqrt{[x'(t)]^2+[y'(t)]^2}\,dt.}

Basically, [x′(t)]2+[y′(t)]2 dt\sqrt{[x'(t)]^2+[y'(t)]^2}\,dt is the length of an infinitesimal chord, obtained from Pythagoras applied to the horizontal displacement x′(t) dtx'(t)\,dt and the vertical displacement y′(t) dty'(t)\,dt; integrating adds all the little chords up.

Example. Find the arc length of one arch of the cycloid

x(t)=r(t−sin⁡t),y(t)=r(1−cos⁡t),0≤t≤2π.x(t)=r(t-\sin t), \qquad y(t)=r(1-\cos t), \qquad 0\leq t\leq2\pi.

(This is the curve traced by a point on the rim of a wheel of radius rr rolling along the xx-axis, and it is closely related to the 'curve of fastest descent' from MATH1131.) Differentiating,

x′(t)=r(1−cos⁡t),y′(t)=rsin⁡t,x'(t)=r(1-\cos t), \qquad y'(t)=r\sin t,

so

[x′(t)]2+[y′(t)]2=r2(1−2cos⁡t+cos⁡2t)+r2sin⁡2t=r2(1−2cos⁡t+1)(since cos⁡2t+sin⁡2t=1)=2r2(1−cos⁡t).\begin{align*} [x'(t)]^2+[y'(t)]^2 &= r^2(1-2\cos t+\cos^2t)+r^2\sin^2t \\ &= r^2(1-2\cos t+1) \quad (\text{since } \cos^2t+\sin^2t=1) \\ &= 2r^2(1-\cos t). \end{align*}

Before substituting this into the formula it is best to express 1−cos⁡t1-\cos t as a square, since we are about to take a square root. The double angle identity 1−cos⁡2θ=2sin⁡2θ1-\cos2\theta=2\sin^2\theta with θ=t/2\theta=t/2 gives

[x′(t)]2+[y′(t)]2=2r2⋅2sin⁡2(t/2)=4r2sin⁡2(t/2).[x'(t)]^2+[y'(t)]^2=2r^2\cdot2\sin^2(t/2)=4r^2\sin^2(t/2).

Now A2=∣A∣\sqrt{A^2}=|A|, not AA — but sin⁡(t/2)\sin(t/2) is positive for 0<t<2π0<t<2\pi, so here the naive square root causes no problems. Hence

ℓ=∫02π4r2sin⁡2(t/2) dt=2r∫02πsin⁡(t/2) dt=2r[−2cos⁡(t/2)]02π=2r(2+2)=8r.\begin{align*} \ell &= \int_0^{2\pi}\sqrt{4r^2\sin^2(t/2)}\,dt \\ &= 2r\int_0^{2\pi}\sin(t/2)\,dt \\ &= 2r\Big[-2\cos(t/2)\Big]_0^{2\pi} \\ &= 2r\big(2+2\big) \\ &= 8r. \end{align*}

Therefore one arch of the cycloid has length 8r8r. Notice how clean that is: the arch spans a horizontal distance of 2πr≈6.28r2\pi r\approx6.28r and the answer is exactly 8r8r, with no π\pi in sight at all.

Example. Use the parametrisation x(t)=rcos⁡tx(t)=r\cos t, y(t)=rsin⁡ty(t)=r\sin t with 0≤t≤2π0\leq t\leq2\pi to find the circumference of a circle of radius rr. What goes wrong with the parametrisation x(t)=rcos⁡2tx(t)=r\cos2t, y(t)=rsin⁡2ty(t)=r\sin2t on the same interval?
For the first parametrisation, x′(t)=−rsin⁡tx'(t)=-r\sin t and y′(t)=rcos⁡ty'(t)=r\cos t, so

ℓ=∫02πr2sin⁡2t+r2cos⁡2t dt=∫02πr dt=2πr,\begin{align*} \ell &= \int_0^{2\pi}\sqrt{r^2\sin^2t+r^2\cos^2t}\,dt \\ &= \int_0^{2\pi}r\,dt \\ &= 2\pi r, \end{align*}

which is the circumference we expect. For the second, x′(t)=−2rsin⁡2tx'(t)=-2r\sin2t and y′(t)=2rcos⁡2ty'(t)=2r\cos2t, and the same computation gives

∫02π4r2sin⁡22t+4r2cos⁡22t dt=∫02π2r dt=4πr,\int_0^{2\pi}\sqrt{4r^2\sin^22t+4r^2\cos^22t}\,dt=\int_0^{2\pi}2r\,dt=4\pi r,

which is twice the circumference. The reason is that as tt runs from 00 to 2π2\pi, the point (rcos⁡2t,rsin⁡2t)(r\cos2t,r\sin2t) goes around the circle twice; the integral faithfully reports the total distance travelled, which is not the length of the curve. To use this parametrisation correctly you would integrate from 00 to π\pi.

For a closed curve (or any curve at all), check that your parametrisation traverses it exactly once before integrating; the formula cannot tell the difference between a long curve and a short curve walked twice.

5.2.3 Arc length for the graph of a function#

Most curves you meet are given as y=f(x)y=f(x) rather than parametrically, but that is just a special case: parametrise the graph of ff on [a,b][a,b] by

x(t)=t,y(t)=f(t),a≤t≤b.x(t)=t, \qquad y(t)=f(t), \qquad a\leq t\leq b.

Then x′(t)=1x'(t)=1 and y′(t)=f′(t)y'(t)=f'(t), so the parametric formula immediately gives ℓ=∫ab1+[f′(t)]2 dt\ell=\int_a^b\sqrt{1+[f'(t)]^2}\,dt. Renaming the variable of integration xx, as one usually does, the arc length ℓ\ell of the graph of ff on [a,b][a,b] is

ℓ=∫ab1+[f′(x)]2 dx.\boxed{\ell=\int_a^b\sqrt{1+[f'(x)]^2}\,dx.}

Notice immediately that ℓ≥b−a\ell\geq b-a, since the integrand is at least 11: the curve is never shorter than the interval it spans, with equality only for a horizontal line. That is the cheapest sanity check available on any arc length answer, and it costs nothing to apply.

Example. Calculate the length of the arc y=x3/2y=x^{3/2} for 0≤x≤10\leq x\leq1.
Here f′(x)=32x1/2f'(x)=\frac{3}{2}x^{1/2}, so [f′(x)]2=9x4[f'(x)]^2=\frac{9x}{4} and

ℓ=∫011+9x4 dx=49⋅23[(1+9x4)3/2]01=827((134)3/2−1)=827(13138−1)=1313−827.\begin{align*} \ell &= \int_0^1\sqrt{1+\frac{9x}{4}}\,dx \\ &= \frac{4}{9}\cdot\frac{2}{3}\left[\left(1+\frac{9x}{4}\right)^{3/2}\right]_0^1 \\ &= \frac{8}{27}\left(\left(\frac{13}{4}\right)^{3/2}-1\right) \\ &= \frac{8}{27}\left(\frac{13\sqrt{13}}{8}-1\right) \\ &= \frac{13\sqrt{13}-8}{27}. \end{align*}

Therefore the arc length is 1313−827≈1.4397\frac{13\sqrt{13}-8}{27}\approx1.4397. The chord from (0,0)(0,0) to (1,1)(1,1) has length 2≈1.4142\sqrt2\approx1.4142, and the curve should be slightly longer than the chord; it is, by about 2%2\%.

Example. Find the arc length of y=13(x2+2)3/2y=\frac{1}{3}(x^2+2)^{3/2} for 0≤x≤30\leq x\leq3.
Differentiating with the chain rule, f′(x)=13⋅32(x2+2)1/2⋅2x=xx2+2f'(x)=\frac{1}{3}\cdot\frac{3}{2}(x^2+2)^{1/2}\cdot2x=x\sqrt{x^2+2}, so

1+[f′(x)]2=1+x2(x2+2)=x4+2x2+1=(x2+1)2.\begin{align*} 1+[f'(x)]^2 &= 1+x^2(x^2+2) \\ &= x^4+2x^2+1 \\ &= (x^2+1)^2. \end{align*}

The square root is now trivial (and x2+1>0x^2+1>0, so no absolute values are needed), giving

ℓ=∫03(x2+1) dx=[x33+x]03=9+3=12.\begin{align*} \ell &= \int_0^3(x^2+1)\,dx \\ &= \left[\frac{x^3}{3}+x\right]_0^3 \\ &= 9+3 \\ &= 12. \end{align*}

Therefore the arc length is exactly 1212 units.

That example is worth pausing on. Arc-length integrands in this course are always engineered so that 1+[f′(x)]21+[f'(x)]^2 collapses to a perfect square, or so that a Pythagorean or hyperbolic identity flattens it; if yours does not collapse, you have almost certainly made an algebra slip — go back and check f′f' before you start hunting for exotic substitutions. The reason is that ∫1+[f′(x)]2 dx\int\sqrt{1+[f'(x)]^2}\,dx is hopeless for almost every ff you could write down: even y=x2y=x^2 leads to ∫1+4x2 dx\int\sqrt{1+4x^2}\,dx, which needs a hyperbolic substitution from the substitution table, and y=sin⁡xy=\sin x leads to an elliptic integral that has no elementary antiderivative at all.

Example. Find the arc length of y=x24−ln⁡x2y=\frac{x^2}{4}-\frac{\ln x}{2} for 1≤x≤21\leq x\leq2.
Differentiating, f′(x)=x2−12xf'(x)=\frac{x}{2}-\frac{1}{2x}, so

1+[f′(x)]2=1+x24−12+14x2=x24+12+14x2=(x2+12x)2,\begin{align*} 1+[f'(x)]^2 &= 1+\frac{x^2}{4}-\frac{1}{2}+\frac{1}{4x^2} \\ &= \frac{x^2}{4}+\frac{1}{2}+\frac{1}{4x^2} \\ &= \left(\frac{x}{2}+\frac{1}{2x}\right)^2, \end{align*}

where the middle term flipped sign from −12-\frac{1}{2} to +12+\frac{1}{2} precisely because of the +1+1. Since x>0x>0 the bracket is positive, so

ℓ=∫12(x2+12x)dx=[x24+ln⁡x2]12=(1+ln⁡22)−14=34+ln⁡22.\begin{align*} \ell &= \int_1^2\left(\frac{x}{2}+\frac{1}{2x}\right)dx \\ &= \left[\frac{x^2}{4}+\frac{\ln x}{2}\right]_1^2 \\ &= \left(1+\frac{\ln2}{2}\right)-\frac{1}{4} \\ &= \frac{3}{4}+\frac{\ln2}{2}. \end{align*}

Therefore the arc length is 34+ln⁡22≈1.0966\frac{3}{4}+\frac{\ln2}{2}\approx1.0966, comfortably more than the span b−a=1b-a=1. This is the standard shape of a hand-made arc length question: f′f' is of the form u−14uu-\frac{1}{4u} or u2−12u\frac{u}{2}-\frac{1}{2u}, so that the cross term in [f′]2[f']^2 is exactly −1-1 and the +1+1 repairs the square.

Example. Find the arc length of the catenary

y=1acosh⁡(ax),x∈[−b,b].y=\frac{1}{a}\cosh(ax), \qquad x\in[-b,b].

A catenary is the shape of a hanging cable, which is why this function keeps reappearing. Here f′(x)=sinh⁡(ax)f'(x)=\sinh(ax), and the hyperbolic Pythagorean identity cosh⁡2t−sinh⁡2t=1\cosh^2t-\sinh^2t=1 does the collapsing for us:

ℓ=∫−bb1+sinh⁡2(ax) dx=∫−bbcosh⁡2(ax) dx=∫−bbcosh⁡(ax) dx(since cosh⁡ is always positive)=2∫0bcosh⁡(ax) dx(since cosh⁡ is even)=2[1asinh⁡(ax)]0b=2asinh⁡(ab)=1a(eab−e−ab).\begin{align*} \ell &= \int_{-b}^{b}\sqrt{1+\sinh^2(ax)}\,dx \\ &= \int_{-b}^{b}\sqrt{\cosh^2(ax)}\,dx \\ &= \int_{-b}^{b}\cosh(ax)\,dx \quad (\text{since } \cosh \text{ is always positive}) \\ &= 2\int_{0}^{b}\cosh(ax)\,dx \quad (\text{since } \cosh \text{ is even}) \\ &= 2\left[\frac{1}{a}\sinh(ax)\right]_0^b \\ &= \frac{2}{a}\sinh(ab) \\ &= \frac{1}{a}\left(e^{ab}-e^{-ab}\right). \end{align*}

Therefore the catenary has arc length 1a(eab−e−ab)\frac{1}{a}\left(e^{ab}-e^{-ab}\right) units. Notice how this is the only function whose arc length integrand simplifies without any factoring at all — the identity 1+sinh⁡2=cosh⁡21+\sinh^2=\cosh^2 is built into it, which is one reason the hyperbolic functions were invented in the first place.

Example. Hence find the length of the suspended cable of Section 5.1, namely y=2cosh⁡(x/2)y=2\cosh(x/2) on [−a,a][-a,a] where a=2cosh⁡−1(3/2)a=2\cosh^{-1}(3/2).
Comparing with the general catenary, 1ccosh⁡(cx)=2cosh⁡(x/2)\frac{1}{c}\cosh(cx)=2\cosh(x/2) forces c=12c=\frac{1}{2}, and the interval endpoint is b=ab=a. So the formula gives

ℓ=2csinh⁡(ca)=4sinh⁡(a2)=4sinh⁡(cosh⁡−1(3/2))=4×52=25,\begin{align*} \ell &= \frac{2}{c}\sinh(ca) \\ &= 4\sinh\left(\frac{a}{2}\right) \\ &= 4\sinh\big(\cosh^{-1}(3/2)\big) \\ &= 4\times\frac{\sqrt5}{2} \\ &= 2\sqrt5, \end{align*}

reusing the value sinh⁡(cosh⁡−1(3/2))=52\sinh(\cosh^{-1}(3/2))=\frac{\sqrt5}{2} computed earlier. Therefore the cable is 25≈4.472\sqrt5\approx4.47 metres long. As a sanity check, the poles are 2a=4cosh⁡−1(3/2)≈3.852a=4\cosh^{-1}(3/2)\approx3.85 metres apart, so the cable is longer than the gap it spans — as it must be, since it sags.

5.2.4 Arc length for a polar curve#

Suppose a curve is described in polar coordinates by r=f(θ)r=f(\theta) for θ0≤θ≤θ1\theta_0\leq\theta\leq\theta_1. Since

x=rcos⁡θ=f(θ)cos⁡θandy=rsin⁡θ=f(θ)sin⁡θ,x=r\cos\theta=f(\theta)\cos\theta \qquad\text{and}\qquad y=r\sin\theta=f(\theta)\sin\theta,

this is a parametrisation of the curve, with θ\theta as the parameter, so there is nothing new to derive — just some algebra to do. By the product rule,

x′(θ)=f′(θ)cos⁡θ−f(θ)sin⁡θ,y′(θ)=f′(θ)sin⁡θ+f(θ)cos⁡θ,x'(\theta)=f'(\theta)\cos\theta-f(\theta)\sin\theta, \qquad y'(\theta)=f'(\theta)\sin\theta+f(\theta)\cos\theta,

and hence

[x′(θ)]2+[y′(θ)]2=[f′(θ)]2cos⁡2θ−2f(θ)f′(θ)sin⁡θcos⁡θ+[f(θ)]2sin⁡2θ+[f′(θ)]2sin⁡2θ+2f(θ)f′(θ)sin⁡θcos⁡θ+[f(θ)]2cos⁡2θ=[f(θ)]2+[f′(θ)]2,\begin{align*} [x'(\theta)]^2+[y'(\theta)]^2 &= [f'(\theta)]^2\cos^2\theta-2f(\theta)f'(\theta)\sin\theta\cos\theta+[f(\theta)]^2\sin^2\theta \\ &\quad+[f'(\theta)]^2\sin^2\theta+2f(\theta)f'(\theta)\sin\theta\cos\theta+[f(\theta)]^2\cos^2\theta \\ &= [f(\theta)]^2+[f'(\theta)]^2, \end{align*}

where the cross terms cancelled and sin⁡2θ+cos⁡2θ=1\sin^2\theta+\cos^2\theta=1 mopped up the rest. Writing f(θ)f(\theta) as rr and f′(θ)f'(\theta) as drdθ\frac{dr}{d\theta}, the arc length ℓ\ell of a polar curve is

ℓ=∫θ0θ1r2+(drdθ)2 dθ.\boxed{\ell=\int_{\theta_0}^{\theta_1}\sqrt{r^2+\left(\frac{dr}{d\theta}\right)^2}\,d\theta.}

Basically, moving through an angle dθd\theta carries you r dθr\,d\theta around the circle and drdθdθ\frac{dr}{d\theta}d\theta outwards, and these two displacements are perpendicular, so Pythagoras gives the little bit of arc. The r2r^2 term is not optional and it is not (drdθ)2\left(\frac{dr}{d\theta}\right)^2 in disguise; forgetting it is the classic polar arc length mistake (it would say that a circle r=1r=1, on which drdθ=0\frac{dr}{d\theta}=0, has length zero).

Example. Find the length of the cardioid r=1+cos⁡θr=1+\cos\theta.
The full cardioid is traced once as θ\theta runs from 00 to 2π2\pi. Here drdθ=−sin⁡θ\frac{dr}{d\theta}=-\sin\theta, so

r2+(drdθ)2=(1+cos⁡θ)2+sin⁡2θ=1+2cos⁡θ+cos⁡2θ+sin⁡2θ=2+2cos⁡θ=4cos⁡2(θ/2),\begin{align*} r^2+\left(\frac{dr}{d\theta}\right)^2 &= (1+\cos\theta)^2+\sin^2\theta \\ &= 1+2\cos\theta+\cos^2\theta+\sin^2\theta \\ &= 2+2\cos\theta \\ &= 4\cos^2(\theta/2), \end{align*}

using the double angle identity 1+cos⁡2α=2cos⁡2α1+\cos2\alpha=2\cos^2\alpha with α=θ/2\alpha=\theta/2 in the last line. Now be careful with the square root: 4cos⁡2(θ/2)=2∣cos⁡(θ/2)∣\sqrt{4\cos^2(\theta/2)}=2|\cos(\theta/2)|, and cos⁡(θ/2)\cos(\theta/2) is positive on (0,π)(0,\pi) but negative on (π,2π)(\pi,2\pi). The cardioid is symmetric about the xx-axis, so the cleanest fix is to integrate over the top half and double:

ℓ=2∫0π2cos⁡(θ/2) dθ=4[2sin⁡(θ/2)]0π=8(1−0)=8.\begin{align*} \ell &= 2\int_0^{\pi}2\cos(\theta/2)\,d\theta \\ &= 4\Big[2\sin(\theta/2)\Big]_0^{\pi} \\ &= 8(1-0) \\ &= 8. \end{align*}

Therefore the cardioid has length 88 units. If you drop the absolute value and integrate 2cos⁡(θ/2)2\cos(\theta/2) straight from 00 to 2π2\pi you get 00, which is obviously nonsense; whenever a half-angle appears inside a square root, check the sign on the second half of the range.

Example. Find the length of the spiral r=eθr=e^{\theta} for 0≤θ≤2π0\leq\theta\leq2\pi.
Here drdθ=eθ=r\frac{dr}{d\theta}=e^{\theta}=r, so the integrand is e2θ+e2θ=2 eθ\sqrt{e^{2\theta}+e^{2\theta}}=\sqrt2\,e^{\theta} and

ℓ=2∫02πeθ dθ=2[eθ]02π=2(e2π−1).\begin{align*} \ell &= \sqrt2\int_0^{2\pi}e^{\theta}\,d\theta \\ &= \sqrt2\Big[e^{\theta}\Big]_0^{2\pi} \\ &= \sqrt2\left(e^{2\pi}-1\right). \end{align*}

Therefore the length is 2(e2π−1)≈755.9\sqrt2\left(e^{2\pi}-1\right)\approx755.9 units — the exponential spiral grows fast enough that a single turn is already enormous.

Example. The spiral r=e−θ/10r=e^{-\theta/10}, where θ≥0\theta\geq0, winds infinitely many times around the origin. Is its total arc length finite? Explain.
Compute the arc length of the segment 0≤θ≤θ10\leq\theta\leq\theta_1 first, and take a limit afterwards. Here drdθ=−110e−θ/10\frac{dr}{d\theta}=-\frac{1}{10}e^{-\theta/10}, so

ℓ(θ1)=∫0θ1(e−θ/10)2+(−110e−θ/10)2 dθ=∫0θ1e−θ/101+1100 dθ=10110∫0θ1e−θ/10 dθ=10110[−10e−θ/10]0θ1=101(1−e−θ1/10).\begin{align*} \ell(\theta_1) &= \int_0^{\theta_1}\sqrt{\left(e^{-\theta/10}\right)^2+\left(-\frac{1}{10}e^{-\theta/10}\right)^2}\,d\theta \\ &= \int_0^{\theta_1}e^{-\theta/10}\sqrt{1+\frac{1}{100}}\,d\theta \\ &= \frac{\sqrt{101}}{10}\int_0^{\theta_1}e^{-\theta/10}\,d\theta \\ &= \frac{\sqrt{101}}{10}\Big[-10e^{-\theta/10}\Big]_0^{\theta_1} \\ &= \sqrt{101}\left(1-e^{-\theta_1/10}\right). \end{align*}

As θ1→∞\theta_1\to\infty we have e−θ1/10→0e^{-\theta_1/10}\to0, so ℓ(θ1)→101\ell(\theta_1)\to\sqrt{101}. Therefore the total arc length is finite and equals 101≈10.05\sqrt{101}\approx10.05 units. This is a curve of infinite extent (in the sense that it never stops winding) with finite length; the windings shrink geometrically, and a geometric series converges. Keep this example in mind — the very last example of the chapter is its evil twin.

5.3 The speed of a moving particle#

In Chapter 4 of the MATH1131 calculus notes we discussed the speed of a particle moving along a straight line. Now we can handle a particle moving along any curve in the plane, because we can measure the length of that curve.

Suppose a particle PP moves in the plane with position (x(t),y(t))(x(t),y(t)) at time tt. The distance s(t)s(t) travelled from time zero to time tt is exactly the arc length of the path traversed in that time,

s(t)=∫0t[x′(u)]2+[y′(u)]2 du.s(t)=\int_0^t\sqrt{[x'(u)]^2+[y'(u)]^2}\,du.

By definition the speed of PP is the rate of change of distance travelled with respect to time, so by the fundamental theorem of calculus the speed v(t)v(t) is

v(t)=s′(t)=[x′(t)]2+[y′(t)]2.\boxed{v(t)=s'(t)=\sqrt{[x'(t)]^2+[y'(t)]^2}.}

If we bundle the position into a vector r(t)=(x(t),y(t))\mathbf{r}(t)=(x(t),y(t)), then r′(t)=(x′(t),y′(t))\mathbf{r}'(t)=(x'(t),y'(t)) is the velocity vector and the formula just says v(t)=∣r′(t)∣v(t)=|\mathbf{r}'(t)|: speed is the length of the velocity. Basically, arc length and speed are the same computation read in two directions — the arc length integrand is the speed, and distance travelled is the integral of speed:

dsdt=v(t)ands=∫v(t) dt.\frac{ds}{dt}=v(t) \qquad\text{and}\qquad s=\int v(t)\,dt.

Example. A stone is thrown horizontally from the deck of the Sydney Harbour Bridge at 2020 metres per second. Its position tt seconds after the throw is

x(t)=20t,y(t)=50−5t2,0≤t≤10,x(t)=20t, \qquad y(t)=50-5t^2, \qquad 0\leq t\leq\sqrt{10},

where y(t)y(t) is the height above the water. Find the speed of the stone an instant before it hits the water.
Differentiating, x′(t)=20x'(t)=20 and y′(t)=−10ty'(t)=-10t, so the speed at time tt is

v(t)=202+100t2=400+100t2.v(t)=\sqrt{20^2+100t^2}=\sqrt{400+100t^2}.

The stone hits the water when y(t)=0y(t)=0, that is when 5t2=505t^2=50, i.e. t=10t=\sqrt{10}. Hence

lim⁡t→(10)−v(t)=400+100×10=1400=1014≈37.42.\lim_{t\to(\sqrt{10})^-}v(t)=\sqrt{400+100\times10}=\sqrt{1400}=10\sqrt{14}\approx37.42.

Therefore the stone is travelling at about 37.4237.42 metres per second just before impact. Notice that the horizontal component of the velocity is a constant 2020 throughout — gravity only ever changes the vertical component — so the speed never drops below 2020.

Example. A projectile is fired from an elevated cannon. Its horizontal distance xx (in metres) from the cannon and height yy (in metres) above the ground, exactly tt seconds after firing, are

x(t)=40t,y(t)=−5t2+40t+45,0≤t≤t1,x(t)=40t, \qquad y(t)=-5t^2+40t+45, \qquad 0\leq t\leq t_1,

where t1t_1 is the time of impact.
(a) Find t1t_1.
(b) Find the speed of the projectile immediately prior to impact.
(c) What was the average height of the projectile above the ground during this period?

(a) Impact happens when y=0y=0:

−5t2+40t+45=0t2−8t−9=0(t−9)(t+1)=0,\begin{align*} -5t^2+40t+45 &= 0 \\ t^2-8t-9 &= 0 \\ (t-9)(t+1) &= 0, \end{align*}

so t=9t=9 or t=−1t=-1; discarding the negative time, t1=9t_1=9 seconds.

(b) Here x′(t)=40x'(t)=40 and y′(t)=−10t+40y'(t)=-10t+40, so at t=9t=9 we have y′(9)=−50y'(9)=-50 and

v(9)=402+(−50)2=1600+2500=4100=1041≈64.03.v(9)=\sqrt{40^2+(-50)^2}=\sqrt{1600+2500}=\sqrt{4100}=10\sqrt{41}\approx64.03.

Therefore the projectile is travelling at about 64.0364.03 metres per second immediately prior to impact.

(c) This is a job for the average value formula applied to yy on [0,9][0,9]:

yˉ=19∫09(−5t2+40t+45)dt=19[−5t33+20t2+45t]09=19(−1215+1620+405)=8109=90.\begin{align*} \bar{y} &= \frac{1}{9}\int_0^9\left(-5t^2+40t+45\right)dt \\ &= \frac{1}{9}\left[-\frac{5t^3}{3}+20t^2+45t\right]_0^9 \\ &= \frac{1}{9}\left(-1215+1620+405\right) \\ &= \frac{810}{9} \\ &= 90. \end{align*}

Therefore the average height of the projectile was 9090 metres. Note carefully that this is the average height with respect to time, which is what the phrasing asks for; averaging with respect to horizontal distance xx would be a different integral (though here it happens to agree, since xx is a linear function of tt).
For completeness, the distance actually flown is ∫091600+(40−10t)2 dt≈427.53\int_0^9\sqrt{1600+(40-10t)^2}\,dt\approx427.53 metres, which needs the substitution u=40−10tu=40-10t followed by a hyperbolic substitution from the substitution table; it is marked as extension material precisely because the algebra is grim, and it is a good illustration that a physically natural arc length usually does not simplify.

Since speed is the derivative of distance travelled, the two quantities determine each other, and it is natural to ask for a parametrisation in which the speed is always 11. Such a parametrisation is said to be by arc length: if v(t)=1v(t)=1 for all tt then s(t)=∫0t1 du=ts(t)=\int_0^t1\,du=t, so the parameter is the distance travelled along the curve.

Example. A particle moves on the circle x(t)=3cos⁡2tx(t)=3\cos2t, y(t)=3sin⁡2ty(t)=3\sin2t. Find its speed and the distance travelled by time tt, and hence reparametrise the circle by arc length.
Differentiating, x′(t)=−6sin⁡2tx'(t)=-6\sin2t and y′(t)=6cos⁡2ty'(t)=6\cos2t, so

v(t)=36sin⁡22t+36cos⁡22t=6,v(t)=\sqrt{36\sin^22t+36\cos^22t}=6,

a constant. Hence s(t)=∫0t6 du=6ts(t)=\int_0^t6\,du=6t, i.e. t=s6t=\frac{s}{6}. Substituting this back into the position,

x=3cos⁡(s3),y=3sin⁡(s3),x=3\cos\left(\frac{s}{3}\right), \qquad y=3\sin\left(\frac{s}{3}\right),

and it is easily checked that this new parametrisation has speed sin⁡2(s/3)+cos⁡2(s/3)=1\sqrt{\sin^2(s/3)+\cos^2(s/3)}=1. Therefore the arc length parametrisation of the circle of radius 33 is (3cos⁡s3,3sin⁡s3)\left(3\cos\frac{s}{3},3\sin\frac{s}{3}\right), and running ss from 00 to 6π6\pi traverses it once — recovering the circumference 2π×3=6π2\pi\times3=6\pi for free. In general the arc-length parametrisation of a circle of radius RR is (Rcos⁡sR,Rsin⁡sR)\left(R\cos\frac{s}{R},R\sin\frac{s}{R}\right); the 1R\frac{1}{R} inside is exactly what makes the speed 11 rather than RR.

Example. The position of a particle PP at time t≥0t\geq0 is

x(t)=cos⁡(π2(cos⁡πt−1)),y(t)=−sin⁡(π2(cos⁡πt−1)).x(t)=\cos\left(\frac{\pi}{2}(\cos\pi t-1)\right), \qquad y(t)=-\sin\left(\frac{\pi}{2}(\cos\pi t-1)\right).

(a) Find a formula for the speed v(t)v(t).
(b) What curve does the trajectory trace out, and what is its length?
(c) How far does the particle travel during the time interval [0,3][0,3]?

(a) Write ϕ(t)=π2(cos⁡πt−1)\phi(t)=\frac{\pi}{2}(\cos\pi t-1), so that x=cos⁡ϕx=\cos\phi and y=−sin⁡ϕy=-\sin\phi. By the chain rule,

x′(t)=−sin⁡ϕ⋅ϕ′(t),y′(t)=−cos⁡ϕ⋅ϕ′(t),x'(t)=-\sin\phi\cdot\phi'(t), \qquad y'(t)=-\cos\phi\cdot\phi'(t),

and hence

v(t)=ϕ′(t)2sin⁡2ϕ+ϕ′(t)2cos⁡2ϕ=∣ϕ′(t)∣=∣π2×(−πsin⁡πt)∣=π22∣sin⁡πt∣.\begin{align*} v(t) &= \sqrt{\phi'(t)^2\sin^2\phi+\phi'(t)^2\cos^2\phi} \\ &= |\phi'(t)| \\ &= \left|\frac{\pi}{2}\times(-\pi\sin\pi t)\right| \\ &= \frac{\pi^2}{2}|\sin\pi t|. \end{align*}

Therefore v(t)=π22∣sin⁡πt∣v(t)=\frac{\pi^2}{2}|\sin\pi t|; in particular the particle is momentarily at rest whenever tt is an integer, and fastest when tt is a half-integer.

(b) Since cos⁡πt−1\cos\pi t-1 takes every value in [−2,0][-2,0], the angle ϕ\phi ranges over [−π,0][-\pi,0], and (cos⁡ϕ,−sin⁡ϕ)(\cos\phi,-\sin\phi) then sweeps out the upper half of the unit circle. So the trajectory is the semicircle of centre (0,0)(0,0) and radius 11 lying in the upper half-plane, whose length is half the circumference, namely π\pi.

(c) Distance travelled is the integral of speed, and ∣sin⁡πt∣|\sin\pi t| has period 11 with ∫01∣sin⁡πt∣ dt=[−cos⁡πtπ]01=2π\int_0^1|\sin\pi t|\,dt=\left[-\frac{\cos\pi t}{\pi}\right]_0^1=\frac{2}{\pi}. Hence

∫03v(t) dt=π22×3×2π=3π.\int_0^3v(t)\,dt=\frac{\pi^2}{2}\times3\times\frac{2}{\pi}=3\pi.

Therefore the particle travels 3π3\pi units in the first three seconds. Compare (b) and (c): the curve is only π\pi units long, but the particle covers 3π3\pi units, because it sweeps back and forth along the same semicircle three times. This is the same trap as the doubly-traversed circle of Section 5.2.2, seen from the other side: arc length is a property of the curve, distance travelled is a property of the motion, and only a parametrisation which does not retrace makes them equal.

5.4 Surface area#

Finding the surface area of a general surface in R3\mathbb{R}^3 is genuinely hard, and the machinery needed for it (the partial derivatives of Functions of several variables, together with double integrals) is second year material. In this section we restrict to surfaces obtained by rotating a plane curve about one of the axes, which is exactly the case where a one-variable integral suffices.

Everything rests on one piece of geometry. The frustum of a right circular cone is what you get by slicing the top off a cone parallel to its base.

Note

Fact
Given a frustum of slant height ss and radii rr and RR, the area AA of its curved surface is given by

A=π(r+R)s.A=\pi(r+R)s.

This is proved by cutting the cone open and flattening it into a sector of a circle; it is an exercise in the tutorial problems, and it needs nothing beyond the formula for the area of a sector. Basically, π(r+R)s=2π(r+R2)s\pi(r+R)s=2\pi\left(\frac{r+R}{2}\right)s, i.e. (circumference at the average radius) ×\times (slant height) — the frustum is a rectangle of width 2πrˉ2\pi\bar{r} and height ss, rolled up.

5.4.1 An heuristic derivation for the surface area of a surface of revolution#

Suppose a curve CC has parametrisation

C={(x(t),y(t))∈R2:a≤t≤b}.C=\big\{(x(t),y(t))\in\mathbb{R}^2:a\leq t\leq b\big\}.

We assume that the curve lies in the upper half-plane (more precisely y(t)≥0y(t)\geq0, meeting the xx-axis at only finitely many points if at all), and that CC is simple: if (x(t0),y(t0))=(x(t1),y(t1))(x(t_0),y(t_0))=(x(t_1),y(t_1)) then t0=t1t_0=t_1, i.e. the curve does not cross itself. Rotating CC about the xx-axis produces a surface of revolution, and we want its area.

The derivation runs exactly parallel to the arc length one. Let A(s)A(s) denote the area of the surface formed by rotating the partial curve {(x(t),y(t)):a≤t≤s}\{(x(t),y(t)):a\leq t\leq s\}, and assume AA is differentiable. Fix t∈(a,b)t\in(a,b), let hh be small and nonzero, and consider P(x(t),y(t))P(x(t),y(t)) and Q(x(t+h),y(t+h))Q(x(t+h),y(t+h)). Then A(t+h)−A(t)A(t+h)-A(t) is the area swept out by the arc from PP to QQ, which for small hh is approximately the area swept by the chord PQPQ — and a rotated chord is precisely a frustum, of radii y(t)y(t) and y(t+h)y(t+h) and slant height equal to the length of PQPQ. So the frustum formula gives

A(t+h)−A(t)≈π(y(t+h)+y(t))[x(t+h)−x(t)]2+[y(t+h)−y(t)]2.A(t+h)-A(t)\approx\pi\big(y(t+h)+y(t)\big)\sqrt{[x(t+h)-x(t)]^2+[y(t+h)-y(t)]^2}.

Dividing by hh,

A(t+h)−A(t)h≈π(y(t+h)+y(t))(x(t+h)−x(t)h)2+(y(t+h)−y(t)h)2.\frac{A(t+h)-A(t)}{h}\approx\pi\big(y(t+h)+y(t)\big)\sqrt{\left(\frac{x(t+h)-x(t)}{h}\right)^2+\left(\frac{y(t+h)-y(t)}{h}\right)^2}.

Since yy is differentiable at tt it is continuous there, so y(t+h)→y(t)y(t+h)\to y(t) as h→0h\to0, and the limit is

A′(t)=π(y(t)+y(t))[x′(t)]2+[y′(t)]2=2πy(t)[x′(t)]2+[y′(t)]2.A'(t)=\pi\big(y(t)+y(t)\big)\sqrt{[x'(t)]^2+[y'(t)]^2}=2\pi y(t)\sqrt{[x'(t)]^2+[y'(t)]^2}.

Applying the fundamental theorem of calculus and using A(a)=0A(a)=0 to kill the constant of integration exactly as before,

A(s)=∫as2πy(t)[x′(t)]2+[y′(t)]2 dt,A(s)=\int_a^s2\pi y(t)\sqrt{[x'(t)]^2+[y'(t)]^2}\,dt,

and putting s=bs=b gives the area of the whole surface.

5.4.2 Surface area formulae and examples#

Assume throughout that CC lies in the upper half-plane and is simple, and that all the derivatives written down exist. Rotating about the xx-axis, the area AA of the surface of revolution is

A=∫ab2πy(t)[x′(t)]2+[y′(t)]2 dt\boxed{A=\int_a^b2\pi y(t)\sqrt{[x'(t)]^2+[y'(t)]^2}\,dt}

if CC is given parametrically by (x(t),y(t))(x(t),y(t)) with a≤t≤ba\leq t\leq b;

A=∫ab2πf(x)1+[f′(x)]2 dx\boxed{A=\int_a^b2\pi f(x)\sqrt{1+[f'(x)]^2}\,dx}

if CC is the graph y=f(x)y=f(x) on [a,b][a,b]; and

A=∫θ0θ12πrsin⁡θr2+(drdθ)2 dθ\boxed{A=\int_{\theta_0}^{\theta_1}2\pi r\sin\theta\sqrt{r^2+\left(\frac{dr}{d\theta}\right)^2}\,d\theta}

if CC is the polar curve r=f(θ)r=f(\theta) with θ0≤θ≤θ1\theta_0\leq\theta\leq\theta_1. Only the first of these needs deriving; the other two follow from it via the parametrisations x(t)=tx(t)=t, y(t)=f(t)y(t)=f(t) and x(θ)=f(θ)cos⁡θx(\theta)=f(\theta)\cos\theta, y(θ)=f(θ)sin⁡θy(\theta)=f(\theta)\sin\theta used in Sections 5.2.3 and 5.2.4. Rotating about the yy-axis instead, the parametric formula becomes

A=∫ab2πx(t)[x′(t)]2+[y′(t)]2 dt,\boxed{A=\int_a^b2\pi x(t)\sqrt{[x'(t)]^2+[y'(t)]^2}\,dt,}

with the graph version A=∫ab2πx1+[f′(x)]2 dxA=\int_a^b2\pi x\sqrt{1+[f'(x)]^2}\,dx obtained the same way.

Look at what changed between the last two boxes: the square root is identical in every formula — it is always the arc length element, i.e. the speed — and the only thing that varies is the factor in front. That factor is 2π×(distance from the point on the curve to the axis of rotation)2\pi\times(\text{distance from the point on the curve to the axis of rotation}): it is yy when you rotate about the xx-axis and xx when you rotate about the yy-axis. Writing 2πy2\pi y out of habit while rotating about the yy-axis is the single most common error in this entire topic. The reason is geometric rather than algebraic: the point traces a circle whose radius is its distance from the axis, and 2π×2\pi\times radius is the circumference of that circle.

Finally, note that these formulae give the area of the surface of revolution only. If the question asks for the surface area of a solid, you must also add the areas of any circular caps at the ends.

Example. Find the surface area of a sphere of radius rr, by rotating the semicircle y=r2−x2y=\sqrt{r^2-x^2}, −r≤x≤r-r\leq x\leq r, about the xx-axis.
Differentiating with the chain rule,

f′(x)=−2x2r2−x2=−xr2−x2,f'(x)=\frac{-2x}{2\sqrt{r^2-x^2}}=\frac{-x}{\sqrt{r^2-x^2}},

so that

1+[f′(x)]2=1+x2r2−x2=r2−x2+x2r2−x2=r2r2−x2.\begin{align*} 1+[f'(x)]^2 &= 1+\frac{x^2}{r^2-x^2} \\ &= \frac{r^2-x^2+x^2}{r^2-x^2} \\ &= \frac{r^2}{r^2-x^2}. \end{align*}

Hence 1+[f′(x)]2=rr2−x2\sqrt{1+[f'(x)]^2}=\frac{r}{\sqrt{r^2-x^2}}, and the crucial cancellation happens:

A=∫−rr2πr2−x2⋅rr2−x2 dx=∫−rr2πr dx=2πr(r−(−r))=4πr2.\begin{align*} A &= \int_{-r}^{r}2\pi\sqrt{r^2-x^2}\cdot\frac{r}{\sqrt{r^2-x^2}}\,dx \\ &= \int_{-r}^{r}2\pi r\,dx \\ &= 2\pi r\big(r-(-r)\big) \\ &= 4\pi r^2. \end{align*}

Therefore the surface area of a sphere of radius rr is 4πr24\pi r^2, as it should be. The same answer falls out of the parametric formula with x(t)=rcos⁡tx(t)=r\cos t, y(t)=rsin⁡ty(t)=r\sin t, 0≤t≤π0\leq t\leq\pi: there [x′]2+[y′]2=r\sqrt{[x']^2+[y']^2}=r and A=∫0π2πrsin⁡t⋅r dt=2πr2[−cos⁡t]0π=4πr2A=\int_0^\pi2\pi r\sin t\cdot r\,dt=2\pi r^2\big[-\cos t\big]_0^\pi=4\pi r^2. Note the limits 0≤t≤π0\leq t\leq\pi, not 0≤t≤2π0\leq t\leq2\pi: the lower semicircle would sweep out the same sphere a second time and double the answer.

The constant integrand 2πr2\pi r in that calculation is worth staring at, because it says something remarkable. The area of the piece of the sphere lying between the planes x=bx=b and x=cx=c (with −r≤b<c≤r-r\leq b<c\leq r) is

∫bc2πr dx=2πr(c−b),\int_b^c2\pi r\,dx=2\pi r(c-b),

which depends only on the width of the slab and not on where it sits. So a thin slice cut near the equator and an equally thin slice cut near the pole have exactly the same area: near the pole the circles are smaller, but the surface is more steeply inclined, and the two effects cancel exactly. This is sometimes called Archimedes' hat-box theorem, and it is the reason a cylindrical map projection of the globe preserves area.

Example. A solid SS is formed by rotating the curve y=2−xy=\sqrt{2-x}, x∈[0,2]x\in[0,2], about the xx-axis. Find the surface area of the solid, making sure that every face is accounted for.
The solid has two faces: the truncated paraboloid, and the flat circular cap at x=0x=0 where the curve has height 2\sqrt2. For the paraboloid, with f(x)=(2−x)1/2f(x)=(2-x)^{1/2} we get f′(x)=−12(2−x)−1/2f'(x)=-\frac{1}{2}(2-x)^{-1/2}, so

A1=∫022π2−x(1+14(2−x))1/2dx=∫022π((2−x)(1+14(2−x)))1/2dx=∫022π(2−x+14)1/2dx=∫02π9−4x dx=π[(9−4x)3/2−6]02=π(−16+276)=13π3,\begin{align*} A_1 &= \int_0^22\pi\sqrt{2-x}\left(1+\frac{1}{4(2-x)}\right)^{1/2}dx \\ &= \int_0^22\pi\left((2-x)\left(1+\frac{1}{4(2-x)}\right)\right)^{1/2}dx \\ &= \int_0^22\pi\left(2-x+\frac{1}{4}\right)^{1/2}dx \\ &= \int_0^2\pi\sqrt{9-4x}\,dx \\ &= \pi\left[\frac{(9-4x)^{3/2}}{-6}\right]_0^2 \\ &= \pi\left(-\frac{1}{6}+\frac{27}{6}\right) \\ &= \frac{13\pi}{3}, \end{align*}

where in the fourth line the 2−x+14=9−4x4=129−4x\sqrt{2-x+\frac14}=\sqrt{\frac{9-4x}{4}}=\frac{1}{2}\sqrt{9-4x} absorbed the factor of 22. Notice how pulling the outside 2−x\sqrt{2-x} inside the root is what made the integrand elementary; that manoeuvre is standard whenever ff and 1f\frac{1}{f} both appear. The cap is a disc of radius 2\sqrt2, so

A2=π(2)2=2π.A_2=\pi\left(\sqrt2\right)^2=2\pi.

Therefore the total surface area of SS is

A1+A2=13π3+2π=19π3≈19.9A_1+A_2=\frac{13\pi}{3}+2\pi=\frac{19\pi}{3}\approx19.9

square units. Read these questions carefully: "the area of the surface of revolution" and "the surface area of the solid" are different numbers, and the cap is worth marks.

Example. Find the area of the surface generated when y=x2y=x^2, 0≤x≤10\leq x\leq1, is rotated about (a) the yy-axis and (b) the xx-axis.
Both parts share the same square root: f′(x)=2xf'(x)=2x, so 1+[f′(x)]2=1+4x2\sqrt{1+[f'(x)]^2}=\sqrt{1+4x^2}.

(a) Rotating about the yy-axis, the distance to the axis is xx, so with the substitution u=1+4x2u=1+4x^2, du=8x dxdu=8x\,dx (and u=1u=1 when x=0x=0, u=5u=5 when x=1x=1),

A=∫012πx1+4x2 dx=2π8∫15u du=π4⋅23[u3/2]15=π6(55−1).\begin{align*} A &= \int_0^12\pi x\sqrt{1+4x^2}\,dx \\ &= \frac{2\pi}{8}\int_1^5\sqrt{u}\,du \\ &= \frac{\pi}{4}\cdot\frac{2}{3}\Big[u^{3/2}\Big]_1^5 \\ &= \frac{\pi}{6}\left(5\sqrt5-1\right). \end{align*}

Therefore the area is π6(55−1)≈5.33\frac{\pi}{6}\left(5\sqrt5-1\right)\approx5.33 square units.

(b) Rotating about the xx-axis, the distance to the axis is y=x2y=x^2, and the integral ∫012πx21+4x2 dx\int_0^12\pi x^2\sqrt{1+4x^2}\,dx is a completely different (and much nastier) beast: the spare xx that the substitution u=1+4x2u=1+4x^2 needs is no longer available, so it requires x=12tan⁡θx=\frac{1}{2}\tan\theta or x=12sinh⁡θx=\frac{1}{2}\sinh\theta from the substitution table, followed by a reduction formula for ∫sec⁡3θ dθ\int\sec^3\theta\,d\theta. Compare this with the uu-substitution versus trigonometric-substitution discussion in Integration Techniques — one factor of xx is again the entire difference. So the answer to the question 'which axis?' changes not only the number but the technique.

Example. Find the area of the surface formed when one arch of the cycloid x(t)=t−sin⁡tx(t)=t-\sin t, y(t)=1−cos⁡ty(t)=1-\cos t, 0≤t≤2π0\leq t\leq2\pi, is rotated about the xx-axis.
From the cycloid arc length example (with r=1r=1), [x′(t)]2+[y′(t)]2=2sin⁡(t/2)\sqrt{[x'(t)]^2+[y'(t)]^2}=2\sin(t/2), which is non-negative on [0,2π][0,2\pi]. Also y=1−cos⁡t=2sin⁡2(t/2)y=1-\cos t=2\sin^2(t/2), so the parametric formula gives

A=∫02π2π(2sin⁡2(t/2))⋅2sin⁡(t/2) dt=8π∫02πsin⁡3(t/2) dt=16π∫0πsin⁡3u du(u=t/2,  dt=2 du)=16π×43=64π3,\begin{align*} A &= \int_0^{2\pi}2\pi\big(2\sin^2(t/2)\big)\cdot2\sin(t/2)\,dt \\ &= 8\pi\int_0^{2\pi}\sin^3(t/2)\,dt \\ &= 16\pi\int_0^{\pi}\sin^3u\,du \quad (u=t/2,\; dt=2\,du) \\ &= 16\pi\times\frac{4}{3} \\ &= \frac{64\pi}{3}, \end{align*}

where ∫0πsin⁡3u du=43\int_0^\pi\sin^3u\,du=\frac{4}{3} either by the odd-power substitution w=cos⁡uw=\cos u or straight from the reduction formula for sin⁡n\sin^n. Therefore the surface area is 64π3≈67.02\frac{64\pi}{3}\approx67.02 square units. Notice how writing yy in terms of sin⁡(t/2)\sin(t/2) before integrating turned the whole integrand into a single power of sin⁡(t/2)\sin(t/2); if you leave it as (1−cos⁡t)sin⁡(t/2)(1-\cos t)\sin(t/2) you end up doing more work for the same answer.

Example. Find the area of the surface formed when the polar curve r=1+cos⁡θr=1+\cos\theta, 0≤θ≤π0\leq\theta\leq\pi, is rotated about the xx-axis.
From the cardioid arc length example, r2+(drdθ)2=2∣cos⁡(θ/2)∣=2cos⁡(θ/2)\sqrt{r^2+\left(\frac{dr}{d\theta}\right)^2}=2|\cos(\theta/2)|=2\cos(\theta/2) on [0,π][0,\pi]. Also r=1+cos⁡θ=2cos⁡2(θ/2)r=1+\cos\theta=2\cos^2(\theta/2) and sin⁡θ=2sin⁡(θ/2)cos⁡(θ/2)\sin\theta=2\sin(\theta/2)\cos(\theta/2), so the polar formula gives

A=∫0π2π(2cos⁡2(θ/2))(2sin⁡(θ/2)cos⁡(θ/2))(2cos⁡(θ/2)) dθ=16π∫0πcos⁡4(θ/2)sin⁡(θ/2) dθ=32π∫01w4 dw(w=cos⁡(θ/2),  dw=−12sin⁡(θ/2) dθ)=32π×15=32π5.\begin{align*} A &= \int_0^{\pi}2\pi\big(2\cos^2(\theta/2)\big)\big(2\sin(\theta/2)\cos(\theta/2)\big)\big(2\cos(\theta/2)\big)\,d\theta \\ &= 16\pi\int_0^{\pi}\cos^4(\theta/2)\sin(\theta/2)\,d\theta \\ &= 32\pi\int_0^1w^4\,dw \quad \left(w=\cos(\theta/2),\; dw=-\tfrac{1}{2}\sin(\theta/2)\,d\theta\right) \\ &= 32\pi\times\frac{1}{5} \\ &= \frac{32\pi}{5}. \end{align*}

Therefore the surface area is 32π5≈20.1\frac{32\pi}{5}\approx20.1 square units. The extra factor here is 2πrsin⁡θ2\pi r\sin\theta, not 2πr2\pi r: the distance from the point to the xx-axis is the yy-coordinate rsin⁡θr\sin\theta, not the radius rr. That is the same 'distance to the axis' rule as before, wearing polar clothing.

We finish with a simple but genuinely startling example. Recall from high school that the volume VV of the solid formed when the graph of f:[a,b]→[0,∞)f:[a,b]\to[0,\infty) is rotated about the xx-axis is

V=∫abπ[f(x)]2 dx.V=\int_a^b\pi[f(x)]^2\,dx.

Example. (Gabriel's horn.) Let f:[1,∞)→Rf:[1,\infty)\to\mathbb{R} be given by f(x)=1xf(x)=\frac{1}{x}, and rotate its graph about the xx-axis. Find the volume of the resulting solid and the area of the resulting surface.
The surface produced is known as Gabriel's horn, after the biblical figure, or Torricelli's trumpet, after Evangelista Torricelli, a pupil of Galileo. It is infinitely long, so both quantities have to be handled as improper integrals: compute over [1,R][1,R] and let R→∞R\to\infty.
For the volume,

VR=∫1Rπx2 dx=π[−1x]1R=π(1−1R),\begin{align*} V_R &= \int_1^R\frac{\pi}{x^2}\,dx \\ &= \pi\left[-\frac{1}{x}\right]_1^R \\ &= \pi\left(1-\frac{1}{R}\right), \end{align*}

and lim⁡R→∞VR=π\lim_{R\to\infty}V_R=\pi. So the solid has finite volume, exactly π\pi cubic units.
For the surface area, f′(x)=−1x2f'(x)=-\frac{1}{x^2}, so

AR=∫1R2π⋅1x1+1x4 dx.A_R=\int_1^R2\pi\cdot\frac{1}{x}\sqrt{1+\frac{1}{x^4}}\,dx.

Finding an antiderivative for this integrand looks unpleasant, but we do not need one — we only need to know whether the improper integral converges, and for that a comparison suffices. Whenever x≥1x\geq1,

1x1+1x4>1x1+0=1x,\frac{1}{x}\sqrt{1+\frac{1}{x^4}}>\frac{1}{x}\sqrt{1+0}=\frac{1}{x},

and ∫1∞1x dx\int_1^{\infty}\frac{1}{x}\,dx diverges, so by the comparison test for integrals (Chapter 8 of the MATH1131 calculus notes) the integral for ARA_R diverges too. Therefore Gabriel's horn has infinite surface area but bounds a solid of finite volume π\pi.
This is the painter's paradox: to paint the outside of the horn you would need infinitely much paint, yet you can fill the whole horn with π\pi cubic units of paint and thereby coat the inside surface completely — just pour the paint in and tip out whatever is not touching the wall. The resolution is that 'painting' in the first sense means covering the surface with a coat of some fixed positive thickness, which requires infinite volume, while the paint filling the horn thins out to zero thickness as you go down the tube; a mathematical surface has no thickness at all, so the two notions of 'painted' are simply not the same thing.
Compare this with the spiral r=e−θ/10r=e^{-\theta/10} at the end of Section 5.2, which was infinitely wound but had finite length. Infinite extent tells you nothing on its own about whether length, area or volume comes out finite; only the integral knows, and that is the point of the whole chapter.