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

The Approximation of a Solid of Revolution

Most math teachers I've had have been able to break down Calculus into two very broad categories: derivatives and integrals. What is truly amazing, is how much you can do with these two tools. By using integration, it is possible to approximate the shape of a 2-D function that is rotated around an axis. This solid created from the rotation is known as a solid of revolution. To explain this concept, we will take a look at the region bounded by the two functions: \[ f(x) = 2^{.25x} - 1 \] and \[ g(x) = e^{.25x} - 1 \] bounded at the line y = 1. This region is meant to represent a cross section of a small bowl. While it may not perfectly represent this practical object, the approximation will be quite textured, and will provide insight into how the process works. The region bounded by the two functions can be rotated around the y-axis to create a fully solid object. This is easy enough to talk about, but what exactly does this new solid look like? Is...

Solids of Revolution Revisited

Introduction In my previous blog post on solids of revolution, we looked at the object formed by rotating the area between \(f(x) = -\frac{1}{9}x^2+\frac{3}{4} \) and \( g(x)=\frac{1}{2}-\frac{1}{2}e^{-x} \) around the \( x \) axis and bounded by \( x = 0 \) and \( x = 1.5 \). When this solid is approximated using 10 washers, the resulting object looks like this: When I was looking back over the 3D prints I’d created for this course, I noticed that the print for this example was the least interesting of the bunch. Looking at the print now, I feel like the shape is rather uninteresting. The curve I chose has such a gradual slope that each of the washers are fairly similar in size and causes the overall shape to just look like a cylinder. Since calculating the changes in the radiuses of the washers is a big part of the washer method, I don’t think this slowly decreasing curve was the best choice to illustrate the concept. The reason I had done this o...

Finding the Center of Mass of a Toy Boat

Consider two people who visit the gym a substantial amount. One is a girl who loves to lift weights and bench press as much as she possibly can. The other is a guy who focuses much more on his legs, trying to break the world record for squat weight. It just so happens that these two are the same height and have the exact same weight, but the center of their weight is not in the same part of their body. This is because the girl has much more weight in the top half of her body and the boy has more weight in the bottom half. This difference in center of mass is a direct result of the different distributions of mass throughout both of their bodies. Moments and Mass There are two main components to finding the center of mass of an object. The first, unsurprisingly, is the mass of the whole object. In this case of the boat example, the mass will be uniform throughout the entire object. This is ideal a majority of the time as it drastically reduces the difficulty...