# Finding volumes between a quadric surface and a plane.

1 Introduction

In this paper, we give descriptions of finding the volume bounded by a quadric surface and a plane. Unless otherwise is specified, we only discuss the cases where there is a finite volume for the bounded solid. In particular, we will focus on the surfaces of an ellipsoid, a hyperboloid of two sheets and a paraboloid. We usually start with a general introduction of a method and lead to a particular example. In this paper, we give several analytical methods using various coordinate systems where most students learn in a multivariable Calculus class. We often see the needs of change of variables and change of bases to transform an arbitrary quadric surface into a standard form  before applying integration techniques, the methods described here are accessible to those who have learned Linear Algebra and multivariable Calculus. We assume readers are familiar with the Divergence and Stokes' Theorems and we will also describe how we may apply these two principles to find the volumes of bounded solid.

We use the following five coordinate systems:

--initial coordinate system X,in which the quadric and the plane are given;

--canonical coordinate system X', used principal axis in which the quadric is in standard form [1, Table 3.5-2, 3.5-7];

--"natural" coordinate system X", in which the origin is located in the center of a cross section of the quadric and the plane and the axes are oriented along the axes of the section;

--"harmonized" system of spherical coordinate ([rho],[??],[[phi]]) connected with natural coordinate system X". The x-axis for "harmonized" system is of [??] = 0, the y-axis is of

[phi] = 0,[??]= [pi]/2, and z-axis is of [phi] = [pi]/n, [??] = [pi]/2;

--two-dimensional coordinates ([theta], [??]), used for parametric representations of quadrics [1, 3.5-22, 3.5-24,3.5-25].

Each of the volumes is calculated by the following six methods. Method 1. Integration by quadratures using the coordinates X".

Method 2. Numerical integration of the distance from point on quadric to the plane using the coordinates X".

Method 3. Numerical integration of the distance from point on quadric to the origin of "natural" coordinate system X" with using "harmonized" system of spherical coordinate {[rho], [??],[phi]).

Method 4. Numerical integration using the Divergence Theorem in the coordinates ([theta], [phi]).

Method 5. Numerical integration with the use of Stokes' Theorem in the coordinates ([theta], [phi]).

Method 6. Integration with the use of discontinuous integrand in the coordinates ([theta], [phi]).

2 General Description of a Quadric surface and a Plane

Since our main focus is on finding the volume bounded by a quadric surface and a plane. We first consider the equation of the quadric surface in a matrix form:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (1)

where [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] are real numbers, i,j 7=0,1,2,3, and [a.sub.11]>0. Alternatively, we sometimes write a quadric surface in the scalar form of:

f (x ,y ,z) = [a.sub.11][x.sup.2] +[a.sub.22][y.sup.2] +[a.sub.33][z.sup.2]+2[a.sub.12]xy +2[a.sub.13]xz +2[a.sub.23]yz +2[a.sub.10]x +2[a.sub.20]y +2[a.sub.30]z+[a.sub.00] = 0. (2)

Also, we consider the equation of the plane in the forms of

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (3)

where [??]=([N.sub.1],[N.sub.2],[N.sub.3]) is a normal vector to the plane, [N.sub.0] is a constant, [??]= ([n.sub.1] ,[n.sub.2] ,[n.sub.3]), and

