Skip to main content
Logo image

APEX Calculus: for University of Lethbridge

Section 15.5 Surface Area

In Section 7.4 we used definite integrals to compute the arc length of plane curves of the form \(y=f(x)\text{.}\) We later extended these ideas to compute the arc length of plane curves defined by parametric or polar equations.
The natural extension of the concept of “arc length over an interval” to surfaces is “surface area over a region.”
Consider the surface \(z=f(x,y)\) over a region \(R\) in the \(xy\)-plane, shown in Figure 15.5.1.(a). Because of the domed shape of the surface, the surface area will be greater than that of the area of the region \(R\text{.}\) We can find this area using the same basic technique we have used over and over: we’ll make an approximation, then using limits, we’ll refine the approximation to the exact value.
Figure 15.5.1. Developing a method of computing surface area
As done to find the volume under a surface or the mass of a lamina, we subdivide \(R\) into \(n\) subregions. Here we subdivide \(R\) into rectangles, as shown in the figure. One such subregion is outlined in the figure, where the rectangle has dimensions \(\dx_i\) and \(\dy_i\text{,}\) along with its corresponding region on the surface.
In part Figure 15.5.1.(b) of the figure, we zoom in on this portion of the surface. When \(\dx_i\) and \(\dy_i\) are small, the function is approximated well by the tangent plane at any point \((x_i,y_i)\) in this subregion, which is graphed in part Figure 15.5.1.(b). In fact, the tangent plane approximates the function so well that in this figure, it is virtually indistinguishable from the surface itself! Therefore we can approximate the surface area \(S_i\) of this region of the surface with the area \(T_i\) of the corresponding portion of the tangent plane.
This portion of the tangent plane is a parallelogram, defined by sides \(\vec u\) and \(\vec v\text{,}\) as shown. One of the applications of the cross product from Section 12.4 is that the area of this parallelogram is \(\norm{\vec u\times \vec v}\text{.}\) Once we can determine \(\vec u\) and \(\vec v\text{,}\) we can determine the area.
The vector \(\vec u\) is tangent to the surface in the direction of \(x\text{,}\) therefore, from Section 14.4, \(\vec u\) is parallel to \(\la 1,0,f_x(x_i,y_i)\ra\text{.}\) The \(x\)-displacement of \(\vec u\) is \(\dx_i\text{,}\) so we know that \(\vec u = \dx_i\la 1,0,f_x(x_i,y_i)\ra\text{.}\) Similar logic shows that \(\vec v = \dy_i\la 0,1,f_y(x_i,y_i)\ra\text{.}\) Thus:
\begin{align*} \text{surface area \(S_i\)} \amp \approx \,\text{area of \(T_i\)}\\ \amp = \norm{\vec u\times \vec v}\\ \amp = \norm{\dx_i\la 1,0,f_x(x_i,y_i)\ra\times\dy_i\la 0,1,f_y(x_i,y_i)\ra}\\ \amp =\sqrt{1+f_x(x_i,y_i)^2+f_y(x_i,y_i)^2}\dx_i\dy_i\text{.} \end{align*}
Note that \(\dx_i\dy_i = \Delta A_i\text{,}\) the area of the \(i\)th subregion.
Summing up all \(n\) of the approximations to the surface area gives
\begin{equation*} \text{surface area over \(R\)} \approx \sum_{i=1}^n \sqrt{1+f_x(x_i,y_i)^2+f_y(x_i,y_i)^2}\Delta A_i\text{.} \end{equation*}
Once again take a limit as all of the \(\dx_i\) and \(\dy_i\) shrink to 0; this leads to a double integral.

Definition 15.5.2. Surface Area.

Let \(z=f(x,y)\) where \(f_x\) and \(f_y\) are continuous over a closed, bounded region \(R\text{.}\) The surface area \(S\) over \(R\) is
\begin{align*} S \amp = \iint_R \, dS\\ \amp =\iint_R \sqrt{1+f_x(x,y)^2+f_y(x,y)^2}\, dA\text{.} \end{align*}
We test this definition by using it to compute surface areas of known surfaces. We start with a triangle.

Example 15.5.3. Finding the surface area of a plane over a triangle.

