computation of surface area of portion of paraboloid


1 Statement of Problem

Using the result of the previous example we may, for instance, compute the area of the portion of paraboloid which lies below the plane z=5:

A=d2A=x2+3y2<51+4x2+36y2𝑑x𝑑y.

2 Integration over x

Despite its innocuous appearance, this integralDlmfPlanetmath is actually rather difficult to compute. As a first step, we shall rewrite it as a double integral and perform the integration over x:

A = -5/3+5/3-5-3y2+5-3y21+4x2+36y2𝑑x𝑑y
= -5/3+5/3[x21+4x2+36y2+9y2log(2x+1+4x2+36y2)]-5-3y2+5-3y2𝑑y
= -5/3+5/3{(5-3y2)(21+24y2)+9y2log(21+24y2+25-3y221+24y2-25-3y2)}𝑑y

3 Notes about the multi-valued integrand and branches

Before proceeding further, it might be worthwhile to dwell on the subtleties involved in the last step of the calculation. Both the square root and the logarithm function are multiply-valued so we need to make sure we are using the right branch of these functionsMathworldPlanetmath in order to get a correct answer. In particular, when substituting the limits of integration into a multiply-valued expression, we need to make sure whether or not we passed from one branch to another along the path of integration. In this case, it is easy to see that it did not change. For the branch of a function to change, the argumentMathworldPlanetmath of the function must pass through a branch point (more generally, if we allow complex arguments, it may also encircle a branch point) In the case of a square root, the branch point occurs at zero. Now, the argument of the square root in our expression is 1+4x2+36y2. It is clear that, as long as y is not equal to zero, this expression cannot equal zero for any real value of x. Hence, it follows that the branch of the square root to be taken at the upper limitMathworldPlanetmath must be the same branch as was taken at the lower limit. (In our case, this is the positive branch.) This is why the first term works out as it does.

The second term needs a little more discussion because of the logarithmMathworldPlanetmath. Recall that the branch point of the logarithm is also located at zero. Thus we must ask whether the argument can equal zero along the path of integration. Setting it equal to zero, we find

2x+1+4x2+36y2=0, or
4x2=1+4x2+36y2.

(It might be worth noting that the choice of sign of square root turns out to be irrelevant here because it is lost in the process of squaring.)

0=1+36y2

Hence we can only encounter a branch point when y=±6i. For any other value, we will use the same branch of the logarithm at the upper limit as at the lower limit. Hence, we arrived at our result for the second term.

4 A simple change of variables

We may simplify the appearance of our expression by making the change of variables y=t5/3.

A=5/3-1+1{105(1-t2)(1+4021y2)+25t2log(21+40t2+25(1-t2)21+40t2-25(1-t2))}𝑑t

Furthermore, since the integrand is an even functionMathworldPlanetmath of t, we may fold the range of integration on half:

A=25/301{105(1-t2)(1+4021y2)+25t2log(21+40t2+25(1-t2)21+40t2-25(1-t2))}𝑑t

For convenience of reference, let us write this result as A=107I1+505/3I2, where

I1=01(1-t2)(1+4021y2)𝑑t

and

I2=01t2log(21+40t2+25(1-t2)21+40t2-25(1-t2))𝑑t.

5 The integral I1

Since the argument of the square root is a quartic polynomial, the integral I1 is an elliptic integralMathworldPlanetmath. Since the quartic has distinct roots, it is not possible to evaluate this integral in tems of elementary functionsMathworldPlanetmath. The best we can do is to express it in terms of the standard elliptic integrals E and K. To begin, we shall write the square root at the quotient of the square by the root

I1=011+1921t2-4021t4(1-t2)(1+4021t2)𝑑t.

To deal with the t4 in the numerator, we shall use a sneaky trick. First, we integrate by parts:

I1 = 01(1-t2)(1+4021t2)𝑑t
= [t(1-t2)(1+4021t2)]01-01t3821t-16021t32(1-t2)(1+4021t2)𝑑t
= -192101t2dt(1-t2)(1+4021t2)+802101t4dt(1-t2)(1+4021t2).

Next, we form the linear combination of our two expressions for I1 which is designed to make the term with the t4 in the numerator cancel.

I1=23I1+13I1=012+1921t23(1-t2)(1+4021t2)𝑑t.

It is easy to re-express this in terms of complete elliptic integrals

I1 = 016140+1940(1+4021t2)3(1-t2)(1+4021t2)𝑑t
= 6112001dt(1-t2)(1+4021t2)+1940011+4021t21-t2𝑑t
= 13K(-4021)+23E(-4021)

6 The integral I2

We now turn our attention to the integral I2. We will start by removing the radicals from the denominator of the argument of the logarithm.

21+40t2+25(1-t2)21+40t2-25(1-t2)
= 21+40t2+25(1-t2)21+40t2+25(1-t2)×21+40t2+25(1-t2)21+40t2-25(1-t2)
= 41+20t2+8(21+40t2)(1-t2)1+60t2

