Что такое краевые условия дифференциального уравнения

Лекция 4. Краевые задачи для дифференциальных уравнений. Задача Штурма-Лиувилля.

Будем рассматривать однородное линейное уравнение второго порядка

Ly ≡ a2(x)y» + a1(x)y’ + a0(x)y = 0. Его можно записать по-другому:

Что такое краевые условия дифференциального уравнения(15)

Однородное уравнение Ly = 0 и неоднородное Ly = f, как известно, имеют бесконечное множество решений. На практике часто бывает нужно из множества решений выделить только одно. Для этого задают некоторые дополнительные условия. Если это начальные условия у(х0) = уo, y'(xo) = y1, то получают задачу Коши. Если задают дополнительные условия на концах некоторого отрезка, то получают задачу, которая называется краевой задачей. Условия, которые задаются на концах отрезка, называются краевыми условиями. Краевые условия иногда именуют также граничными условиями и тогда говорят о граничной задаче.
Мы будем задавать линейные краевые условия вида

Что такое краевые условия дифференциального уравнения(16)

где α1, α2, β1, β2, A, B — заданные числа, причем по крайней мере одно из чисел α1, α2, и одно из чисел β1, β2, отличны от нуля. Если в (16) хотя бы одно из чисел А и В не равно нулю, то краевые условия называют неоднородными. Если А = В = 0, то условия (16) называются однородными. Краевая задача называется однородной, если рассматривается однородное уравнение (15) Ly = 0 и однородные краевые условия (16). Решением краевой задачи называется такое решение дифференциального уравнения, которое удовлетворяет заданным краевым условиям. Заметим сразу, что однородная краевая задача всегда имеет решение у ≡ 0 (тривиальное решение).

Наряду с уравнением (15) рассмотрим уравнение

Что такое краевые условия дифференциального уравнения(17)

содержащее некоторый числовой параметр λ. Здесь функции р(х), q(x), r(x) действительные, а число λ может быть, вообще говоря, и комплексным. Краевая задача (17), (16) при А = В = 0 является однородной. Поэтому при любых λ она имеет тривиальное решение. Нас будут интересовать такие значения λ, при которых эта задача обладает не только тривиальными решениями.

Задача Штурма-Лиувилля. Найти те значения параметра λ, при которых уравнение (17) имеет нетривиальное решение, удовлетворяюшее однородным краевым условиям (16). В дальнейшем будем ее записывать в виде

<Lλy = 0, l1y = 0, l2y = 0>.
Те значения параметра λ, при которых задача Штурма-Лиувилля имеет ненулевое решение, называются собственными значениями (собственными числами) задачи, а сами эти решения — собственными функциями. Задачу Штурма-Лиувилля называют также задачей на собственные значения. В силу однородности уравнения и краевых условий собственные функции задачи Штурма-Лиувилля определены с точностью до постоянного множителя. Это означает, что если y(х) -собственная функция при некотором значении λ, то произведение Cy(x), где С — произвольная постоянная, также является собственной функцией при том же значении параметра λ. В связи с этим часто в качестве собственной функции рассматривают нормированную функцию у <х), у которой ||у(х)|| = 1. Такая собственная функция определена, по существу, однозначно (с точностью до знака ±). Далее мы подробно изучим наиболее простой случай задачи Штурма-Лиувилля, когда уравнение имеет вид

y» + λy = 0.(18)

Из множества краевых условий вида (16) ограничимся тремя частными случаями:

1) краевые условия первого рода

y(a) = y(b) = 0,(19)

2) краевые условия второго рода

y'(a) = y'(b) = 0,(20)

3) краевые условия третьего рода

Что такое краевые условия дифференциального уравнения(21)

Общая задача Штурма-Лиувилля будет обладать свойствами, очень похожими на свойства в этих простых случаях, если на коэффициенты уравнения (17) наложить дополнительные условия: р(х), q(x), f(x) -непрерывные функции, причем р(х) имеет, кроме того, непрерывную производную на [а, b], р(х) > 0, q(x) ≥ 0.

Основные свойства собственных значений и собственных функций задачи Штурма-Лиувилля.

