Scribd
Upload
64 views4 pages

Volume Under SurfacesEXP - 3B

This document discusses calculating volumes under surfaces using double integrals in MATLAB. It provides: 1) The definitions and formulas for calculating volumes under surfaces using double integrals over a region R. 2) Examples of calculating volumes under various surfaces using MATLAB code and visualizing the surfaces. 3) Exercises for students to calculate additional volumes under surfaces.

Uploaded by

vishwadeep
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
64 views4 pages

Volume Under SurfacesEXP - 3B

This document discusses calculating volumes under surfaces using double integrals in MATLAB. It provides: 1) The definitions and formulas for calculating volumes under surfaces using double integrals over a region R. 2) Examples of calculating volumes under various surfaces using MATLAB code and visualizing the surfaces. 3) Exercises for students to calculate additional volumes under surfaces.

Uploaded by

vishwadeep
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 4

Department of Mathematics

School of Advanced Sciences


BMAT 101P – Calculus (MATLAB)
Experiment 3–B
Evaluating Volume under surfaces

In this experiment, we consider a continuous function f such that f ( x, y)  0 for all


( x, y) in a region R in the xy  plane, then the volume of the solid region that lies above R
and below the graph of f is defined as the double integral V   f ( x, y )dA , where R is the
R

region bounded by the curves y  1 ( x ) and y   2 ( x) between x  a and x  b .

In this case, the inner integration is with respect to y and outer integration is with respect to
x . Hence

b 2 ( x )
V   f ( x, y)dA    f ( x, y)dy]dx
[
R x  a y 1 ( x )

MATLAB Syntax

int(int(f, y, phi1, phi2),x,a,b) where y is the inner variable, x is the outer


variable.

When R is a region bounded by the curves x   1 ( y ) and x   2 ( y ) between y  c and


y  d , i.e., the inner integration is with respect to x and outer integration is with respect to
d  2 ( y)
y . Then V   f ( x, y)dA   [  f ( x, y)dx]dy
R y c x  1 ( y )

MATLAB Syntax

int(int(f, x, psi1, psi2),y,c,d) where x is the inner variable, y is the outer


variable.

Supporting files required:

To visualize the surfaces two additional m-files viz., viewSolid.m, viewSolidone.m


are required. These files are to be included in the current working directory before execution.
Students are advised to upload these files (viewSolid.m and viewSolidone.m) to
their MATLAB drive. These supporting files should not be edited.

Download viewSolid.m from the following link:


https://drive.google.com/file/d/1qEsq7VCgrmI60GI-C0bMY6kl8yRWBznk/view?usp=sharing

Download viewSolidone.m from the following link:


https://drive.google.com/file/d/1H1cJOfJArmUujQNVxeSeuGfBJejtbGgJ/view?usp=sharing

Department of Mathematics, Vellore Institute of Technology, Vellore. Page 1


Syntax for visualization of the surfaces:

viewSolid(z,0+0*x+0*y,f,y,phi1,phi2,x,a,b)

viewSolidone (z,0+0*x+0*y,f,x,psi1,psi2,y,c,d)

x y
2 x
Example. 1 To find 
1 x/2
4
dydx .

syms x y z
int(int((x+y)/4,y,x/2,x),x,1,2)
viewSolid(z,0+0*x+0*y,(x+y)/4,y,x/2,x,x,1,2)

Output
ans =
49/96

In this figure the required volume is above


the plane z=0 (shown in red) and above the
x y
surface z  (shown in green) .
4

Example. 2 To find the volume of the prism whose base is the triangle in the xy  plane
bounded by the x  axis and the lines y  x and x  1 and whose top lies in the plane
z  f ( x, y)  3  x  y . The limits of integration here are y  0 to 1 while x  y to 1 .
1 1
Hence  (3  x  y)dA    (3  x  y)dxdy
R 0 y

syms x y z
int(int(3-x-y,x,y,1),y,0,1)
viewSolidone(z,0+0*x+0*y,3-x-y,x,y,1,y,0,1)

Output:
ans =
1

In this figure the triangular region on the


xy plane is shown in red, while the plane
surface z=3–x–y above the xy plane is
shown in green.

Department of Mathematics, Vellore Institute of Technology, Vellore. Page 2


2 2x
Example 3 Evaluate the integral   (4 x  2)dydx by changing the order of integration.
0 x2

As per the given limits of integration x  0 to 2 while y  x 2 to 2 x .

MATLAB Code:
syms x y z
int(int((4*x+2),y,x^2,2*x),x,0,2)
viewSolid(z,0+0*x+0*y, 4*x+2,y,x^2,2*x,x,0,2)

Output
ans =

By changing the order of integration, the limits are


y
y  0 to 4 while x  to y .
2

int(int(4*x+2,x,y/2,sqrt(y)),y,0,4)

ans =

Example 4: Evaluate  ( x  3 y
2
)dA where R  ( x, y ) | 0  x  2, 1  y  2
R

clc
clear all
syms x y z
viewSolid(z, 0+0*x+0*y, x-3*y^2,y,1+0*x, 2+0*x,x,0,2)
int(int(x-3*y^2,y,1,2),x,0,2)

Output:

>> ans
-12

In this figure the required volume is below the


plane z  0 (shown in red) and above the
surface z  x  3y 2 (shown in blue). The reason
why the answer is negative is that the surface
z  x  3y 2 is below z  0 for the given
domain of integration.

Department of Mathematics, Vellore Institute of Technology, Vellore. Page 3


Example 5: Evaluate  y sin( xy) dA where R  [1, 2]  [0,  ]
R

MATLAB Code:

syms x y z
viewSolidone(z, 0+0*x+0*y, y*sin(x*y),x,1+0*y, 2+0*y,y,0,pi)
int(int(y*sin(x*y),x,1,2),y,0,pi)

Output:
>> ans
0

For a function f ( x, y) that takes on both


positive and negative values ∬ is a
difference of volumes V1-V2, V1 is the volume
above R and below the graph of f and V2 is the
volume below R and above the graph. The
integral in this example is 0 means V1=V2

Example 6: Find the volume of the solid that lies under the paraboloid z  x 2  y 2 and above
the region D in the xy-plane bounded by the lines y  2 x and the parabola y  x 2

MATLAB Code:

syms x y z
viewSolidone(z, 0+0*x+0*y, x^2+y^2,x,y/2,sqrt(y),y,0,4)
int(int(x^2+y^2,x,y/2,sqrt(y)),y,0,4)

Output:

>> ans
216/35

Exercise:

1. Find the volume of the solid S that is bounded by the elliptic paraboloid
x 2  2 y 2  z  16 , the planes x  2 and y  2 , and the three coordinate planes.
2. Evaluate  sin x cos y dA where R  [0,  / 2]  [0,  / 2]
R

Department of Mathematics, Vellore Institute of Technology, Vellore. Page 4

You might also like

pFad - Phonifier reborn

Pfad - The Proxy pFad of © 2024 Garber Painting. All rights reserved.

Note: This service is not intended for secure transactions such as banking, social media, email, or purchasing. Use at your own risk. We assume no liability whatsoever for broken pages.


Alternative Proxies:

Alternative Proxy

pFad Proxy

pFad v3 Proxy

pFad v4 Proxy