Multiple Integral Article Index for
Multiple
Shopping
Integral
Website Links For
Multiple
 

Information About

Multiple Integral




The multiple integral is a kind of ''definite'' ''' Integral ''' extended to Functions of more than one real Variable (i.e. f(x,y)\,\! or f(x,y,z)\,\!).


MULTIPLE INTEGRALS ARE NOT THE SAME AS ITERATED INTEGRALS

It is easy to confuse the concepts of ''multiple integral'' and Iterated Integral , especially since the same notation is often used for either concept. The notation

:\int_0^1\int_0^1 f(x,y)\,dy\,dx

in some cases means an iterated integral rather than a true double integral. In an iterated integral, the outer integral

:\int_0^1 \cdots \, dx

is the integral with respect to ''x'' of the following function of ''x'':

:g(x)=\int_0^1 f(x,y)\,dy.

A double integral, on the other hand is defined with respect to area in the ''xy''-plane. If the double integral exists, then it is equal to either of the two iterated integrals (either "''dy dx''" or "''dx dy''"; see Fubini's Theorem ) and one often computes it by computing either of the iterated integrals. But sometimes the two iterated integrals exist when the double integral does not, and in some such cases the two iterated integrals are different numbers, i.e., one has

:\int_0^1\int_0^1 f(x,y)\,dy\,dx
eq \int_0^1\int_0^1 f(x,y)\,dx\,dy.

For an elementary example (doable by the methods of first-year Calculus ), see Examples Of Fubini's Theorem . This is an instance of rearrangement of a Conditionally Convergent integral.

The notation

:\int_{[0,1] imes[0,1]} f(x,y)\,dx\,dy

may be used if one wishes to be emphatic about intending a double integral rather than an iterated integral.


MULTIPLE INTEGRALS

If, conceptually, the definite integral for functions of one variable represent the Area of the region on the graph between the curve defined by the function and the ''x''-axis, that for functions of two variables ( Double Integral ) represent the measure of the space between the surface defined by the function and the plane which contains its Domain , so they describe not an area but a Volume of a particular solid called Cylindroid ; you obtain the same value if you consider the ''triple integrals'' (functions of three variables) calculated with the constant ''f''(''x'', ''y'', ''z'') = 1. If the number of variables is higher you will calculate " Hypervolume s" (volumes of solid of more than three dimensions) that you cannot graph.

In that example the volume of the Parallelepiped of sides 4x6x5 is obtainable by two ways:
  • by the Double Integral \iint_D 5 \ dx\, dy of the function ''f''(''x'', ''y'') = 5 calculated in the "bi-dimensional interval" ''D'' (region contained in the ''xy''-plane)

  • by the triple integral \iiint_\mathrm{parallelepiped} 1 \, dx\, dy\, dz of the constant function 1 calculated on the "three-dimensional" interval coinciding with the same parallelepiped; in that case the volume is considered as the "sum" of all the Infinitesimal elements that compose the domain.


Because is impossible to calculate the Antiderivative of a function of more than one variable, ''indefinite'' multiple integrals do not exist so they are all ''definite'' integrals.


SOME PRACTICAL APPLICATIONS


These integrals are used in a lot of Physics applications.

In Mechanics the Moment Of Inertia is calculated as a volume integral (that is a triple integral) of the Density weighed with the square of the distance from the axis:

:I_z = \int_V^. ho r^2\, dV.

In Electromagnetism the Maxwell's Equations can be written with multiple integrals to calculate the total magnetic and electric fields. In the example the Electric Field produced by a Distribution Of Charges is obtained by a ''triple integral'' of a vector function:



:and obtain the same value.


Normal domains on R3

Is immediate the extension of these formulas to the triple integrals.

''T'' is a domain perpendicular to the ''xy''-plane respect to the α (''x'',''y'',''z'') and β(''x'',''y'',''z'') functions. Then:

:\iiint_T f(x,y,z) \ dx\, dy\, dz = \iint_D dx\, dy \int_{\alpha (x,y,z)}^{\beta (x,y,z)} f(x,y,z) \, dz

(this definition is the same for the other five normality cases on R3).


Change of variables


Often, because the limits of integration are not easily interchangeable (without normality or with complex formulas to integrate) one makes a ''change of variables'' to rewrite the integral on more a "comfortable" region, that is described in simpler way through formulas. To do that, the function must be changed to the new coordinates.

:Example (1-a):
::The function is f(x, y) = (x-1)^2 +\sqrt y;
::if you adopt this substitution x' = x-1, \ y= y' \, \! therefore x = x' + 1, \ y=y' \,\!
::you obtain new function f(x,y)_2 = x' ^2 +\sqrt y.

There exist three main "kinds" of changes of variable (one in R2, two in R3), however is possible to hand with this method so as to to operate the substitution that more thinks good.


Polar coordinates

See Also: coordinates (mathematics)



In R2 if the domain as a circular "symmetry" (if it describes a Circular Crown ) and the function has some "particular" characteristics you can apply the ''passate to polar coordinates'' (see the example in the picture) which means that the generic points ''P(x,y)'' in cartesian coordinates switch to their respective points in polar coordinates. That allows to change the "shape" of the domain and simplify the operations.

The fundamental relation to make the transformation is the follow:

:f(x,y) ightarrow f( ho \ \cos \phi, ho \ \sin \phi ).

Example (2-a):

:The function is f(x,y) = x + y\,\!

:and applying the transformation you get

::f( ho, \phi) = ho \cos \phi + ho \sin \phi = ho \ (\cos \phi + \sin \phi ).

Example (2-b):
:The function is f(x,y) = x^2 + y^2\,\!
:In this case you get

::f( ho, \phi) = ho^2 (\cos \phi^2 + \sin \phi^2) = ho^2\,\!

:using the Pythagorean Trigonometric Identity (very useful to simplify this operations).

The transformation of the domain is made by defining the radius' crown's length and the amplitude of the described angle to define the ρ, φ intervals starting from ''x'', ''y''.

Example (2-c):
:The domain is D = x^2 + y^2 \le 4\,\!, that is a circumference of radius 2; it's evident that the described angle is the Circle Angle , so φ varies from 0 to 2π, while the crown's radius varies from 0 to 2 (the crown with the inside radius null is just a circle).

Example (2-d):
:The domain is D = \{ x^2 + y^2 \le 9, \ x^2 + y^2 \ge 4, \ y \ge 0 \}, that is the circular crown in the semiplane of positive ''y'' (please see the picture in the example); you note that φ describes a Plane Angle while ρ varies from 2 to 3. Therefore the transformed domain will be the following Rectangle :

::T = \{ 2 \le ho \le 3, \ 0 \le \phi \le \pi \}.

The Jacobian Determinant of that transformation is the following:

: rac{\partial (x,y)}{\partial ( ho, \phi)} =
\begin{vmatrix}
\cos \phi & - ho \sin \phi \
\sin \phi & ho \cos \phi
\end{vmatrix} = ho


which has been got by inserting the partial derivatives of ''x'' = ρ cos(φ), ''y'' = ρ sin(φ) in the first column respect to ρ and in the second respect to φ, so the ''dx dy'' differentials in this transformation becomes ρ ''d''ρ ''d''φ.

Once transformed the function and evaluated the domain, it's possible to define the formula for the change of variables in polar coordinates:

:\iint_D f(x,y) \ dx\, dy = \iint_T f( ho \cos \phi, ho \sin \phi) ho \, d ho\, d \phi.

Please note that φ is valid in the interval while ρ, because it is a measure of a length, can only have positive values.

Example (2-e):
:The function is f(x,y) = x\,\! and as the domain the same in 2-d example.
:From the previous analysis of ''D'' we know the intervals of ρ (from 2 to 3) and of φ (from 0 to π). Now let's change the function:

::f(x,y) = x \longrightarrow f( ho,\phi) = ho \ \cos \phi.

:finally let's apply the integration formula:

::\iint_D x \, dx\, dy = \iint_T ho \cos \phi \ ho \, d ho\, d\phi.

:Once known the intervals you have

::\int_0^{\pi} \int_2^3 ho^2 \cos \phi \ d \phi \ d ho = \int_0^{\pi} \cos \phi \ d \phi \left[ rac{ ho^3}{3} ight]_2^3 = \left \sin \phi ight _0^{\pi} \ \left(9 - rac{8}{3} ight) = 0.


Cylindrical coordinates


In R3 the integration on domains with a circular base can be made by the ''passage in cylindrical coordinates''; the transformation of the function is made by the following relation:
f(x,y,z) ightarrow f( ho \ \cos \phi, ho \ \sin \phi, z)


The domain's transformation is easy because graphically only the shape of the base varies while the three-dimensional development follows that of the starting region.

Example (3-a):

:The region is D = \{ x^2 + y^2 \le 9, \ x^2 + y^2 \ge 4, \ 0 \le z \le 5 \} (that is the "tube" whose base is the circular crown of the 2-d example and whose height is 5); if you apply the transformation you'll get this region: T = \{ 2 \le ho \le 3, \ 0 \le \phi \le \pi, \ 0 \le z \le 5 \} (that is the parallelepiped whose base is the rectangle in 2-d example and whose height is 5).

Because the ''z'' component is unvaried during the transformation, the ''dx dy dz'' differentials vary as in the passage in polar coordinates, therefore they become ''ρ dρ dφ dz''.

Finally, it is possible to apply the final formula for the passage in cylindrical coordinates:

:\iiint_D f(x,y,z) \, dx\, dy\, dz = \iiint_T f( ho \cos \phi, ho \sin \phi, z) ho \, d ho\, d\phi\, dz.

This method is convenient in case of cylindrical or conical domains or in regions where is easy to individuate the ''z'' interval and even transform the circular base and the function.

Example (3-b):
:The function is f(x,y,z) = x^2 + y^2 + z\,\! and as integration domain this Cylinder : D = \{ x^2 + y^2 \le 9, \ -5 \le z \le 5 \}.
:The transformation of ''D'' in cylindrical coordinates is the following:

::T = \{ 0 \le ho \le 3, \ 0 \le \phi \le 2 \pi, \ -5 \le z \le 5 \}.

:while the function becomes

::f( ho \ \cos \phi, ho \ \sin \phi, z) = ho^2 + z\,\!


:Finally you can apply the integration's formula:

::\iiint_D (x^2 + y^2 +z) \, dx\, dy\, dz = \iiint_T ( ho^2 + z) ho \, d ho\, d\phi\, dz;

:developing the formula you have

::\int_{-5}^5 dz \int_0^{2 \pi} d\phi \int_0^3 ( ho^3 + ho z )\, d ho = 2 \pi \int_{-5}^5 \left[ rac{ ho^4}{4} + rac{ ho^2 z}{2} ight]_0^3 \, dz

::= 2 \pi \int_{-5}^5 \left( rac{81}{4} + rac{9}{2} z ight)\, dz = \cdots = 2 \pi \left( rac{405}{2} + 225 ight).


Spherical coordinates


In R3 some domains have a spherical symmetry, so it's possible to determinate the coordinates of every point of the integration's region by two angles and one distance. It's possible to use therefore the ''passage in spherical coordinates''; the function is transformed by this relation:
f(x,y,z) \longrightarrow f( ho \sin heta \cos \phi, ho \sin heta \sin \phi, ho \cos heta)\,\!


Please note that points on ''z'' axis don't have a precise characterization in spherical coordinates, so ''θ'' can varies from 0 to π .

The better integration's domain for this passage is obviously the sphere.

Example (4-a):

:The domain is D = x^2 + y^2 + z^2 \le 16 (sphere with radius 4 and center in the origin); applying the transformation you get this region: T = \{ 0 \le ho \le 4, \ 0 \le \phi \le 2 \pi, \ 0 \le heta \le \pi \}.

:The Jacobian determinant of this transformation is the following:

:: rac{\partial (x,y,z)}{\partial ( ho, heta, \phi)} =
\begin{vmatrix}
\sin heta \cos \phi & ho \cos heta \cos \phi & - ho \sin heta \sin \phi \
\sin heta \sin \phi & ho \cos heta \sin \phi & ho \sin heta \cos \phi \
\cos heta & - ho \sin heta & 0
\end{vmatrix} = ho^2 \sin heta


:The ''dx dy dz'' differentials therefore are transformed to ρ2 sin(θ) ''d''ρ ''d''θ ''d''φ.

:Finally you obtain the final integration formula:

::\iiint_D f(x,y,z) \, dx\, dy\, dz = \iiint_T f( ho \sin heta \cos \phi, ho \sin heta \sin \phi, ho \cos heta) ho^2 \sin heta \, d ho\, d heta\, d\phi.

:It's better to use this method in case of spherical domains and in case of functions that you can easily simplify by the first fundamental relation of trigonometry extended in '''R'''3 (please see the 4-b example); in the other cases it's better to use the passage in cylindrical coordinates (please see the 4-c example).

Example (4-b):

D


:Its transformation is very easy:

::f( ho \sin heta \cos \phi, ho \sin heta \sin \phi, ho \cos heta) = ho^2,\,

:while we know the intervals of the transformed region ''T'' from ''D'':


::(0 \le ho \le 4, \ 0 \le \phi \le 2 \pi, \ 0 \le heta \le \pi).\,

:Let's apply therefore the integration's formula:

::\iiint_D (x^2 + y^2 +z^2) \, dx\, dy\, dz = \iiint_T ho^2 \ ho^2 \sin heta \, d ho\, d heta\, d\phi,

:and developing we get

::\iiint_T ho^4 \sin heta \, d ho\, d heta\, d\phi = \int_0^{\pi} \sin heta \,d heta \int_0^4 ho^4 d ho \int_0^{2 \pi} d\phi = 2 \pi \int_0^{\pi} \sin heta \left[ rac{ ho^5}{5} ight]_0^4 \, d heta

::= 2 \pi \left[ rac{ ho^5}{5} ight]_0^4 \left \cos heta ight _0^{\pi} = 4 \pi \cdot rac{1024}{5} = rac{4096 \pi}{5}.

Example (4-c):
:The domain ''D'' is the ball with center in the origin and radius ''3a'' (D = x^2 + y^2 + z^2 \le 9a^2 \,\!) and f(x,y,z) = x^2 + y^2\,\! is the function to integrate.

:Looking at the domain seems convenient adopting the passage in spherical coordinates, in fact are immediate the intervals of the variables that delimit the new ''T'' region:

::0 \le ho \le 3a, \ 0 \le \phi \le 2 \pi, \ 0 \le heta \le \pi.\,

:However applying the transformation we get

::f(x,y,z) = x^2 + y^2 \longrightarrow ho^2 \sin^2 heta \cos^2 \phi + ho^2 \sin^2 heta \sin^2 \phi = ho^2 \sin^2 heta.

:Applying the formula for integration it would be obtained:

::\iiint_T ho^2 \sin^2 heta ho^2 \sin heta \, d ho\, d heta\, d\phi = \iiint_T ho^4 \sin^3 heta \, d ho\, d heta\, d\phi

:very hard to solve. This problem will be solved by using the passage in cylindrical coordinates. The new ''T'' intervals are

::0 \le ho \le 3a, \ 0 \le \phi \le 2 \pi, \ - \sqrt{9a^2 - ho^2} \le z \le \sqrt{9a^2 - ho^2};

:the ''z'' interval has been obtained by dividing the ball in two semispheres simply by solving the Inequality from the formula of ''D'' (and then directly transforming ''x2 + y2'' in ''ρ2''). The new function is simply ''ρ2''. Therefore let's apply the integration formula

::\iiint_T ho^2 ho \ d ho d \phi dz.

:Then we get

::\int_0^{2 \pi} d\phi \int_0^{3a} ho^3 d ho \int_{- \sqrt{9a^2 - ho^2} }^{\sqrt{9 a^2 - ho^2} }\, dz = 2 \pi \int_0^{3a} 2 ho^3 \sqrt{9 a^2 - ho^2} \, d ho.

:Now let's apply the transformation

::9 a^2 - ho^2 = t\,\! \longrightarrow dt = -2 ho\, d ho \longrightarrow d ho = rac{d t}{- 2 ho}\,\!

:(the new intervals become 0, 3a \longrightarrow 9 a^2, 0). We get

::- 2 \pi \int_{9 a^2}^{0} ho^2 \sqrt{t}\, dt

:because ho^2 = 9 a^2 - t\,\!, we get

::-2 \pi \int_{9 a^2}^0 (9 a^2 - t) \sqrt{t}\, dt,

:after inverting the integration's bounds and multiplying the terms between parenthesis is possible to decompose the integral in two parts that can be directly solved:

::2 \pi \left \int_0^{9 a^2} 9 a^2 \sqrt{t} \, dt - \int_0^{9 a^2} t \sqrt{t} \, dt ight = 2 \pi \left a^2 rac{2}{3} t^{ rac{3}{2} } - rac{2}{5} t^{ rac{5}{2}} ight _0^{9 a^2}

::= 2 \cdot 27 \pi a^5 ( 6 - rac{2}{5} ) = 54 \pi rac{28}{5} a^5 = rac{1512 \pi}{5} a^5.

:Thanks to the passage in cylindrical coordinates it was possible to reduce the triple integral to an easier one-variable integral.

See also the differential volume entry in Nabla In Cylindrical And Spherical Coordinates .


EXAMPLE OF MATHEMATICAL APPLICATIONS - CALCULATIONS OF VOLUME


Thanks to the methods previously described are possible to demonstrate the value of the volume of some solid ones.

::\mathrm{Volume} = \int_0^{2 \pi } d \phi \int_0^r h ho \ d ho = h 2 \pi \left[ rac{ ho^2}{2 } ight]_0^R = \pi R^2 h


::\mathrm{Volume} = \int_0^{2 \pi }\, d \phi \int_0^{ \pi } \sin heta\, d heta \int_0^R ho^2\, d ho = 2 \pi \int_0^{ \pi } \sin heta rac{R^3}{3 }\, d heta = rac{2}{3 } \pi R^3 \cos heta _0^{ \pi } = rac{4}{3 } \pi R^3.

:: \mathrm{Volume} = \int_0^\ell dx \int_0^{\ell-x }\, dy \int_0^{\ell-x-y }\, dz = \int_0^\ell dx \int_0^{\ell-x } (\ell - x - y)\, dy

:: = \int_0^\ell (\ell^2 - 2\ell x + x^2 - rac{ (\ell-x)^2 }{2 })\, dx = \ell^3 - \ell \ell^2 + rac{\ell^3}{3 } - \left[ rac{\ell^2}{2 } - \ell x + rac{x^2}{2 } ight]_0^\ell =

:: = rac{\ell^3}{3 } - rac{\ell^3}{6 } = rac{\ell^3}{6}

: ''Verification'': Volume = base area × height/3 = rac{\ell^2}{2 } \cdot \ell/3 = rac{\ell^3}{6}.


MULTIPLE IMPROPER INTEGRAL

In case of unbounded domains or unbounded integrands near the frontier of the dominion you have the double Improper Integral or the triple improper integral.


SEE ALSO



BIBLIOGRAPHY

''Sources verified on: Robert A. Adams - Calcolo differenziale 2, Funzioni di più variabili ISBN 8840810242''