Лемма. Определитель Вронского двух собственных функций задачи Штурма-Лиувилля на концах отрезка [а, b] равен нулю.

Доказательство. Напомним, что определителем Вронского функций у = y1(x) и у = у2(x) называется определитель вида

Что такое краевые условия дифференциального уравнения

Рассмотрим однородные краевые условия общего вида (16). Пусть у1(x) и у2(x) — две любые собственные функции. Это означает, что в точке x = а выполняются равенства

Что такое краевые условия дифференциального уравнения

Числа α1, и α2 не могут одновременно равняться нулю. Значит, алгебраическая система двух однородных уравнений с двумя неизвестными имеет ненулевое решение. Это возможно только в том случае, когда определитель этой системы равен нулю:

Что такое краевые условия дифференциального уравнения

Этот определитель совпадает с определителем Вронского в точке x = а, то есть W(a) = 0.

Аналогичные рассуждения, проведенные для точки x = b, показывают, что W(b) = 0.

Свойство 1. Две собственные функции задачи Штурма-Лиувилля, соответствующие одному и тому же собственному значению λ, линейно зависимые.

Доказательство. Так как собственные функции являются решениями одного и того же однородного уравнения (17) (по условию число λ одно), то в случае их линейной независимости определитель Вронского не равен нулю ни в одной точке отрезка [а, b]. Это противоречит только что доказанной лемме. Следовательно, y1(x) и у2(x) — линейно зависимые функции.

Свойство 2. Две собственные функции у1(x) и у2(x), соответствующие различным собственным значениям λ1 и λ21 ≠ λ2), на отрезке [а, b] ортогональны.

Доказательство этого свойства проведем для собственных функций такой задачи, в которой уравнение имеет вид (18). Составим определитель Вронского функций у1 и у2 и продифференцируем его:

Что такое краевые условия дифференциального уравнения

Так как у1 и у2 — решения уравнения (18) при λ = λ1 и λ = λ2, соответственно, то получим

Что такое краевые условия дифференциального уравнения

Проинтегрируем по отрезку [а, b] левую и правую части полученного равенства. С учетом леммы будем иметь

Что такое краевые условия дифференциального уравненияПо условию λ1 — λ2 ≠0, следовательно Что такое краевые условия дифференциального уравнения
Функции y1(x) 0 и у2(х) 0, поэтому Что такое краевые условия дифференциального уравнения
Значит, y1(x) и у2(х) на отрезке [а, b] ортогональны.

Если уравнение, входящее в задачу Штурма-Лиувилля, имеет вид (17), где r(х) > 0 и r(x) 1, то под ортогональностью функций в этом случае подразумевают ортогональность с весом r(х): две функции y1(x) и у2(х) ортогональны на отрезке [а, b] с весом r(x), если

Что такое краевые условия дифференциального уравнения

Под нормой функции ||у(x)|| в этом случае также подразумевают весовую норму:

Что такое краевые условия дифференциального уравнения

Свойство 3. Собственные функции, соответствующие различным собственным значениям, образуют линейно независимую систему функций.

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

Свойство 4. Собственные значения задачи Штурма-Лиувилля действительные.

Доказательство. Предположим, что задача Штурма-Лиувилля <Lλy = 0, l1y = 0, l2y = 0> имеет комплексное собственное значение λ = α + βi,β ≠ 0. Пусть ему соответствует собственная функция у(х) (вообще говоря, тоже комплекснозначная). Так как все коэффициенты уравнения и краевых условий имеют действительные значения, то

Что такое краевые условия дифференциального уравнения

Здесь черта означает переход к комплексно сопряженному выражению. В нашем случае

Что такое краевые условия дифференциального уравнения

Значит число Что такое краевые условия дифференциального уравнениятакже является собственным значением той же задачи Штурма-Лиувилля и ему соответствует собственная функция Что такое краевые условия дифференциального уравнения. Так как в силу свойства 2 функции y(x) и Что такое краевые условия дифференциального уравненияортогональны на [а, b], то

Что такое краевые условия дифференциального уравнения