Let \(f(x,y) = 4-x-2y\text{,}\) and let \(R\) be the region in the plane bounded by \(x=0\text{,}\) \(y=0\) and \(y=2-x/2\text{,}\) as shown in Figure 15.5.4. Find the surface area of \(f\) over \(R\text{.}\)
Figure 15.5.4. Finding the area of a triangle in space in Example 15.5.3
Solution.
We follow Definition 15.5.2. We start by noting that \(f_x(x,y) = -1\) and \(f_y(x,y) = -2\text{.}\) To define \(R\text{,}\) we use bounds \(0\leq y\leq 2-x/2\) and \(0\leq x\leq 4\text{.}\) Therefore
\begin{align*} S \amp = \iint_R\, dS\\ \amp = \int_0^4\int_0^{2-x/2} \sqrt{1+(-1)^2+(-2)^2}\, dy\, dx\\ \amp = \int_0^4 \sqrt{6}\left(2-\frac x2\right)\, dx\\ \amp = 4\sqrt{6}\text{.} \end{align*}
Because the surface is a triangle, we can figure out the area using geometry. Considering the base of the triangle to be the side in the \(xy\)-plane, we find the length of the base to be \(\sqrt{20}\text{.}\) We can find the height using our knowledge of vectors: let \(\vec u\) be the side in the \(xz\)-plane and let \(\vec v\) be the side in the \(xy\)-plane. The height is then \(\norm {\vec u - \proj{u}{v}} = 4\sqrt{6/5}\text{.}\) Geometry states that the area is thus
\begin{equation*} \frac 12\cdot4\sqrt{6/5}\cdot\sqrt{20} = 4\sqrt{6}\text{.} \end{equation*}
We affirm the validity of our formula.
It is “common knowledge” that the surface area of a sphere of radius \(r\) is \(4\pi r^2\text{.}\) We confirm this in the following example, which involves using our formula with polar coordinates.

Example 15.5.5. The surface area of a sphere.

Find the surface area of the sphere with radius \(a\) centered at the origin, whose top hemisphere has equation \(f(x,y)=\sqrt{a^2-x^2-y^2}\text{.}\)
Solution.
We start by computing partial derivatives and find
\begin{equation*} f_x(x,y) = \frac{-x}{\sqrt{a^2-x^2-y^2}} \text{ and } f_y(x,y) = \frac{-y}{\sqrt{a^2-x^2-y^2}}\text{.} \end{equation*}
As our function \(f\) only defines the top upper hemisphere of the sphere, we double our surface area result to get the total area:
\begin{align*} S \amp = 2\iint_R \sqrt{1+ f_x(x,y)^2+f_y(x,y)^2}\, dA\\ \amp = 2\iint_R \sqrt{1+ \frac{x^2+y^2}{a^2-x^2-y^2}}\, dA\text{.} \end{align*}
The region \(R\) that we are integrating over is bounded by the circle, centered at the origin, with radius \(a\text{:}\) \(x^2+y^2=a^2\text{.}\) Because of this region, we are likely to have greater success with our integration by converting to polar coordinates. Using the substitutions \(x=r\cos(\theta)\text{,}\) \(y=r\sin(\theta)\text{,}\) \(dA = r\, dr\, d\theta\) and bounds \(0\leq\theta\leq2\pi\) and \(0\leq r\leq a\text{,}\) we have:
\begin{align} S \amp = 2\int_0^{2\pi}\int_0^a \sqrt{1+\frac{r^2\cos^2(\theta) +r^2\sin^2(\theta) }{a^2-r^2\cos^2(\theta) -r^2\sin^2(\theta) }}\, r\, dr\, d\theta\notag\\ \amp =2\int_0^{2\pi}\int_0^ar\sqrt{1+\frac{r^2}{a^2-r^2}}\, dr\, d\theta\notag\\ \amp =2\int_0^{2\pi}\int_0^ar\sqrt{\frac{a^2}{a^2-r^2}}\, dr\, d\theta.\tag{15.5.1}\\ \end{align}
Apply substitution \(u=a^2-r^2\) and integrate the inner integral, giving
\begin{align} \amp = 2\int_0^{2\pi} a^2\, d\theta\notag\\ \amp = 4\pi a^2\text{.}\notag \end{align}
Our work confirms our previous formula.

Example 15.5.6. Finding the surface area of a cone.

