Формула Грина 🔄

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

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

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

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

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

Джордж Грин — один из самых необычных персонажей в истории математики XIX века. Он родился в 1793 году в английском городе Ноттингем в семье пекаря, получил всего около года формального школьного образования и с юности работал на мельнице своего отца — сначала помогая молоть зерно, а позже управляя ею самостоятельно. Никакого университета, никакого научного окружения, никаких учителей математики — и тем не менее именно этот человек, самоучка по книгам, взятым в местной библиотеке, к тридцати годам самостоятельно освоил математический анализ на уровне, недоступном многим профессиональным учёным того времени.

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

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

Точная формулировка и условия применимости

Интуиция: граница знает всё о внутренности

Давай разберёмся, в чём именно состоит фокус формулы Грина, на интуитивном уровне, прежде чем переходить к строгой формулировке. Представь плоскую область $D$, ограниченную замкнутой кривой $\partial D$ — например, озеро причудливой формы на карте, а его берег и есть контур $\partial D$. Пусть в каждой точке плоскости задано векторное поле $\vec F=(P,Q)$ — скажем, поле скоростей течения воды. Циркуляция поля по берегу озера (криволинейный интеграл $\oint_{\partial D}P\,dx+Q\,dy$) — это суммарный эффект того, насколько вода «закручивается» вдоль всего периметра.

Ключевая идея Грина: этот эффект на границе можно полностью восстановить, если знать, насколько поле «закручивается» в каждой отдельной точке внутри озера, а затем просуммировать (проинтегрировать) эти локальные завихрения по всей площади. Величина $\dfrac{\partial Q}{\partial x}-\dfrac{\partial P}{\partial y}$ как раз и измеряет локальное завихрение поля в конкретной точке — это плоский аналог ротора, с которым ты уже познакомился в прошлом уроке про теорию поля. Формула Грина утверждает: суммарная циркуляция по границе равна сумме всех локальных завихрений внутри области. Это очень похоже на то, как суммарный ток через замкнутый контур в электротехнике равен сумме токов, «рождённых» внутри контура каждым отдельным источником, — общий мотив «граница = сумма внутренних вкладов» встречается в физике постоянно.

Определение

Теорема (формула Грина). Пусть $D$ — ограниченная односвязная область на плоскости (то есть область без «дырок»), а её граница $\partial D$ — кусочно-гладкая замкнутая кривая, положительно ориентированная (обход происходит против часовой стрелки, так что область $D$ всё время остаётся слева от направления движения). Пусть функции $P(x,y)$ и $Q(x,y)$ непрерывно дифференцируемы в некоторой открытой области, содержащей $D$ вместе с её границей. Тогда справедливо равенство:

$$\oint_{\partial D} P\,dx+Q\,dy=\iint_D\left(\frac{\partial Q}{\partial x}-\frac{\partial P}{\partial y}\right)dx\,dy$$

Условия применимости здесь не формальность, а строгая необходимость: контур должен быть обойдён именно в положительном направлении (иначе обе части формулы отличаются знаком), область должна быть односвязной (иначе нужно отдельно учитывать ориентацию внутренних «дырок» — мы разберём это ниже), а функции $P$ и $Q$ обязаны быть непрерывно дифференцируемыми во всей замкнутой области целиком, без единого исключения — даже одна особая точка внутри контура, где производные не определены, способна полностью сломать равенство.

Обрати внимание на прямую параллель с прошлым уроком: выражение $\dfrac{\partial Q}{\partial x}-\dfrac{\partial P}{\partial y}$ — это ровно тот же критерий, который определяет потенциальность поля. Если это выражение тождественно равно нулю всюду внутри области, правая часть формулы Грина автоматически обращается в ноль, а значит и циркуляция по любому замкнутому контуру внутри такой области равна нулю. Формула Грина, таким образом, не просто вычислительный трюк — она даёт строгое объяснение того, почему потенциальные поля не совершают суммарной работы по замкнутым траекториям.

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

Пример 1 (лёгкий). Проверить формулу Грина для поля $\vec F=(P,Q)=(xy,\,x^2)$ по границе единичного квадрата $D=[0,1]\times[0,1]$, обходимой против часовой стрелки.

Сначала правая часть. Находим $\dfrac{\partial Q}{\partial x}=2x$, $\dfrac{\partial P}{\partial y}=x$, разность равна $x$:

Теперь левая часть — напрямую по четырём сторонам квадрата. На стороне $C_1$ ($y=0$, $x:0\to1$, $dy=0$): вклад $\int_0^1 x\cdot0\,dx=0$. На стороне $C_2$ ($x=1$, $y:0\to1$, $dx=0$): вклад $\int_0^1 1^2\,dy=1$. На стороне $C_3$ ($y=1$, $x:1\to0$, $dy=0$): вклад $\int_1^0 x\cdot1\,dx=-\dfrac12$. На стороне $C_4$ ($x=0$, $y:1\to0$, $dx=0$): вклад $\int_1^0 0^2\,dy=0$. Сумма: $0+1-\dfrac12+0=\dfrac12$.

Ответ: обе части равны $\dfrac12$ — формула подтверждена.

