Диференціальні рівняння

Диференці́альні рівня́ння[1][2][3] — рівняння, що встановлюють залежність між незалежними змінними, числами (параметрами), невідомими функціями та їхніми похідними. Невідома функція може бути як скалярною, так і векторною.

Візуалізація повітряного потоку з рівняння Нав'є-Стокса

Візуалізація теплообміну у корпусі насоса, отримана шляхом розв'язування рівняння теплопровідності

Ісаак Ньютон

Ґотфрід Лейбніц

Такі залежності віднаходяться в різних областях знань: у механіці, фізиці, хімії, біології, економіці та ін. Диференціальні рівняння широко використовуються на практиці, зокрема для опису перехідних процесів, коливань, теплопровідності, деформації балок і пластин, поширення електричного струму у провіднику тощо.

У застосуваннях математики часто виникають задачі, в яких залежність одного параметра від іншого є невідомою, але можливо записати вираз для швидкості зміни одного параметра відносно іншого (похідної). У цьому випадку задача зводиться до знаходження функції за її похідною відносно з деяких інших виразів.

Диференціальні рівняння, або теорія диференціальних рівнянь — розділ математики, який розглядає теорію та способи розв'язування диференціальних рівнянь.

Основні поняття і означення

У випадку одного аргументу диференціальне рівняння називається звичайним; у випадку декількох аргументів — диференціальним рівнянням з частинними похідними. Складнішими є інтегро-диференціальні рівняння.

Порядком диференціального рівняння називається найвищий порядок похідної, що входить до рівняння.

Степенем диференціального рівняння називається найвищий степінь, до якого піднесено похідну найбільшого порядку n, що входить до рівняння.

Розв'язком диференціального рівняння порядку n називається функція, що має похідні, до n-ного порядку включно на деякому інтервалі, підставлення якої у рівняння перетворює його у тотожність. Якщо рівняння має розв'язок, то не один, а нескінченну множину; розв'язок може залежати не лише від аргументу, але також від однієї або декількох довільних сталих чи функцій. Якщо розв'язок рівняння отримано у формі неявної функції, то його називають інтегралом  рівняння.

Початковими умовами або граничними умовами називаються додаткові умови, що накладаються на функцію при розв'язку конкретної задачі, що приводить до диференціального рівняння. За цих умов розв'язок може виявитись єдиним. Розв'язок рівняння, що залежить від довільних сталих, кількість яких дорівнює порядку рівняння і які можуть бути підібраними так, щоб задовольнити будь-яким початковим та граничним умовам, що допускають єдиний розв'язок, називається загальним розв'язком. Частинним розв'язком диференціального рівняння називається будь-який розв'язок, що може бути отриманий із загального при визначених числових значеннях довільних сталих. Довільні сталі, що входять в загальний розв'язок, визначаються з початкових або граничних умов.

Диференціальне рівняння називається інтегровним в квадратурах, якщо задачу знаходження усіх розв'язків можна звести до обчислення скінченного числа інтегралів від відомих функцій і простих алгебраїчних операцій. Через те, що багато рівнянь не можуть бути виражені через прості функції, тому деякі, рішення, що часто зустрічаються в таких задачах, отримали власні назви, були досліджені їх значення і взаємозв'язок, і тепер вони входять у число спеціальних функцій.

Спочатку диференціальні рівняння виникли із задач механіки, в яких брали участь координати тіл, їхні швидкості та прискорення, розглянуті як функції від часу, пізніше вони знайшли застосування практично в усіх розділах фізики — такі основні для своїх областей рівняння як рівняння Максвелла в електродинаміці, рівняння Ейнштейна у загальній теорії відносності та рівняння Шредінгера у квантовій механіці є диференціальними. Багато моделей з інших наук, таких як біологія, хімія і економіка також описуються різноманітними диференціальними рівняннями.

Для багатьох з цих рівнянь, в тому числі практично важливих, наприклад, рівняння Нав'є-Стокса, допоки що не знайдено розв'язку в загальному вигляді. Проте в реальних задачах за допомогою чисельних методів можна знайти їх рішення з будь-якою необхідною точністю.

Історія

Диференціальні рівняння винайдені Ньютоном (1642—1727). Ньютон вважав цей свій винахід настільки важливим, що зашифрував його у вигляді анаграми, смисл якої в сучасних термінах можна вільно передати так: «закони природи виражаються диференціальними рівняннями».

