For this development, we will use a parametrized surface given by
\begin{equation*}
\vr(s,t) = x(s,t) \vi + y(s,t) \vj + z(s,t) \vk
\end{equation*}
with bounds on our parameters given by \(a \leq s \leq b\) and \(c \leq t \leq d\text{.}\) In this case, we are using a rectangular domain in the \(st\)-plane but in general, this same approach will work for non-rectangular domains. The partial derivatives of this parameterization,
\begin{align*}
\vr_s(s,t) \amp = x_s(s,t) \vi + y_s(s,t) \vj + z_s(s,t) \vk\\
\vr_t(s,t) \amp = x_t(s,t) \vi + y_t(s,t) \vj + z_t(s,t) \vk
\end{align*}
will give vectors tangent to the surface in the direction of changes in \(s\) and \(t\text{,}\) respectively.
We will break our section of the curved surface into steps in
\(s\) and
\(t\text{.}\) In particular, we partition the interval of
\(s\)-values
\([a,b]\) into
\(m\) subintervals of length
\(\Delta s = \frac{b-a}{m}\) and let
\(s_0\text{,}\) \(s_1\text{,}\) \(\ldots\text{,}\) \(s_m\) be the endpoints of these subintervals, where
\(a = s_0\lt s_1\lt s_2 \lt \cdots \lt s_m = b\text{.}\) We also partition the interval of
\(t\)-values
\([c,d]\) into
\(n\) subintervals of equal length
\(\Delta t = \frac{d-c}{n}\) and let
\(t_0\text{,}\) \(t_1\text{,}\) \(\ldots\text{,}\) \(t_n\) be the endpoints of these subintervals, where
\(c = t_0\lt t_1\lt t_2 \lt \cdots \lt t_n = d\text{.}\)
These subintervals partition the rectangle
\(R = [a,b] \times [c,d]\) in
\(st\)-coordinates into
\(mn\) sub-rectangles
\(R_{ij}\) with opposite vertices
\((s_{i-1},t_{j-1})\) and
\((s_i, t_j)\) for
\(i\) between
\(1\) and
\(m\) and
\(j\) between
\(1\) and
\(n\text{.}\) These rectangles all have equal area
\(\Delta A = \Delta s \cdot \Delta t\text{.}\)
Each of these rectangles from the
\(st\)-domain will split our curved surface into a
\(m\) by
\(n\) grid of pieces, but these pieces of the curved surface will not all be the same size (in surface area). Looking at the surface shown in
FigureΒ 13.10.9, we can see the traces given for five equally spaced values of
\(s\) and
\(t\) breaks our curved surface into a four by four grid of pieces. Note that these pieces of the curved surface do not have the same surface area. So we will need to estimate the surface area for each particular piece of our grid using our tools from the parameterization.
Each of the pieces of the curved surface correspond to increasing
\(s\) by a small amount
\(\Delta s\) or increasing
\(t\) by a small amount
\(\Delta t\) from the point
\((s_{i-1},t_{j-1})\) in the
\(st\)-parameter plane. We want to find the vectors corresponding to each of the sides of the orange parallelogram in
FigureΒ 13.10.9. In particular, the vector
\(\vr_t \Delta t\) will be the vector that corresponds to the side of the parallelogram along a constant value of
\(s\) (shown as a yellow vector in
FigureΒ 13.10.9) and the vector
\(\vr_s \Delta s\) will be the vector that corresponds to the side of the parallelogram along a constant value of
\(t\) (shown as a red vector in
FigureΒ 13.10.9).
The area of the parallelogram will then be \(S_{ij}=\vecmag{(\vr_s \Delta s) \times (\vr_t \Delta t)}\) evaluated at the point given by \((s_{i-1},t_{j-1})\text{.}\) This will be the approximation of the surface area for a piece of our grid of the curved surface. We can simplify this magnitude of a cross product to be \(S_{ij}=\vecmag{\vr_s \times \vr_t} \Delta s \Delta t\text{,}\) which means that our approximation over all of these pieces of the grid will sum to
\begin{equation}
SA\approx \sum_{i=1}^m \sum_{j=1}^n S_{ij}= \sum_{i=1}^m \sum_{j=1}^n \vecmag{\vr_s(s_{i-1},t_{j-1}) \times \vr_t(s_{i-1},t_{j-1})} \Delta s \Delta t\text{.}\tag{13.10.3}
\end{equation}
This has satisfied both the first and second steps of the classic calculus approach because we have approximated the surface area and quantified how that approximation will change with a smaller scale of steps (in
\(s\) and
\(t\)). You can use
FigureΒ 13.10.12 to geometrically see how increasing the number of pieces will give successively better approximations for the surface area of each piece.
Note that this approximation takes the form of a Riemann sum corresponding to the setup of a double integral. In particular, if we take the limit as
\(\Delta s\) and
\(\Delta t\) go to zero, our Riemann sum will be come a double integral of
\(\vecmag{\vr_s \times \vr_t}\) over the rectangular region of integration with
\(a \leq s \leq b\) and
\(c \leq t \leq d\text{,}\) which gives the following formula for calculating the surface area.