Отсюда следует, что у(x) ≡ 0 на [а, b]. Значит ни одно комплексное число λ не может быть собственным значением.

Свойство 5. Пусть коэффициенты уравнения (17) удовлетворяют условиям: р(х), q(x), r(x) — непрерывные функции и, кроме того, р(х) имеет непрерывную производную на [а, b], р(х) > 0, q(x) > 0, r(х) > 0. Тогда задача Штурма-Лиувилля <Lλ y = 0, l1 y = 0, l2 y = 0> имеет бесконечное число собственных значений λ 1, λ2, . λn, . Если краевые условия имеют вид (19) или (20), или (21), то собственные значения соответствующей задачи Штурма-Лиувилля удовлетворяют неравенствам

Что такое краевые условия дифференциального уравнения

Теорема Стеклова.Всякая непрерывная функция f(x), удовлетворяющая однородным краевым условиям : l1f = 0 и l2f = 0 , и имеющая непрерывные производные до второго порядка на отрезке [а, b], разлагается на этом отрезке в сходящийся ряд Фурье по собственным функциям yn(х) задачи Штурма-Лиувилля <Lλ y = 0, l1 y = 0, l2 y = 0> :

Что такое краевые условия дифференциального уравнения

где коэффициенты Фурье Сn вычисляются по формулам:

Что такое краевые условия дифференциального уравнения

Эта теорема применяется при решении уравнений математической физики методом Фурье.

Решение задач Штурма-Лиувилля

Вначале рассмотрим уравнение (18) y» + λy = 0. и краевые условия первого рода (19) y(a) = y(b) = 0. Для удобства будем считать, что a = 0 и b = l > 0. К такой задаче можно всегда свести данную задачу, если сделать замену переменной x’ = x — a, при этом вид уравнения не изменится.

Вид общего решения уравнения (18) зависит от значений параметра λ. Разберем три случая: 1) λ 0. В первом случае обозначим λ = — k 2 . Тогда характеристическое уравнение r 2 — k 2 = 0 будет иметь действительные различные корни r1 = k, r2 = — k: Поэтому, общее решение дифференциального уравнения запишется в виде y = C1e kx + C2e -kx . Подставим краевые условия в общее решение и получим

Что такое краевые условия дифференциального уравненияОпределитель этой системы равен Что такое краевые условия дифференциального уравнения

Следовательно, система имеет только нулевое (тривиальное) решение C1 = C2 = 0. Значит, при λ 2 и получим характеристическое уравнение r 2 + k 2 = 0. Оно имеет комплексные корни r1 = ki и r2 = -ki и общее решение дифференциального уравнения в этом случае запишется в виде y = C1cos kx + C2sin kx. Подставим краевые условия в общее решение:

Что такое краевые условия дифференциального уравнения(22)

Для того, чтобы эта система имела нетривиальные решения, необходимо и достаточно, чтобы sin kl = 0. Следовательно kl = πn, то есть Что такое краевые условия дифференциального уравненияТак как Что такое краевые условия дифференциального уравнениято можно ограничиться только положительными значениями n = 1, 2, . . Таким образом, собственные значения данной задачи имеют вид Что такое краевые условия дифференциального уравненияПри этих значениях алгебраическая система (22) имеет решения:C1 = 0, C2 — любое действительное число. Подставим эти значения в общее решение дифференциального уравнения и получим собственные функции задачи Что такое краевые условия дифференциального уравнения

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

Что такое краевые условия дифференциального уравнения

По тому же алгоритму решаются задачи Штурма-Лиувилля следующего вида:

Что такое краевые условия дифференциального уравнения(23)

и

Что такое краевые условия дифференциального уравнения(24)

Эти задачи так же, как и предыдущая, при λ 0 не имеют собственных значений. В случае λ > 0 общее решение уравнения записывается в виде y = C1cos kx +C2sin kx, где Что такое краевые условия дифференциального уравненияПосле подстановки у в краевые условия, получим:

а) для задачи (23)

Что такое краевые условия дифференциального уравнения

Что такое краевые условия дифференциального уравнения