[n.sub.i] = [N.sub.i] / [square root of [N.sup.2.sub.1], [N.sup.2.sub.3, [N.sup.2.sub.3]

2.1. Transformation of the quadrics equations. We are interested in those solids bounded by quadric surfaces and a plane produce finite volumes. The quadric surface we discuss in this paper will be an ellipsoid, an elliptic paraboloid, a hyperboloid of two sheets or a cone; however, we will leave it to reader to explore the case when quadric surface is a cone. We classify the quadric surfaces by using the following scenarios:

Case 1. [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] and

where [[??].sub.3] is the normalize eigenvector corresponds to [[lambda].sub.3] = 0. Then the quadric is an elliptic paraboloid. Indeed, we have

[[lambda].sub.1] = [I+[square root of [I.sup.2]-4j]]/2, [[lambda].sub.2] = [I-[square root of [I.sup.2]-4J]]/2 [[lambda].sub.3] =0, (4)

and the corresponding eigenvectors to [[lambda].sub.1] and [[lambda].sub.2], in case of three different eigenvalues, are

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (5)

respectively. If we have a repeated eigenvalue, we use (5) for the eigenvector [[??].sub.1] associated with the single eigenvalue, and the vector [[??].sub.2] can be choose to be any vector perpendicular to [[??].sub.1]. If one of these vectors is 0, we rename corresponding [lambda] as [[lambda].sub.2]. We remind the readers to make necessary renumbering of vectors in various special cases. We normalize [[??].sub.1] and write [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] for i = 1,2. We set [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] and form the matrix R=([[??].sub.1], [[??].sub.1]2, [[??].sub.1]). If we set

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

then the substitution [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] will transform the original quadric into

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (6)

We note that the signs of [[lambda].sub.1] and [[lambda].sub.2] are the same since J=[[lambda].sub.1] [[lambda].sub.2]>0, which gives us an alternative way of categorizing the quadric to be an elliptic paraboloid. We denote [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. If the sign of [??]-[[??].sub.3] coincides with the sign of [[lambda].sub.1], we change the direction of the axis z' to the opposite direction. As a result, we obtain the equation of the paraboloid in the form of [(x').sup.2]/[a.sup.2] + [(y').sup.2]/[b.sup.2] = z'. In matrix form we obtain the following equation: a b

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] where [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], and [a'.sub.00]=0. (7)

Case 2. Suppose [DELTA][not equal to 0] and the signs of [[lambda].sub.1], [[lambda].sub.2] and [[lambda].sub.3] are the same and they are different from the sign of a [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. Then such a quadric corresponds to an ellipsoid. In particular, the substitution [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] will transform the original ellipsoid into

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (8)

We denote [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. As a result, we

[(x').sup.2]/[a.sup.2] + [(y').sup.2]/[b.sup.2] + [(z').sup.2]/[c.sup.2]= 1

obtain the equation in the form of [(x').sup.2]/[a.sup.2] + [(y').sup.2]/[b.sup.2] + [(z').sup.2]/[c.sup.2] = 1. In matrix form we obtain the following equation:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (9)

Case 3. Suppose [lambda] [not equal to 0] and not all of the signs of [[lambda].sub.1], [[lambda].sub.2] and [[lambda].sub.3] are the same, among the numbers [[lambda].sub.1], [[lambda].sub.2] and [[lambda].sub.3] and [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] , one number differs in sign from the other three numbers. Such quadric corresponds to a hyperboloid of two sheets. We assume the sign of [[lambda].sub.1] differs from the signs of [[lambda].sub.2], [[lambda].sub.3], and [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. We name this situation as sub-case 1. The substitution [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] , and [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] will transform equation (8) into the form of [[xi].sup.2.sub.1] [[xi].sup.2.sub.2]-[[xi].sup.2.sub.3]= 1; that is, to the canonical equation of the hyperboloid of two sheets. In matrix form we reach the following equation:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (10)

Exercise. Obtain the same results in other sub-cases within case 3.

The other cases beyond cases 1,2 and 3 will not produce finite volumes for the solids, which we will not discuss here.

2.2. Equation for the curve of intersection - natural coordinates. Let us show that the intersection curve of the quadric and the plane is an ellipse. The equation of the quadric surface, we use in the form of (1), is equivalent to (7), (9) and (10). The equation of the plane will be written in the form of [??]x[??]'+[m.sub.0]=0, and [parallel][??][parallel] = 1.

2.2.0. We adopt the coordinate system in which the equation of the cross section has the forms of x" = 0 and [a.sub.23]y"z" = 0. Let the origin of the coordinate system x",y",z", [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] be at point [[??].sub.0] which will be defined later, the basis in this coordinate system

contains vectors [??],[??], and [??]X[??]. The corresponding transformation is [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]). Vector [??] lies on the plane of the cross section, and it will be defined explicitly later.

The substitution in the equation (1) yields

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (11)

and the substitution in the plane yields [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. The element of the matrix [r.sup.T] A'r with the (2,3) position, second row and third column of [r.sup.T]A'r, is 0 or [([r.sup.T]A'r).sub.23] =[[??].sup.T]A'[[??]X[??]] = 0, since this determines the coefficient of the product y"z".

We find the basis vector [??] by solving the following simultaneous equations:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (12)

Each of the roots of the quadratic equation (12) corresponds to one of the axes of the quadric surface, we may use either solution for further investigation. These solutions corresponds to [??] and [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. We get the components of [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] by using the following equations:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

[U.sub.x] = 1+[m.sub.2][U.sub.y]/[m.sub.1], and [U.sub.z] = 1/[m.sub.3]. We next normalize the vector [??].

With information from the vector J, it allows us to define the coefficient of [(y").sup.2] to be [l.sub.2] = [([r.sup.T]A'r).sup.22] =[[??].sup.T] A'[??] and the coefficient of [(z").sup.2] to be [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. We define the origin of coordinates [[??].sub.0] to be on the plane of the cross section [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. We make the coefficients of y" and of z" equal to zero. That is [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], and [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. Then [[??].sub.0] defines the center for the ellipse of the cross section. The constant term of the equation (11) in the coordinates x",y" ,z" is equal to [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. The semi-axes of the ellipse are respectively

[a.sub.e] = [square root of [-a.sub.0]"/[l.sub.2] and be = [b.sub.e] = [square root of [-a.sub.0]"/[l.sub.3]. (13)

Hence, the cross section is an ellipse if and only if the signs of the numbers [l.sub.2] and [l.sub.3] are the same and differ from the sign of [a.sub.0]".

Let us write the equations defining the intersection curve in each of the cases. For convenience, we use the point [??] to denote the point where the tangent plane to the quadric is parallel to the given plane; obtaining such [??] is possible due to . We have one such point for the paraboloid. We have two such points for the ellipsoid and the hyperboloid of two sheets, and the points are symmetrical about the origin. Simple calculations yield the following results.

2.2.1. We use the substitution [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] for the paraboloid in the equation (2) and yields

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (14)

The center of the ellipse [[??].sub.0] and the point [??] have the following coordinates respectively:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (15)

The semi-axes of the ellipse are [a.sub.e] and [b.sub.e] respectively: [a.sub.e] = [square root of [-a.sub.0]"/ [l.sub.2] and be = [b.sub.e] = [square root of [-a.sub.0]"/[l.sub.3] (16)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (16)

The area of the ellipse is [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (17)

The coordinates of any point of the quadric in the solid of intersection have the following standard vector form:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], where t [member of] [0, [Z.sub.A] - [Z.sub.t]]. (18)

2.2.2. For the ellipsoid, the substitution [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] in the equation (2) gives

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (19)

The center of the ellipse [[??].sub.0] and the point [??] have the following coordinates:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII],

[x'.sub.T]= [m.sub.1][a.sup.2]/w, [z'.sub.T]= [m.sub.2][b.sup.2]/w, [z'.sub.T]= [m.sub.3][c.sup.2]/w, and [??]'= [??]. (20)

The semi-axes of the ellipse are [a.sub.e] and [b.sub.e] respectively: [a.sub.e] = [a.sub.e] = [square root of [-a.sub.0]"/ [l.sub.2] and be = [b.sub.e] = [square root of [-a.sub.0]"/[l.sub.3],

where [-a".sub.0] = 1-[m.sup.0.sub.2]/[w.sup.2],

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (21)

The area of the ellipse is

[S.sub.e] ([[mu].sub.0]) = [pi][a.sub.e][b.sub.e] = [-a.sub.0]"[pi]/[square root of [l.sub.2][l.sub.3] = [pi] abc ([m.sup.2.sub.1] + [m.sup.2.sub.2] + [m.sup.2.sub.3]/w (1-[[mu].sup.2.sub.0], where [[mu].sub.0] = [m.sub.0]/w. (22)

The coordinates of any point of the quadric in the solid of intersection have the following standard vector form:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (23)

2.2.3. For the hyperboloid of two sheets, if the sign of [[lambda].sub.1] differs from the signs of [[lambda].sub.2], [[lambda].sub.3] and [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], then the substitution [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] in the equation (2) gives equation (19).

The center of the ellipse, [[??].sub.0] and the point T have the following coordinates:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (24)

The semi-axes of the ellipse are [a.sub.e] and [b.sub.e] respectively: [a.sub.e] = [square root of [-a.sub.0]"/[-l.sub.2] and be = [b.sub.e] = [square root of [-a.sub.0]"/[-l.sub.3], where

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (25)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (26)

The coordinates of any point of the quadric in the solid of intersection have standard vector form:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], and [mu] [member of] [[[mu].sub.0]-1]. (27)

The other cases will not produce finite cross section, which we will not discuss here.

Remarks. We make the following observations in view of the equations (22-27):

1. A set the of the cross sections of the given quadric by parallel planes is a family of the similar ellipses.

2. We may say that:

a) a diameter of the paraboloid is a beam, parallel to its axis and beginning at point [??];

b) a diameter of the ellipsoid is a segment passing through its center and connecting points [??]' and

T;

c) a diameter of the hyperboloid of two sheets is a beam, passing through its center and outgoing into the infinity from points [??]' and [??] and lying on the straight line [??][??]'.

The diameter contains the centers of the ellipses obtained from the intersection of the quadric and the family of planes parallel to the given plane, which touches the quadric at the point T. The Figure 1 illustrates the statements 1 and 2. It shows the solids cut off by the plane from the paraboloid (Fig. 1a), the ellipsoid (Fig. 1b) and the hyperboloid of two sheets (Fig.1c). By varying the parameter [m.sub.0], we construct family of parallel planes and observe the center of the elliptical cross-section slides in the direction along the diameter.

2.3. Equation for the curve of intersection-classic coordinates

In the description of the quadric surfaces, we use the canonical variables X' as initial variables and rewrite the equation in the internal variables for the surface.

2.3.1. For the paraboloid described by the equations (7) and (14), we make the substitution

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (28)

We transform the plane equation into [t.sup.2]-2[t.sub.s]t + [m.sub.0]/[m.sub.3], [t.sub.s]=-cos[??] + [y'.sub.A]/b sin [??](29)

The solution of this equation is t = [t.sub.s] [+ or -] [square root of [t.sup.2.sub.s]- [m.sub.0]/[m.sub.3]. If the vertex of the paraboloid [X.sub.0] belongs to the intersection solid, then we have a unique positive solution t = [t.sub.s] + [square root of [t.sup.2.sub.s]- [m.sub.0]/[m.sub.3]for any value of [??]. If [X.sub.0] does not belong to the solid of intersection, then there are two positive solutions t > 0, which exist only if [??] satisfying the condition of [t.sup.2.sub.s] [greater than or equal to] [m.sub.0]/[m.sub.3].

2.3.2. For the ellipsoid described by the equations (8) and (19) we make the following substitution

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (30)

and the intersection curve between the ellipsoid and the plane is

(a[m.sub.1]cos [??] + b[m.sub.2] sin[??]) sin[theta] +c[m.sub.3] cos[theta] + [m.sub.0] = 0. (31)

Under the condition [parallel][m.sub.3]c[parallel]>[parallel][m.sub.0][parallel], the above equation has a unique positive solution for [theta] when [??][member of] [0,2[pi]]

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] and [??]([phi]) = [pi]/2 arccos, [t.sub.s]/ [square root of [c.sup.2][m.sup.2.sub.3] - arccos -[m.sub.0]/ [square root of [c.sup.2][m.sup.2.sub.3], where [t.sub.s]=a[m.sub.1]cos[??]+b[m.sub.2]sin[??] (32)

2.3.3. For the hyperboloid of two sheets described by the equations (10) and (19), we make the substitution

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (33)

The curve of intersection between the hyperboloid and the plane is

(b[m.sub.2]cos[??]+c[m.sub.3]sin[??]sinht +a[m.sub.1]cosht+[m.sub.0] = 0. (34)

Under the condition -[m.sub.0]/[m.sub.1] >a, the equation (34) has a unique positive solution for t for any [phi] [member of] [0,2 [pi]]: t=ln l+[tau]/ l-[tau], where [tau] = [t.sub.s] - [square roots of [t.sup.2.sub.s] + [m.sup.2.sub.0] - [a.sup.2][m.sup.2.sub.1] / [m.sub.0]-[m.sub.1]a

and [t.sup.s] = b[m.sub.2] cos[??]+c[m.sub.3]sin [??]. (35)

3 Examples

We consider the following three randomly selected solids. Let us first consider a paraboloid whose vertex belongs to the solid of intersection.

Example 1 We consider the paraboloid given by the equation of [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], where

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] or we may write the paraboloid in the form of f (x,y,z)=5[x.sup.2]+8xy +6[y.sup.2] +14xz + 2yz + 243/14 [z.sup.2] - 4x - 6y +14z -7=0, and consider the plane in the form of g(x,y,z) = x - 3y + 2z - 2=0. Let us name the solid bounded by the quadric surface and the plane [S.sub.1]. We note that the eigenvalues for the paraboloid are 397+[square root of 36145]/28, 397-[square root of 36145]/28 and 0, respectively. The eigenvectors for the paraboloid are [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], respectively. If we write f (x,y,z) in the canonical form by using the canonical variable X', we get [(x'/2).sup.2] + [(x'/2).sup.2] = z' [approximately equal to] 0.46373, and b [approximately equal to] 0.78120. The equation of the plane when using the canonical variable X' is [m.sub.1]x' +[m.sub.2]y'+[m.sub.3]z'+[m.sub.0] = 0 where [m.sub.1] [approximately equal to] 0.4485, [m.sub.2] [approximately equal to] -0.77026, [m.sub.3] [approximately equal to] -0.4533, and [m.sub.0] [approximately equal to] - 2.51798. The center of the elliptical cross section is at the point

[[??].sub.0] =(21217/6241, - 4919 /6241 - 11746)/6241). We remark the calculations above are obtained from Maple (see ,).

Next we consider a solid formed by an ellipsoid and a plane as follows: Example 2 Consider the ellipsoid of the form [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], where

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], or we may write the surface as

f (x,y,z)= [x.sup.2] - 2xy + 5[y.sup.2] +4yz + 2[z.sup.2] - 4x + 6y - 6z + 5=0, and consider the semi-space g(x,y,z) = 3x - 4y +3z - 21 [less than or equal to] 0. We note that the eigenvectors for the ellipsoid are [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], respectively. The center of the ellipsoid is at the point [[??].sub.0] =(0, - 2, 3.5), the semi-axes for the ellipsoid are a [approximately equal to] 1.3667, b [approximately equal to] 2.8982, and c [approximately equal to] 4.9226, respectively. The equation for the ellipsoid is [(x'/a).sup.2] + [(y'/a).sup.2] + [(z'/a).sup.2]-1 =0, and for the plane is [m.sub.1]x'+[m.sub.2]y'+[m.sub.3]z'+[m.sub.0]=0, where [m.sub.1] [approximately equal to] 0.4777,[m.sub.2] [approximately equal to] 0.8654,[m.sub.3] [approximately equal to]- 0.1514, and [m.sub.0] [approximately equal to] 0.4287. Let us name the solid of this intersection [S.sub.2]. The center of the cross section of the solid, an ellipse, is at the point [[??].sub.0] =(5/43, - 9/49, - 168/43) * The calculations above are obtained from Maple (see ,).

Let us consider the following solid formed by a hyperboloid of two sheets and a plane as follows:

Example 3 We consider the hyperboloid of two sheets of the form [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], where

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] or we may write the surface as f (x,y,z)=5[x.sup.2]-2xy-3[y.sup.2]-8xz-6[z.sup.2]+4yz-4x-2y-14z-1=0 and consider the plane of g(x,y,z)=9x-7y-x+19=0. The eigenvectors are [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] respectively. The center of the quadric is at the point [??]0 =(-3/5, -29/35, -73/70). If the variable X' is used, the equations for the hyperboloid of two sheets becomes [(x'/a).sup.2] - [(y'/b).sup.2] - [(z'/c).sup.2] -1 = 0, where the semi axes for the hyperboloid are a [approximately equal to] 0.84365, b [approximately equal to] 1.31066, and c [approximately equal to] 0.77172, respectively. The plane can be written as [m.sub.1] x' + [m.sub.2]y' + [m.sub.3]z' + [m.sub.0] = 0 where [m.sub.1] [approximately equal to] -0.8619, [m.sub.2] [approximately equal to] 0.3925, [m.sub.3] [approximately equal to] 0.3210, and [m.sub.0] [approximately equal to] 1.7861. Let us name the solid of such intersection [S.sub.3]. The center of the elliptical cross section is at the point [[??].sub.0] = (-5343/796, -4789/796, 144/199). The calculations above are obtained from Maple (see ,).

4 Calculation methods

4.1 Volume and Natural Coordinates

Let the solid be bounded by the given quadric and the plane of [m.sub.1]x' + [m.sub.2]y' + [m.sub.3]z' + [m.sub.0] = 0. We consider the set of cross sections of the solid by parallel planes of m [member of] ([m.sub.T],[m.sub.0]], where [m.sub.T] = [-m.sub.1][x'.sub.T]-[m.sub.2][y'.sub.T]-[m.sub.3][z'.sub.T]. The set contains similar ellipses of the area [S.sub.e](m), whose centers [??] are located on the segment [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. The point T is excluded. In this point the ellipse degenerates into a point. The volume of the investigated solid is easy to find by integrating the area with respect to the coordinate z'. It is noted that for a paraboloid we see [m.sub.3] [not equal to 0], and yet [m.sub.3] = 0 for other two solids. We perform integration using different coordinate when [m.sub.3] [not equal to 0], which does not affect the result. The volume of the solid bounded by the quadric and the plane is:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (36)

4.1.1 Volume of the solid bounded by the paraboloid and the plane. We consider a paraboloid of the form [(x').sup.2]/[a.sup.2] + [(y').sup.2]/[b.sup.2] = z' and a plane of the form [m.sub.1] x' + [m.sub.2]y' + [m.sub.3]z' + [m.sub.0] [=.sub.0]. It follows from equation (17) that the area of the cross section is [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. Furthermore, we see from

equation (36) that:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (37)

An Exercise. Find the volume of [S.sub.1] by using the method mentioned in this section.

Answer: The volume of the solid [S.sub.1], calculated by using formula (37), is approximately 20.80907893464 (see Maple worksheet in ).

4.1.2. Let the solid be bounded by an ellipsoid of [(x').sup.2]/[a.sup.2] + [(y').sup.2]/[b.sup.2] + (z')2/[c.sup.2] = 1 and a plane of an equation

[m.sub.1]x' + [m.sub.2] y' + [m.sub.3]z' + [m.sub.0] = 0. Then [S.sub.e] (m) = [pi] abc [square root of [m.sup.2.sub.1] + [m.sup.2.sub.2] + [m.sup.2.sub.3] / w (1- [m.sup.2]/w). The volume of the solid bounded by the ellipsoid and given plane is:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (38)

The volume of the solid [S.sub.2], calculated by formula (38), is approximately 50.496501830226. We note the values of V([??])=0 and V([??]')=3/4 [pi]abc. [??]([??]') are poles in which polar planes are parallel to given plane ([1, 3,5-8,c]) (see Maple worksheet in ).

4.1.3. Let the solid be bounded by a hyperboloid of [(x').sup.2]/[a.sup.2] + [(y').sup.2]/[b.sup.2] + (z')2/[c.sup.2] = 1 and a plane [m.sub.1]x' + [m.sub.2] y' + [m.sub.3]z' + [m.sub.0] = 0. Then

V = [pi] abc/2 (2+3[mu] - [[mu].sup.3]), [mu] = - [parallel][m.sub.0][parallel]/[square root of w] <0. (39) 3 V w

The volume of the solid [S.sub.3] calculated using formula (39) is approximately 46.9249144931, [mu] [approximately equal to] -3.96665310075 (see Maple worksheet in ).

4.1.4 Mapping Ellipsoid to Sphere

To calculate the volume of the two pieces of an ellipsoid intersected by a plane, we may transform an ellipsoid to a sphere by stretching. Specifically, we turn the ellipsoid into the sphere by stretching along perpendicular axes, which coincide with the principal axes of the ellipsoid. The secant plane by stretching along perpendicular axes goes into the plane. Let us consider the distance from the secant plane to the center of the sphere l, and the ratio between l and the radius of the sphere [rho], say

k = 1/[rho]. Then the ratio of the volume for the segments of the sphere is equal to

K = 2-3k + [k.sup.3] / 2+ 3k-[k.sub.3]. (40)

It is clear that the ratio K under compression or stretching will stay the same. If each of the principal axes of the ellipsoid is stretched with the coefficient [square root of [lambda].sub.1], or x'=x'[square root of [lambda].sub.1], y'=y'[square root of [lambda].sub.2], and z'=z' [square root of [lambda].sub.1], then the ellipsoid becomes a sphere with the radius [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] and the volume [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. Subsequently, the volume of the ellipsoid is [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], where [DELTA]= [[lambda].sub.1], [[lambda].sub.2], [[lambda].sub.3] (41)

To find k, we use the equation of the plane of the form (19 [m.sub.1]x' + [m.sub.2] y' + [m.sub.3]z' + [m.sub.0] = 0, where [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. As a result of stretching we get

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (42)

Therefore, the distance from the plane to the center of the sphere; that is the distance from the plane (42) to the origin is

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (43)

An Exercise. Find the volume of the solid S2 by using the method mentioned in this section.

Answer: We calculate the volume of the ellipsoid [S.sub.2] by using formula (41), and it is approximately 81.67805252196. The volumes of two ellipsoidal pieces calculated by formulas (40) and (43) are approximately 50.496501830226 and 31.181550691739. Since the solid [S.sub.2] contains the center of the ellipsoid, the volume of [S.sub.2] should be more than half of the ellipsoid, and hence the volume of the solid [S.sub.2] is approximately 50.496501830226 (see Maple worksheet in ).

4.2 Volume and Spherical Coordinates

Our goal now is to find the volume of the solid of the intersection between the quadric surface and the plane by using spatial spherical coordinates.

Let point [[??].sub.0] be located on the given plane. We use coordinate system X" with the origin at [[??].sub.0] and axis x" perpendicular to the plane. For considered solids [S.sub.1],[S.sub.2], and [S.sub.3], it is convenient to choose a coordinate system with origin at the center of the elliptical cross section and the axes oriented along the ellipse axes. In general, the point [[??].sub.0] can be chosen arbitrary. We recall that the equation (1) of the quadric is

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], the equation of the plane (2) is x" = 0. We use the following spherical coordinate:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (44)

We note the Jacobian J = [[rho].sup.2] sin [??], and by substituting into the quadric equation, we reach a quadratic equation with respect to p([??],[phi]). In other words, we have

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (45)

where [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. We note that the paraboloid and the ellipsoid have the value of [k.sub.0] > 0 and unique positive solution of the equation (45) exists as follows:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (46)

In the case of the hyperboloid of two sheets, we observe that [k.sub.0] < 0 and we obtain two positive solutions of the quadric equation (45) which correspond to two different sheets of the hyperboloid. We need the solution that corresponds to the desired bounded volume.

If the basis vector of the axis x" is directed toward the surface, the volume is

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (47)

Otherwise, in the integral over [??], the upper limit should be changed to - [pi]/2.

The choice of point [[??].sub.0] is not unique. It is sufficient that point is located on the surface of the intersection solid. Its position is indicated by the sign of [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. We may assume that the calculation accuracy increases if point is located near the center of the intersection curve of the quadric surface and the plane.

We make calculations using Maple student[simpson] with 256 points (see Maple worksheet in ). We use 15 digits in all calculations. We obtain that the volume of the solid [S.sub.1] is approximately 20.8090792. The volumes for the solid [S.sub.2] and [S.sub.3] are approximately 50.4965018300 and 46.924952, respectively.

The method has an interesting trick in the numerical solution. If point [[??].sub.0] is on the intersecting plane, but outside the solid (outside the intersecting ellipse) and we use standard Maple solving command solve, Maple return a complex solution for [rho] and may obtain a complex number as the result for the volume calculation. This is surprising for an inexperienced user. But the correct result is equal to zero and also may be calculated using Maple. The authors apply the following trick to perform the calculations in cases where complex solutions occur: we add the inequalities that are identically satisfied (for example x+1111 > 0) into the simultaneous equations which we solve using Maple. After this Maple calculations with solve command, the answers often become correct.

Users should be noted that there is a significant difference between a theoretical solution and computational one, when Maple command solve is used.

4.3 Volume and Divergence Theorem

We first discuss a general strategy of how we may use the Divergence theorem to find the volume of the solid bounded by a quadric surface and a plane. We assume the followings:

(1) the standard form equation of the quadric f (x',y',z') = 0 uses the canonical coordinate system X' and principal axis,

(2) the equation of the plane is [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], where [??] = ([m.sub.1],[m.sub.2],[m.sub.3]) is a unit vector,

(3) the transformation of parametric representations to canonical coordinate system is either (t, [phi]-X' [1, 3.5-25] or ([??],[phi])-X' [1, 3.5-22, 3.5-24].

(4) the equation of the intersection curve uses the coordinates t ([phi]) or [theta]([phi]).

In this subsection, we discuss how we may find the volumes of the regions bounded by a quadric and a plane by using the following equation from Divergence theorem:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (48)

where [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. We define a vector field [??] such that [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] where [C.sub.0] [not equal to 0]. It follows from [nabla][??][C.sub.0] that [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], which is the volume of the bounded solid multiplied by [C.sub.0].

On the other hand, since [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], and the vector [??] (on the plane) is collinear with [??], we see that [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. Therefore the volume of the solid of intersection is the surface integral of [??] over the given quadric. In other words we have

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (49)

Next we attempt to find the simplest form for [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

Suppose we write [??] in the form of

[??] =([c.sub.11]x' + [c.sub.12]y' + [c.sub.13]z' ,[c.sub.21]x' + [c2.sub.2] y' + [c.sub.23]z' ,[c.sub.31]x'+[c.sub.32]y' + [c.sub.33]z'). (50)

We get [nabla][??]=[c.sub.11]+[c.sub.22]+[c.sub.33]=[C.sub.0], [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], where [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (51)

We have 3 equations for 9 variables and one inequality [C.sub.0] [not equal to] 0. We consider the following scenarios.

Case 1. If the surface is a paraboloid, by using the change of variables equation (28), we get [??] [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

with the components

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (52)

We proceed to find [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], where we use [[??].sub.p] for the vector field [??] on the paraboloid. It follows from the equations (51) and (52) that [[??].sub.p] x [??] is a polynomial of t. We observe that the right side of equation (49) has nine terms of [c.sub.ij]. We can find seven terms of [c.sub.ij] to be equal to zero.

Furthermore, we find [C.sub.0] = 1 and [[??].sub.p]=(-[m.sub.3]/[m.sub.1] z', 0 ,z'), and note that

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] Consequently, we get

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (53)

Where t([??]) is the solution to the equation (29) for the curve of intersection.

Case 2. In the case of an ellipsoid, we use the change of variables equation (30), to get [??] with the components

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (54)

Let [[??].sub.e] be the vector field [??] on the ellipsoid. We find [C.sub.0]=2 [m.sup.2.sub.1][a.sup.2]-[m.sup.2.sub.2] [b.sup.2] -[m.sup.2.sub.3][c.sup.2],

[T.sub.e] =(-[m.sup.2.sub.3][c.sup.2],x '-[m.sub.1][a.sup.2] ([m.sub.2]y'+[m.sub.3]z'),[m.sub.1] ([m.sub.2][b.sub.2] x' + [m.sub.1][a.sup.2]y'),[m.sub.1]([m.sup.3][c.sup.2]x' + [m.sub.1] [a.sup.2]z ')).

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. We get

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (55)

where [??] ([phi]) is the solution to the equation (32) for the curve of intersection.

Case 3. In the case of a hyperboloid of two sheets, we use the change of variables equation (33), to get [??] with the components

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (56)

Let [[??].sub.h] be the vector field [??] on the hyperboloid. We obtain [C.sub.0]=2[m.sup.2.sub.1][a.sup.2]+[m.sup.2.sub.2][b.sup.2]+ [m.sup.2.sub.3][c.sup.2],

[T.sub.e] =(-[m.sup.2.sub.3][c.sup.2],x '-[m.sub.1][a.sup.2] ([m.sub.2]y'+[m.sub.3]z'),[m.sub.1] ([-m.sub.2][b.sub.2] x' + [m.sub.1][a.sup.2]y'),[-m.sub.1]([m.sup.3][c.sup.2]x' + [-m.sub.1] [a.sup.2]z'))

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. We get

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (57)

We now apply the equations (53), (55), and (57) to find the volume bounded by the quadric and the plane as follows:

Example 4 We consider the scenario of the volume discussed in Example 1. We substitute the solution (29) for the intersection curve in the equation (53). We choose the unique positive solution, taking into account the signs of the numbers [n.sub.0] and [n.sub.0] [n.sub.3]/[??]x[[??].sub.3]. As a result of calculations we obtain the volume of the solid [S.sub.1] to be about 20.80907893464 (see Maple worksheet in ).

Example 5 We would like to find the volume bounded by the ellipsoid and the plane. We consider the scenario of the volume discussed in Example 2. We substitute the solution (32) for the intersection curve in the equation (53). As a result of the calculations, we find the volume of the solid [S.sub.2] is approximately 50.496501830226. The calculation runs faster if we use [[??].sub.e] and [??](([phi]). In this case Maple is able to perform the symbolic integration command int. We use different possible Maple int subcommands such as method = _d01akc and method = Ncrule. They lead to the same value of the volume but takes different time for calculations. (See Maple worksheet in ).

Example 6 We consider the scenario of the volume discussed in Example 3. We substitute the solution (35) for the intersection curve in the equation (57). As a result of calculations we find the volume of the solid [S.sub.3] to be about 46.9249144928 (see Maple worksheet in ).

We observe from the calculations made in this section 4.3 that the vector field [??] can be any vector field satisfying [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] In the first calculation we use the simplest form for the vector field

[??]=(0, [n.sub.3]y'-[n.sub.3]z',- [n.sub.2]y' +(l-[n.sub.3])z'). (58)

In this case, we do not need to load the Maple package, student[simpson] for calculation. We note that the computation time will be much less when using Maple command, int, after obtaining [[??].sub.p], [[??].sub.e] , and [[??].sub.h] individually.

4.4 Volume and the Stokes' Theorem-Part 1

Our focus now is to apply the Stokes' theorem to find the volume of the solid bounded by the quadric surface and the plane. Assume that we know the canonical equation of the quadric, the equation of the plane [??] [??]'+[m.sub.0]=0, where [??] = ([m.sub.1], [m.sub.2], [m.sub.3]) is a unit vector, and the equation of the curve of intersection. We consider the equation of [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], from Stokes' theorem (59)

where [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. We define a vector field [??] satisfying [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] For example, [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. Then [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] and we note that

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (60)

On the other hand, since [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] we have [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. In view of the equation (59) and plan using the y'- component of the integral, we obtain

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (61)

We now apply the equation (61) to find the volume bounded by the quadric and the plane as follows:

Example 7 We consider the scenario of the volume discussed in Example 1. As a result of the calculations we have the volume of the solid [S.sub.1] to be about 20.80907893464 in both cases.(See Maple worksheet in )

Example 8 We consider the scenario of the volume discussed in Example 2. With the change of variables using the spherical coordinates we get

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (62)

Similarly, we may use the z' -- component for computing the integral for equation (60). We substitute [??] and get

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (63)

Moreover, we have

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (64)

With the help of Maple, we perform numerical calculations when both components are used. We obtain the volume of [S.sub.2] to be about 50.496501830226 when y' -- component is used for computation and [S.sub.2] to be about 50.496501830226 when z'-- component is used for computation. (See Maple worksheet in ).

Example 9 We consider the scenario of the volume bounded by the hyperboloid and the plane discussed in Example 3. We follow the calculation technique from Example 5, and with the help of Maple, we perform numerical calculations for both components and obtain the volume of the solid S 3 to be about 46.9249144925 for the first case and 46.9249144927for the second case. (See Maple worksheet in ).

4.5. Calculation of the volume using Stokes' Theorem Part 2

The idea of using Stokes' theorem for calculating the volume of a solid is described as follows. We choose a vector field [??] satisfying [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], For example, we set

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (65)

We observe that

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII],

and for the z" component we get

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (66)

We see that the left-hand side of the equation (66) is equal to the volume of the solid piece to which x" <0. In particular, we have the following:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII],

and the right side of (66) contains the surface integral over the entire surface of the quadric surface.

We now apply the equations (65) and (66) to find the volume bounded by the quadric and the plane as follows:

Example 10 We consider the scenario of the volume discussed in Example 2. We get the volume of the solid [S.sub.2], a piece of the ellipsoid with x" <0, by using the formula of

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

As a result of calculation, we obtain the volume of the solid [S.sub.2] to be approximately 50.496500. The same method can be applied for calculating the volume of the other piece of the ellipsoid, which is about 31.181599 (See Maple worksheet in )

4.6. Calculation of the volume using the distance to the plane

Here we apply an integration method for calculating the volume of the intersecting solid by using the distance from the quadric to the intersecting plane. However, the approach described here needs to be addressed separately for each solid [S.sub.1],[S.sub.2], and [S.sub.3], which we describe as follows:

4.6.1. For the considered solids [S.sub.1], [S.sub.2], and [S.sub.3], it is convenient to choose a coordinate system X" so that the equation of the quadric looks like [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], and the equation of the plane is of the form of x = 0. The cross section of the solids is the ellipse with the semi-axes [a.sub.e] = [square root of [a.sub.0]"/ [l.sub.y"], and [b.sub.e] = [square root of [a.sub.0]"/ [l.sub.z"], respectively. We perform the integration as follows:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (67)

For the solid [S.sub.3] the whole solid is inside an elliptical cylinder, whose generator is perpendicular to the plane, see Figure 4 for demonstration. Therefore the calculated integral (67) coincides with the volume of the solid of intersection. As the result of calculations, we see the volume of the solid [S.sub.3] to be approximately 46.9249144958. (See Maple worksheet in ).

For the solids [S.sub.1] and [S.sub.2], the result of the calculations by the formula (67) is different from the required volume. The given integral allows us to find only the volume of the piece of the solid inside the elliptical cylinder, whose generator is perpendicular to the plane.

We describe an alternative way of calculating the volume for the bounded solid by using the distance from the solid to the intersecting plane as follows:

4.6.2. Let an arbitrary surface f ([??]') = 0 be cut by a plane of the form [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], where [parallel][??][parallel] = 1. We would like to find the volume of the solid bounded by the surface and the plane.

We assume the distance from point [??]' to the plane, in other words, the length of the perpendicular dropped from point [??]' onto the plane, is equal to [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. The projection of the elementary area [??] of the surface f ([??]')=0 on the plane is equal to [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. Therefore, the volume bounded by the surface and the plane can be calculated by using the following formula

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (68)

The calculations have been performed in the same way as in Section 4.4. As a result we obtained that the volume of the solid [S.sub.1] is approximately 20.809078934644. Similarly, [S.sub.2] and [S.sub.3] can be calculated and they are approximately 50.496501830224 and 46.9249144928, respectively. (See Maple worksheet in ). The resulting formula (67) is similar to a similar one obtained in  (see Theorem 12 in ).
```

5 Comparing The Results of Calculations

Paraboloid     Ellipsoid    Hyperboloid
[S.sub.1]     [S.sub.2]     [S.sub.3]

4.1    20.80907893   50.49650183   46.92491449
4.2    20.8090792    50.49650183    46.924952
4.3    20.80907893   50.49650183   46.92491449
4.4    20.80907893   50.49650183   46.92491449
20.80907893   50.49650183   46.92491449
4.5         -          50.4965          -
4.6.1        -             -        46.9249145
4.5.2   20.80907893   50.49650183   46.92491449
```

6 Conclusion

In this paper, we summarize several methods of calculating the volume of a solid bounded by a quadric surface and a plane. Through the combination use of the CAS, Maple, (see ) and the geometric software, GInMA (see ), we are able not only to perform complex computations algebraically but also visualize if our algebraic solutions coincide with our theories analytically. Through advanced and evolving technological tools, students and researchers are able to explore more challenging real-life problems.

References

 G.A.Corn, and T.M.Corn. Mathematical Handbook for Scientists and Engineers. Definitions, Theorems, and Formulas for Reference and Review. DOVER PUBLICATIONS, INC. Mineola, New York.

 W.-C.Yang, and M.-L.Lo. Finding Signed Areas and Volumes Inspired by Technology., Electronic Journal of Mathematics and Technology (eJMT), ISSN 1933-2823, Issue 2, Vol. 2, 2008.

 W.-C.Yang, and V.Shelomovskii. Mean Value Theorems in Higher Dimensions and Their Applications., Electronic Journal of Mathematics and Technology (eJMT), Issue 1, Vol.6, 2012.  X. Shengxiang, W.-C.Yang, and V.Shelomovskii. Computing Signed Areas and Volumes with Maple., Electronic Journal of Mathematics and Technology (eJMT), Issue 2, Vol.6, 2012.

Software Packages

 [GInMA] GInMA, 2012, Nosulya, S., Shelomovskii, D. and Shelomovskii, V. http://deoma-cmd.ru/en/Products/Geometry/GInMA.aspx

 [Maple] A product of Maplesoft, http://www.maplesoft.com/.

Supplemental Electronic Materials

All GInMA supplemental materials that accompany this paper can be used from figures. Install the GInMA software from the website and click on the picture.

 Shelomovskii, V., Maple worksheet for Examples 1-3. This worksheet may take up to 250 seconds to run.

 Shelomovskii, V., Maple worksheet for section 4.2.

 Shelomovskii, V., Maple worksheet for sections 4.1,4.3,4.4,4.6.2. This worksheet may take up to 250 seconds to run.

 Shelomovskii, V., Maple worksheet for section 4.5. This worksheet may take up to 520 seconds to run.

 Shelomovskii, V., Maple worksheet for section 4.6.1.

Russia

E-mail: vvsss@rambler.ru

Wei-Chi Yang

Department of Mathematica and Statistica

Author: Printer friendly Cite/link Email Feedback Shelomovskii, Vladimir V.; Yang, Wei-Chi Electronic Journal of Mathematics and Technology Report 1USA Feb 1, 2013 8465 Software tools for visualizing multivalued functions. Optimization of the courses in geometry by the usage of Dynamic Geometry Software Sam. Planes (Geometry) Surfaces (Geometry) Theorems (Mathematics) Volume (Cubic content) Volume (Cubic measurement)