For other surfaces, the basic idea is to cut the surface up into pieces, each of which is a graph in one of the three directions. This will be described here for a parametric surface, or a collection of parametric surface pieces.
As noted in
Subsection 6.6.4,
Oriented Surfaces, a smooth level surface
\(G(x,y,z) = 0\) with unit normal
\(\N = \del G/\|\del G\|\) defined everywhere can be described as a collection of parametric surface pieces, (indeed, graphs of functions), so the result works there too.
Writing the normal as
\(\N = \vector{n_x, n_y, n_z}\text{,}\) The Implicit Function Theorem in 3D can be used to show that if
\(n_z \neq 0\) at a point of the surface
\(\vec{r}(u,v) = \vector{x(u,v), y(u,v), z(u,v)}\text{,}\) one can solve nearby for the parameters
\(u\) and
\(v\) in terms of
\(x\) and
\(y\text{,}\) and thus get
\(z\) in terms of
\(x\) and
\(y\text{:}\) this part of the surface is a graph
\(z = g(x,y)\text{,}\) as handled in the first part of this proof.
Similarly, having either \(n_x \neq 0\) or \(n_y \neq 0\) at a point allows the surface to be descrbed as a graph in one of the other two directions.
Since \(\N \neq 0\) anywhere, at least one of these three conditions holds at every point on the surface, so every part of the surface can be covered by at least one of these types of graph.
This huge collection of "graph pieces" covers the surface many times over; then one can cut the collection and the domains down to cover each point once, except for overlap along the edges of their domains.
Then the sum of the surface integrals over each piece is the total surface integral, while the sum of the boundary integrals over the pieces is
the integral around the boundary of the whole surface, plus
various path integrals along the internal edges produced by cutting.
These internal edge path integrals come in pairs going in opposite directions (along the "opposite sides" of each cut), and so cancel out, leaving just that outside boundary integral, as needed.