Пример 2 (средний, нарушение условий применимости). Рассмотрим вихревое поле с особенностью в начале координат $\vec F=(P,Q)=\left(-\dfrac{y}{x^2+y^2},\ \dfrac{x}{x^2+y^2}\right)$ и окружность радиуса $R$ с центром в начале координат. Посмотрим, что произойдёт, если применить формулу Грина, не проверив условия.

Прямое вычисление циркуляции: параметризуем $x=R\cos t$, $y=R\sin t$, $t\in[0,2\pi]$, тогда $dx=-R\sin t\,dt$, $dy=R\cos t\,dt$, $P=-\dfrac{\sin t}{R}$, $Q=\dfrac{\cos t}{R}$. Подставляем:

Заметь: результат $2\pi$ не зависит от радиуса $R$ вообще. Теперь попробуем формально применить формулу Грина, не глядя на условия: прямым дифференцированием можно проверить, что $\dfrac{\partial Q}{\partial x}=\dfrac{\partial P}{\partial y}$ всюду, где поле определено (то есть везде, кроме точки $(0,0)$), а значит формально правая часть должна была бы дать ноль. Но прямое вычисление даёт $2\pi\neq0$ — противоречие.

Ответ: формула Грина здесь неприменима в исходном виде, потому что поле не определено (и не дифференцируемо) в точке $(0,0)$, которая лежит внутри контура, — нарушено ключевое условие непрерывной дифференцируемости во всей замкнутой области. Вся «недостача» в $2\pi$ сосредоточена ровно в этой одной особой точке — мы вернёмся к этому наблюдению в разделе про упрощение сложных контуров, потому что именно оно и лежит в основе одного из самых мощных приёмов применения формулы.

Пример 3 (сложный, машинное обучение). Пусть функция потерь модели с двумя весами задаётся гармоническим многочленом $L(w_1,w_2)=w_1^3-3w_1w_2^2$. Антиградиентное поле оптимизатора равно $\vec F=-\nabla L=(-3w_1^2+3w_2^2,\ 6w_1w_2)$. Проверить, что циркуляция этого поля по границе единичного квадрата в пространстве весов равна нулю, используя формулу Грина, — и объяснить, почему этот результат верен для любой дважды непрерывно дифференцируемой функции потерь.

Находим $P=-3w_1^2+3w_2^2$, $Q=6w_1w_2$. Тогда $\dfrac{\partial Q}{\partial w_1}=6w_2$ и $\dfrac{\partial P}{\partial w_2}=6w_2$ — эти выражения совпадают тождественно, значит разность равна нулю в каждой точке плоскости, а не только внутри квадрата:

Ответ: циркуляция антиградиентного поля по границе квадрата равна нулю. Причина этого результата глубже, чем конкретный вид функции $L$: для антиградиента $\vec F=-\nabla L=(-\partial L/\partial w_1,\,-\partial L/\partial w_2)$ выражение $\dfrac{\partial Q}{\partial w_1}-\dfrac{\partial P}{\partial w_2}$ равно $-\dfrac{\partial^2L}{\partial w_1\partial w_2}+\dfrac{\partial^2L}{\partial w_2\partial w_1}$, а эти две смешанные производные обязаны совпадать для любой дважды непрерывно дифференцируемой функции — по теореме Шварца, знакомой тебе из курса частных производных. Значит формула Грина автоматически доказывает общий факт: циркуляция чисто градиентного (антиградиентного) поля по любому замкнутому контуру в пространстве весов всегда равна нулю, каким бы причудливым ни был сам контур и какой бы конкретной ни была функция потерь.

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

Точная формулировка формулы Грина — это не педантичная формальность, а страховка от ошибок, которые на практике встречаются постоянно. Пропустишь проверку особой точки внутри контура (как в примере 2) — получишь неверный ноль вместо правильного $2\pi$. Перепутаешь направление обхода — получишь результат с противоположным знаком. В машинном обучении критерий $\partial Q/\partial x=\partial P/\partial y$, который формула Грина связывает с нулевой циркуляцией, — это ровно тот тест, которым можно быстро проверить, ведёт ли себя эффективное «поле обновлений» оптимизатора как честный градиентный спуск или нет, не вычисляя интегралы по каждому конкретному пути вручную.

Идея доказательства: от простой области до всей мозаики

Интуиция: мозаика с сокращающимися рёбрами

Давай разберёмся, откуда вообще берётся формула Грина — не зазубрим её, а увидим, почему она обязана быть верной. Представь плоскую область, выложенную, как мозаика, из множества маленьких кусочков — например, из маленьких квадратиков или треугольников, плотно прилегающих друг к другу. Пусть для каждого отдельного кусочка мозаики формула Грина уже доказана (мы разберём это чуть ниже для одного простого кусочка). Что произойдёт, если сложить циркуляции по границам всех кусочков вместе?

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

Определение

Утверждение (доказательство для элементарной области). Пусть область $D$ — «элементарная», то есть одновременно допускает описание $\varphi_1(x)\le y\le\varphi_2(x)$ для $x\in[a,b]$ (ограничена снизу и сверху графиками функций от $x$) и аналогичное описание с ролями $x$ и $y$ поменянными местами. Тогда для такой области справедливо:

$$\oint_{\partial D}P\,dx=-\iint_D\frac{\partial P}{\partial y}\,dx\,dy,\qquad \oint_{\partial D}Q\,dy=\iint_D\frac{\partial Q}{\partial x}\,dx\,dy$$

Докажем первое равенство. По формуле Ньютона-Лейбница, применённой к внутреннему интегрированию по $y$:

$$\iint_D\frac{\partial P}{\partial y}\,dy\,dx=\int_a^b\Bigl[P(x,\varphi_2(x))-P(x,\varphi_1(x))\Bigr]dx$$

С другой стороны, при обходе границы против часовой стрелки нижняя кривая $y=\varphi_1(x)$ проходится слева направо ($x$ от $a$ до $b$), а верхняя $y=\varphi_2(x)$ — справа налево ($x$ от $b$ до $a$), поэтому

$$\oint_{\partial D}P\,dx=\int_a^b P(x,\varphi_1(x))\,dx+\int_b^a P(x,\varphi_2(x))\,dx=\int_a^b P(x,\varphi_1(x))\,dx-\int_a^b P(x,\varphi_2(x))\,dx$$

Это выражение — в точности минус то, что получилось для двойного интеграла, что и требовалось доказать. Второе равенство доказывается совершенно аналогично, только через описание области по $x$-правильности. Складывая оба равенства, получаем полную формулу Грина для элементарной области. Для произвольной области нужное разбиение на элементарные кусочки и последующее сокращение внутренних рёбер, описанное выше, завершают доказательство в общем случае.

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

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

Пример 1 (лёгкий). Проверить равенство $\oint_{\partial D}P\,dx=-\iint_D\dfrac{\partial P}{\partial y}\,dx\,dy$ для $P=y^2$ на треугольнике $D$ с вершинами $(0,0)$, $(1,0)$, $(0,1)$.

Правая часть: $\dfrac{\partial P}{\partial y}=2y$, и

значит правая часть равна $-\dfrac13$.

Теперь левая часть по трём сторонам треугольника (обход против часовой стрелки: $(0,0)\to(1,0)\to(0,1)\to(0,0)$). На стороне $(0,0)\to(1,0)$: $y=0$, значит $P=y^2=0$, вклад $0$. На стороне $(1,0)\to(0,1)$: параметризуем $x=1-t$, $y=t$, $t\in[0,1]$, $dx=-dt$, $P=t^2$, вклад $\int_0^1t^2\cdot(-1)\,dt=-\dfrac13$. На стороне $(0,1)\to(0,0)$: $x=0$ постоянно, $dx=0$, вклад $0$. Сумма: $0-\dfrac13+0=-\dfrac13$.

Ответ: обе части равны $-\dfrac13$ — равенство подтверждено.

Пример 2 (средний). Проверить равенство $\oint_{\partial D}Q\,dy=\iint_D\dfrac{\partial Q}{\partial x}\,dx\,dy$ для $Q=x^2$ на области $D$ между кривыми $y=x^2$ и $y=x$ при $x\in[0,1]$.

Правая часть: $\dfrac{\partial Q}{\partial x}=2x$, область описывается как $x^2\le y\le x$ для $x\in[0,1]$:

Левая часть: граница состоит из нижней кривой $y=x^2$ (проходится слева направо, $x:0\to1$) и верхней кривой $y=x$ (проходится справа налево, $x:1\to0$) — именно такой обход даёт положительную ориентацию, поскольку область лежит выше нижней кривой и ниже верхней. На нижней кривой: $x=t$, $y=t^2$, $dy=2t\,dt$, $Q=t^2$, вклад $\int_0^1 t^2\cdot2t\,dt=\dfrac12$. На верхней кривой: $x=t$, $y=t$, $t:1\to0$, $dy=dt$, $Q=t^2$, вклад $\int_1^0t^2\,dt=-\dfrac13$. Сумма: $\dfrac12-\dfrac13=\dfrac16$.

Ответ: обе части равны $\dfrac16$ — равенство подтверждено.

Пример 3 (сложный, топология данных). Проиллюстрируем идею «сокращения внутренних рёбер» на составной области — квадрате $[0,1]\times[0,1]$, разбитом диагональю $y=x$ на два треугольника $T_1$ (вершины $(0,0)$, $(1,0)$, $(1,1)$) и $T_2$ (вершины $(0,0)$, $(1,1)$, $(0,1)$). Возьмём $P=0$, $Q=x$ и вычислим $\oint Q\,dy$ по каждому треугольнику отдельно, а затем сложим результаты.

Для $T_1$ (обход $(0,0)\to(1,0)\to(1,1)\to(0,0)$): сторона $(0,0)\to(1,0)$ даёт вклад $0$ ($y=0$ постоянно, $dy=0$); сторона $(1,0)\to(1,1)$ даёт $\int_0^1 1\,dy=1$; диагональная сторона $(1,1)\to(0,0)$ (параметризация $x=t=y$, $t:1\to0$) даёт $\int_1^0 t\,dt=-\dfrac12$. Сумма по $T_1$: $0+1-\dfrac12=\dfrac12$.

