Поверхностные интегралы 🌐

В прошлом уроке ты научился интегрировать вдоль произвольной кривой — одномерного объекта, изогнутого в двумерном или трёхмерном пространстве. Но мир не ограничивается кривыми линиями: купол собора, обшивка самолёта, поверхность озера, мембрана барабана, кожа трёхмерной модели, полученная сканером глубины, — всё это двумерные поверхности, вложенные в трёхмерное пространство. И к ним возникают ровно те же самые вопросы, что и к кривым, только теперь в двух измерениях сразу: какова масса этой оболочки, если её плотность меняется от точки к точке? Сколько жидкости или излучения протекает сквозь эту мембрану за единицу времени? Обычный двойной интеграл по плоской области здесь не годится напрямую — поверхность может быть изогнутой, наклонной, вообще не лежать ни в одной координатной плоскости. Нужен новый инструмент — поверхностный интеграл, — который сегодня мы построим шаг за шагом.

Как и в случае с кривыми, поверхностных интегралов тоже два, и они принципиально разной природы. Интеграл первого рода, $\iint_S f\,dS$, интегрирует скалярную функцию (например, плотность материала) по площади поверхности, и ему совершенно всё равно, с какой стороны на поверхность смотреть, — результат один и тот же. Интеграл второго рода, поток векторного поля через поверхность, интегрирует уже не скаляр, а проекцию векторного поля на нормаль к поверхности, и здесь выбор стороны — ориентации — решающим образом влияет на знак ответа. Это прямая параллель с прошлым уроком: там направление обхода кривой определяло знак интеграла второго рода, здесь эту роль играет выбор нормали.

А вот и обещанная связь с прикладным анализом данных. Представь трёхмерную модель объекта, восстановленную по облаку точек лидара или по серии фотографий в задаче 3D-реконструкции сцены. Такая модель — это, по сути, дискретная поверхность (полигональная сетка из треугольников), и очень многие задачи компьютерного зрения сводятся именно к поверхностным интегралам, только записанным на дискретном языке треугольников вместо непрерывного языка $dS$. Хочешь узнать суммарную яркость, зафиксированную камерой по всей видимой стороне объекта? Это интеграл первого рода с весовой функцией «яркость». Хочешь понять, сколько света (в физическом смысле — энергии фотонов) приходит на точку поверхности со всех направлений полусферы над ней — базовая задача расчёта освещённости в фотореалистичном воссоздании изображений и в алгоритмах восстановления материалов по фотографиям? Это в точности поток векторного поля через полусферу направлений, то есть поверхностный интеграл второго рода. А когда в графовых нейросетях, работающих прямо на поверхностях (поверхностные свёрточные сети, геодезические свёртки по полигональной сетке), признаки соседних вершин усредняются с весами, пропорциональными площади прилегающего кусочка сетки, — это дискретная аппроксимация того самого элемента площади $dS$, который мы сегодня строго выведем.

План на сегодня такой. Сначала разберёмся с интегралом первого рода — массой и площадью искривлённой оболочки — и параллельно освоим параметризацию поверхности как способ её описания. Затем выведем универсальную формулу элемента площади $dS$ через частные производные — самый технический, но и самый важный кусок урока, без которого ни один поверхностный интеграл посчитать не получится. После этого перейдём к интегралу второго рода: потоку векторного поля и понятию ориентации поверхности. И завершим тем, как оба типа интеграла на практике сводятся к обычному двойному интегралу по плоской области параметров — тем самым инструментом, которым ты пользуешься уже не первый месяц.

История: откуда это взялось?

Понятие потока через поверхность родилось из гидродинамики и электродинамики рубежа XVIII–XIX веков. Леонард Эйлер, разрабатывая уравнения движения жидкости, естественным образом пришёл к вопросу: сколько жидкости протекает сквозь произвольно изогнутую воображаемую перегородку внутри потока за единицу времени? Ответ требовал суммирования (интегрирования) нормальной составляющей скорости жидкости по всей площади этой перегородки — ровно то, что мы сегодня называем потоком векторного поля через поверхность. Чуть позже похожая идея понадобилась Шарлю Кулону и его последователям для описания электрического поля: поток напряжённости через замкнутую поверхность оказался напрямую связан с зарядом внутри неё — эта связь легла в основу закона Гаусса, к строгому доказательству которого мы вернёмся в следующем уроке про теорию поля.

Строгий математический аппарат поверхностных интегралов оформился в первой половине XIX века в работах Карла Фридриха Гаусса и русского математика Михаила Васильевича Остроградского. Остроградский в 1828 году (независимо и, по мнению многих историков науки, даже раньше западных коллег) доказал теорему, связывающую поток векторного поля через замкнутую поверхность с интегралом от дивергенции этого поля по объёму, ограниченному поверхностью, — сегодня эта теорема носит имя Остроградского-Гаусса и станет главным героем следующего урока. Но прежде чем формулировать такую теорему, нужно было сначала строго определить, что вообще такое поверхностный интеграл, — эту фундаментальную работу проделал сам Гаусс, введя понятие элемента площади поверхности через параметрическое представление и векторное произведение касательных векторов, которым мы пользуемся и сегодня, спустя два столетия почти без изменений.

Отдельного упоминания заслуживает вопрос ориентируемости поверхностей — казалось бы, чисто техническая деталь оказалась одним из самых глубоких открытий топологии XIX века. В 1858 году немецкий математик Август Мёбиус и независимо от него Иоганн Листинг обнаружили поверхности, у которых попросту нет двух сторон, — знаменитый лист Мёбиуса, полоску бумаги, перекрученную на пол-оборота перед склейкой. На такой поверхности невозможно последовательно, без разрывов, выбрать нормаль так, чтобы обойти всю поверхность и вернуться в исходную точку с той же самой нормалью, с которой начал. Это открытие показало, что понятие ориентации поверхности — не формальность, а содержательное свойство, которое может как выполняться, так и не выполняться, и без него определение потока векторного поля через поверхность попросту повисает в воздухе.

Поверхностный интеграл первого рода и параметризация поверхности

Интуиция: оболочка переменной плотности

Давай разберёмся с самого простого вопроса. Представь тонкую металлическую оболочку — например, купол или обшивку корпуса — произвольной формы $S$ в трёхмерном пространстве. Плотность материала (масса на единицу площади) в разных точках оболочки разная: где-то материал толще и тяжелее, где-то — тоньше. Обозначим эту поверхностную плотность через $f(x,y,z)$. Как найти полную массу оболочки?