Substituting this simplified expression back into the integral and using the fact that the logarithm of a quotient equals the differencePlanetmathPlanetmath of logarithms, we obtain the following:

I2 = 01y2log(21+40t2+26(1-t2)21+40t2-25(1-t2))𝑑t
= 01t2log(41+20t2+8(21+40t2)(1-t2))𝑑t-01t2log(1+60t2)𝑑t

For convenience, let us define the notation I2=I3-I4, where

I3=01t2log(41+20t2+8(21+40t2)(1-t2))𝑑t

and

I4=01t2log(1+60t2)𝑑t.

7 The integral I4

Since it can be evaluated in terms of elementary functions rather easily, let us now dispose of the integral I4. Begin by rescaling the variable:

I4=01t2log(1+60t2)𝑑t=60-3/2060t2log(t2+1)𝑑t

We may turn this into the integration of a rational function by integrating by parts. Then we may make a change of variable s=t2.

060t2log(t2+1)𝑑t = [13t3log(t2+1)]060-13060t3t2+1𝑑t
= [13t3log(t2+1)]060-16060s2s+1𝑑s
= [13t3log(t2+1)]060-16[s22-s+log(s+1)]060
= 13603/2log61-300+10-16log61

Hence, we have

I4=13log61-1510800log61-290

8 Using integration by parts to eliminate a logarithm

The integral I3 may be simplified by using integration by parts to get rid of the logarithm in the integrand:

I3 = 01t2log(41+20t2+8(21+40t2)(1-t2))𝑑t
= 1301log(41+20t2+8(21+40t2)(1-t2))d(t3)
= 13[t3log(41+20t2+8(21+40t2)(1-t2))]01-
1301t340t+4(38t-160t3)/(21+40t2)(1-t2)41+20t2+8(21+40t2)(1-t2)𝑑t
= 13log61-1301t440+(152-640t2)/(21+40t2)(1-t2)41+20t2+8(21+40t2)(1-t2)𝑑t

9 Simplifying a fraction

As we did once earlier, we may simplify the fraction by eliminating the radicals from the denominator.

40+(152-640t2)/(21+40t2)(1-t2)41+20t2+8(21+40t2)(1-t2)
= 1(21+40t2)(1-t2)×152-640t2+40(21+40t2)(1-t2)41+20t2+8(21+40t2)(1-t2)
= 1(21+40t2)(1-t2)×152-640t2+40(21+40t2)(1-t2)41+20t2+8(21+40t2)(1-t2)×
  41+20t2-8(21+40t2)(1-t2)41+20t2-8(21+40t2)(1-t2)
= -488-29280t2-1280t4+(-1052+5920t2)(21+40t2)(1-t2)(337-1052t2-2160t4)(21+40t2)(1-t2)

This may be split into three terms like so:

-1052+5920t2337-1052t2-2160t4 + 1627(21+40t2)(1-t2)
- 18568-773728t227(337-1052t2-2160t4)(21+40t2)(1-t2)

Substituting this back into the integral, we obtain

I3 = 13log3-1301-t2+t2(1+2t2)(1-t2)(1+2t2)(1-t2)𝑑t
= 13log3-1301t2𝑑t+1301t2dt(1+2t2)(1-t2)

As for the first integral, it is elementary. As for the second integral, we shall express it in terms of standard elliptic integrals just as we did with I1 earlier.

I3 = 13log3-19+1301-1/2+(1+2t2)/2(1+2t2)(1-t2)𝑑t
= 13log3-19-1601dt(1+2t2)(1-t2)+16011+2t21-t2𝑑t
= 13log3-19-16K(-2)+16E(-2)

10 Endgame

Having done all the work, all that remains is to collect the pieces and arrive at an expression for the area of our paraboloid. Recall that

A=205/3I1+305/3I2

where

  • I1=13K(-2)+23E(-2),

  • I2=I3-I4,

  • I4=-29, and, most recently,

  • I3=13log3-19-16K(-2)+16E(-2).

Combining these equations, we find that

I2=13log3+19-16K(-2)+16E(-2),

hence

A=5915K(-2)+25915E(-2)+10315log3-10915.

11 Closing Remarks

In many ways, the complexities of the present calculation are typical of what happens when one attempts to compute even the simplest of integrals with respect to area explicitly. Because of the square root which appears in the fomulae for reducing integrals with respect to surface area, one is faced with the evaluation of a double integral involving a square root. This can easily lead to elliptic integrals and even more complicated entities such as hyperelliptic integrals.

http://planetmath.org/node/6660main entry http://planetmath.org/node/6672previous example

Title computation of surface area of portion of paraboloid
Canonical name ComputationOfSurfaceAreaOfPortionOfParaboloid
Date of creation 2013-03-22 14:58:22
Last modified on 2013-03-22 14:58:22
Owner rspuzio (6075)
Last modified by rspuzio (6075)
Numerical id 61
Author rspuzio (6075)
Entry type Example
Classification msc 28A75