The general formula for a right cone with height \(h\) and base radius \(a\) is
\begin{equation*} \ds f(x,y) = h-\frac{h}a\sqrt{x^2+y^2}\text{,} \end{equation*}
shown in Figure 15.5.7. Find the surface area of this cone.
Figure 15.5.7. Finding the surface area of a cone in Example 15.5.6
Solution.
We begin by computing partial derivatives.
\begin{equation*} f_x(x,y) = -\frac{xh}{a\sqrt{x^2+y^2}} \text{ and } f_y(x,y)=-\frac{yh}{a\sqrt{x^2+y^2}}\text{.} \end{equation*}
Since we are integrating over the disk bounded by \(x^2+y^2=a^2\text{,}\) we again use polar coordinates. Using the standard substitutions, our integrand becomes
\begin{equation*} \sqrt{1+ \left(\frac{hr\cos(\theta) }{a\sqrt{r^2}}\right)^2 + \left(\frac{hr\sin(\theta) }{a\sqrt{r^2}}\right)^2}\text{.} \end{equation*}
This may look intimidating at first, but there are lots of simple simplifications to be done. It amazingly reduces to just
\begin{equation*} \sqrt{1+\frac{h^2}{a^2}} = \frac1a\sqrt{a^2+h^2}\text{.} \end{equation*}
Our polar bounds are \(0\leq\theta\leq2\pi\) and \(0\leq r\leq a\text{.}\) Thus
\begin{align*} S \amp = \int_0^{2\pi}\int_0^ar\frac1a\sqrt{a^2+h^2}\, dr\, d\theta\\ \amp = \int_0^{2\pi} \left.\left(\frac12r^2\frac1a\sqrt{a^2+h^2}\right)\right|_0^ad\theta\\ \amp = \int_0^{2\pi} \frac12a\sqrt{a^2+h^2} \, d\theta\\ \amp = \pi a\sqrt{a^2+h^2}\text{.} \end{align*}
This matches the formula found in the back of this text.

Example 15.5.8. Finding surface area over a region.

Find the area of the surface \(f(x,y) = x^2-3y+3\) over the region \(R\) bounded by \(-x\leq y\leq x\text{,}\) \(0\leq x\leq 4\text{,}\) as pictured in Figure 15.5.9.
Figure 15.5.9. Graphing the surface in Example 15.5.8
Solution.
It is straightforward to compute \(f_x(x,y) = 2x\) and \(f_y(x,y) = -3\text{.}\) Thus the surface area is described by the double integral
\begin{equation*} \iint_R \sqrt{1+(2x)^2+(-3)^2}\, dA = \iint_R \sqrt{10+4x^2}\, dA\text{.} \end{equation*}
As with integrals describing arc length, double integrals describing surface area are in general hard to evaluate directly because of the square-root. This particular integral can be easily evaluated, though, with judicious choice of our order of integration.
Integrating with order \(dx\, dy\) requires us to evaluate \(\int \sqrt{10+4x^2}\, dx\text{.}\) This can be done, though it involves Integration By Parts and \(\sinh^{-1}(x)\text{.}\) Integrating with order \(dy\, dx\) has as its first integral \(\int \sqrt{10+4x^2}\, dy\text{,}\) which is easy to evaluate: it is simply \(y\sqrt{10+4x^2}+C\text{.}\) So we proceed with the order \(dy\, dx\text{;}\) the bounds are already given in the statement of the problem.
\begin{align*} \iint_R\sqrt{10+4x^2}\, dA \amp = \int_0^4\int_{-x}^x\sqrt{10+4x^2}\, dy \, dx\\ \amp = \int_0^4\left.\big(y\sqrt{10+4x^2}\big)\right|_{-x}^x dx\\ \amp =\int_0^4\big(2x\sqrt{10+4x^2}\big)\, dx.\\ \end{align*}
Apply substitution with \(u = 10+4x^2\text{:}\)
\begin{align*} \amp = \left.\left(\frac16\big(10+4x^2\big)^{3/2}\right)\right|_0^4\\ \amp = \frac13\big(37\sqrt{74}-5\sqrt{10}\big) \approx 100.825\,\text{units}^2\text{.} \end{align*}
So while the region \(R\) over which we integrate has an area of 16 square units, the surface has a much greater area as its \(z\)-values change dramatically over \(R\text{.}\)
In practice, technology helps greatly in the evaluation of such integrals. High powered computer algebra systems can compute integrals that are difficult, or at least time consuming, by hand, and can at the least produce very accurate approximations with numerical methods. In general, just knowing how to set up the proper integrals brings one very close to being able to compute the needed value. Most of the work is actually done in just describing the region \(R\) in terms of polar or rectangular coordinates. Once this is done, technology can usually provide a good answer.
We have learned how to integrate integrals; that is, we have learned to evaluate double integrals. In the next section, we learn how to integrate double integrals — that is, we learn to evaluate triple integrals, along with learning some uses for this operation.