Действуй ровно так же, как при выводе любого интеграла: разбей поверхность $S$ на очень много крошечных кусочков площадью $\Delta S_i$ каждый. На каждом таком бесконечно малом кусочке плотность практически постоянна, поэтому масса кусочка приближённо равна $f(x_i,y_i,z_i)\cdot\Delta S_i$. Просуммируй эти вклады по всем кусочкам и устреми размер разбиения к нулю — получишь поверхностный интеграл первого рода. Идея дословно та же, что и с криволинейным интегралом первого рода в прошлом уроке, только теперь суммирование идёт не по длине кривой, а по площади поверхности.

Определение

Определение (поверхностный интеграл первого рода). Пусть $S$ — гладкая поверхность в пространстве, а функция $f(x,y,z)$ определена во всех её точках. Разобьём $S$ на $n$ малых кусочков площадью $\Delta S_1,\Delta S_2,\dots,\Delta S_n$, выберем на каждом кусочке произвольную точку $(x_i,y_i,z_i)$ и составим интегральную сумму $\sum_{i=1}^n f(x_i,y_i,z_i)\Delta S_i$. Поверхностным интегралом первого рода (по площади поверхности) называется предел этой суммы при неограниченном измельчении разбиения:

$$\iint_S f(x,y,z)\,dS = \lim_{\max\Delta S_i\to0}\sum_{i=1}^n f(x_i,y_i,z_i)\Delta S_i$$

Ключевое свойство: этот интеграл не зависит от выбора стороны поверхности — площадь кусочка $\Delta S_i$ всегда положительна независимо от того, с какой стороны на поверхность смотреть. Если $S^-$ обозначает ту же поверхность с противоположной ориентацией, то $\iint_{S^-}f\,dS=\iint_S f\,dS$.

Как и в случае с кривыми, частные случаи этого интеграла тебе уже интуитивно знакомы. При $f(x,y,z)\equiv1$ интеграл $\iint_S 1\,dS$ даёт в точности площадь самой поверхности $S$. А если $f$ — поверхностная плотность материала, интеграл даёт полную массу оболочки $m=\iint_S f\,dS$.

Параметризация поверхности

Чтобы вообще вычислять такой интеграл, нужен способ описать поверхность математически. Кривую в прошлом уроке мы описывали одним параметром $t$: точка на кривой — это функция одной переменной времени. Поверхность — объект двумерный, поэтому ей нужны уже два параметра. Представь навигационный трек спутникового позиционирования, но теперь не пешехода, а сканирующего лазера: он записывает координаты $(x,y,z)$ каждой точки поверхности объекта, перебирая её как бы по строкам и столбцам — по двум независимым направлениям сканирования $u$ и $v$. Это и есть параметризация поверхности:

где $D$ — некоторая плоская область изменения параметров. Самый частый частный случай — явный график функции $z=g(x,y)$, где роль параметров $u,v$ играют сами $x$ и $y$: $\vec r(x,y)=(x,\,y,\,g(x,y))$. Но бывают и поверхности, которые невозможно задать явным графиком (например, полную сферу нельзя выразить как единую функцию $z=g(x,y)$, поскольку каждой паре $(x,y)$ внутри экватора соответствуют сразу две точки сферы — верхняя и нижняя), и тогда параметризация через два независимых угла или две криволинейные координаты — единственный практичный способ описания.

Разбор примеров

Пример 1 (лёгкий). Найти площадь части плоскости $z=2x+2y$ над квадратом $D=[0,1]\times[0,1]$.

Это явный график $z=g(x,y)=2x+2y$. Чтобы найти площадь, нам понадобится элемент $dS$ — формулу для него мы строго выведем в следующем разделе, а пока используем готовый результат $dS=\sqrt{1+g_x^2+g_y^2}\,dx\,dy$, знакомый тебе, возможно, из курса физики. Здесь $g_x=2$, $g_y=2$, значит $dS=\sqrt{1+4+4}\,dx\,dy=3\,dx\,dy$.

Ответ: площадь равна $3$.

Пример 2 (средний). Найти массу верхней полусферы $x^2+y^2+z^2=R^2$, $z\ge0$, с поверхностной плотностью $\sigma(x,y,z)=z$.

Параметризуем сферу сферическими координатами: $x=R\sin\varphi\cos\theta$, $y=R\sin\varphi\sin\theta$, $z=R\cos\varphi$, где $\varphi\in[0,\pi/2]$ (верхняя половина) и $\theta\in[0,2\pi]$. Для сферы радиуса $R$ элемент площади (мы выведем это в следующем разделе) равен $dS=R^2\sin\varphi\,d\varphi\,d\theta$. Плотность на поверхности: $\sigma=z=R\cos\varphi$.

Внутренний интеграл: $\int_0^{\pi/2}\sin\varphi\cos\varphi\,d\varphi=\left[\dfrac{\sin^2\varphi}2\right]_0^{\pi/2}=\dfrac12$. Тогда $m=R^3\cdot2\pi\cdot\dfrac12=\pi R^3$.

Ответ: $m=\pi R^3$. Заметь, насколько чистый получился результат — совпадение не случайно: для сферы плотность $z$ и элемент площади согласованно зависят от угла $\varphi$ так, что все «лишние» коэффициенты аккуратно сокращаются.

Пример 3 (сложный, машинное обучение). Часть параболоида $z=x^2+y^2$ над кругом $x^2+y^2\le1$ моделирует куполообразную поверхность объекта после 3D-реконструкции (типичная форма выпуклой «чаши» на карте глубины). На поверхности задана интенсивность признака (например, уверенность модели в оценке глубины в данной точке) $f(x,y,z)=x^2+y^2$ — интенсивность растёт от центра купола к его краю. Найти суммарное («агрегированное») значение этого признака по всей поверхности купола.

Элемент площади для параболоида: $g=x^2+y^2$, $g_x=2x$, $g_y=2y$, $dS=\sqrt{1+4x^2+4y^2}\,dx\,dy$. Переходим к полярным координатам $x=r\cos\theta$, $y=r\sin\theta$, $r\in[0,1]$, $\theta\in[0,2\pi]$: $dS=\sqrt{1+4r^2}\,r\,dr\,d\theta$, а $f=x^2+y^2=r^2$.