Для $T_2$ (обход $(0,0)\to(1,1)\to(0,1)\to(0,0)$): диагональная сторона $(0,0)\to(1,1)$ ($x=t=y$, $t:0\to1$) даёт $\int_0^1 t\,dt=\dfrac12$; сторона $(1,1)\to(0,1)$ даёт $0$ ($y=1$ постоянно); сторона $(0,1)\to(0,0)$ даёт $0$ ($Q=x=0$ на ней). Сумма по $T_2$: $\dfrac12+0+0=\dfrac12$.

Заметь: диагональ входит в сумму дважды — один раз с вкладом $-\dfrac12$ (со стороны $T_1$, обход от $(1,1)$ к $(0,0)$) и один раз с вкладом $+\dfrac12$ (со стороны $T_2$, обход от $(0,0)$ к $(1,1)$) — эти два вклада полностью взаимно уничтожаются, как и предсказывает идея доказательства.

Ответ: $T_1$ даёт $\dfrac12$, $T_2$ даёт $\dfrac12$, сумма равна $1$. Проверим напрямую по границе всего квадрата: только правая сторона ($x=1$, $y:0\to1$) даёт ненулевой вклад $\int_0^1 1\,dy=1$, все остальные стороны дают $0$ — итого тоже $1$, полное совпадение. Эта же самая идея — разбить сложный объект на простые кусочки (симплексы), а затем сложить локальные вклады с автоматическим сокращением общих внутренних граней — лежит в основе конструкции, которую топологи анализа данных называют цепным комплексом: облако точек триангулируется на треугольники, тетраэдры и симплексы более высоких размерностей, и оператор границы, отслеживающий, как эти симплексы стыкуются друг с другом, оказывается дискретным родственником ровно того же самого механизма сокращения, который мы только что наблюдали на конкретных числах.

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

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

Площадь плоской фигуры как частный случай формулы Грина

Интуиция: измерить площадь, ни разу не заглянув внутрь

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

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

Определение

Теорема (площадь через формулу Грина). Если в формуле Грина выбрать $P$ и $Q$ так, чтобы $\dfrac{\partial Q}{\partial x}-\dfrac{\partial P}{\partial y}\equiv1$, правая часть формулы становится равна площади области $D$. Есть несколько удобных стандартных вариантов такого выбора, дающих один и тот же результат:

$$S=\oint_{\partial D}x\,dy=-\oint_{\partial D}y\,dx=\frac12\oint_{\partial D}\bigl(x\,dy-y\,dx\bigr)$$

Третий вариант (симметричная форма) обычно предпочтительнее на практике, потому что он менее чувствителен к выбору начала координат и часто даёт более простые вычисления после параметризации контура.

Это очень мощный результат: площадь области — величина сугубо «внутренняя» по своему смыслу — полностью восстанавливается по данным, собранным исключительно на границе. Для параметрического контура $x=x(t)$, $y=y(t)$, $t\in[\alpha,\beta]$, формула площади превращается в конкретный определённый интеграл:

А для кривой, заданной в полярных координатах уравнением $r=r(\theta)$, при подстановке $x=r\cos\theta$, $y=r\sin\theta$ эта же формула сворачивается в хорошо знакомую тебе из курса определённого интеграла формулу площади криволинейного сектора $S=\dfrac12\int r^2(\theta)\,d\theta$ — то есть привычная формула площади в полярных координатах на самом деле является частным случаем формулы Грина.

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

Пример 1 (лёгкий). Найти площадь эллипса $\dfrac{x^2}{a^2}+\dfrac{y^2}{b^2}=1$ через формулу Грина.

Параметризуем эллипс: $x=a\cos t$, $y=b\sin t$, $t\in[0,2\pi]$. Тогда $x'(t)=-a\sin t$, $y'(t)=b\cos t$.

Ответ: $S=\pi ab$ — узнаваемая школьная формула площади эллипса, теперь строго выведенная.

Пример 2 (средний). Найти площадь фигуры, ограниченной астроидой $x=a\cos^3t$, $y=a\sin^3t$, $t\in[0,2\pi]$.

Находим производные: $x'(t)=-3a\cos^2t\sin t$, $y'(t)=3a\sin^2t\cos t$. Подставляем в $x\,dy-y\,dx$:

Используем $\sin^2t\cos^2t=\dfrac14\sin^2(2t)=\dfrac18\bigl(1-\cos4t\bigr)$:

Ответ: $S=\dfrac{3\pi a^2}8$ — классический результат для площади астроиды.

Пример 3 (сложный, машинное обучение). В компьютерном зрении и анализе данных объект, найденный алгоритмом сегментации, часто представлен не аналитической формулой, а просто списком координат вершин многоугольного контура. Найти площадь такого сегментированного объекта с вершинами (в порядке обхода против часовой стрелки) $(0,0)$, $(3,0)$, $(4,2)$, $(2,4)$, $(0,3)$, используя дискретный частный случай формулы Грина — так называемую формулу площади многоугольника (формулу шнурков).