Для того, чтобы эти системы уравнений имели нетривиальные решения, необходимо и достаточно, чтобы coskl = 0. Следовательно, Что такое краевые условия дифференциального уравнениято есть Что такое краевые условия дифференциального уравненияОтрицательные значения n можно не рассматривать, так как Что такое краевые условия дифференциального уравненияТаким образом, собственные значения у этих задач одинаковые Что такое краевые условия дифференциального уравнения

Собственные функции задачи (23) имеют вид Что такое краевые условия дифференциального уравненияА у задачи (24) они другие:

Что такое краевые условия дифференциального уравнения

Некоторые отличия возникают при решении задачи Штурма-Лиувилля в случае краевых условий второго рода

y» + λy = 0, y'(0) = y'(l) = 0.(25)

Рассуждениями, аналогичными тем, которые проводились для краевых условий первого рода, можно показать, что задача (25) при λ 0. В этом случае, общее решение уравнения имеет вид y = C1cos kx + C2sin kx, Что такое краевые условия дифференциального уравненияНайдем производную этой функции и подставим в нее краевые условия (25):

Что такое краевые условия дифференциального уравнения

Эта алгебраическая система имеет нетривиальные решения тогда и только тогда, когда, sinkl = 0 то есть kl = πn или Что такое краевые условия дифференциального уравненияТаким образом, числа Что такое краевые условия дифференциального уравнениятакже являются собственными значениями задачи. Собственные функции при этих значениях имеют вид Что такое краевые условия дифференциального уравнения. Окончательно, задача (25) имеет собственные значения Что такое краевые условия дифференциального уравненияи собственные функции Что такое краевые условия дифференциального уравнения

Для задачи Штурма-Лиувилля с краевыми условиями третьего рода (21) уже не удается получить собственные значения в явном виде. В качестве примера рассмотрим одну такую задачу, когда Что такое краевые условия дифференциального уравнения

y» + λy = 0, y'(0) = y(0), y'(l) = 0.(26)

При Что такое краевые условия дифференциального уравнениязадача (26) не имеет собственных значений и собственных функций. Доказательство этого проводится так же, как и для краевых условий первого рода. При λ > 0 общее решение уравнения записывается в виде y = C1coskx + C2sinkx, где Что такое краевые условия дифференциального уравнения. После дифференцирования этой функции и подстановки её производной и самой функции в краевые условия (26) будем иметь:

Что такое краевые условия дифференциального уравненияили

Что такое краевые условия дифференциального уравнения(27)

Получившаяся алгебраическая система будет иметь нетривиальные решения только в том случае, когда

coskl — ksinkl = 0 или

ctgkl = k(28)

Уравнение (28) является трансцендентным уравнением относительно k. Оно не решается в явном виде. Однако, построив графики левой и правой частей уравнения (28), видно, что оно имеет бесконечно много решений (см. рис.13). Обозначим корни уравнения (28) через rn, n = 1,2, . . Тогда Что такое краевые условия дифференциального уравненияпри Что такое краевые условия дифференциального уравнения

Что такое краевые условия дифференциального уравнения
Рис.13

Численными методами можно найти приближенные значения rn. Из системы (27) при k = rn получим C1n = rnC2n , где C2n -произвольные постоянные. При этих значениях постоянных решения дифференциального уравнения будут иметь вид

Они являются собственными функциями краевой задачи (26) с собственными значениями Что такое краевые условия дифференциального уравнения

Видео:Задача Коши ➜ Частное решение линейного однородного дифференциального уравненияСкачать

Задача Коши ➜ Частное решение линейного однородного дифференциального уравнения

Понятие о краевых задачах для дифференциальных уравнений

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

Например, в задаче о движении материальной точки массы m под действием заданной силы Что такое краевые условия дифференциального уравнениячасто требуется найти закон движения, если в начальный момент t = t0 точка находилась в положении, характеризуемом радиус-вектором Что такое краевые условия дифференциального уравнения, а в момент t = t1 должна попасть в точку Что такое краевые условия дифференциального уравнения. Задача сводится к интегрированию дифференциального уравнения движения