Сделаем замену $u=1+4r^2$, $du=8r\,dr$, тогда $r^2=\dfrac{u-1}4$ и $r^3\,dr=r^2\cdot r\,dr=\dfrac{u-1}4\cdot\dfrac{du}8=\dfrac{u-1}{32}\,du$:

При $r=1$, $u=5$: $\dfrac25\cdot5^{5/2}-\dfrac23\cdot5^{3/2}=\dfrac25\cdot25\sqrt5-\dfrac23\cdot5\sqrt5=10\sqrt5-\dfrac{10\sqrt5}3=\dfrac{20\sqrt5}3$. При $r=0$, $u=1$: $\dfrac25-\dfrac23=-\dfrac4{15}$. Разность, делённая на $32$:

Умножаем на $2\pi$:

Ответ: $\dfrac{\pi(25\sqrt5+1)}{60}\approx2{,}979$. Такая агрегированная величина — это ровно та операция, которую выполняют графовые нейросети, работающие прямо на поверхности 3D-объекта (поверхностные и геодезические свёрточные сети): каждый узел сетки вносит вклад, пропорциональный площади прилегающего к нему кусочка поверхности, и суммарная «энергия» признака по всему объекту — это в точности дискретный аналог интеграла первого рода.

Почему это важно

Поверхностный интеграл первого рода — естественное расширение двойного интеграла на случай, когда «плоская область интегрирования» сама изогнута в пространстве. Он отвечает на вопросы вида «сколько всего» — масса, площадь, суммарная интенсивность признака, — и ему в принципе всё равно, с какой стороны на поверхность смотреть. Это делает его удобным инструментом всякий раз, когда речь идёт об агрегировании скалярной величины по искривлённому объекту, будь то физическая оболочка или трёхмерная модель, восстановленная алгоритмом компьютерного зрения. Но чтобы реально вычислить любой из примеров выше, нам пришлось воспользоваться готовой формулой для $dS$ — настало время вывести её строго.

Элемент площади поверхности: векторное произведение частных производных

Интуиция: параллелограмм из касательных векторов

Давай разберёмся, откуда вообще берётся формула для $dS$. Возьми параметризацию $\vec r(u,v)=(x(u,v),y(u,v),z(u,v))$ и зафиксируй малое приращение параметров $\Delta u$ и $\Delta v$ вокруг некоторой точки $(u_0,v_0)$. Если менять только $u$, точка на поверхности сдвинется приблизительно на вектор $\vec r_u\,\Delta u$, где $\vec r_u=\dfrac{\partial\vec r}{\partial u}$ — частная производная по $u$, вектор, касательный к поверхности вдоль «линии постоянного $v$». Аналогично, если менять только $v$, точка сдвинется на вектор $\vec r_v\,\Delta v$. Малый криволинейный кусочек поверхности, соответствующий прямоугольнику $\Delta u\times\Delta v$ в плоскости параметров, приближённо превращается в плоский параллелограмм, натянутый на векторы $\vec r_u\,\Delta u$ и $\vec r_v\,\Delta v$.

А площадь параллелограмма, натянутого на два вектора, — это, как ты помнишь из курса линейной алгебры и аналитической геометрии, в точности модуль их векторного произведения. Значит:

Устремляя $\Delta u,\Delta v\to0$, получаем точную формулу для элемента площади.

Определение

Формула элемента площади поверхности. Пусть поверхность $S$ задана параметрически: $\vec r(u,v)=(x(u,v),y(u,v),z(u,v))$, $(u,v)\in D$, с непрерывно дифференцируемыми компонентами. Введём касательные векторы $\vec r_u=(x_u,y_u,z_u)$ и $\vec r_v=(x_v,y_v,z_v)$. Тогда элемент площади поверхности равен

$$dS = \bigl|\vec r_u\times\vec r_v\bigr|\,du\,dv$$

Вектор $\vec n_0=\vec r_u\times\vec r_v$ (ненормированный) при этом перпендикулярен поверхности в каждой точке — он и задаёт направление нормали, о роли которой мы поговорим в следующем разделе. В координатной форме через компоненты (якобианы пар координат по $u,v$):

$$\vec r_u\times\vec r_v=\left(\frac{\partial(y,z)}{\partial(u,v)},\ \frac{\partial(z,x)}{\partial(u,v)},\ \frac{\partial(x,y)}{\partial(u,v)}\right)$$

Частный, но крайне важный случай — явный график $z=g(x,y)$, где роль параметров играют сами $x,y$: $\vec r(x,y)=(x,y,g(x,y))$. Тогда $\vec r_x=(1,0,g_x)$, $\vec r_y=(0,1,g_y)$, и их векторное произведение

Модуль этого вектора равен $\sqrt{g_x^2+g_y^2+1}$, откуда получаем уже знакомую тебе по прошлым примерам формулу:

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

Разбор примеров

Пример 1 (лёгкий). Найти вектор нормали (ненормированный) и элемент площади $dS$ для плоскости $z=3x-2y+5$.

Это явный график с $g_x=3$, $g_y=-2$. Вектор нормали $\vec r_x\times\vec r_y=(-3,\,2,\,1)$. Модуль: $\sqrt{9+4+1}=\sqrt{14}$.

Ответ: $\vec n_0=(-3,2,1)$, $dS=\sqrt{14}\,dx\,dy$.

Пример 2 (средний). Вывести элемент площади для боковой поверхности цилиндра $x^2+y^2=R^2$, $z\in[0,h]$, и найти площадь боковой поверхности при $R=3$, $h=4$.

Параметризуем цилиндр углом и высотой: $x=R\cos\theta$, $y=R\sin\theta$, $z=z$, $(\theta,z)\in[0,2\pi]\times[0,h]$. Тогда $\vec r_\theta=(-R\sin\theta,\,R\cos\theta,\,0)$, $\vec r_z=(0,0,1)$.

Модуль: $\sqrt{R^2\cos^2\theta+R^2\sin^2\theta}=R$. Значит $dS=R\,d\theta\,dz$ — постоянный элемент площади (что логично: боковая поверхность цилиндра разворачивается в прямоугольник без искажений). Площадь при $R=3$, $h=4$:

Ответ: $dS=R\,d\theta\,dz$, площадь $=24\pi\approx75{,}40$.

