Кратный интеграл — Википедия

Кратный интеграл — определённый интеграл, взятый от ${\displaystyle n>1}$ переменных; например:

${\displaystyle \underbrace {\int \cdots \int } _{n}f(x_{1},\ldots ,x_{n})dx_{1}\cdots dx_{n}}$.

В классическом анализе чаще всего используются двойной интеграл➤ и тройной интеграл➤ — интегралы от двух и трёх переменных соответственно.

Этот раздел должен быть полностью переписан. На странице обсуждения могут быть пояснения.

В общем случае кратный интеграл определяется для функции ${\displaystyle f\colon B\to \mathbb {R} }$ на жорданово измеримом множестве ${\displaystyle B\subset \mathbb {R} ^{n}}$ с использованием понятия разбиения ${\displaystyle \sigma }$ — набора попарно непересекающихся подмножеств ${\displaystyle \sigma =\left\{{{U}_{i}}\subset B\right\}}$, которые в объединение дают всё ${\displaystyle B}$. Мелкостью ${\displaystyle \left|\sigma \right|}$ измерения называется наибольший диаметр множеств ${\displaystyle U_{i}\in \sigma }$:

${\displaystyle \left|\sigma \right|=\max \left\{\operatorname {diam} \left({{U}_{i}}\right)\right\}}$.

Разбиение называется конечным, если является конечным множеством, и измеримым, если все его элементы — измеримые (в данном случае — по Жордану) множества.

Кратным (${\displaystyle n}$-кратным) интегралом функции ${\displaystyle f}$ на ${\displaystyle B}$ называется число ${\displaystyle I}$ (если оно существует), такое, что при любой ${\displaystyle \varepsilon }$-окрестности числа ${\displaystyle I}$ всегда найдётся такое разбиение множества ${\displaystyle B}$ и набор промежуточных точек, что сумма произведений значения функции в промежуточной точке разбиения на меру разбиения будет попадать в эту окрестность:

${\displaystyle \forall \varepsilon >0\;\;\exists \delta >0}$ : ${\displaystyle \sigma =\{U_{i}\}_{i=1}^{m}}$ : ${\displaystyle \left|\sigma \right|<\delta \;\forall \xi _{i}\in U_{i}\;\;\left|{\sum _{i=1}^{m}f(\xi _{i})\mu (U_{i})-I}\right|<\varepsilon }$.

