Skip to main content

Double Integration over Rectangles over a Non-Rectangular Region

Double Integrals over Rectangles over a Non-Rectangular Region

Previously we have discussed methods of approximating the volume of a solid with regard to solids of known cross section and solids that are symmetric about an axis, via shells and washers. We are now discussing ways to approximate the volumes of more irregular solids—i.e. those without known cross section and or are not symmetric about an axis. We will first discuss the solid we will be approximating, followed by the method of approximation and calculating an actual volume by double integrals, and conclude with actual calculations of volume for the solid in question.

The solid that we will be considering in bounded on the \(xy\)-plane by the function \(y = \left| x \right| \) and \( y = 5 \); which results in the region shown below:

Vertically, or with respect to the \(z\)-axis, the solid is bounded by \( f(x,y) = 2y - x^2 \) which is shown in the following 3D image:

Methods

When approximating the volume of a solid of unknown cross section that is not symmetric about an axis one method is to divide the region in which the solid exists into a region of many rectangles, then by setting the height of the rectangle as an output of the function that bounds the top of our solid. In this way we can find the volume of this rectangular prism which in turn approximates solid's volume within that particular rectangular region. So, if we take one rectangle that totally encompasses the region where our solid lives we get a very rough approximation of its volume, more like an upper bound. For example, any solid object placed inside of a box must have volume less than or equal to that of the box. Now as we divide the region into more and more rectangles and subsequent rectangular prisms our approximation becomes more and more accurate as there is less "extra" volume in each prism that is not also in our solid; and if we increase these prisms infinitely we get an infinitely more accurate approximation. This infinitely more accurate approximation is actually the correct volume of said solid that is given by double integration as it essentially yields the volume of a very small prism, that could more precisely be referred to as a line, over every point that lies within the bounded region that makes up the base of our solid.

Another consideration is that if the height of these prisms are an output of the function that bounds the top of our solid, which actual point gives the output that we will use as our height? And simply you can use any point within the rectangle that forms the base of the rectangular prism as the one that gives the height. I have chosen the upper-righthand vertex for my approximation, and would recommend that you choose one of the four vertices, the center point, or the midpoint of one of the sides.

I decided to use 1x1 rectangles, which results in 30 rectangles of which 23 have heights, and thus volumes, that are of interest to use. The approximation of the solid is shown below:

Calculating the actual volume of our solids via integration requires the use of double integrals but is really rather simple. We define the region on the \(xy\)-plane within which our solid lives as \(R\) where \(x\) is bounded by \(a\) and \(b\) and \(y\) is bounded by \(g_1(x)\) and \(g_2(x)\) where \(g_1(x)\) is the lower-bound on \(R\) and \(g_2(x)\) is the upper-bound of \(R\) and \(f(x,y)\) is the function that bounds the top of our solid; or \[ R = { (x,y) \mid a \leq x \leq b, g_1(x) \leq y \leq g_2(x) } \\ V = \int_{R} \int f(x,y) dA = \int_{a}^b \int_{g_1(x)}^{g_2(x)} f(x,y) dy dx. \]

Calculations and Solutions

As stated above, I divided the region our solid lives in into 30 1x1 rectangles, of these 30 rectangles 23 of them had heights greater than zero, those with heights less than or equal to zero were ignored as they lay below or right on the plane and were thus not part of our volumes. Those rectangles and their sum, which is our approximate volume are shown in the table below:

Note that the \(z\) coordinates highlighted in green denote rectangles that actually contribute towards our approximations volume. As you can see from the table our approximation yielded a volume of 106 units\(^2\). Now we can see how good our approximation is by preforming the double integration— describes in the "Methods" section. This integration is setup can completed as follows: \[ V = \int_{-5}^5 \int_{\left| x \right|}^5 2y-x^2 dydx \\ = \int_{-5}^5 (y^2-x^2y) \mid_{\left| x \right|}^5 dx \\ \approx 62.5. \] Obviously, our approximation and volume via the integral are not that similar, in fact our actual volume from our double integration is only about 59% of our approximated volume. So, now let's see if using more rectangles gets us a more accurate approximate volume. These values and subsequent volume can be found in the table below, also with the green marking rectangles we care about:

Now our approximation has gone the other way and is actually smaller than our actual volume. This is due to volume that is lost when the upper-righthand vertex lies inside the top curve such that some volume is lost, but now our approximate volume is 85% of our actual volume which is a fairly precise approximation.

Comments

Popular posts from this blog

Do Over: Double Integrals over Regions

Introduction Over the semester we've looked at many topics and created 3D models. For this we are going to revisit an old topic, double integrals over a region. In this we found the volume of a surface in the xyz-plane bounded by two curves. From the many topics I chose to revisit this topic. I have a couple reason to why I chose to redo this. First, the model did not print correctly. The print added spaces between the rectangular prisms. Another reason was that I think the surface and curves did not represent the topic entirely. The surface I chose just increased between the curves. Improvements When making the model on Onshape there were no spaces between the rectangles, which can be seen on the right. However, when printing this spaces were being added. The second issue was with the surface I chose which was \(f(x,y)=xy+x\). This function only increased over the two curves I chose \begin{align*} f(x) &= \sqrt{x} ...

The Septoil Knot

Knots are a very interesting topic and a field that has not quite been fully discovered, so mathematicians are still discovering new ideas and invariances about knots even today. While it may seem like knots are a simple skill you learn at camp, they actually have a lot of mathematical properties and in this blog post we are going to look at just a few. By mathematical definition, a knot is a closed curve in three dimensional space that does not intersect itself. Since we are working with three dimensional space and you are reading this on a two dimensional screen, we need a way to look at knots in two dimensions and that is where knot projections come in. A knot projection is simply a picture of a knot in two dimensions and where a knot crosses itself in the projection is simply a crossing of that projection. The number of crossings of a knot is the smallest number of crossings among all projections of a knot. Since a knot is not necessarily solid, one...

Knot 10-84

Introduction In mathematics, a knot is simply a closed loop. The simplest version of this is the unknot, which is a just a closed circle (imagine a ponytail holder). Knots, however, quickly become more complicated than this more basic example. This post will examine a particular knot (knot 10-84) and a few of its knot invariants. Crossing Number Knots are often defined by their crossing number, which is the number of times the knot’s strands cross each other. As indicated in its name, knot 10-84 is a 10 crossing knot. In order to visualize the knot, we can look at its knot projection, in which the knot is represented by a line segment broken only at its undercrossings: Tricolorability Now that we’ve looked at knot crossings, we will examine a potential property of knots: tricolorability. In order to understand tricolorability, it is first important to know that one strand of a knot is defined as an unbroken line segment in the knot p...