Пример 3 (сложный, машинное обучение). Вывести элемент площади для сферы радиуса $R$ в параметризации $x=R\sin\varphi\cos\theta$, $y=R\sin\varphi\sin\theta$, $z=R\cos\varphi$ (использовавшейся в примере 2 предыдущего раздела без доказательства), а затем найти, какую долю всей площади сферы покрывает конус обзора датчика (например, широкоугольной камеры или лидара) с половинным углом раскрыва $60°=\pi/3$.

Находим касательные векторы: $\vec r_\varphi=(R\cos\varphi\cos\theta,\,R\cos\varphi\sin\theta,\,-R\sin\varphi)$, $\vec r_\theta=(-R\sin\varphi\sin\theta,\,R\sin\varphi\cos\theta,\,0)$. Вычисляя векторное произведение покомпонентно, получаем

Последний множитель в скобках — единичный вектор (сумма квадратов его компонент равна $\sin^2\varphi+\cos^2\varphi=1$), поэтому модуль всего произведения равен просто $R^2\sin\varphi$ (при $\varphi\in[0,\pi]$ синус неотрицателен). Значит

— это и есть формула, которой мы уже несколько раз пользовались без доказательства. Теперь применим её: площадь конической «шапочки» сферы единичного радиуса ($R=1$) при $\varphi\in[0,\pi/3]$, $\theta\in[0,2\pi]$:

Площадь всей сферы единичного радиуса равна $4\pi$. Доля покрытия:

Ответ: $dS=R^2\sin\varphi\,d\varphi\,d\theta$; конус обзора с половинным углом $60°$ покрывает ровно четверть сферы. Похожий расчёт используется на практике при проектировании массивов датчиков (например, лидаров кругового обзора или камер для панорамной 3D-реконструкции сцены): зная угол обзора одного сенсора, легко посчитать, сколько таких сенсоров нужно расставить, чтобы без пробелов покрыть всю сферу направлений вокруг объекта.

Почему это важно

Формула $dS=|\vec r_u\times\vec r_v|\,du\,dv$ — это фундамент, на котором держится вообще любой численный расчёт поверхностного интеграла. Без неё интеграл первого рода остаётся красивой, но невычислимой абстракцией. С ней — любая гладкая поверхность, какой бы причудливой ни была её форма, сводится к паре чисел на входе (значения $u$ и $v$) и одной формуле для локального «растяжения» площади. Ровно та же логика, кстати, лежит в основе того, как в компьютерной графике и обработке трёхмерных сеток вычисляют площадь треугольника полигональной сетки по его вершинам — площадь треугольника равна половине модуля векторного произведения двух его сторон, что является простейшим дискретным аналогом только что выведенной формулы.

Поверхностный интеграл второго рода: ориентация и поток векторного поля

Интуиция: сколько поля протекает сквозь сетку

Теперь представь другую картину: через поверхность $S$ «протекает» векторное поле $\vec F(x,y,z)=(P,Q,R)$ — это может быть реальный поток жидкости, электрическое поле, поле скоростей ветра или, как мы увидим в примерах, поле яркости в задаче рендеринга. Вопрос: сколько «вещества» этого поля проходит сквозь поверхность $S$ за единицу времени (или, если поле стационарно, — какой суммарный поток через неё)?

Ключевая идея: сквозь поверхность «проходит» не всё поле целиком, а только та его часть, что направлена перпендикулярно поверхности, — составляющая вдоль нормали. Составляющая поля, направленная вдоль самой поверхности («по касательной»), никуда сквозь неё не протекает — она просто скользит вдоль плёнки, не пересекая её. Поэтому на каждом крошечном кусочке площадью $\Delta S$ вклад в поток равен скалярному произведению поля на единичную нормаль, умноженному на площадь: $\vec F\cdot\vec n\,\Delta S$. Суммируя (интегрируя) эти вклады по всей поверхности, получаем поверхностный интеграл второго рода.

Ориентация поверхности

Но у любой точки гладкой поверхности есть не одна, а сразу две возможные единичные нормали, направленные в противоположные стороны, — и заранее непонятно, какую из них выбрать. Выбор одной из двух нормалей, непрерывно (без скачков) продолжаемой по всей поверхности, называется ориентацией поверхности. Для замкнутой поверхности (например, сферы) обычно выбирают внешнюю нормаль (наружу от тела) — она задаёт так называемую положительную ориентацию. Для незамкнутого куска поверхности (например, графика функции) чаще выбирают нормаль «вверх» — с положительной $z$-компонентой.

Важная тонкость: не каждая поверхность вообще допускает такой непрерывный выбор нормали. Классический контрпример — лист Мёбиуса: если начать двигать нормаль вдоль его срединной линии, обойти весь лист и вернуться в исходную точку, нормаль окажется развёрнутой в противоположную сторону. Такие поверхности называются неориентируемыми, и поток векторного поля через них попросту не определён. К счастью, почти все поверхности, с которыми ты столкнёшься на практике (сферы, плоскости, графики функций, цилиндры, поверхности 3D-моделей без «перекруток»), ориентируемы, и выбор нормали для них — вопрос соглашения, а не принципиальной невозможности.

Определение

Определение (поверхностный интеграл второго рода, или поток). Пусть $S$ — гладкая ориентированная поверхность с выбранной единичной нормалью $\vec n(x,y,z)$, а $\vec F=(P,Q,R)$ — векторное поле, определённое на $S$. Потоком поля $\vec F$ через поверхность $S$ называется поверхностный интеграл первого рода от скалярного произведения:

$$\Phi=\iint_S \vec F\cdot\vec n\,dS=\iint_S \bigl(P\,dy\,dz+Q\,dz\,dx+R\,dx\,dy\bigr)$$

где запись $P\,dy\,dz+Q\,dz\,dx+R\,dx\,dy$ — традиционная координатная форма записи того же самого интеграла (проекции ориентированной площадки на координатные плоскости).

Ключевое свойство: при смене ориентации на противоположную ($\vec n\to-\vec n$) поток меняет знак:

$$\iint_{S^-}\vec F\cdot\vec n\,dS=-\iint_S \vec F\cdot\vec n\,dS$$

Обрати внимание на структурную параллель с прошлым уроком: там смена направления обхода кривой меняла знак криволинейного интеграла второго рода, здесь смена стороны поверхности меняет знак поверхностного интеграла второго рода. И там, и там причина одна и та же — интеграл первого рода интегрирует «безориентационную» величину (длину или площадь, всегда положительную), а интеграл второго рода интегрирует величину, которая по своей физической природе обязана менять знак при развороте — работу против направления движения или поток против направления нормали.