Для многоугольника с вершинами $(x_i,y_i)$ дискретный аналог интеграла $\dfrac12\oint(x\,dy-y\,dx)$ превращается в сумму по рёбрам:

Вычисляем каждое слагаемое по очереди. От $(0,0)$ до $(3,0)$: $0\cdot0-3\cdot0=0$. От $(3,0)$ до $(4,2)$: $3\cdot2-4\cdot0=6$. От $(4,2)$ до $(2,4)$: $4\cdot4-2\cdot2=12$. От $(2,4)$ до $(0,3)$: $2\cdot3-0\cdot4=6$. От $(0,3)$ обратно до $(0,0)$: $0\cdot0-0\cdot3=0$. Сумма: $0+6+12+6+0=24$.

Ответ: $S=12$. Именно эта формула (в компьютерной графике и обработке изображений её называют формулой шнурков или формулой Гаусса) стоит за функциями вычисления площади сегментированных объектов практически во всех библиотеках компьютерного зрения — например, за вычислением площади под ROC-кривой (метрика AUC), которая на практике часто приближается ломаной по нескольким контрольным точкам и вычисляется как площадь многоугольника между этой ломаной и осью абсцисс — тем же самым приёмом, который мы только что применили к сегментированному объекту.

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

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

Упрощение сложных контуров: переход к простой области

Интуиция: замена сложного пути простым

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

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

Определение

Следствие (случай многосвязной области). Пусть область $D$ ограничена внешним контуром $C_1$ (обход против часовой стрелки) и содержит внутри исключённую «дырку», ограниченную внутренним контуром $C_2$ (для сохранения положительной ориентации всей кольцевой области $C_2$ обходится по часовой стрелке — так, чтобы сама область $D$ всё время оставалась слева от направления движения по любой части границы). Если $P$ и $Q$ непрерывно дифференцируемы всюду в $D$ (то есть все особенности поля вынесены внутрь исключённой «дырки»), формула Грина принимает вид:

$$\oint_{C_1}P\,dx+Q\,dy+\oint_{C_2}P\,dx+Q\,dy=\iint_D\left(\frac{\partial Q}{\partial x}-\frac{\partial P}{\partial y}\right)dx\,dy$$

В частности, если правая часть равна нулю (поле потенциально всюду в кольцевой области), циркуляция по сложному внешнему контуру $C_1$ автоматически равна циркуляции по простому внутреннему контуру, обойдённому в том же (против часовой стрелки) направлении — вне зависимости от того, насколько причудливую форму имеет $C_1$.

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

Пример 1 (лёгкий). Вычислить $\oint_L P\,dx+Q\,dy$ для поля $P=0$, $Q=xy$ по границе верхнего полукруга радиуса $2$ (сложный контур — дуга плюс диаметр), перейдя к двойному интегралу по полукругу вместо прямого обхода дуги.

Находим $\dfrac{\partial Q}{\partial x}-\dfrac{\partial P}{\partial y}=y-0=y$. Переходим к полярным координатам: область — верхний полукруг радиуса $2$, то есть $\theta\in[0,\pi]$, $r\in[0,2]$, а $y=r\sin\theta$:

Ответ: $\dfrac{16}3\approx5{,}333$. Для проверки посчитаем и напрямую: на диаметре ($y=0$) вклад нулевой, поскольку $P=0$ и $dy=0$; на дуге ($x=2\cos\theta$, $y=2\sin\theta$, $\theta:0\to\pi$, $dy=2\cos\theta\,d\theta$, $Q=xy=4\sin\theta\cos\theta$) вклад равен $\int_0^\pi4\sin\theta\cos\theta\cdot2\cos\theta\,d\theta=8\int_0^\pi\sin\theta\cos^2\theta\,d\theta=8\cdot\dfrac23=\dfrac{16}3$ — оба способа дают одинаковый ответ, но переход к двойному интегралу оказался заметно короче прямого обхода границы.

Пример 2 (средний). Найти циркуляцию вихревого поля $\vec F=\left(-\dfrac{y}{x^2+y^2},\ \dfrac{x}{x^2+y^2}\right)$ по границе квадрата $[-2,2]\times[-2,2]$, не вычисляя интеграл напрямую по четырём сторонам квадрата.

Квадрат содержит особую точку $(0,0)$, поэтому формулу Грина нельзя применить прямо к квадрату. Вырежем из квадрата маленький круг радиуса $\varepsilon$ вокруг начала координат — в получившейся кольцевой области поле всюду непрерывно дифференцируемо, а разность $\dfrac{\partial Q}{\partial x}-\dfrac{\partial P}{\partial y}$ тождественно равна нулю (это стандартное свойство данного поля, проверяемое прямым дифференцированием). Значит правая часть формулы Грина для кольца равна нулю:

А циркуляция по окружности любого радиуса вокруг начала координат (против часовой стрелки) была вычислена в примере 2 первого раздела и равна $2\pi$ вне зависимости от радиуса.

Ответ: циркуляция по границе квадрата равна $2\pi$ — вместо утомительного вычисления по четырём сторонам квадрата (с обратными тангенсами и логарифмами) мы свели задачу к уже известному результату для простой окружности.

