The definite integral of a function of a single variable is defined to be the limit over all applicable partitions of a Riemann sum. In imitation of this, the double integral of a function of two variables is defined to be the (bivariate) limit of an appropriate Riemann sum.
In two dimensions, where the domain of integration is rectangular in shape and oriented so the sides are parallel to the coordinate axes, the horizontal and vertical edges are partitioned so that a grid of sub-rectangles covers the domain. The function defined over the domain is evaluated at one point in each sub-rectangle, and this value is multiplied by the area of the corresponding sub-rectangle. The Riemann sum is the sum of all the volumes of the rectangular "blocks" so determined. The limit as the number of blocks becomes infinite is, under the appropriate conditions, taken as the meaning of the double integral.
Figure 5.1.2(a) shows the surface graphed over the rectangle . It is the volume beneath this surface that is to be computed by the double integral.
|
Figure 5.1.2(a) Surface
|
|
|
|
>
|
use plots,plottools in
module()
local p1,p2,p3,p4,p5,i,j;
p1:=rectangle([-1,1/2],[1,-1/2],style=line,thickness=3);
p2:=pointplot([seq(seq([-.8+.4*i,-.4+.2*j],i=0..4),j=0..4)],symbol=solidcircle,symbolsize=15,color=red);
p3:=plot([-.3,-.1,.1,.3],x=-1..1,color=green);
p4:=implicitplot([x=-.6,x=-.2,x=.2,x=.6],x=-1..1,y=-1/2..1/2,color=green);
p5:=display(p1,p2,p3,p4,scaling=constrained,labels=[x,y],tickmarks=[3,3],axes=frame);
print(p5);
end module:
end use:
|
|
Figure 5.1.2(b) Grid with midpoint evaluation
|
|
|
|
|
Figure 5.1.2(c) Volume-approximating blocks
|
|
|
|
|
|
Figure 5.1.2(b) shows a grid for which a Riemann sum is constructed. The red dot at the center of each sub-rectangle is the point at which the height of the surface is sampled. Figure 5.1.2(c) shows the "blocks" whose volumes are analogous to the rectangles that approximate area in a Riemann sum for a single integral.
Computationally, the content of Figures 5.1.2(a - c) can be implemented once the coordinates of the "red dots" in Figure 5.1.2(b) are specified as . The sum of the volumes of the 25 blocks so determined is
Passage to the limit as the density of the grid increases requires discretizing so that and , in which case the red-dot midpoints become
and the approximating sum becomes
which, in the limit as and become infinite, approaches as a limiting value.
A simpler Riemann sum can be constructed if the evaluation point for each sub-rectangle is taken as the "lower-left" corner where the generic coordinate would be and the approximating sum would be
Again, the limiting value as and become infinite is . Finally, if the evaluation point in each sub-rectangle is taken as the "lower-right" corner where the generic coordinate would still be , the approximating sum would be
with limiting value .