Вычисление через параметризацию

Если поверхность задана параметризацией $\vec r(u,v)$, а векторное произведение $\vec r_u\times\vec r_v$ указывает в ту же сторону, что и выбранная ориентация $\vec n$, то поток вычисляется без явного вычисления модуля нормали — они взаимно сокращаются:

Если же $\vec r_u\times\vec r_v$ направлен противоположно выбранной ориентации, результат нужно взять со знаком минус (или, что то же самое, поменять местами $u$ и $v$ в параметризации). Для явного графика $z=g(x,y)$ с ориентацией «нормаль вверх» подставляем найденный ранее вектор $(-g_x,-g_y,1)$ и получаем удобную рабочую формулу:

Разбор примеров

Пример 1 (лёгкий). Найти поток постоянного поля $\vec F=(1,2,3)$ через часть плоскости $z=2x+2y$ над квадратом $D=[0,1]\times[0,1]$, ориентация — нормаль вверх.

Здесь $g=2x+2y$, $g_x=2$, $g_y=2$. По формуле для явного графика: $P=1$, $Q=2$, $R=3$.

Ответ: $\Phi=-3$. Отрицательный знак означает, что выбранное поле в среднем направлено противоположно нормали «вверх» на этой наклонной плоскости — сам факт, что плоскость наклонена, а не горизонтальна, влияет на знак, даже несмотря на то что $z$-компонента поля положительна.

Пример 2 (средний). Найти поток радиального поля $\vec F=(x,y,z)$ через сферу радиуса $R=3$ с внешней нормалью.

На сфере единичная внешняя нормаль равна $\vec n=\dfrac{(x,y,z)}{R}$ (радиус-вектор, делённый на длину). Скалярное произведение на поверхности сферы:

— оказывается постоянным во всех точках сферы! Это значит, что интеграл превращается в произведение постоянной на площадь сферы:

Ответ: $\Phi=108\pi\approx339{,}29$. Этот приём — заметить, что скалярное произведение постоянно на всей сфере благодаря симметрии, — многократно упростит тебе жизнь и в следующем уроке про теорию поля: он же лежит в основе закона Гаусса.

Пример 3 (сложный, машинное обучение). В компьютерной графике и обратном рендеринге (восстановлении материалов и освещения по фотографиям — важная задача 3D-реконструкции сцен) освещённость точки поверхности вычисляется интегрированием входящего излучения по полусфере направлений над этой точкой, с весом, зависящим от косинуса угла падения (закон Ламберта). Рассмотрим упрощённую модель: единичная верхняя полусфера $x^2+y^2+z^2=1$, $z\ge0$, направленность излучения задаётся полем $\vec F=(0,0,z)$ (энергия пропорциональна высоте направления над горизонтом), ориентация — внешняя нормаль. Найти суммарный поток («полную энергию, приходящую со всех направлений верхней полусферы»).

На единичной сфере внешняя нормаль $\vec n=(x,y,z)$. Скалярное произведение: $\vec F\cdot\vec n=(0,0,z)\cdot(x,y,z)=z^2$.

Параметризуем сферу как в предыдущем разделе: $z=\cos\varphi$, $dS=\sin\varphi\,d\varphi\,d\theta$, $\varphi\in[0,\pi/2]$, $\theta\in[0,2\pi]$:

Ответ: $\Phi=\dfrac{2\pi}3\approx2{,}094$. Это, по сути, дискретизированный вариант классического интеграла Ламберта по полусфере, который лежит в основе расчёта диффузного освещения в компьютерной графике и в нейросетевых алгоритмах обратного рендеринга, восстанавливающих материалы объекта по фотографиям: чтобы понять, сколько света в среднем «накапливает» точка поверхности с матовым (ламбертовским) отражением, нужно проинтегрировать входящую яркость с весом косинуса угла по всей полусфере направлений над этой точкой — в точности та операция, которую мы только что проделали.

Почему это важно

Поток векторного поля через поверхность — это язык, на котором строго формулируются вопросы вида «сколько чего-то проходит сквозь границу». В физике это поток жидкости, электрического или магнитного поля. В компьютерном зрении и графике — это поток излучения через поверхность объекта, лежащий в основе расчёта освещённости и восстановления материалов по изображениям. А как ты увидел на примере с радиальным полем через сферу, симметрия поверхности часто позволяет свести громоздкий интеграл к произведению константы на площадь — приём, который в следующем уроке превратится в один из самых красивых результатов всего курса, связывающий поток через замкнутую поверхность с поведением поля внутри неё.

Сведение обоих интегралов к двойному интегралу по параметрам

Интуиция: поверхность как плоская область, «натянутая» на пространство

Мы уже, по сути, применяли этот метод в каждом примере выше — настало время собрать его в единую схему и потренироваться на разных типах поверхностей. Идея та же, что и с параметризацией кривой в прошлом уроке: если поверхность описана как $\vec r(u,v)$, $(u,v)\in D$, то любой поверхностный интеграл, каким бы пугающим он ни казался в абстрактном виде, превращается в обычный двойной интеграл по плоской области $D$ — тот самый, которым ты уверенно пользуешься начиная с середины курса. Вся «геометрическая сложность» изогнутой поверхности при этом упаковывается в один множитель: модуль векторного произведения $|\vec r_u\times\vec r_v|$ для интеграла первого рода или сам вектор $\vec r_u\times\vec r_v$ для интеграла второго рода.

Определение

Сводные формулы вычисления поверхностных интегралов. Пусть поверхность $S$ задана параметризацией $\vec r(u,v)=(x(u,v),y(u,v),z(u,v))$, $(u,v)\in D$. Тогда:

Интеграл первого рода:

$$\iint_S f(x,y,z)\,dS=\iint_D f\bigl(\vec r(u,v)\bigr)\,\bigl|\vec r_u\times\vec r_v\bigr|\,du\,dv$$

Интеграл второго рода (поток), при согласовании $\vec r_u\times\vec r_v$ с выбранной ориентацией:

$$\iint_S \vec F\cdot d\vec S=\iint_D \vec F\bigl(\vec r(u,v)\bigr)\cdot\bigl(\vec r_u\times\vec r_v\bigr)\,du\,dv$$

Для явного графика $z=g(x,y)$ с ориентацией «вверх» частные формулы:

$$dS=\sqrt{1+g_x^2+g_y^2}\,dx\,dy, \qquad \iint_S(P\,dy\,dz+Q\,dz\,dx+R\,dx\,dy)=\iint_D\bigl[-P\,g_x-Q\,g_y+R\bigr]dx\,dy$$

Обрати внимание: как и с параметризацией кривой, сама параметризация одной и той же поверхности не единственна — можно «просканировать» ту же поверхность разными параметрами (например, разными скоростями изменения $u$ и $v$), и результат интеграла первого рода от этого не изменится (для интеграла второго рода важно лишь сохранить согласованность направления $\vec r_u\times\vec r_v$ с ориентацией). Это прямое следствие формулы замены переменных в кратном интеграле, о которой ты уже знаешь из темы двойных и тройных интегралов.

Разбор примеров

Пример 1 (лёгкий). Найти площадь боковой поверхности конуса $z=\sqrt{x^2+y^2}$ над кругом $x^2+y^2\le4$.

Здесь $g=\sqrt{x^2+y^2}$, значит $g_x=\dfrac x{\sqrt{x^2+y^2}}$, $g_y=\dfrac y{\sqrt{x^2+y^2}}$, и $g_x^2+g_y^2=\dfrac{x^2+y^2}{x^2+y^2}=1$ — тождественно, в любой точке конуса (кроме вершины). Значит $dS=\sqrt2\,dx\,dy$ — постоянный множитель, не зависящий от радиуса. Площадь:

Ответ: $4\sqrt2\,\pi\approx17{,}77$.

Пример 2 (средний). Найти массу части боковой поверхности цилиндра $x^2+y^2=4$, $0\le z\le3$, с плотностью $\sigma(x,y,z)=z$.

Из примера 2 предыдущего раздела для цилиндра радиуса $R=2$ элемент площади равен $dS=2\,d\theta\,dz$. Плотность на поверхности $\sigma=z$.

Ответ: $18\pi\approx56{,}55$.

Пример 3 (сложный, машинное обучение). Функция потерь двух весов вблизи седловой точки описывается моделью гиперболического параболоида — «седла» $z=w_1w_2$ (классическая локальная форма ландшафта функции потерь рядом с седловой точкой, где одно направление — «долина», а перпендикулярное — «хребет»). На этой поверхности задано поле $\vec F=(w_1,\,w_2,\,z)$ — радиус-вектор в расширенном пространстве «веса + значение функции потерь». Найти поток этого поля через кусок седла над квадратом $D=[0,1]\times[0,1]$, ориентация вверх.

Здесь $g=w_1w_2$, $g_{w_1}=w_2$, $g_{w_2}=w_1$. Обозначим $P=w_1$, $Q=w_2$, $R=z=w_1w_2$. По формуле для явного графика:

Поскольку область $D=[0,1]\times[0,1]$ — прямоугольник, интеграл распадается на произведение:

Ответ: $\Phi=-\dfrac14=-0{,}25$. Отрицательный знак здесь содержателен: он говорит о том, что вблизи седловой точки радиальное поле в среднем направлено «внутрь» относительно выбранной нормали седла — иначе говоря, седловая структура ландшафта устроена так, что она «затягивает» поле с одной комбинации направлений весов сильнее, чем «выталкивает» с другой. Именно из-за такой асимметрии седловые точки — не локальные минимумы и не локальные максимумы, а места, где оптимизатор может надолго застревать, двигаясь очень медленно вдоль «хребта».

Почему это важно

Сводные формулы этого раздела — это финальный мост, превращающий поверхностный интеграл из красивой геометрической идеи в конкретную вычислительную процедуру: выбрал параметризацию, посчитал $\vec r_u\times\vec r_v$ (или его частный случай $(-g_x,-g_y,1)$ для явного графика), подставил в двойной интеграл — и дальше работают все приёмы, которые ты знаешь с курса кратных интегралов: полярные координаты, замена переменных, разбиение области на простые куски. Владение этой техникой окончательно закрывает вопрос «как вообще считать интегралы по искривлённым объектам» — а в следующем уроке ты увидишь, что для замкнутых поверхностей есть ещё более мощный трюк, позволяющий вообще избежать прямого вычисления потока, — теорему Остроградского-Гаусса.

Практика: 30 заданий

Базовые задания (1–10)

Задание 1: Найти площадь части плоскости $z=x+y$ над квадратом $D=[0,2]\times[0,2]$.

Задание 2: Вычислить $\iint_S 1\,dS$, где $S$ — часть плоскости $2x+2y+z=4$ над треугольником с вершинами $(0,0)$, $(1,0)$, $(0,1)$.

Задание 3: Найти массу горизонтального диска $x^2+y^2\le1$ (лежащего в плоскости $z=1$) с поверхностной плотностью $\sigma(x,y,z)=x^2+y^2$ (модель распределения точек лидара, плотность которых растёт к краю круга сканирования).

Задание 4: Найти площадь боковой поверхности цилиндра $x^2+y^2=9$, $0\le z\le4$.

Задание 5: Найти площадь части конуса $z=\sqrt{x^2+y^2}$ над кругом $x^2+y^2\le1$.

Задание 6: Вычислить $\iint_S z\,dS$, где $S$ — часть плоскости $z=x$ над квадратом $[0,1]\times[0,1]$.

Задание 7: Найти поток постоянного поля излучения $\vec F=(0,0,5)$ (модель равномерного вертикального потока фотонов) через горизонтальный плоский датчик — диск $x^2+y^2\le4$ в плоскости $z=0$, ориентация вверх.

Задание 8: Найти поток поля $\vec F=(1,0,0)$ через часть вертикальной плоскости $x=1$, $0\le y\le2$, $0\le z\le3$, ориентация — нормаль $(1,0,0)$.

Задание 9: Найти элемент площади $dS$ и ненормированный вектор нормали для поверхности $z=x^2+y^2$ в точке $(1,1,2)$.

Задание 10: Найти площадь верхней полусферы $x^2+y^2+z^2=4$, $z\ge0$.

Средние задания (11–20)

Задание 11: Найти площадь части параболоида $z=x^2+y^2$ над кругом $x^2+y^2\le1$.

Задание 12: Найти поток радиального поля $\vec F=(x,y,z)$ через сферу радиуса $R=5$ с внешней нормалью.

Задание 13: Вычислить $\iint_S xy\,dS$, где $S$ — часть плоскости $x+y+z=1$ в первом октанте.

Задание 14: Найти массу боковой поверхности конуса $z=\sqrt{x^2+y^2}$, $0\le z\le2$, с плотностью $\sigma(x,y,z)=z$.