Основним аналітичним досягненням Ньютона було розкладання всіляких функцій в ступеневі ряди (сенс другої, довгої анаграми Ньютона в тому, що для вирішення будь-якого рівняння потрібно підставити в рівняння ряд і прирівняти члени однакового степеня). Особливе значення мала тут відкрита ним формула бінома Ньютона (зрозуміло, не тільки з цілими показниками, для яких формулу знав, наприклад, Вієт (1540—1603), але і, що особливо важливе, з дробовими і негативними показниками). Ньютон розклав в «ряди Тейлора» всі основні елементарні функції (раціональні, радикали, тригонометричні, експоненту і логарифм). Це, разом з складеною ним таблицею первісних (яка перейшла в майже незмінному вигляді в сучасні підручники аналізу), дозволяло йому, за його словами, порівнювати площі будь-яких фігур «за половину чверті години».

Ньютон указував, що коефіцієнти його рядів пропорційні послідовним похідним функції, але не зупинявся на цьому детально, оскільки він справедливо вважав, що всі обчислення в аналізі зручніше проводити не за допомогою кратних диференціювань, а шляхом обчислення перших членів ряду. Для Ньютона зв'язок між коефіцієнтами ряду і похідними був скоріше засобом обчислення похідних, чим засобом складання ряду. Одним з найважливіших досягнень Ньютона є його теорія сонячної системи, викладена в «Математичних принципах натуральної філософії» («Principia») без допомоги математичного аналізу. Зазвичай вважають, що Ньютон відкрив за допомогою свого аналізу закон всесвітнього тяжіння. Насправді Ньютону (1680) належить лише доказ еліптичності орбіт в полі тяжіння за законом зворотних квадратів: сам цей закон був вказаний Ньютону Гуком (1635—1703) і, мабуть, вгадувався ще декількома вченими.

З величезного числа робіт XVIII століття з диференціальних рівнянь виділяються роботи Ейлера (1707—1783) і Лагранжа (1736—1813). У цих роботах була передусім розвинена теорія малих коливань, а отже — теорія лінійних систем диференціальних рівнянь; попутно виникли основні поняття лінійної алгебри (власні числа і вектори в n-мірному випадку). Характеристичне рівняння лінійного оператора довго називали секулярним, оскільки саме з такого рівняння визначаються секулярні (вікові, тобто повільні в порівнянні з річним рухом) збурення планетних орбіт згідно з теорією малих коливань Лагранжа. Услід за Ньютоном Лаплас і Лагранж, а пізніше Гаус (1777—1855) розвивають також методи теорії збуджень.

Коли була доведена нерозв'язність алгебраїчних рівнянь в радикалах, Жозеф Ліувілль (1809—1882) побудував аналогічну теорію для диференціальних рівнянь, встановивши неможливість рішення низки рівнянь (зокрема таких класичних, як лінійні рівняння другого порядку) в елементарних функціях і квадратурі. Пізніше Софус Лі (1842—1899), аналізуючи питання про інтегрування рівнянь в квадратурі, прийшов до необхідності детально досліджувати групи дифеоморфізмів (що отримали згодом ім'я груп Лі) — так з теорії диференціальних рівнянь виникла одна з найплідніших областей сучасної математики, подальший розвиток якої був тісно пов'язаний зовсім з іншими питаннями (алгебри Лі ще раніше розглядали Сімеон-Дені Пуассон (1781—1840) і, особливо, Карл Густав Якоб Якобі (1804—1851)).

Новий етап розвитку теорії диференціальних рівнянь починається з робіт Анрі Пуанкаре (1854—1912), створена ним «якісна теорія диференціальних рівнянь» разом з теорією функцій комплексних змінних привела до заснування сучасної топології. Якісна теорія диференціальних рівнянь, або, як тепер її частіше називають, теорія динамічних систем, зараз розвивається найактивніше і має найважливіші застосування теорії диференціальних рівнянь в природознавстві.

Звичайні диференціальні рівняння