Пример 3 (сложный, машинное обучение). Функция потерь модели с седловой точкой задана как $L(w_1,w_2)=w_1^2-w_2^2$. Антиградиентное поле $\vec F=-\nabla L=(-2w_1,\ 2w_2)$. Найти циркуляцию этого поля вдоль сложного «звездообразного» контура в пространстве весов, заданного в полярных координатах уравнением $r(\theta)=1+0{,}3\cos5\theta$, $\theta\in[0,2\pi]$, — не вычисляя интеграл по этому контуру напрямую.

Находим $\dfrac{\partial Q}{\partial w_1}=0$ (здесь $Q=2w_2$ не зависит от $w_1$) и $\dfrac{\partial P}{\partial w_2}=0$ (здесь $P=-2w_1$ не зависит от $w_2$). Разность тождественно равна нулю на всей плоскости — седловая точка $(0,0)$ не создаёт никакой особенности поля, оно остаётся гладким везде.

Ответ: циркуляция равна $0$, причём этот ответ получен без единого вычисления интеграла по самому звездообразному контуру — прямая параметризация такой кривой потребовала бы интегрирования тригонометрических выражений высокой степени с аргументом $5\theta$, что заняло бы страницу вычислений. Содержательный вывод для машинного обучения: если эффективное поле обновлений весов остаётся чистым (анти)градиентом даже вблизи седловой точки, суммарная циркуляция вдоль любой, даже самой замысловатой траектории в пространстве весов равна нулю — а значит никакая петля оптимизатора вокруг седла сама по себе не порождает «даровой» дополнительной работы.

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

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

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

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

Задание 1: Вычислить $\oint_{\partial D}x^2\,dx+y^2\,dy$ по границе единичного квадрата $[0,1]\times[0,1]$ через формулу Грина.

Задание 2: Вычислить $\oint_{\partial D}2xy\,dx+x^2\,dy$ по границе треугольника с вершинами $(0,0)$, $(1,0)$, $(1,1)$.

Задание 3: Найти площадь круга радиуса $R=3$ через формулу Грина $S=\oint x\,dy$.

Задание 4: Вычислить $\oint_{\partial D}y\,dx-x\,dy$ по границе окружности единичного радиуса с центром в начале координат через формулу Грина.

Задание 5: Вычислить $\oint_{\partial D}3x\,dy$ по границе прямоугольника $[0,2]\times[0,3]$ через формулу Грина.

Задание 6: Найти площадь прямоугольника с вершинами $(0,0)$, $(4,0)$, $(4,2)$, $(0,2)$ через формулу $S=\dfrac12\oint(x\,dy-y\,dx)$.

Задание 7: Вычислить $\oint_{\partial D}e^x\,dx+x\,dy$ по границе единичного круга через формулу Грина.

Задание 8: Проверить потенциальность поля $P=2xy^3$, $Q=3x^2y^2$ и найти его циркуляцию по границе треугольника $(0,0)$, $(2,0)$, $(0,2)$.

Задание 9: Вычислить $\oint x\,dy$ по границе верхнего полукруга радиуса $2$ (дуга и диаметр) и сверить результат с площадью полукруга.

Задание 10: Найти площадь области, ограниченной параболой $y=x^2$ и прямой $y=1$ (от $x=-1$ до $x=1$), через формулу Грина, и сравнить с обычным определённым интегралом.

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

Задание 11: Вычислить $\oint_{\partial D}x^2y\,dx+xy^2\,dy$ по границе единичного круга через переход к полярным координатам в двойном интеграле.

Задание 12: Найти работу поля $\vec F=(y^2,x^2)$ вдоль границы треугольника $(0,0)$, $(1,0)$, $(1,1)$ через формулу Грина.

Задание 13: Проверить потенциальность поля $\vec F=(\cos y,\,-x\sin y)$ и найти работу вдоль любого пути от $(0,0)$ до $(\pi,0)$.

Задание 14: Найти площадь эллипса $\dfrac{x^2}9+\dfrac{y^2}4=1$ через формулу Грина.

Задание 15: Найти площадь фигуры, ограниченной одной аркой циклоиды $x=a(t-\sin t)$, $y=a(1-\cos t)$, $t\in[0,2\pi]$, и отрезком оси $x$.

Задание 16: Вычислить циркуляцию поля $\vec F=(x^2-y^2,\,2xy)$ вдоль окружности радиуса $R$ с центром в начале координат.

Задание 17: Найти площадь фигуры, ограниченной кардиоидой $r=1+\cos\theta$, через формулу площади в полярных координатах (частный случай формулы Грина).

Задание 18: Проверить, зависит ли работа поля $\vec F=(3x^2y,\,x^3+y)$ от пути, вычислив его циркуляцию по границе единичного квадрата через формулу Грина.

Задание 19 (машинное обучение): Функция потерь $L(w_1,w_2)=w_1^4+w_2^2$. Проверить потенциальность антиградиентного поля $\vec F=-\nabla L$ через критерий формулы Грина и вычислить его циркуляцию по границе прямоугольника $[-1,1]\times[-1,1]$ в пространстве весов.