Задание 15: Найти поток поля $\vec F=(0,0,z)$ через часть горизонтальной плоскости $z=3$ над квадратом $[0,2]\times[0,2]$, ориентация вверх.

Задание 16: Вычислить $\iint_S(x^2+y^2)\,dS$, где $S$ — боковая поверхность цилиндра $x^2+y^2=1$, $0\le z\le2$.

Задание 17: Найти поток поля $\vec F=(x,y,0)$ через часть плоскости $z=x+y$ над квадратом $D=[0,1]\times[0,1]$, ориентация вверх.

Задание 18: Найти площадь пояса сферы $x^2+y^2+z^2=9$ между углами $\varphi=\pi/6$ и $\varphi=\pi/3$.

Задание 19: Найти массу части плоскости $x+2y+2z=6$ в первом октанте с плотностью $\sigma(x,y,z)=z$.

Задание 20: Найти поток поля $\vec F=(x,y,2z)$ через верхнюю полусферу $x^2+y^2+z^2=1$, $z\ge0$, с внешней нормалью.

Продвинутые задания (21–30)

Задание 21 (машинное обучение): Датчик (лидар) с конусом обзора, половинный угол которого $45°=\pi/4$, установлен в центре сферы радиуса $R=10$ (модель охвата сканирования). Найти площадь покрываемой конусом «шапочки» сферы и долю от общей площади сферы.

Задание 22 (машинное обучение): На боковой поверхности цилиндра $x^2+y^2=4$, $0\le z\le1$ (цилиндрическая развёртка отсканированного объекта) задан признак $f(x,y,z)=x^2$. Найти суммарное («агрегированное») значение признака по поверхности.

Задание 23 (машинное обучение): Найти поток «закрученного» поля $\vec F=(y,-x,z)$ через часть параболоида $z=4-x^2-y^2$ над кругом $x^2+y^2\le4$, ориентация вверх.

Задание 24: Найти массу сферы $x^2+y^2+z^2=1$ с плотностью $\sigma=z^2$.

Задание 25 (машинное обучение): Найти поток «радиального затухающего» поля $\vec F=\dfrac{\vec r}{|\vec r|^3}$ (модель влияния точечного признака в пространстве эмбеддингов, убывающего как в законе Кулона) через сферу радиуса $R=2$ с центром в начале координат, внешняя нормаль.

Задание 26: Найти площадь части боковой поверхности цилиндра $y^2+z^2=4$ (ось вдоль $Ox$), $z\ge0$, $0\le x\le3$.

Задание 27 (машинное обучение): Найти поток постоянного вертикального поля $\vec F=(0,0,1)$ через график «ландшафта внимания» $z=\sin x\cos y$ над квадратом $D=[0,\pi]\times[0,\pi]$, ориентация вверх.

Задание 28: Найти площадь боковой поверхности конуса $z=2-\sqrt{x^2+y^2}$ над кругом $x^2+y^2\le9$.

Задание 29 (машинное обучение): Найти поток радиального поля $\vec F=(x,y,z)$ через боковую поверхность цилиндра $x^2+y^2=1$, $-1\le z\le1$ (без верхней и нижней крышек), внешняя нормаль (модель оценки суммарного «выхода» признака через боковую стенку цилиндрического рецептивного поля свёрточной сети).

Задание 30 (машинное обучение): Ландшафт функции потерь вблизи седловой точки моделируется поверхностью $z=w_1w_2$ (седло) над кругом $w_1^2+w_2^2\le1$. Найти площадь этого куска поверхности — геометрическую характеристику «остроты» седла в окрестности критической точки.

Частые ошибки

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

• Путают формулу $dS$ для явного графика с общей формулой через параметризацию. Формула $dS=\sqrt{1+g_x^2+g_y^2}\,dx\,dy$ работает только для явного графика $z=g(x,y)$. Для произвольной параметризации нужно честно вычислять $\vec r_u\times\vec r_v$ — попытка «подогнать» под знакомую формулу сферу или тор в общем случае приведёт к неверному ответу.

• Забывают, что интеграл первого рода не зависит от ориентации, а второго — зависит. Это частая путаница именно потому, что оба интеграла выглядят похоже. Если в задаче не сказано «поток» или «векторное поле», а стоит скалярная функция $f$, ориентация вообще не имеет значения; как только речь заходит о потоке — выбор нормали становится обязательной частью условия задачи.

• Не проверяют согласованность направления $\vec r_u\times\vec r_v$ с требуемой ориентацией. При вычислении потока через параметризацию легко забыть, что векторное произведение может «случайно» смотреть в противоположную сторону от требуемой нормали — тогда весь ответ получится с обратным знаком. Всегда проверяй знак хотя бы в одной характерной точке (например, для сферы — смотрит ли $\vec r_u\times\vec r_v$ наружу или внутрь при конкретных значениях углов).

• Путают элемент площади $dS$ сферы в разных параметризациях по $\theta$ и $\varphi$ с полярным элементом площади плоскости $r\,dr\,d\theta$. Это разные формулы: у сферы стоит $R^2\sin\varphi$, а не $r$, и путать их — источник частых ошибок при переходе между темами.

• Забывают про множитель площади при переходе от «плотности» к «массе». Иногда подставляют плотность $\sigma(x,y,z)$ прямо как ответ, забывая домножить на $dS$ и проинтегрировать по всей области — а ведь смысл интеграла именно в накоплении вклада по всей поверхности, а не в вычислении плотности в одной точке.

• Неправильно определяют пределы интегрирования по параметрам для «половинок» и «долек» поверхностей. Например, для верхней полусферы забывают ограничить $\varphi\in[0,\pi/2]$ вместо полного $[0,\pi]$, что задаёт всю сферу, а не только её верхнюю часть — всегда рисуй мысленную (или настоящую) картинку области параметров $D$ перед вычислением.

Главное запомнить

• Поверхностный интеграл первого рода $\iint_S f\,dS$ интегрирует скалярную функцию по площади поверхности и не зависит от выбора стороны (ориентации) поверхности.

• Поверхностный интеграл второго рода (поток) $\iint_S\vec F\cdot\vec n\,dS$ интегрирует нормальную составляющую векторного поля и обязательно требует выбора ориентации — при её смене знак интеграла меняется на противоположный.