Звичайні диференціальні рівняння — це рівняння виду $F(t,x,x',x'',...,x^{(n)})=0$ , де $x=x(t)$ — невідома функція (можливо, вектор-функція; в такому випадку часто говорять про систему диференціальних рівнянь), що залежить від змінної часу $t$ , штрих означає диференціювання по $t$ . Число $n$ називається порядком диференціального рівняння.

Розв'язком (або рішенням) диференціального рівняння називається функція, що диференціюється n разів, і задовольняє рівнянню в усіх точках своєї області визначення. Зазвичай існує ціла множина таких функцій, і для вибору однієї з них на розв'язок потрібно накласти додаткові умови: наприклад, вимагати, щоб рішення приймало в певній точці певне значення.

Основні завдання і результати теорії диференціальних рівнянь: існування і єдиність рішення різних задач для ЗДР, методи розв'язання простих ЗДР, якісне дослідження рішень ЗДР без знаходження їхнього явного вигляду.

Диференціальні рівняння в частинних похідних

Диференціальні рівняння в частинних похідних — це рівняння, що містять невідомі функції від декількох змінних та їх частинних похідних.

Загальний вид таких рівнянь можна представити у вигляді:

F\left(x_{1},x_{2},\dots ,x_{m},z,{\frac {\partial z}{\partial x_{1}}},{\frac {\partial z}{\partial x_{2}}},\dots ,{\frac {\partial z}{\partial x_{m}}},{\frac {\partial ^{2}z}{\partial x_{1}^{2}}},{\frac {\partial ^{2}z}{\partial x_{1}\partial x_{2}}},{\frac {\partial ^{2}z}{\partial x_{2}^{2}}},\dots ,{\frac {\partial ^{n}z}{\partial x_{m}^{n}}}\right)=0

,

де $x_{1},x_{2},\dots ,x_{m}$ — незалежні змінні, а $z\!$ — функція цих змінних.

Лінійні та нелінійні диференціальні рівняння

Як звичайні диференціальні рівняння, так і рівняння у частинних похідних можна поділити на лінійні та нелінійні.

Лінійні диференціальні рівняння

Диференціальне рівняння є лінійним, якщо невідома функція і її похідні входять у рівняння лише у першому степені (й не перемножаються одна з одною). Для таких рівнянь розв'язки утворюють афінний підпростір простору функцій. Теорія лінійних диференціальних рівнянь розвинена значно глибше, ніж теорія нелінійних рівнянь. Загальний вигляд лінійного диференціального рівняння $n$ -го порядку:

p_{n}(x)y^{(n)}(x)+p_{n-1}(x)y^{(n-1)}(x)+\cdots +p_{0}(x)y(x)=r(x),

де $p_{i}(x)$ — відомі функції незалежної змінної, що називаються коефіцієнтами рівняння. Функція $r(x)$ у правій частині називається вільним членом (єдиний доданок, що є незалежним від невідомої функції). Важливим частковим класом лінійних рівнянь є лінійні диференціальні рівняння із сталими коефіцієнтами.

Підкласом лінійних рівнянь є однорідні диференціальні рівняння — рівняння, що не мають вільного члена: $r(x)=0$ . Для однорідних диференціальних рівнянь виконується принцип суперпозиції: лінійна комбінація часткових розв'язків такого рівняння також буде його розв'язком. Усі інші лінійні диференціальні рівняння називаються неоднорідними диференціальними рівняннями.

Нелінійні диференціальні рівняння

Нелінійне диференціальне рівняння — це рівняння, в якому невідомою величиною є деяка функція та у диференціальне рівняння входить не лише вона, але й різні її похідні в нелінійному виді. Розрізняють звичайні нелінійні диференціальні рівняння і нелінійні диференціальні рівняння в частинних похідних.

Нелінійні диференціальні рівняння виникли із задач нелінійної механіки, в яких фігурували координати тіл, їх швидкості та прискорення, розглянуті як функції від часу.

Нелінійні диференціальні рівняння у загальному випадку не мають розроблених методів розв'язування, крім деяких часткових випадків. В деяких випадках (із застосуванням тих чи інших наближень) вони можуть бути зведені до лінійних. Наприклад, лінійне рівняння гармонічного осцилятора ${\frac {d^{2}y}{dx^{2}}}+\omega ^{2}y=0$ може розглядатись як наближення нелінійного рівняння математичного маятника ${\frac {d^{2}y}{dx^{2}}}+\omega ^{2}\sin y=0$ для випадку малих амплітуд, коли $y\approx \sin y$ .

Існування розв'язків

Розв'язування диференціальних рівнянь не схоже на розв'язування алгебраїчних рівнянь. Рішення часто є неочевидним і складним, і навіть такі питання як існування рішення, або його унікальність є напрочуд цікавими.

Для рівнянь першого порядку, теорема Пеано дає набір достатніх умов, при яких у рівняння є принаймні один розв'язок. Теорема формулюється так:

Нехай функція $f(t,x)$ неперервна в деякій області $|t-t_{0}|\leqslant a,|x-x_{0}|\leqslant b$ і $\beta$ - максимум $|f(t,x)|$ у цій області. Якщо $h=min{\mathcal {f}}(a,{\frac {b}{\beta }}{\mathcal {g}})$ , то на відрізку $[t_{0}-h,t_{0}+h]$ існує принаймні одне рішення рівняння ${\frac {dx}{dt}}=f(t,x)$ , що задовільняє умові $x(t_{0})=x_{0}$ .

Для рівнянь більш високих порядків існує наступна теорема:

Для рівняння

y^{(n)}+y^{(n-1)}f_{1}(x)+\cdots +f_{n-1}(x)y'+f_{n}(x)y=g(x)

Тоді для будь-якого x₀ з інтервалу [a, b], на якому функції f_n, …, f₀ і g(x) є визначенними, і для будь яких С₁, С₂ … С_n, що задають граничні умови

y(x_{0})=C_{1},y'(x_{0})=C_{2},\cdots ,y^{(n)}(x_{0})=C_{n}

Існує одне і тільки одне рішення, і воно також визначене на цьому інтервалі.[4]

Точні розв'язки

Деякі диференціальні рівняння мають розв'язки, що можна подати точною формулою. Такі класи рівнянь подані нижче.

В таблиці, H(x), Z(x), H(y), Z(y), чи H(x,y), Z(x,y) — довільні інтегровні функції від x чи y (або від обидвох параметрів), a A, B, C, I, L, N, M — константи. В загальному A, B, C, I, L, є дійсними числами, а N, M, P та Q можуть бути комплексними. Диференціальні рівняння подані в альтернативній формі, що дозволяє їх розв'язати методом інтегрування.

	Диференціальні рівняння	Загальний розв'язок
1	${\frac {\mathrm {d} y}{\mathrm {d} x}}=F(x)\,\!$ $\mathrm {d} y=F(x)\,\mathrm {d} x\,\!$	$y=\int F(x)\,\mathrm {d} x\,\!$
2	${\frac {\mathrm {d} y}{\mathrm {d} x}}=F(y)\,\!$ $\mathrm {d} y=F(y)\,\mathrm {d} x\,\!$	$x=\int {\frac {\mathrm {d} y}{F(y)}}\,\!$
3	$H(y){\frac {\mathrm {d} y}{\mathrm {d} x}}+Z(x)=0\,\!$ $H(y)\,\mathrm {d} y+Z(x)\,\mathrm {d} x=0\,\!$	$\int H(y)\,{\mathrm {d} y}+\int Z(x)\,{\mathrm {d} x}=C\,\!$
4	${\frac {\mathrm {d} y}{\mathrm {d} x}}+H(x)y+Z(x)=0\,\!$ $\mathrm {d} y+H(x)y\,\mathrm {d} x+Z(x)\,\mathrm {d} x=0\,\!$	$y=-e^{-\int H(x)\,\mathrm {d} x}\int e^{\int H(x)\mathrm {d} x}Z(x)\,{\mathrm {d} x}\,\!$
5	${\frac {\mathrm {d} y}{\mathrm {d} x}}=F\left({\frac {y}{x}}\right)\,\!$	$\ln Cx=\int {\frac {\mathrm {d} r}{F(r)-r}}\,\!$ розв'язком може бути неявна фунція від x та y, отримана обчисленням наведеного інтегралу використовуючи заміну змінних $r=y/x\,\!$
6	${\frac {\mathrm {d} ^{2}y}{\mathrm {d} x^{2}}}=F(y)\,\!$	$x=\pm \int {\frac {\mathrm {d} y}{\sqrt {2\int F(y)\,\mathrm {d} y+C_{1}}}}+C_{2}\,\!$
7	$H(x,y){\frac {\mathrm {d} y}{\mathrm {d} x}}+Z(x,y)=0\,\!$ $H(x,y)\,{\mathrm {d} y}+Z(x,y)\,{\mathrm {d} x}=0\,\!$	Якщо ДР є точним, тобто ${\frac {\partial H}{\partial x}}={\frac {\partial Z}{\partial y}}\,\!$ тоді розв'язок задається формулою: $F(x,y)=\int \left[H(x,y)\,\mathrm {d} y+Z(x,y)\,\mathrm {d} x\right]+\gamma (y)+\chi (x)=C\,\!$ де $\gamma (y)\,\!$ та $\chi (x)\,\!$ — певні функції, залежні від інтегралів, що дозволяють коректно визначити функцію $F(x,y)\,\!$ hold. Якщо рівняння не є точним, з функцій H(x,y) та Z(x,y) можна визначити інтегральний множник, після домноження рівняння на який воно розв'язується аналогічно до точного.
8	${\frac {\mathrm {d} ^{2}y}{\mathrm {d} x^{2}}}+I{\frac {\mathrm {d} y}{\mathrm {d} x}}+Ly=0\,\!$	Якщо $I^{2}>4L\,\!$ тоді $y=Ne^{\left(-I+{\sqrt {I^{2}-4L}}\right){\frac {x}{2}}}+Me^{-\left(I+{\sqrt {I^{2}-4L}}\right){\frac {x}{2}}}\,\!$ Якщо $I^{2}=4L\,\!$ тоді $y=(Ax+B)e^{-Ix/2}\,\!$ Якщо $I^{2}<4L\,\!$ тоді $y=e^{-I{\frac {x}{2}}}\left[P\sin {\left({\sqrt {\left\|I^{2}-4L\right\|}}{\frac {x}{2}}\right)}+Q\cos {\left({\sqrt {\left\|I^{2}-4L\right\|}}{\frac {x}{2}}\right)}\right]\,\!$
9	$\sum _{\alpha =1}^{d}I_{\alpha }{\frac {\mathrm {d} ^{\alpha }y}{\mathrm {d} x^{\alpha }}}=0\,\!$	$\sum _{\alpha =1}^{d}A_{\alpha }e^{B_{\alpha }x}=0\,\!$ де $B_{\alpha }\,\!$ — d розв'зки поліному степеня d: $\prod _{\alpha =1}^{d}\left(B-B_{\alpha }\right)=0\,\!$

Зауважте, що 3 і 4 є частковими випадками 7, вони досить поширені і презентовані для повноти.

Також 8 рівняння є частковим випадком 9, але 8 досить поширена форма рівнянь, особливо у простих фізичних та інженерних задачах.

Приклади

Другий закон Ньютона можна записати у формі диференціального рівняння

m{\frac {d^{2}x}{dt^{2}}}=F(x,t)

,

де $m$ — маса тіла, $x$ — його координата, $F(x,t)$ — сила, діюча на тіло з координатою $x$ у момент часу $t$ . Його розв'язком є траєкторія руху тіла під дією вказаної сили.

Коливання струни описується рівнянням

{\frac {\partial {}^{2}u}{\partial t^{2}}}=a^{2}{\frac {\partial ^{2}u}{\partial x^{2}}}

,

де $u=u(x,t)$ — відхилення струни в точці з координатою $x$ у момент часу $t$ , параметр $a$ задає властивості струни.

Диференціальне рівняння прогину пластини під дією рівномірно розподіленого навантаження $q$ :

{\frac {\partial ^{4}w}{\partial x^{4}}}+{\frac {\partial ^{4}w}{\partial x^{2}\partial y^{2}}}+{\frac {\partial ^{4}w}{\partial y^{4}}}={\frac {q}{D}}

,

де $w$ — вертикальні прогини пластини, $D$ — циліндрична жорсткість пластини при згині.

Програмне забезпечення

Існує програмне забезпечення, яке може розв'язувати диференціальні рівняння:

ExpressionsinBar
Maple:[5] dsolve
SageMath
Xcas:[6] desolve(y'=k*y,y)

Див. також

Асимптотично еквівалентні системи

Примітки

Диференціальні рівняння
ЗВИЧАЙНІ ДИФЕРЕНЦІАЛЬНІ РІВНЯННЯ ПЕРШОГО ПОРЯДКУ
Диференціальні рівняння для інформатиків
Виленкин Н. Я. и др. Дифференциальные уравнения
dsolve - Maple Programming Help. www.maplesoft.com. Процитовано 12 травня 2020.
Symbolic algebra and Mathematics with Xcas.

Література

Українською

Бугрій О.М.; Процах Н.П.; Бугрій Н.В. (2011 р.). Основи диференціальних рівнянь: теорія, приклади та задачі : Навчальний посібник. Львів. ISBN 978-966-2645-01-9.
Диференціальні рівняння : навч. посіб. / Головатий Ю. Д., Кирилич В. М., Лавренюк С. П.; Львів. нац. ун-т ім. Івана Франка. - Львів : ЛНУ ім. І. Франка, 2011. - 468 с. - ISBN 978-966-613-859-3
Диференціальні рівняння : навч. посіб. / [Каленюк П. І. та ін.] ; Нац. ун-т "Львів. політехніка". - Львів : Вид-во Львів. політехніки, 2014. - 378 с. - ISBN 978-617-607-564-6
Диференціальні рівняння : навч. посіб. / [Л. С. Тесленко та ін.] ; Миколаїв. нац. ун-т ім. В. О. Сухомлинського. - Миколаїв : Іліон, 2013. - 336 с. : рис. - ISBN 978-617-534-141-4
Диференціальні рівняння : навч. посіб. / Т. П. Гой, О. В. Махней ; Прикарпат. нац. ун-т ім. Василя Стефаника. - Івано-Франківськ : Сімик, 2012. - 351 с. - ISBN 978-966-8067-90-7
Диференціальні та інтегральні рівняння : підручник / Кривошея С.А., Перестюк М.О., Бурим В.М. - К. : Либідь, 2004. - 407 с.: рис. - ISBN 966-06-0348-7
Лінійні динамічні системи і звичайні диференціальні рівняння : навч. посібник / П. М. Гащук. - Львів : Українські технології, 2002. - 607 с.: рис. - ISBN 966-666-024-5
Самойленко А. М.; Перестюк М. О.; Парасюк I.О. (2003 р.). Диференціальні рівняння. Київ: Либідь. ISBN 966-06-0249-9. Архів оригіналу за 17 червня 2014. Процитовано 2 грудня 2015.

Іншими мовами

Teschl, Gerald (2012). Ordinary Differential Equations and Dynamical Systems. Graduate Studies in Mathematics 140. American Mathematical Society. ISBN 978-0-8218-8328-0. (англ.)
Tenenbaum, Morris; Pollard, Harry (1985). Ordinary Differential Equations. Dover. ISBN 0-486-64940-7. (англ.)
Hartman, Philip (2002) [1964]. Ordinary differential equations. Classics in Applied Mathematics 38. Philadelphia: Society for Industrial and Applied Mathematics. ISBN 978-0-89871-510-1. (англ.)
Strauss, Walter A. (2008). Partial Differential Equations: An Introduction (вид. 2nd). John Wiley & Sons. ISBN 978-0470054567. (англ.)
Pontryagin, Lev (1962). Ordinary Differential Equations. Adiwes International Series in Mathematics. Pergamon Press. (англ.)
Л. С. Понтрягин. Обыкновенные дифференциальные уравнения — 4-е изд. — Москва, 1974. (рос.) (Підручник удостоєний державної премії СРСР)
В. И. Арнольд. Обыкновенные дифференциальные уравнения — 4-е изд. — Ижевск: Удм. гос. университет, 2000. — 367 с. ISBN 5-89806-029-4 (рос.)

Посилання

Диференційні рівняння // Вища математика в прикладах і задачах / Клепко В.Ю., Голець В.Л.. — 2-ге видання. — К. : Центр учбової літератури, 2009. — С. 468. — 594 с.
Gabriel Nagy (2016), Ordinary Differential Equations, Michigan State University. (англ.)
ДИФЕРЕНЦІА́ЛЬНИХ РІВНЯ́НЬ ТЕО́РІЯ (ДРТ) //ЕСУ (укр.)

This article is issued from Wikipedia. The text is licensed under Creative Commons - Attribution - Sharealike. Additional terms may apply for the media files.

[1] Диференціальні рівняння

[2] ЗВИЧАЙНІ ДИФЕРЕНЦІАЛЬНІ РІВНЯННЯ ПЕРШОГО ПОРЯДКУ

[3] Диференціальні рівняння для інформатиків

[4] Виленкин Н. Я. и др. Дифференциальные уравнения

[5] solve - Maple Programming Help. www.maplesoft.com. Процитовано 12 травня 2020.

[6] Symbolic algebra and Mathematics with Xcas.