In this chapter we point calculus at four geometric problems:
finding the average height of a cable, or the average temperature over a time interval;
finding the length of a curve;
finding the speed of a particle that travels along a curve in the plane; and
finding the surface area of certain solids.
Every one of them is attacked in exactly the same way: chop the object into n 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→∞. 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.
Suppose a cable is suspended between two poles, so that each end of the cable is 3 metres above the ground and the lowest point of the cable is 2 metres above the ground. What is the average height of the cable above the ground? Clearly the answer is somewhere between 2 and 3 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)=c1cosh(cx)
over some interval [a,b], where c 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 f on [a,b]?
Averaging finitely many numbers is easy, but a function takes uncountably many values on [a,b], so we have to say what we even mean. Divide [a,b] into n subintervals of equal length, sample the height of the graph in the kth subinterval by choosing a point ck in it, and take the ordinary arithmetic mean of the sampled heights:
an=n1(f(c1)+f(c2)+⋯+f(cn))=n1k=1∑nf(ck).
As n increases, an 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−a,
Notice what the last sum is: f(ck) is the height and nb−a is the width of the kth rectangle in the usual picture, so it is a Riemann sum for f on [a,b]. Hence, provided f is Riemann integrable,
Definition
Suppose that f is integrable on a closed interval [a,b]. Then the average valuefˉ of f on [a,b] is defined by the formula
fˉ=b−a1∫abf(x)dx.
Rearranging this formula gives a second way to read it: fˉ is the unique constant with
(b−a)fˉ=∫abf(x)dx.
Geometrically, fˉ is the unique y-value for which the rectangle of width b−a and height fˉ has the same area as the region under the graph of f. 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 over [0,2π] is 0, not because 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 on [0,π].
By the definition,
Therefore the average value of sin on [0,π] is π2≈0.6366. As a sanity check, the graph sits between 0 and 1 and spends most of its time up near the top, so an average a bit under 32 is exactly what you would expect; a common wrong guess is 21, 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,
where the x-axis runs along the ground, the y-axis passes through the vertex of the curve, and a=2cosh−1(3/2). Find, to the nearest centimetre, the average height of the cable above the ground.
Set f(x)=2cosh(x/2). The interval has length a−(−a)=2a, so
fˉ=2a1∫−aa2cosh(x/2)dx=2a2∫0a2cosh(x/2)dx(since cosh is even)=a1[4sinh(x/2)]0a=a4sinh(a/2)=cosh−1(3/2)2sinh(cosh−1(3/2)).
The remaining sinh of an inverse cosh is unwound with the identity cosh2t−sinh2t=1:
where the positive root is correct because cosh−1(3/2)>0 and sinh is positive there. Hence
fˉ=cosh−1(3/2)2×25=cosh−1(3/2)5≈2.32.
Therefore the average height of the cable above the ground is about 2.32 metres, which does indeed sit between 2 and 3 as predicted. Note that the fact that the x-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), measured in degrees Celsius t hours after noon, is given by
T(t)=25+2t−3t2.
Find the average temperature between noon and 5 p.m.
Noon is t=0 and 5 p.m. is t=5, so
Therefore the average temperature between noon and 5 p.m. is 9245=2792≈27.2∘C. Do not average the endpoint values T(0)=25 and T(5)=2632 and call it a day; that gives 25.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 f is continuous on [a,b]. Then there is a number c in (a,b) such that
∫abf(t)dt=f(c)(b−a).
Proof. Define F:[a,b]→R by
F(x)=∫axf(t)dt.
By the fundamental theorem of calculus, F is continuous on [a,b], differentiable on (a,b), and F′(x)=f(x). So F satisfies the hypotheses of the ordinary mean value theorem, which supplies a c∈(a,b) with
b−aF(b)−F(a)=F′(c).
But F(a)=0, F(b)=∫abf(t)dt and F′(c)=f(c), so this reads
b−a1∫abf(t)dt=f(c),
which is the claim (in the form fˉ=f(c)). ■
Basically, the mean value theorem for integrals says the horizontal line y=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 0 on [0,1) and 1 at x=1 has average value 0... but a function equal to 0 on [0,1) and 2 on [1,2] has average 1 on [0,2] and never takes the value 1 anywhere.
Example. Verify the mean value theorem for integrals for f(x)=x2 on [0,3] by finding the point c explicitly.
The function is continuous, and its average value is
fˉ=31∫03x2dx=31[3x3]03=31×9=3.
So we need c∈(0,3) with c2=3, i.e. c=3≈1.732 (rejecting c=−3, which is not in the interval). Therefore the theorem holds with c=3. Notice that c is not the midpoint 1.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.
Suppose P0(x0,y0) and P1(x1,y1) are two points in R2. The distance between them is
dist(P0,P1)=(x1−x0)2+(y1−y0)2,
and if P0 and P1 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 C is a curve in R2 which can be expressed in parametric form as
C={(x(t),y(t))∈R2:a≤t≤b},
where x and y are differentiable functions of t. We also assume the parametrisation is chosen so that the moving point (x(t),y(t)) does not retrace its steps, either forwards or backwards. For a≤s≤b, let ℓ(s) denote the arc length of the partial curve
Cs={(x(t),y(t))∈R2:a≤t≤s},
so that ℓ(a)=0 and ℓ(b) is the number we want.
Take a small piece of the curve and approximate it by a secant. Fix t with a<t<b, let h be a small nonzero real number, and consider the two points P(x(t),y(t)) and Q(x(t+h),y(t+h)). The length of the arc from P to Q is ℓ(t+h)−ℓ(t), and it is approximately the length of the chord PQ, which the distance formula gives as
ℓ(t+h)−ℓ(t)≈[x(t+h)−x(t)]2+[y(t+h)−y(t)]2.
Dividing both sides by h (and pushing the h inside the square root, where it becomes h2),
hℓ(t+h)−ℓ(t)≈(hx(t+h)−x(t))2+(hy(t+h)−y(t))2.
The approximation improves as h shrinks, so assuming that ℓ is a differentiable function, letting h→0 gives
ℓ′(t)=[x′(t)]2+[y′(t)]2.
By the fundamental theorem of calculus,
ℓ(s)=∫as[x′(t)]2+[y′(t)]2dt+K
for some constant K, and putting s=a gives 0=ℓ(a)=0+K, so K=0. Taking s=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: x and y must be differentiable, the parametrisation must trace the curve once and only once, and ℓ itself is assumed differentiable.
Collecting the derivation above: if a curve C is given parametrically by
C={(x(t),y(t))∈R2:a≤t≤b},
where x and y are differentiable functions of t, then its arc length ℓ is
ℓ=∫ab[x′(t)]2+[y′(t)]2dt.
Basically, [x′(t)]2+[y′(t)]2dt is the length of an infinitesimal chord, obtained from Pythagoras applied to the horizontal displacement x′(t)dt and the vertical displacement y′(t)dt; integrating adds all the little chords up.
Example. Find the arc length of one arch of the cycloid
x(t)=r(t−sint),y(t)=r(1−cost),0≤t≤2π.
(This is the curve traced by a point on the rim of a wheel of radius r rolling along the x-axis, and it is closely related to the 'curve of fastest descent' from MATH1131.) Differentiating,
Before substituting this into the formula it is best to express 1−cost as a square, since we are about to take a square root. The double angle identity 1−cos2θ=2sin2θ with θ=t/2 gives
[x′(t)]2+[y′(t)]2=2r2⋅2sin2(t/2)=4r2sin2(t/2).
Now A2=∣A∣, not A — but sin(t/2) is positive for 0<t<2π, so here the naive square root causes no problems. Hence
Therefore one arch of the cycloid has length 8r. Notice how clean that is: the arch spans a horizontal distance of 2πr≈6.28r and the answer is exactly 8r, with no π in sight at all.
Example. Use the parametrisation x(t)=rcost, y(t)=rsint with 0≤t≤2π to find the circumference of a circle of radius r. What goes wrong with the parametrisation x(t)=rcos2t, y(t)=rsin2t on the same interval?
For the first parametrisation, x′(t)=−rsint and y′(t)=rcost, so
ℓ=∫02πr2sin2t+r2cos2tdt=∫02πrdt=2πr,
which is the circumference we expect. For the second, x′(t)=−2rsin2t and y′(t)=2rcos2t, and the same computation gives
∫02π4r2sin22t+4r2cos22tdt=∫02π2rdt=4πr,
which is twice the circumference. The reason is that as t runs from 0 to 2π, the point (rcos2t,rsin2t) 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 0 to π.
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.
Most curves you meet are given as y=f(x) rather than parametrically, but that is just a special case: parametrise the graph of f on [a,b] by
x(t)=t,y(t)=f(t),a≤t≤b.
Then x′(t)=1 and y′(t)=f′(t), so the parametric formula immediately gives ℓ=∫ab1+[f′(t)]2dt. Renaming the variable of integration x, as one usually does, the arc length ℓ of the graph of f on [a,b] is
ℓ=∫ab1+[f′(x)]2dx.
Notice immediately that ℓ≥b−a, since the integrand is at least 1: 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/2 for 0≤x≤1.
Here f′(x)=23x1/2, so [f′(x)]2=49x and
Therefore the arc length is 271313−8≈1.4397. The chord from (0,0) to (1,1) has length 2≈1.4142, and the curve should be slightly longer than the chord; it is, by about 2%.
Example. Find the arc length of y=31(x2+2)3/2 for 0≤x≤3.
Differentiating with the chain rule, f′(x)=31⋅23(x2+2)1/2⋅2x=xx2+2, so
1+[f′(x)]2=1+x2(x2+2)=x4+2x2+1=(x2+1)2.
The square root is now trivial (and x2+1>0, so no absolute values are needed), giving
ℓ=∫03(x2+1)dx=[3x3+x]03=9+3=12.
Therefore the arc length is exactly 12 units.
That example is worth pausing on. Arc-length integrands in this course are always engineered so that 1+[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′ before you start hunting for exotic substitutions. The reason is that ∫1+[f′(x)]2dx is hopeless for almost every f you could write down: even y=x2 leads to ∫1+4x2dx, which needs a hyperbolic substitution from the substitution table, and y=sinx leads to an elliptic integral that has no elementary antiderivative at all.
Example. Find the arc length of y=4x2−2lnx for 1≤x≤2.
Differentiating, f′(x)=2x−2x1, so
Therefore the arc length is 43+2ln2≈1.0966, comfortably more than the span b−a=1. This is the standard shape of a hand-made arc length question: f′ is of the form u−4u1 or 2u−2u1, so that the cross term in [f′]2 is exactly −1 and the +1 repairs the square.
Example. Find the arc length of the catenary
y=a1cosh(ax),x∈[−b,b].
A catenary is the shape of a hanging cable, which is why this function keeps reappearing. Here f′(x)=sinh(ax), and the hyperbolic Pythagorean identity cosh2t−sinh2t=1 does the collapsing for us:
ℓ=∫−bb1+sinh2(ax)dx=∫−bbcosh2(ax)dx=∫−bbcosh(ax)dx(since cosh is always positive)=2∫0bcosh(ax)dx(since cosh is even)=2[a1sinh(ax)]0b=a2sinh(ab)=a1(eab−e−ab).
Therefore the catenary has arc length a1(eab−e−ab) units. Notice how this is the only function whose arc length integrand simplifies without any factoring at all — the identity 1+sinh2=cosh2 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) on [−a,a] where a=2cosh−1(3/2).
Comparing with the general catenary, c1cosh(cx)=2cosh(x/2) forces c=21, and the interval endpoint is b=a. So the formula gives
reusing the value sinh(cosh−1(3/2))=25 computed earlier. Therefore the cable is 25≈4.47 metres long. As a sanity check, the poles are 2a=4cosh−1(3/2)≈3.85 metres apart, so the cable is longer than the gap it spans — as it must be, since it sags.
where the cross terms cancelled and sin2θ+cos2θ=1 mopped up the rest. Writing f(θ) as r and f′(θ) as dθdr, the arc length ℓ of a polar curve is
ℓ=∫θ0θ1r2+(dθdr)2dθ.
Basically, moving through an angle dθ carries you rdθ around the circle and dθdrdθ outwards, and these two displacements are perpendicular, so Pythagoras gives the little bit of arc. The r2 term is not optional and it is not (dθdr)2 in disguise; forgetting it is the classic polar arc length mistake (it would say that a circle r=1, on which dθdr=0, has length zero).
Example. Find the length of the cardioid r=1+cosθ.
The full cardioid is traced once as θ runs from 0 to 2π. Here dθdr=−sinθ, so
using the double angle identity 1+cos2α=2cos2α with α=θ/2 in the last line. Now be careful with the square root: 4cos2(θ/2)=2∣cos(θ/2)∣, and cos(θ/2) is positive on (0,π) but negative on (π,2π). The cardioid is symmetric about the x-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.
Therefore the cardioid has length 8 units. If you drop the absolute value and integrate 2cos(θ/2) straight from 0 to 2π you get 0, 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θ for 0≤θ≤2π.
Here dθdr=eθ=r, so the integrand is e2θ+e2θ=2eθ and
ℓ=2∫02πeθdθ=2[eθ]02π=2(e2π−1).
Therefore the length is 2(e2π−1)≈755.9 units — the exponential spiral grows fast enough that a single turn is already enormous.
Example. The spiral r=e−θ/10, where θ≥0, winds infinitely many times around the origin. Is its total arc length finite? Explain.
Compute the arc length of the segment 0≤θ≤θ1 first, and take a limit afterwards. Here dθdr=−101e−θ/10, so
As θ1→∞ we have e−θ1/10→0, so ℓ(θ1)→101. Therefore the total arc length is finite and equals 101≈10.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.
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 P moves in the plane with position (x(t),y(t)) at time t. The distance s(t) travelled from time zero to time t is exactly the arc length of the path traversed in that time,
s(t)=∫0t[x′(u)]2+[y′(u)]2du.
By definition the speed of P is the rate of change of distance travelled with respect to time, so by the fundamental theorem of calculus the speed v(t) is
v(t)=s′(t)=[x′(t)]2+[y′(t)]2.
If we bundle the position into a vector r(t)=(x(t),y(t)), then r′(t)=(x′(t),y′(t)) is the velocity vector and the formula just says v(t)=∣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:
dtds=v(t)ands=∫v(t)dt.
Example. A stone is thrown horizontally from the deck of the Sydney Harbour Bridge at 20 metres per second. Its position t seconds after the throw is
x(t)=20t,y(t)=50−5t2,0≤t≤10,
where y(t) is the height above the water. Find the speed of the stone an instant before it hits the water.
Differentiating, x′(t)=20 and y′(t)=−10t, so the speed at time t is
v(t)=202+100t2=400+100t2.
The stone hits the water when y(t)=0, that is when 5t2=50, i.e. t=10. Hence
t→(10)−limv(t)=400+100×10=1400=1014≈37.42.
Therefore the stone is travelling at about 37.42 metres per second just before impact. Notice that the horizontal component of the velocity is a constant 20 throughout — gravity only ever changes the vertical component — so the speed never drops below 20.
Example. A projectile is fired from an elevated cannon. Its horizontal distance x (in metres) from the cannon and height y (in metres) above the ground, exactly t seconds after firing, are
x(t)=40t,y(t)=−5t2+40t+45,0≤t≤t1,
where t1 is the time of impact.
(a) Find t1.
(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=0:
−5t2+40t+45t2−8t−9(t−9)(t+1)=0=0=0,
so t=9 or t=−1; discarding the negative time, t1=9 seconds.
(b) Here x′(t)=40 and y′(t)=−10t+40, so at t=9 we have y′(9)=−50 and
v(9)=402+(−50)2=1600+2500=4100=1041≈64.03.
Therefore the projectile is travelling at about 64.03 metres per second immediately prior to impact.
Therefore the average height of the projectile was 90 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 x would be a different integral (though here it happens to agree, since x is a linear function of t).
For completeness, the distance actually flown is ∫091600+(40−10t)2dt≈427.53 metres, which needs the substitution u=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 1. Such a parametrisation is said to be by arc length: if v(t)=1 for all t then s(t)=∫0t1du=t, so the parameter is the distance travelled along the curve.
Example. A particle moves on the circle x(t)=3cos2t, y(t)=3sin2t. Find its speed and the distance travelled by time t, and hence reparametrise the circle by arc length.
Differentiating, x′(t)=−6sin2t and y′(t)=6cos2t, so
v(t)=36sin22t+36cos22t=6,
a constant. Hence s(t)=∫0t6du=6t, i.e. t=6s. Substituting this back into the position,
x=3cos(3s),y=3sin(3s),
and it is easily checked that this new parametrisation has speed sin2(s/3)+cos2(s/3)=1. Therefore the arc length parametrisation of the circle of radius 3 is (3cos3s,3sin3s), and running s from 0 to 6π traverses it once — recovering the circumference 2π×3=6π for free. In general the arc-length parametrisation of a circle of radius R is (RcosRs,RsinRs); the R1 inside is exactly what makes the speed 1 rather than R.
Example. The position of a particle P at time t≥0 is
x(t)=cos(2π(cosπt−1)),y(t)=−sin(2π(cosπt−1)).
(a) Find a formula for the speed 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]?
(a) Write ϕ(t)=2π(cosπt−1), so that x=cosϕ and y=−sinϕ. By the chain rule,
Therefore v(t)=2π2∣sinπt∣; in particular the particle is momentarily at rest whenever t is an integer, and fastest when t is a half-integer.
(b) Since cosπt−1 takes every value in [−2,0], the angle ϕ ranges over [−π,0], and (cosϕ,−sinϕ) then sweeps out the upper half of the unit circle. So the trajectory is the semicircle of centre (0,0) and radius 1 lying in the upper half-plane, whose length is half the circumference, namely π.
(c) Distance travelled is the integral of speed, and ∣sinπt∣ has period 1 with ∫01∣sinπt∣dt=[−πcosπt]01=π2. Hence
∫03v(t)dt=2π2×3×π2=3π.
Therefore the particle travels 3π units in the first three seconds. Compare (b) and (c): the curve is only π units long, but the particle covers 3π 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.
Finding the surface area of a general surface in R3 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 s and radii r and R, the area A of its curved surface is given by
A=π(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π(2r+R)s, i.e. (circumference at the average radius) × (slant height) — the frustum is a rectangle of width 2πrˉ and height s, rolled up.
5.4.1 An heuristic derivation for the surface area of a surface of revolution#
Suppose a curve C has parametrisation
C={(x(t),y(t))∈R2:a≤t≤b}.
We assume that the curve lies in the upper half-plane (more precisely y(t)≥0, meeting the x-axis at only finitely many points if at all), and that C is simple: if (x(t0),y(t0))=(x(t1),y(t1)) then t0=t1, i.e. the curve does not cross itself. Rotating C about the x-axis produces a surface of revolution, and we want its area.
The derivation runs exactly parallel to the arc length one. Let A(s) denote the area of the surface formed by rotating the partial curve {(x(t),y(t)):a≤t≤s}, and assume A is differentiable. Fix t∈(a,b), let h be small and nonzero, and consider P(x(t),y(t)) and Q(x(t+h),y(t+h)). Then A(t+h)−A(t) is the area swept out by the arc from P to Q, which for small h is approximately the area swept by the chordPQ — and a rotated chord is precisely a frustum, of radii y(t) and y(t+h) and slant height equal to the length of PQ. So the frustum formula gives
Assume throughout that C lies in the upper half-plane and is simple, and that all the derivatives written down exist. Rotating about the x-axis, the area A of the surface of revolution is
A=∫ab2πy(t)[x′(t)]2+[y′(t)]2dt
if C is given parametrically by (x(t),y(t)) with a≤t≤b;
A=∫ab2πf(x)1+[f′(x)]2dx
if C is the graph y=f(x) on [a,b]; and
A=∫θ0θ12πrsinθr2+(dθdr)2dθ
if C is the polar curve r=f(θ) with θ0≤θ≤θ1. Only the first of these needs deriving; the other two follow from it via the parametrisations x(t)=t, y(t)=f(t) and x(θ)=f(θ)cosθ, y(θ)=f(θ)sinθ used in Sections 5.2.3 and 5.2.4. Rotating about the y-axis instead, the parametric formula becomes
A=∫ab2πx(t)[x′(t)]2+[y′(t)]2dt,
with the graph version A=∫ab2πx1+[f′(x)]2dx 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): it is y when you rotate about the x-axis and x when you rotate about the y-axis. Writing 2πy out of habit while rotating about the y-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π× 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 r, by rotating the semicircle y=r2−x2, −r≤x≤r, about the x-axis.
Differentiating with the chain rule,
f′(x)=2r2−x2−2x=r2−x2−x,
so that
1+[f′(x)]2=1+r2−x2x2=r2−x2r2−x2+x2=r2−x2r2.
Hence 1+[f′(x)]2=r2−x2r, and the crucial cancellation happens:
Therefore the surface area of a sphere of radius r is 4πr2, as it should be. The same answer falls out of the parametric formula with x(t)=rcost, y(t)=rsint, 0≤t≤π: there [x′]2+[y′]2=r and A=∫0π2πrsint⋅rdt=2πr2[−cost]0π=4πr2. Note the limits 0≤t≤π, not 0≤t≤2π: the lower semicircle would sweep out the same sphere a second time and double the answer.
The constant integrand 2π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=b and x=c (with −r≤b<c≤r) is
∫bc2πrdx=2π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 S is formed by rotating the curve y=2−x, x∈[0,2], about the x-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=0 where the curve has height 2. For the paraboloid, with f(x)=(2−x)1/2 we get f′(x)=−21(2−x)−1/2, so
where in the fourth line the 2−x+41=49−4x=219−4x absorbed the factor of 2. Notice how pulling the outside 2−xinside the root is what made the integrand elementary; that manoeuvre is standard whenever f and f1 both appear. The cap is a disc of radius 2, so
A2=π(2)2=2π.
Therefore the total surface area of S is
A1+A2=313π+2π=319π≈19.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=x2, 0≤x≤1, is rotated about (a) the y-axis and (b) the x-axis.
Both parts share the same square root: f′(x)=2x, so 1+[f′(x)]2=1+4x2.
(a) Rotating about the y-axis, the distance to the axis is x, so with the substitution u=1+4x2, du=8xdx (and u=1 when x=0, u=5 when x=1),
Therefore the area is 6π(55−1)≈5.33 square units.
(b) Rotating about the x-axis, the distance to the axis is y=x2, and the integral ∫012πx21+4x2dx is a completely different (and much nastier) beast: the spare x that the substitution u=1+4x2 needs is no longer available, so it requires x=21tanθ or x=21sinhθ from the substitution table, followed by a reduction formula for ∫sec3θdθ. Compare this with the u-substitution versus trigonometric-substitution discussion in Integration Techniques — one factor of x 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−sint, y(t)=1−cost, 0≤t≤2π, is rotated about the x-axis.
From the cycloid arc length example (with r=1), [x′(t)]2+[y′(t)]2=2sin(t/2), which is non-negative on [0,2π]. Also y=1−cost=2sin2(t/2), so the parametric formula gives
where ∫0πsin3udu=34 either by the odd-power substitution w=cosu or straight from the reduction formula for sinn. Therefore the surface area is 364π≈67.02 square units. Notice how writing y in terms of sin(t/2)before integrating turned the whole integrand into a single power of sin(t/2); if you leave it as (1−cost)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θ, 0≤θ≤π, is rotated about the x-axis.
From the cardioid arc length example, r2+(dθdr)2=2∣cos(θ/2)∣=2cos(θ/2) on [0,π]. Also r=1+cosθ=2cos2(θ/2) and sinθ=2sin(θ/2)cos(θ/2), so the polar formula gives
Therefore the surface area is 532π≈20.1 square units. The extra factor here is 2πrsinθ, not 2πr: the distance from the point to the x-axis is the y-coordinate rsinθ, not the radius r. 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 V of the solid formed when the graph of f:[a,b]→[0,∞) is rotated about the x-axis is
V=∫abπ[f(x)]2dx.
Example. (Gabriel's horn.) Let f:[1,∞)→R be given by f(x)=x1, and rotate its graph about the x-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] and let R→∞.
For the volume,
VR=∫1Rx2πdx=π[−x1]1R=π(1−R1),
and limR→∞VR=π. So the solid has finite volume, exactly π cubic units.
For the surface area, f′(x)=−x21, so
AR=∫1R2π⋅x11+x41dx.
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≥1,
x11+x41>x11+0=x1,
and ∫1∞x1dx diverges, so by the comparison test for integrals (Chapter 8 of the MATH1131 calculus notes) the integral for AR diverges too. Therefore Gabriel's horn has infinite surface area but bounds a solid of finite volume π.
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 π 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−θ/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.