• Элемент площади поверхности выводится через векторное произведение касательных векторов параметризации: $dS=|\vec r_u\times\vec r_v|\,du\,dv$.

• Для явного графика $z=g(x,y)$ формула упрощается до $dS=\sqrt{1+g_x^2+g_y^2}\,dx\,dy$, а ненормированный вектор нормали равен $(-g_x,-g_y,1)$.

• Ориентация поверхности — это непрерывный (без разрывов) выбор одной из двух возможных единичных нормалей в каждой точке; не всякая поверхность ориентируема (классический пример неориентируемой поверхности — лист Мёбиуса).

• Оба типа интеграла сводятся к обычному двойному интегралу по плоской области параметров $D$ — вся геометрическая сложность поверхности «упаковывается» в множитель $|\vec r_u\times\vec r_v|$ или в вектор $\vec r_u\times\vec r_v$.

• При $f\equiv1$ интеграл первого рода даёт площадь поверхности; при $f$, равной плотности, — массу оболочки.

• На поверхностях с высокой симметрией (сфера, цилиндр) скалярное произведение $\vec F\cdot\vec n$ часто оказывается постоянным, и интеграл сводится к произведению константы на площадь — мощный практический приём.

• Поток постоянного вертикального поля через любой график функции равен площади проекции области определения на плоскость — форма самой поверхности при этом не важна.

• Параметризация поверхности — универсальный инструмент, применимый к любой гладкой поверхности, вне зависимости от того, выражается ли она явной функцией $z=g(x,y)$ или нет.

Связь с другими темами курса

Поверхностные интегралы напрямую опираются на параметризацию кривых из прошлого урока (217) — сама логика «превратить геометрический объект в функцию параметров и свести к обычному интегралу» дословно повторяется, только на один параметр больше. Формула элемента площади через векторное произведение частных производных использует аппарат частных производных функций нескольких переменных и векторное произведение из курса линейной алгебры и аналитической геометрии. Вычисление самих двойных интегралов по области параметров опирается на технику кратных интегралов, включая переход к полярным координатам. А впереди, в следующем уроке о теории поля (219), поток через поверхность станет центральным героем теоремы Остроградского-Гаусса, связывающей его с дивергенцией поля внутри объёма, — и ты увидишь, как многие сегодняшние вычисления (например, поток радиального поля через сферу) можно будет получить буквально в одну строчку без всякой параметризации.

Интересные факты

• Лист Мёбиуса — не просто занятная игрушка: он служит стандартным контрпримером в топологии для демонстрации того, что понятие ориентации поверхности нетривиально, и используется как учебная модель в компьютерной графике для тестирования алгоритмов текстурирования, которые обязаны корректно обрабатывать неориентируемые поверхности.

• Формула Остроградского-Гаусса, о которой пойдёт речь в следующем уроке, была независимо переоткрыта несколько раз в XIX веке — помимо Остроградского и Гаусса, схожие результаты получали Джордж Грин и (в частном случае) Жозеф Луи Лагранж, что типично для фундаментальных теорем, «витавших в воздухе» одновременно у нескольких исследователей, работавших над сходными физическими задачами.

• В современной компьютерной графике уравнение визуализации изображения, сформулированное Джеймсом Каджией в 1986 году, — это, по сути, поверхностный интеграл по полусфере направлений в каждой точке видимой поверхности; численное приближение этого интеграла методом Монте-Карло (случайная выборка направлений вместо честного вычисления) лежит в основе трассировки лучей, используемой при фотореалистичном построении изображений и в нейросетевых системах синтеза картинок.

• Площадь поверхности сферы, вычисленная через поверхностный интеграл ($4\pi R^2$), была известна ещё Архимеду более двух тысяч лет назад — он доказал её геометрически, без всякого интегрального исчисления, сравнив площадь сферы с боковой поверхностью описанного вокруг неё цилиндра; поверхностный интеграл лишь дал этому результату строгий и универсальный вычислительный аппарат, применимый к произвольным, а не только идеально симметричным поверхностям.

Лайфхаки и полезные трюки

• Прежде чем параметризовать поверхность, проверь, не является ли она явным графиком $z=g(x,y)$ (или $y=g(x,z)$, или $x=g(y,z)$) — тогда формула для $dS$ и потока сильно упрощается, и векторное произведение можно не вычислять вручную.

• Для сфер и цилиндров держи под рукой готовые формулы $dS=R^2\sin\varphi\,d\varphi\,d\theta$ и $dS=R\,d\theta\,dz$ — они возникают настолько часто, что имеет смысл выучить их наизусть, а не выводить каждый раз заново.

• Если задание просит найти поток через сферу с центром в начале координат, всегда сначала попробуй посчитать $\vec F\cdot\vec n$ на самой поверхности — не исключено, что оно окажется постоянным (как в примерах с радиальными полями), и тогда весь интеграл схлопнется до произведения константы на площадь сферы.

• При вычислении потока через явный график всегда сначала распиши, какая ориентация требуется («вверх» или «вниз»), и только потом подставляй числа в формулу $-Pg_x-Qg_y+R$ — если требуется ориентация «вниз», результат нужно взять с обратным знаком.

• Если поле имеет «вращательную» компоненту вида $(y,-x,0)$ или похожую антисимметричную структуру, проверь сначала, не сокращается ли она в выражении $-Pg_x-Qg_y$ — часто именно так и происходит, и это сильно упрощает вычисления (как в задании 23).

• Для интегралов по симметричным областям (кругам, прямоугольникам с центром в начале координат) заранее прикинь, не занулится ли ответ из-за нечётности подынтегральной функции — это сэкономит время и убережёт от лишних вычислений.

• Всегда делай финальную проверку размерности и порядка величины ответа: площадь части плоскости не может быть меньше площади её проекции на плоскость $Oxy$ (поскольку $\sqrt{1+g_x^2+g_y^2}\ge1$), а масса не может получиться отрицательной, если плотность везде положительна, — такие быстрые проверки ловят большинство арифметических ошибок.

Поверхностные интегралы — это тема, где геометрическая интуиция, аккуратная алгебра частных производных и техника кратных интегралов наконец-то сходятся вместе в одном мощном инструменте. Освоив её, ты готов к следующему, ещё более элегантному шагу — теореме Остроградского-Гаусса, которая свяжет всё, что ты сегодня вычислял вручную через параметризацию, с куда более простым и глубоким взглядом на поведение векторного поля внутри объёма. Не останавливайся — самое красивое ещё впереди.