Задание 20: Найти площадь фигуры, ограниченной астроидой $x=2\cos^3t$, $y=2\sin^3t$ (то есть при $a=2$), используя результат из теоретического раздела $S=\dfrac{3\pi a^2}8$.

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

Задание 21: Найти площадь одной петли розы $r=\sin3\theta$ (при $\theta\in\left[0,\dfrac\pi3\right]$) через формулу площади в полярных координатах.

Задание 22: Для вихревого поля $\vec F=\left(-\dfrac{y}{x^2+y^2},\dfrac{x}{x^2+y^2}\right)$ вычислить циркуляцию по границе квадрата $[1,3]\times[1,3]$, не содержащего особую точку $(0,0)$.

Задание 23: Для того же вихревого поля вычислить циркуляцию по границе кольца между окружностями радиуса $1$ и радиуса $3$ с центром в начале координат (внешняя — против часовой стрелки, внутренняя — по часовой, как того требует корректная ориентация кольцевой области).

Задание 24: Вычислить работу поля $\vec F=(y^3,\,3xy^2)$ от точки $(0,0)$ до точки $(2,4)$ вдоль дуги параболы $y=x^2$, используя потенциальность поля, и сверить с прямым вычислением по параболе.

Задание 25 (машинное обучение): Для седловой функции потерь $L(w_1,w_2)=w_1^2-w_2^2$ и антиградиентного поля $\vec F=-\nabla L$ определить циркуляцию вдоль сложного звездообразного контура вокруг седловой точки, не вычисляя интеграл напрямую.

Задание 26: Найти площадь фигуры, ограниченной параболой $y=4-x^2$ ($x\in[-2,2]$) и отрезком оси $x$, через формулу Грина, и сравнить с обычным определённым интегралом.

Задание 27: Найти площадь фигуры, ограниченной аркой синусоиды $y=\sin x$ ($x\in[0,\pi]$) и отрезком оси $x$, через формулу Грина.

Задание 28 (машинное обучение): Алгоритм сегментации изображения вернул замкнутый контур объекта с вершинами (в порядке обхода против часовой стрелки) $(0,0)$, $(5,0)$, $(6,3)$, $(3,5)$, $(0,4)$. Найти площадь этого объекта по формуле площади многоугольника.

Задание 29: Найти площадь фигуры, ограниченной лемнискатой Бернулли $r^2=2a^2\cos2\theta$ (обе петли), через формулу площади в полярных координатах.

Задание 30 (машинное обучение): Эффективное поле обновления весов при оптимизации с моментом содержит вихревую добавку $\vec F=(w_2,\,-w_1)$. Проверить потенциальность и найти циркуляцию этого поля по эллиптическому контуру $\dfrac{w_1^2}4+w_2^2=1$ (полуоси $a=2$, $b=1$) в пространстве весов.

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

Формула Грина обманчиво проста по внешнему виду, и именно поэтому с ней особенно легко совершить содержательную, а не просто арифметическую ошибку.

Неправильная ориентация контура. Это самая распространённая и самая коварная ошибка. Формула Грина требует обхода контура строго против часовой стрелки (область должна оставаться слева от направления движения). Если обойти контур в противоположном направлении, весь ответ поменяет знак — и, что хуже всего, такая ошибка не вызывает никакого явного противоречия в вычислениях, ты просто получишь неверный, но правдоподобно выглядящий результат. Прежде чем подставлять пределы интегрирования, всегда мысленно проверь: если встать в точку контура лицом по направлению движения, область должна оказаться слева от тебя.

Игнорирование особых точек внутри контура. Формула Грина требует, чтобы $P$ и $Q$ были непрерывно дифференцируемы во всей замкнутой области, а не «почти везде». Одна-единственная точка, где поле не определено или не дифференцируемо (как в примерах с вихревым полем $\vec F=(-y/(x^2+y^2),x/(x^2+y^2))$), способна дать разницу между «правильными» нулём и настоящим ненулевым ответом. Всегда проверяй область определения $P$ и $Q$ на предмет знаменателей, обращающихся в ноль, или других источников недифференцируемости внутри контура, прежде чем применять формулу напрямую.

Забытая неодносвязность области. Если область содержит «дырку» (например, кольцо), формулу Грина нельзя применить к одному только внешнему контуру — обязательно нужно учитывать внутреннюю границу, причём обойдённую в противоположном направлении относительно внешней. Пропуск этого шага — частая причина неверных ответов в задачах с кольцевыми областями.

Путаница знака в частном выражении ротора. Формула Грина использует именно $\dfrac{\partial Q}{\partial x}-\dfrac{\partial P}{\partial y}$, а не $\dfrac{\partial P}{\partial y}-\dfrac{\partial Q}{\partial x}$. Перепутать порядок вычитания — значит получить ответ с противоположным знаком, причём внешне вычисления будут выглядеть совершенно корректными.

Неверная параметризация контура при вычислении площади. При использовании формулы $S=\dfrac12\oint(x\,dy-y\,dx)$ легко забыть, что контур должен быть обойдён именно в положительном направлении. Если параметр естественным образом возрастает в направлении по часовой стрелке (как иногда бывает при стандартной параметризации некоторых кривых, включая циклоиду, что мы явно увидели в одном из примеров), результат интегрирования получится отрицательным — и правильный ответ нужно брать по модулю, а не пугаться «неправильного» знака.

Смешение интеграла первого и второго рода при применении формулы. Формула Грина работает исключительно для криволинейного интеграла второго рода $\oint P\,dx+Q\,dy$ — тот интеграл, что зависит от направления обхода. Попытка применить её к интегралу первого рода $\int f\,ds$ (по длине дуги, не зависящему от направления) не имеет смысла в принципе — эти два типа интегралов, как ты помнишь из прошлого урока, устроены принципиально по-разному.

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

• Формула Грина связывает криволинейный интеграл по замкнутому контуру с двойным интегралом по ограниченной им области: $\oint_{\partial D}P\,dx+Q\,dy=\iint_D\left(\dfrac{\partial Q}{\partial x}-\dfrac{\partial P}{\partial y}\right)dx\,dy$.

• Три обязательных условия применимости: контур положительно ориентирован (против часовой стрелки, область слева), область односвязна (без исключённых внутренних точек или «дырок»), функции $P$ и $Q$ непрерывно дифференцируемы во всей замкнутой области целиком.

• Выражение $\dfrac{\partial Q}{\partial x}-\dfrac{\partial P}{\partial y}$ — это плоский аналог ротора поля; если оно тождественно равно нулю, поле потенциально, и циркуляция по любому замкнутому контуру равна нулю.

• Доказательство строится по принципу «разбей на простые кусочки, докажи для одного, сложи — внутренние границы взаимно сократятся».

• Площадь плоской фигуры вычисляется как частный случай формулы Грина: $S=\oint x\,dy=-\oint y\,dx=\dfrac12\oint(x\,dy-y\,dx)$.

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

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

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

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

• Циркуляция чисто градиентного (антиградиентного) поля по любому замкнутому контуру в пространстве весов модели всегда равна нулю — это прямое следствие формулы Грина и теоремы о равенстве смешанных производных.

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

Формула Грина опирается напрямую на криволинейные интегралы второго рода и понятие потенциальности поля из предыдущего урока — критерий $\partial P/\partial y=\partial Q/\partial x$ был введён там же, но именно формула Грина даёт ему строгое обоснование через сведение к двойному интегралу. Она также использует двойные интегралы и переход к полярным координатам, изученные несколькими уроками ранее, — без уверенного владения этой техникой вычисление площадей через формулу Грина было бы затруднительно.

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

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

Джордж Грин почти никогда не покидал окрестности своего родного Ноттингема и до тридцати пяти лет вообще не публиковал ни одной научной работы — вся его знаменитая формула была выведена самоучкой без учителей, без научного руководителя и без доступа к современным ему передовым журналам, только по книгам из скромной местной библиотеки.

Работа Грина 1828 года была издана настолько малым тиражом (по разным оценкам, около полусотни экземпляров) и настолько скромно распространена, что практически двадцать лет оставалась почти неизвестной научному сообществу — её случайно обнаружил и оценил по достоинству Уильям Томсон, будущий лорд Кельвин, уже после смерти автора.

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

Формула площади многоугольника (формула шнурков), которую мы использовали в нескольких примерах и заданиях для вычисления площади по координатам вершин, — это точный дискретный аналог непрерывной формулы Грина, и она до сих пор ежедневно используется в геоинформационных системах, компьютерной графике и алгоритмах компьютерного зрения.

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

Прежде чем применять формулу Грина, всегда сначала проверь особые точки: найди, где знаменатели функций $P$ и $Q$ обращаются в ноль (или где иначе нарушается дифференцируемость), и убедись, что таких точек нет внутри контура. Это займёт двадцать секунд, а сэкономит от ошибки, которую потом придётся долго искать.

Если разность $\dfrac{\partial Q}{\partial x}-\dfrac{\partial P}{\partial y}$ получается тождественно равной константе или простому выражению, почти всегда выгоднее перейти к двойному интегралу, а не вычислять циркуляцию напрямую по контуру, особенно если контур состоит из нескольких кусков с разными параметризациями.

Для вычисления площади через контурный интеграл выбирай симметричную форму $S=\dfrac12\oint(x\,dy-y\,dx)$ — она часто даёт более простые промежуточные выражения после подстановки параметризации, чем несимметричные варианты $\oint x\,dy$ или $-\oint y\,dx$.

Если контур задан в полярных координатах, сразу переходи к формуле $S=\dfrac12\int r^2(\theta)\,d\theta$ — не переводи параметризацию обратно в декартовы координаты, это лишний и ненужный шаг.

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

Когда нужно проверить потенциальность поля быстрее, чем искать сам потенциал, вычисли именно $\partial Q/\partial x-\partial P/\partial y$ — если это выражение равно нулю, ты сразу знаешь, что циркуляция по любому замкнутому контуру равна нулю, и можешь пропустить утомительный поиск функции-потенциала, если она тебе не нужна для самой задачи.

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