Exercises Exercises

Terms and Concepts

1.
“Surface area” is analogous to what previously studied concept?
2.
To approximate the area of a small portion of a surface, we computed the area of its plane.
3.
We interpret \(\ds \iint_R\, dS\) as “sum up lots of little .”
4.
Why is it important to know how to set up a double integral to compute surface area, even if the resulting integral is hard to evaluate?
5.
Why do \(z=f(x,y)\) and \(z=g(x,y)=f(x,y)+h\text{,}\) for some real number \(h\text{,}\) have the same surface area over a region \(R\text{?}\)
6.
Let \(f(x,y)\) be a function defined over a region \(R\) and let \(g(x,y)=2f(x,y)\text{.}\) Why is the surface area of \(z=g(x,y)\) over \(R\) not twice the surface area of \(z=f(x,y)\) over \(R\text{?}\)

Problems

Exercise Group.
In the following exercises, set up the iterated integral that computes the surface area of the graph of the given function over the region \(R\text{.}\)
7.
\(f(x,y) = \sin(x) \cos(y)\text{;}\)\(R\) is the rectangle with bounds \(0\leq x\leq 2\pi\text{,}\) \(0\leq y\leq2\pi\text{.}\)
8.
\(\ds f(x,y) = \frac{1}{x^2+y^2+1}\text{;}\)\(R\) is bounded by the circle \(x^2+y^2=9\text{.}\)
9.
\(\ds f(x,y) = x^2-y^2\text{;}\)\(R\) is the rectangle with opposite corners \((-1,-1)\) and \((1,1)\text{.}\)
10.
\(\ds f(x,y) = \frac{1}{e^{x^2}+1}\text{;}\)\(R\) is the rectangle bounded by
\(-5\leq x\leq 5\) and \(0\leq y\leq 1\text{.}\)
Exercise Group.
In the following exercises, find the area of the given surface over the region \(R\text{.}\)
11.
\(z = 3x-7y+2\text{;}\) \(R\) is the rectangle with opposite corners \((-1,0)\) and \((1,3)\text{.}\)
12.
\(z = 2x+2y+2\text{;}\) \(R\) is the triangle with corners \((0,0)\text{,}\) \((1,0)\) and \((0,1)\text{.}\)
13.
\(z = x^2+y^2+10\text{;}\) \(R\) is bounded by the circle \(x^2+y^2=16\text{.}\)
14.
\(z = -2x+4y^2+7\) over \(R\text{,}\) the triangle bounded by \(y=-x\text{,}\) \(y=x\text{,}\) \(0\leq y\leq 1\text{.}\)
15.
\(z = x^2+y\) over \(R\text{,}\) the triangle bounded by \(y=2x\text{,}\) \(y=0\) and \(x=2\text{.}\)
16.
\(z = \frac23x^{3/2}+2y^{3/2}\) over \(R\text{,}\) the rectangle with opposite corners \((0,0)\) and \((1,1)\text{.}\)
17.
\(z = 10-2\sqrt{x^2+y^2}\) over the region \(R\) bounded by the circle \(x^2+y^2=25\text{.}\) (This is the cone with height 10 and base radius 5; be sure to compare your result with the known formula.)
18.
Find the surface area of the sphere with radius 5 by doubling the surface area of \(f(x,y) = \sqrt{25-x^2-y^2}\) over the region \(R\) bounded by the circle \(x^2+y^2=25\text{.}\) (Be sure to compare your result with the known formula.)
19.
Find the surface area of the ellipse formed by restricting the plane \(f(x,y) = cx+dy+h\) to the region \(R\) bounded by the circle \(x^2+y^2=1\text{,}\) where \(c\text{,}\) \(d\) and \(h\) are some constants. Your answer should be given in terms of \(c\) and \(d\text{;}\) why does the value of \(h\) not matter?