]>
Home | 18.013A | Chapter 25 |
||
|
When a surface is described in parametric form, that is, you are given functions and , and a region of values for and , you can use a spreadsheet to obtain the flux integral of any spreadsheet definable vector function (that is the dot product of with the normal to the surface) over it.
You can do so by exploiting the ability of the spreadsheet to copy a single instruction to a large array of locations in two dimensions, row and column.
Let us suppose you want to integrate over first and then , and you will integrate from to , and from to . Suppose you want to break each range into pieces.
The plan for this is as follows: you set up arrays of values, values and then and values and then and values, and finally an array for the integral; each of these arrays will be by in size, but can be constructed by constructing at most 3 or 4 entries. The rest is all copying.
Set up an array of values, assigning to the first row that you use in the first column, and in the next row beneath it, and beneath that. Copy that into the next 98 rows. Then set the next column entry to "= entry to left" ; copy that into the next 98 columns and rows.
This should give a single constant value in each row, but different ones in each column.
Next set up two columns that compute and for each value.
Set up an array of values by putting in the first column, (use dollar signs on the ) in the next column, and in the next column.
Then copy this down 98 rows, and copy the third column 98 rows across.
Set up the and and arrays again of the same size, by putting in the top left entry of the array, and copying it into the entire array. This can be done similarly for the other 5 arrays.
To compute the integral you can add up the contribution from each "pair of intervals" and , that is
To compute this for every 2x2 rectangle in the arrays again you need only compute it for the top left one and copy again to a 100 by 100 array.
This can be done in one step or once for each term. The integral will be the sum of this over the entire array.
What do you do to get or or for a pair of intervals?
Suppose the interval corners are and .
Then the analogue of the trapezoid rule in two dimensions gives
Though it is messy to create this last array, it really only has to be done in one square and copied, and also only has to be done once, and can be reused with different integrals, if the other arrays have the same start points.
Notice that to change the domain for and you need only change the first two arrays; to change the parametrization you need only change the and arrays, and to change the vector whose flux integral you are computing you need only change its three arrays.
You can get integrals over subsets of the surface by taking different sums over the last array.
Here is an example
Notice that changing or each require changing only one entry, and then copying it into the entire corresponding array. (Or wherever appropriate. In the case of and you need copy nothing, and for and that is a column.)
It is quite easy to extrapolate these results: you can change to 50 and sum over the top left quarter of the last array, and change to 25 and sum over the top left quarter of that array.
The results for this integral are 132.2450414 for , 132.1941481 for , and 131.9909671 for .
The difference between these is roughly 4 times smaller between 50 and 100 than it is between 25 and 50. Multiplying the better by and subtracting the worse multiplied by gives 132.262006 and 132.261875. Extrapolating assuming that the remaining error goes down by a factor of 16 gives an estimate of 132.2620145, which is probably accurate to 6 decimal places.
This can be verified by enlarging the arrays and trying , something which requires only copying, and forming one new sum.
Exercises:
25.6 Create the spreadsheet just described and verify or disprove the claims above.
25.7 What is the best estimate you can get for the integral described?
25.8 Change and to go from 0 to and , make and equal what they are in spherical coordinates, and let be the vector . Find the integral by this approach to some reasonable accuracy.
|