(здесь ${\displaystyle \mu (U_{i})}$ — мера множества ${\displaystyle U_{i}}$.

Это определение можно сформулировать в другой форме с использованием интегральных сумм: для данного разбиения ${\displaystyle \sigma =\{U_{i}\}_{i=1}^{m}}$ и множества точек ${\displaystyle \xi =\{\xi _{i}\in U_{i}\}}$ рассматривается интегральная сумма:

${\displaystyle \zeta (f,\sigma ,\xi )=\sum _{i=1}^{m}f(\xi _{i})\mu (U_{i})}$,

тогда кратным интегралом функции ${\displaystyle f\colon B\to \mathbb {R} }$ называют предел:

${\displaystyle I=\lim _{\left|\sigma \right|\to 0}\zeta (f,\sigma ,\xi )}$,

если он существует. Предел берётся по множеству всех последовательностей разбиений, с мелкостью стремящейся к 0. Разумеется, это определение отличается от предыдущего, по сути, лишь используемым языком.

В случае ${\displaystyle n=1}$ кратный интеграл совпадает с интегралом Римана.

В современном анализе чаще всего используется векторная форма записи^[1]:

${\displaystyle \int \limits _{G}{f\left(X\right)dX}=I}$.

Иногда ставят знак интеграла ${\displaystyle n}$ раз, записывают функцию и ${\displaystyle n}$ дифференциалов:

${\displaystyle \underbrace {\int {\cdots }\int } _{n}f({x_{1}},\ldots ,{x_{n}})d{x_{1}}\cdots d{x_{n}}=I}$.

Для двойного и тройного интегралов используются также обозначения ${\displaystyle \iint }$ и ${\displaystyle \iiint }$ соответственно.

Многократное использование знака характерно для учебной литературы, в современных научных работах используется редко.

Если функция непрерывна на измеримом по Жордану компакте, то она интегрируема на нем. Как следствие, неограниченная функция на множестве может быть не интегрируемой, даже если она непрерывна. Например, функция ${\displaystyle y=1/x}$ не интегрируема на интервале ${\displaystyle \left(0;1\right)}$.

Если функция определена на измеримом по Жордану множестве, у которого существуют сколь угодно мелкие разбиения, для которых данная функция неограничена на объединении всех их элементов положительной меры, то эта функция неинтегрируема на этом множестве.

Критерий Дарбу: если существуют верхний ${\displaystyle I^{*}}$ и нижний ${\displaystyle I_{*}}$ интегралы Дарбу функции на ${\displaystyle G}$, тогда, если верхний и нижний интегралы Дарбу равны, то данная функция интегрируема на ${\displaystyle G}$, причём:

${\displaystyle I^{*}=I_{*}=\int \limits _{G}{f\left(X\right)dX}}$.

Критерий Лебега: если ${\displaystyle G}$ — измеримое по Жордану множество, то ${\displaystyle f}$ интегрируема на ${\displaystyle G}$ если она ограничена на ${\displaystyle G}$ и непрерывна на ${\displaystyle G\setminus E}$, где множество ${\displaystyle E}$ имеет меру Лебега нуль.

Кратные интегралы линейны по функции: если ${\displaystyle G}$ измеримо, функции ${\displaystyle f}$ и ${\displaystyle g}$ интегрируемы на ${\displaystyle G}$, то:

${\displaystyle \forall \lambda ,\mu \in \mathbb {R} :~\int \limits _{G}{\left(\lambda f+\mu g\right)dX}=\lambda \int \limits _{G}{fdX}+\mu \int \limits _{G}{gdX}}$.

Аддитивность по множеству интегрирования позволяет раздельно проинтегрировать по частям областей: если ${\displaystyle {G}_{1}}$ и ${\displaystyle G_{2}}$ измеримы, ${\displaystyle G_{1}\cap G_{2}=\varnothing }$ и ${\displaystyle G_{1}\cup G_{2}=G}$, а функция ${\displaystyle f(X)}$ определена и интегрируема на каждом из множеств ${\displaystyle G_{1}}$ и ${\displaystyle G_{2}}$, то интеграл по ${\displaystyle G}$ существует и равен сумме интегралов по ${\displaystyle G_{1}}$ и ${\displaystyle G_{2}}$:

${\displaystyle \int _{G}f(X)dX=\int _{G_{1}}f(X)dX+\int _{G_{2}}f(X)dX}$.

Кратные интегралы монотонны по функции: если ${\displaystyle G}$ измеримо, функции ${\displaystyle f}$ и ${\displaystyle g}$ интегрируемы на ${\displaystyle G}$, причём ${\displaystyle \forall X\in G\colon f\left(X\right)\leqslant g\left(X\right)}$, то:

${\displaystyle \int _{G}f(X)dX\leqslant \int _{G}g(X)dX}$.

Интегральное неравенство треугольника (следствие монотонности):

${\displaystyle \left|\int _{G}f(X)dX\right|\leqslant \int _{G}\left|f(X)\right|dX}$

Интегральная теорема о среднем: если ${\displaystyle G}$ — компакт, функция ${\displaystyle f(X)}$ непрерывна и интегрируема на ${\displaystyle G}$, тогда:

${\displaystyle \exists Y\in G\colon \int _{G}f(X)dX=f(Y)\mu (G)}$.

Постоянная функция ${\displaystyle f(X)=c}$ интегрируема на любом измеримом множестве ${\displaystyle G}$, причём:

${\displaystyle \int _{G}f(X)dX=c\cdot \mu (G)}$.

Как следствие, ${\displaystyle \ \int _{G}dX=\mu (G)}$.

Пусть ${\displaystyle D\subset {{\mathbb {R} }^{d-1}}}$ — измеримое множество, ${\displaystyle G=\left\{\left({{x}_{1}},\ldots ,{{x}_{n}}\right):\left({{x}_{1}},\ldots ,{{x}_{n-1}}\right)\in D;\varphi \left({{x}_{1}},\ldots ,{{x}_{n-1}}\right)\leq {{x}_{n}}\leq \psi \left({{x}_{1}},\ldots ,{{x}_{n-1}}\right)\right\}}$ — также измеримое множество, ${\displaystyle f\left(X\right)}$ определена и интегрируема на ${\displaystyle \ G}$. Тогда

• ${\displaystyle \int \limits _{\ \varphi \left({{x}_{1}},\ldots ,{{x}_{n-1}}\right)}^{\psi \left({{x}_{1}},\ldots ,{{x}_{n-1}}\right)}{f\left({{x}_{1}},\ldots ,{{x}_{n}}\right)d{{x}_{n}}}\equiv I\left({{x}_{1}},\ldots ,{{x}_{n-1}}\right)}$ существует всюду на ${\displaystyle D}$, кроме множества ${\displaystyle D_{0}}$ Лебеговой меры нуль (${\displaystyle D_{0}}$ может быть пустым);
• существует ${\displaystyle \int \limits _{D}{\tilde {I}}\left({{x}_{1}},\ldots ,{{x}_{n-1}}\right)d{{x}_{1}}\ldots d{{x}_{n-1}}\equiv \int \limits _{D}{\left[\int \limits _{\,\varphi \left({{x}_{1}},\ldots ,{{x}_{n-1}}\right)}^{\psi \left({{x}_{1}},\ldots ,{{x}_{n-1}}\right)}{f\left({{x}_{1}},\ldots ,{{x}_{n}}\right)d{{x}_{n}}}\right]d{{x}_{1}}\ldots d{{x}_{n-1}}}}$, где

${\displaystyle {\tilde {I}}\left({{x}_{1}},\ldots ,{{x}_{n-1}}\right)\equiv {\begin{cases}I\left({{x}_{1}},\ldots ,{{x}_{n-1}}\right),&({{x}_{1}},\ldots ,{{x}_{n-1}})\in D\backslash D_{0}\\\qquad \quad 0&({{x}_{1}},\ldots ,{{x}_{n-1}})\in D_{0},\end{cases}}}$

называемый повторным интегралом от функции ${\displaystyle f\left({{x}_{1}},\ldots ,{{x}_{n-1}}\right)}$ по множеству ${\displaystyle G}$;

• ${\displaystyle \int \limits _{G}{f\left({{x}_{1}},\ldots ,{{x}_{n}}\right)d{{x}_{1}}\ldots d{{x}_{n}}}=\int \limits _{D}{\left[\int \limits _{\,\varphi \left({{x}_{1}},\ldots ,{{x}_{n-1}}\right)}^{\psi \left({{x}_{1}},\ldots ,{{x}_{n-1}}\right)}{f\left({{x}_{1}},\ldots ,{{x}_{n}}\right)d{{x}_{n}}}\right]d{{x}_{1}}\ldots d{{x}_{n-1}}}}$.

Любой ${\displaystyle n}$-мерный интеграл можно свести к ${\displaystyle n}$ одномерным.

Пусть задано биективное отображение ${\displaystyle {{\mathbb {R} }^{d}}\leftrightarrow {{\mathbb {R} }^{d}}}$, переводящее область ${\displaystyle \ {{D}'}}$ в ${\displaystyle \ D}$:

${\displaystyle \left\{{\begin{aligned}{{t}_{1}}={{\psi }_{1}}\left({{x}_{1}},\ldots ,{{x}_{n}}\right)\\{{t}_{2}}={{\psi }_{2}}\left({{x}_{1}},\ldots ,{{x}_{n}}\right)\\\ldots \\{{t}_{n}}={{\psi }_{n}}\left({{x}_{1}},\ldots ,{{x}_{n}}\right)\\\end{aligned}}\right.}$,

где ${\displaystyle t}$ — «старые» координаты, а ${\displaystyle x}$ — «новые» координаты. Пусть далее функции, задающие отображение, имеют в области ${\displaystyle \ {{D}'}}$ непрерывные частные производные первого порядка, а также ограниченный и отличный от нуля якобиан

${\displaystyle {\frac {D\left(t\right)}{D\left(x\right)}}={\frac {D\left({{t}_{1}},\ldots ,{{t}_{n}}\right)}{D\left({{x}_{1}},\ldots ,{{x}_{n}}\right)}}}$.

Тогда при условии существования интеграла

${\displaystyle \int \limits _{D}{f\left(T\right)dT}=\int {\int \limits _{D}{\ldots \int {f\left({{t}_{1}},\ldots ,{{t}_{n}}\right)d{{t}_{1}}\ldots d{{t}_{n}}}}}}$

справедлива формула замены переменных:

${\displaystyle \int {\int \limits _{D}{\ldots \int {f\left({{t}_{1}},\ldots ,{{t}_{n}}\right)d{{t}_{1}}\ldots d{{t}_{n}}}}}=\int {\int \limits _{{D}'}{\ldots \int {f\left({{\psi }_{1}}\left({{x}_{1}},\ldots ,{{x}_{n}}\right),\ldots ,{{\psi }_{n}}\left({{x}_{1}},\ldots ,{{x}_{n}}\right)\right)\left|{\frac {D\left({{t}_{1}},\ldots ,{{t}_{n}}\right)}{D\left({{x}_{1}},\ldots ,{{x}_{n}}\right)}}\right|d{{x}_{1}}\ldots d{{x}_{n}}}}}}$

Если область интегрирования симметрична относительно начала координат по крайней мере для одной из переменных интегрирования и подынтегральная функция нечётна по этой переменной, интеграл равен нулю, поскольку интегралы по двум половинкам области интегрирования имеют одно и то же абсолютное значение, но противоположные знаки. Если подынтегральная функция чётна по этой переменной, интеграл равен удвоенному интегралу по одной из половинок области интегрирования, поскольку интегралы по каждой из половинок равны.

Например, если функция ${\displaystyle f(x,y)=2\sin(x)-3y^{3}+5}$ интегрируется по области:

${\displaystyle T=\left\{(x,y)\in \mathbf {R} ^{2}\ :\ x^{2}+y^{2}\leq 1\right\}}$ —

кругу радиуса 1 с центром в начале координат, то используя свойство линейности, интеграл можно разложить на три части:

${\displaystyle \iint _{T}(2\sin x-3y^{3}+5)\,dx\,dy=\iint _{T}2\sin x\,dx\,dy-\iint _{T}3y^{3}\,dx\,dy+\iint _{T}5\,dx\,dy}$

${\displaystyle 2\sin x}$ и ${\displaystyle 3y^{3}}$ являются нечётными функциями и, кроме того, очевидно, что диск ${\displaystyle T}$ симметричен как относительно оси ${\displaystyle x}$, так и по оси ${\displaystyle y}$. Таким образом, вклад в конечный результат даёт только константа 5.

В примере для интегрирования функции ${\displaystyle f(x,y,z)=x\exp(y^{2}+z^{2})}$ по сфере радиуса 2 с центром в начале координат:

${\displaystyle T=\left\{(x,y,z)\in \mathbb {R} ^{3}\ :\ x^{2}+y^{2}+z^{2}\leq 4\right\}}$

достаточно проинтегрировать по оси ${\displaystyle x}$, чтобы показать, что интеграл равен 0, поскольку по этой переменной функция нечётна.

Двойным интегралом называют кратный интеграл с ${\displaystyle n=2}$:

${\displaystyle \iint \limits _{D}{f\left(P\right)d\sigma }}$.

Здесь ${\displaystyle \ d\sigma }$ — элемент площади в рассматриваемых координатах.

В прямоугольных координатах: ${\displaystyle \iint \limits _{D}{f\left(x,y\right)dxdy}}$, где ${\displaystyle \ dxdy}$ — элемент площади в прямоугольных координатах.

Пусть функция ${\displaystyle f\left(x,y\right)}$ принимает в области ${\displaystyle D}$ только положительные значения. Тогда двойной интеграл ${\displaystyle \iint \limits _{D}{f\left(x,y\right)d\sigma }::}$ численно равен объёму ${\displaystyle \ V}$ вертикального цилиндрического тела, построенного на основании ${\displaystyle \ D}$ и ограниченного сверху соответствующим куском поверхности ${\displaystyle z=f\left(x,y\right)}$.

В некоторых случаях двойной интеграл проще считать не в прямоугольных, а в полярных координатах, так как при этом может произойти существенное упрощение вида области интегрирования и всего процесса интегрирования в целом.

Для выражения через полярные координаты используется следующая замена переменных:

${\displaystyle x=r\cos \varphi }$, ${\displaystyle y=r\sin \varphi }$.

Модуль якобиана отображения равен ${\displaystyle r}$. Таким образом:

${\displaystyle \iint \limits _{D}{f\left(x,y\right)dxdy}=\iint \limits _{{D}'}{g\left(r,\varphi \right)rdrd\varphi },\quad }$ где ${\displaystyle g\left(r,\varphi \right)=f\left(r\cos \varphi ,r\sin \varphi \right)}$.

Здесь ${\displaystyle rdrd\varphi }$ является элементом площади в полярных координатах.

Наименование величины	Общее выражение	Прямоугольные координаты	Полярные координаты
Площадь плоской фигуры	${\displaystyle S=\iint \limits _{G}{d\sigma }}$	${\displaystyle \iint \limits _{G}{dxdy}}$	${\displaystyle \iint \limits _{G}{rdrd\varphi }}$
Масса тонкой плоской пластинки плотностью ${\displaystyle \mu }$	${\displaystyle m=\iint \limits _{G}{\mu \left(\sigma \right)d\sigma }}$	${\displaystyle \iint \limits _{G}{\mu \left(x,y\right)dxdy}}$	${\displaystyle \iint \limits _{G}{\mu \left(r,\varphi \right)rdrd\varphi }}$
Площадь куска поверхности${\displaystyle ^{1)}}$	${\displaystyle S=\iint \limits _{G}{\frac {d\sigma }{\cos \gamma }}}$	${\displaystyle \iint \limits _{G}{{\sqrt {1+{{\left({\frac {\partial z}{\partial x}}\right)}^{2}}+{{\left({\frac {\partial z}{\partial y}}\right)}^{2}}}}dxdy}}$	${\displaystyle \iint \limits _{G}{{\sqrt {{{r}^{2}}+{{r}^{2}}{{\left({\frac {\partial z}{\partial r}}\right)}^{2}}+{{\left({\frac {\partial z}{\partial \varphi }}\right)}^{2}}}}drd\varphi }}$
Объем цилиндрического тела, стоящего на плоскости ${\displaystyle XOY}$	${\displaystyle V=\iint \limits _{G}{zd\sigma }}$	${\displaystyle \iint \limits _{G}{zdxdy}}$	${\displaystyle \iint \limits _{G}{zrdrd\varphi }}$
Момент инерции плоской фигуры${\displaystyle ^{2)}}$ относительно оси ${\displaystyle O{{Z}^{3)}}}$	${\displaystyle {{I}_{z}}=\iint \limits _{G}{{{r}^{2}}d\sigma }}$	${\displaystyle \iint \limits _{G}{\left({{x}^{2}}+{{y}^{2}}\right)dxdy}}$	${\displaystyle \iint \limits _{G}{{{r}^{3}}drd\varphi }}$
Момент инерции плоской фигуры${\displaystyle ^{2)}}$ относительно оси ${\displaystyle O{{X}^{3)}}}$	${\displaystyle {{I}_{z}}=\iint \limits _{G}{{{y}^{2}}d\sigma }}$	${\displaystyle \iint \limits _{G}{{{y}^{2}}dxdy}}$	${\displaystyle \iint \limits _{G}{{{r}^{3}}{{\sin }^{2}}\varphi drd\varphi }}$
Координаты центра масс однородной пластинки${\displaystyle ^{3)}}$	${\displaystyle {{x}_{c}}={\frac {1}{S}}{\iint \limits _{G}{xd\sigma }}}$ ${\displaystyle {{y}_{c}}={\frac {1}{S}}{\iint \limits _{G}{yd\sigma }}}$	${\displaystyle {\begin{aligned}&{\frac {1}{S}}{\iint \limits _{G}{xdxdy}}\\&{\frac {1}{S}}{\iint \limits _{G}{ydxdy}}\\\end{aligned}}}$	${\displaystyle {\begin{aligned}&{\frac {1}{S}}{\iint \limits _{G}{{{r}^{2}}\cos \varphi drd\varphi }}\\&{\frac {1}{S}}{\iint \limits _{G}{{{r}^{2}}\sin \varphi drd\varphi }}\\\end{aligned}}}$
Примечания	1) Область ${\displaystyle G}$ — проекция на плоскость ${\displaystyle XOY}$; в каждую точку области проектируется только одна точка поверхности; ${\displaystyle \gamma }$ — угол между касательной плоскостью и плоскостью ${\displaystyle XOY}$. 2) Совмещенной с плоскостью ${\displaystyle XOY}$. 3) Или, что то же, относительно центра О.

Тройным интегралом называют кратный интеграл с ${\displaystyle \ d=3}$:

${\displaystyle \iiint \limits _{D}{f\left(P\right)dV}}$

где ${\displaystyle \ dV}$ — элемент объёма в рассматриваемых координатах.

В прямоугольных координатах тройной интеграл имеет следующий вид:

${\displaystyle \iiint \limits _{D}f(x,y,z)\,dx\,dy\,dz}$

где ${\displaystyle \ dxdydz}$ — элемент объёма в прямоугольных координатах.

Аналогично в некоторых случаях тройной интеграл проще считать не в прямоугольных, а в цилиндрических координатах. Применим теорему о замене переменных. Соответствующее переходу преобразование имеет вид:

${\displaystyle \left\{{\begin{aligned}&x=r\cos \varphi \\&y=r\sin \varphi \\&z=h\\\end{aligned}}\right.}$

Модуль якобиана отображения равен ${\displaystyle \ r}$. Таким образом получаем, что

${\displaystyle \iiint \limits _{D}{f\left(x,y,z\right)dxdydz}=\iiint \limits _{{D}'}{f\left(r,\varphi ,h\right)rdrd\varphi dh}}$

где ${\displaystyle rdrd\varphi dh}$ — элемент объёма в цилиндрических координатах.

Кроме цилиндрических можно также переходить и в сферические координаты. Применим теорему о замене переменных. Соответствующее переходу преобразование имеет вид:

${\displaystyle \left\{{\begin{aligned}&x=r\sin \theta \cos \varphi \\&y=r\sin \theta \sin \varphi \\&z=r\cos \theta \\\end{aligned}}\right.}$

Модуль якобиана отображения равен ${\displaystyle \ {{r}^{2}}\sin \theta }$. Таким образом получаем, что

${\displaystyle \iiint \limits _{D}{f\left(x,y,z\right)dxdydz}=\iiint \limits _{{D}''}{f\left(r,\varphi ,\theta \right){{r}^{2}}\sin \theta drd\varphi d\theta }}$

где ${\displaystyle {{r}^{2}}\sin \theta drd\varphi d\theta }$ — элемент объёма в сферических координатах.

Наименование величины	Общее выражение	Прямоугольные координаты	Цилиндрические координаты	Сферические координаты
Объем тела	${\displaystyle V=\iiint \limits _{G}{dV}}$	${\displaystyle \iiint \limits _{G}{dxdydz}}$	${\displaystyle \iiint \limits _{G}{rdrd\varphi dh}}$	${\displaystyle \iiint \limits _{G}{{{\rho }^{2}}\sin \theta d\rho d\varphi d\theta }}$
Момент инерции геометрического тела относительно оси ${\displaystyle OZ}$	${\displaystyle {{I}_{z}}=\iiint \limits _{G}{{{r}^{2}}dV}}$	${\displaystyle \iiint \limits _{G}{\left({{x}^{2}}+{{y}^{2}}\right)dxdydz}}$	${\displaystyle \iiint \limits _{G}{{{r}^{3}}drd\varphi dh}}$	${\displaystyle \iiint \limits _{G}{{{\rho }^{4}}{{\sin }^{3}}\theta d\rho d\varphi d\theta }}$
Масса физического тела с плотностью ${\displaystyle \mu }$	${\displaystyle m=\iiint \limits _{G}{\mu dV}}$	${\displaystyle \iiint \limits _{G}{\mu dxdydz}}$	${\displaystyle \iiint \limits _{G}{\mu rdrd\varphi dh}}$	${\displaystyle \iiint \limits _{G}{\mu {{\rho }^{2}}\sin \theta d\rho d\varphi d\theta }}$
Координаты центра масс однородного тела	${\displaystyle {\begin{aligned}&{{x}_{c}}={\frac {1}{V}}{\iiint \limits _{G}{xdV}}\\&{{y}_{c}}={\frac {1}{V}}{\iiint \limits _{G}{ydV}}\\&{{z}_{c}}={\frac {1}{V}}{\iiint \limits _{G}{zdV}}\\\end{aligned}}}$	${\displaystyle {\begin{aligned}&{\frac {1}{V}}{\iiint \limits _{G}{xdxdydz}}\\&{\frac {1}{V}}{\iiint \limits _{G}{ydxdydz}}\\&{\frac {1}{V}}{\iiint \limits _{G}{zdxdydz}}\\\end{aligned}}}$	—	—

• Выгодский, М. Я. Дифференцирование и интегрирование функций нескольких аргументов // Справочник по высшей математике. — М.: Астрель, АСТ, 2005. — 991 с. — 10 000 экз. — ISBN 5-17-012238-1, 5-271-03651-0.
• Ильин, В. А., Позняк, Э. Г. Глава 2. Двойные и n-кратные интегралы // Основы математического анализа. — 4. — М.: ФИЗМАТЛИТ, 2001. — Т. 2. — 464 с. — (Курс высшей математики и математической физики). — 5000 экз. — ISBN 5-9221-0131-5.
• Кудрявцев, Л. Д. Глава 6. Интегральное исчисление функций многих переменных // Курс математического анализа. — М.: Высшая школа, 1981. — Т. 2. — 584 с.
• Будак, Б. М., Фомин С. В. Кратные интегралы и ряды. — М.: Наука, 1967. — 608 с.