All published worksheets from http://sagenb.org
Image: ubuntu2004
Find the flux of the vector field through the cylinder ,
We use the general form of the cylinder: , so and
We will integrate over the surfaces: Part 1: Cylinder + Part 2: Bottom + Part3: Top. Then we will check our answer with the Part 4: Divergence Theorem (aka Gauss-Ostrogradsky).
Surface 2 is the top of the cylinder S. We want
As always, we ran the SOLVER first (below) and then used that information to get the graphs here.
Graph 1: We use the parameterization of S that we get below to graph the surface. We split it up so we could "see" the orientation.
Graph 2: Our integrand from below is: . We translate this to:
Click and drag to rotate. Notice that the arrows are not all the same length. We scale the vectors by vs.
As always, to get the arrows, we use the parameterization of S to get the "base point" and then use "base point" plus ( Integrand) to get the "end point".
Preparation for using the SOLVER.
We parameterize the surface with 2 parameters and find a unit vector normal and pointing out from the surface.
Work particular to this problem:
== Remember that we parameterize a disk with and . We set .
== Remember we have and .
== So:
== By examination, we see: (That is, the vector pointing straight up.)
SOLVER
Step 1: We define , and . Changes with problem; input your parametrization of Sf and your vector n.
Step 2: The program finds the partial derivatives.
Step 3: The program checks whether our parameterization is "orientation-preserving" or not. Check that we get 1 or -1.
If we get 0 here, we need to change the point, e.g. (u=2.0,v=2.0). Use "real" values with decimal points. (If we get something other than +1, -1 or 0, we have made an error in step 1.)
Step 4: The program defines a standard sage function for "change of variables" for a 2-variable parametrization (surface).
The program changes the variables of .
Step 5: Then program finds the mixed product and multiplies it by the orientation from step 3.
Step 6: Then program computes the integral (flux). Changes with problem - we must put in our intervals of integration.
So our flux over the cylinder is: . For general constants a and b, we have
Now we must do the integral over the bottom of the cylinder, but this is easy as we need only change z=b with z=0 and . Then we will check our answer with the divergence (Gauss') theorem.