Что такое краевые условия дифференциального уравнения

с краевыми условиями Что такое краевые условия дифференциального уравнения.

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

Рассмотрим краевые задачи для линейных уравнений второго порядка

Что такое краевые условия дифференциального уравнения

Что такое краевые условия дифференциального уравнения

краевые условия сводятся к нулевым условиям z(x0)=z(x1)=0, причем линейность дифференциального уравнения не нарушается. Умножим это уравнение на exp(∫p1(x)dx), получим

Что такое краевые условия дифференциального уравнения

где p(x)=exp(∫p1(x)dx). Без ограничения общности изучение исходной краевой задачи можем заменить изучением следующей краевой задачи:

Что такое краевые условия дифференциального уравнения(1Что такое краевые условия дифференциального уравнения)

с граничными условиями

Что такое краевые условия дифференциального уравнения(2Что такое краевые условия дифференциального уравнения)

Пусть функция f(x) равна нулю на всем отрезке [x0, x1], за исключением ε-окрестности точки x=s, s — ε < x < s + ε, причем

Что такое краевые условия дифференциального уравнения

Обозначим эту функцию f(x)=fε(x, s).

Определение. Функцией Грина G(x, s) краевой задачи (1Что такое краевые условия дифференциального уравнения), (2Что такое краевые условия дифференциального уравнения) называется функция, удовлетворяющая следующим условиям:

  • 1) G(x, s) непрерывна по x при фиксированном s при x0 ≤ x ≤ x1, x0 ≤ s ≤ x1;
  • 2) G(x, s) является решением соответствующего однородного уравнения

Что такое краевые условия дифференциального уравнения

на всем отрезке [x0, x1], за исключением точки x = s;
3) G(x, s) удовлетворяет граничным условиям:

Что такое краевые условия дифференциального уравнения

  • 4) В точке x = s производная G′x(x, s) должна иметь разрыв первого рода со скачком 1/p(s).
  • Что такое краевые условия дифференциального уравнения

    является решением уравнения (1Что такое краевые условия дифференциального уравнения). Краевые условия (2Что такое краевые условия дифференциального уравнения) удовлетворяются по 2 свойству. Метод построения функции Грина. Рассмотрим решение y1(x) уравнения

    Что такое краевые условия дифференциального уравнения

    определяемое начальными условиями

    Что такое краевые условия дифференциального уравнения

    Это решение, вообще говоря, не удовлетворяет второму граничному условию y(x1)=0. Предполагаем, что не существует нетривиального решения y(x) этого однородного уравнения, удовлетворяющего нулевым граничным условиям y(x0)=y(x1)=0.

    Очевидно, что решения C1y1(x), где C1 — произвольная постоянная, также удовлетворяет граничному условию y(x0)=0. Аналогично находим нетривиальное решение y2(x) однородного уравнения, удовлетворяющее второму граничному условию y2(x1)=0. Этому же условию удовлетворяют все решения C2y2(x), где C2 — произвольная постоянная.

    Функцию Грина ищем в виде

    Что такое краевые условия дифференциального уравнения

    причем постоянные C1, C2 выберем так, чтобы были выполнены свойства 1) и 4), то есть, чтобы функция G(x, s) была непрерывна по x при фиксированном s, и, в частности, непрерывна в точке x=s:

    Что такое краевые условия дифференциального уравнения

    и чтобы G′x(x, s) в точке x=s имела скачок l/p(s):

    Что такое краевые условия дифференциального уравнения

    В силу предположения, что y1(x1) ≤ 0, решения y1(x) и y2(x) линейно независимы, так как все линейно зависимые от y1(x) решения имеют вид C1y1(x) и, следовательно, при C1 ≠ 0 не обращаются в нуль в точке x1, в которой обращается в нуль решение y2(x). Следовательно, определитель системы

    Что такое краевые условия дифференциального уравнения

    являющийся определителем Вронского: W(y1(x), y2(x))=W(x) в точке x=s, отличен от нуля и постоянные C1 и C2, удовлетворяющие этой системе, легко определяются:

    Видео:Дифференциальные уравнения | задача Штурма - Лиувилля | классические краевые задачи | 1Скачать

    Дифференциальные уравнения | задача Штурма - Лиувилля | классические краевые задачи | 1

    Начальные и граничные условия

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

    Введение

    Итак, для использования численных методов при решении дифференциального уравнения необходимо дополнительные условия. Если искомая функция(концентрация, температура и т.д) является функцией времени u=u(t), то требуются начальные условия, которые являются значением этой функции в момент времени, принятый за начальный:

    Если начальная функция также зависит и от пространственных координат u=u(t,x), то начальное условие характеризуют ее распределение в пространстве в начальный момент времени:

    В последнем случае помимо начальных условий требуются еще и граничные условия, которые имеют значения функции u(t,x) на границе изучаемой системы для любого момента времени. Причем, если искомая функция зависит от нескольких пространственных координат, то необходимо задавать граничные условия по каждой из них.

    Небольшой пример

    Например для следующего уравнения:

    • начальное условие
    • 2 граничных условия по координате
    • 1 граничное условие по координате
    • 2 граничных условия по координате

    Сразу же возникает вопрос, почему именно так? Так вот, порядок производной определяет количество граничных условий для переменной. Как вы заметили, по y присутствует только первая производная, поэтому и одно граничное условие.

    Классификация граничных условий

    Для лучшего понимания рассмотрим классификацию на примере уравнения:

    будет изменятся от до , соответственно при , будет левая граница, а при , будет правая.

      Граничные условия 1-ого рода

    Записываются следующим образом:

    — функции, зависящие от , как пример:

    Граничные условия 2-ого рода

    Здесь вместо самих функций используются их первые производные.

    Граничные условия 3-ого рода

    Смешанные граничные условия

    В этом случае левое и правое граничные условия могут быть разных родов:

    Заключение

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

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

    💡 Видео

    Частное решение дифференциального уравнения. 11 класс.Скачать

    Частное решение дифференциального уравнения. 11 класс.

    Дифференциальные уравнения, 1 урок, Дифференциальные уравнения. Основные понятияСкачать

    Дифференциальные уравнения, 1 урок, Дифференциальные уравнения. Основные понятия

    1. Что такое дифференциальное уравнение?Скачать

    1. Что такое дифференциальное уравнение?

    Линейное неоднородное дифференциальное уравнение второго порядка с постоянными коэффициентамиСкачать

    Линейное неоднородное дифференциальное уравнение второго порядка с постоянными коэффициентами

    Дифференциальные уравнения. 11 класс.Скачать

    Дифференциальные уравнения. 11 класс.

    Общее и частное решение дифференциального уравненияСкачать

    Общее и частное решение дифференциального уравнения

    Краевые Задачи для Обыкновенных Дифференциальных УравненийСкачать

    Краевые Задачи для Обыкновенных Дифференциальных Уравнений

    Краевая задача.Функция Грина.Дифференциальное ур.Скачать

    Краевая задача.Функция Грина.Дифференциальное ур.

    Откуда появляются дифференциальные уравнения и как их решатьСкачать

    Откуда появляются дифференциальные уравнения и как их решать

    Геометрический смысл дифференциального уравненияСкачать

    Геометрический смысл дифференциального уравнения

    18+ Математика без Ху!ни. Дифференциальные уравнения.Скачать

    18+ Математика без Ху!ни. Дифференциальные уравнения.

    Операционное исчисление. Решить неоднородное дифференциальное уравнение 2 порядкаСкачать

    Операционное исчисление. Решить неоднородное дифференциальное уравнение 2 порядка

    13. Как решить дифференциальное уравнение первого порядка?Скачать

    13. Как решить дифференциальное уравнение первого порядка?

    Основные понятия дифференциальных уравнений от bezbotvyСкачать

    Основные понятия дифференциальных уравнений от bezbotvy

    Задача Коши для дифференциальных уравненийСкачать

    Задача Коши для дифференциальных уравнений

    Решение физических задач с помощью дифференциальных уравненийСкачать

    Решение  физических задач с помощью дифференциальных уравнений
    Поделиться или сохранить к себе: