CC BY
13
Онлайн-доступ к журналу: http: / / mathizv.isu.ru
Серия «Математика»
2019. Т. 30. С. 114-124
УДК 517.957 MSG 35L75
DOI https://doi.org/10.26516/1997-7670.2019.30.114
Многомерные точные решения системы нелинейных уравнений типа Буссинеска*
А. А. Косов1, Э. И. Семенов1, В. В. Тирских2
1 Институт динамики систем и теории управления имени, В. М. Матросова СО РАН, Иркутск, Российская Федерация
2 Иркутский государственный университет путей сообщения, Иркутск, Российская Федерация
Аннотация. Изучается система двух нелинейных уравнений в частных производных четвертого порядка. Правые части системы уравнений содержат многомерные аналоги уравнения Буссинеска, выражаемые через двукратные операторы Лапласа и квадраты градиентов искомых функций, а также линейные функции взаимосвязи. Такого рода уравнения, близкие к уравнениям Навье — Стокса, встречаются в задачах гидродинамики. Предлагается искать решение в виде анзаца, содержащего квадратичную зависимость от пространственных переменных и произвольные функции от времени. Использование предложенного анзаца позволяет декомпозировать процесс отыскания компонент решения зависящих от пространственных переменных и от времени. Для отыскания зависимости от пространственных переменных необходимо решать алгебраическую систему матричных, векторных и скалярного уравнения. Найдено общее решение этой системы уравнений в параметрическом виде. Для отыскания компонент решения исходной системы, зависящих от времени, возникает система нелинейных обыкновенных дифференциальных уравнений. Эта система сведена к одному уравнению четвертого порядка, для которого найдены частные решения. Приводится ряд примеров построенных точных решений исходной системы уравнений типа Буссинеска, в том числе выражаемые через функции Якоби по времени и анизотропные по пространственным переменным.
Ключевые слова: нелинейная система, нелинейные уравнения типа Буссинеска, редукция, точные решения.
* Работа выполнена при частичной финансовой поддержке РФФИ (проект № 1908-00746)
1. Введение
Некоторые задачи гидродинамики сводятся к нелинейному уравнению Буссинеска [1;6;8] ии = д2 (ßuxx + ßu2 +еи). Это уравнение описывает волны, которые могут перемещаться как влево, так и вправо. Хорошо известно [1] его решение в виде уединенной волны (солитона),
которое имеет вид
^ /-
u(x,t) =--7Г-,-Г) ш = су е — 4рсг,
ß COS^ (X — Lüt + Со)
где с ф 0, Со — произвольные постоянные. Многомерный аналог уравнения Буссинеска можно записать как
utt = А (ßAu + ßu2 + ей) . (1.1)
Цель нашей работы — построение точных многомерных решений системы двух нелинейных уравнений типа Буссинеска с линейными взаимосвязями:
Г utt = А (ДД-и + ß\u2 + е\и) + aiv, , .
\vtt = A (ß2Av + ß2v2 + e2v) + a2u.
Здесь и = u(x,i), v = v(x,i), t € R+, x e Rra, n € N, n > 2, А -n - мерный оператор Лапласа; V — оператор набла взятия градиента; ßi ф 0, ßi ф 0, £i, оц ф 0, г = 1, 2 — произвольные параметры.
Точные решения системы (1.2) будем отыскивать методом обобщенного разделения переменных [9; 11]
u(x,i)=^i(i)W(x)+pi(i), v(*,t)=Mt)W(*)+Mt), (1.3)
W(x) = ^(Ах,х) + (В,х)+С, (1.4)
где ipi(t), ipi(t), г = 1,2 — неизвестные функции времени; ненулевая числовая симметрическая матрица А размера пхп, постоянный вектор В ё Ги константа С € R подлежат определению. Здесь и далее (•, •) — скалярное произведение в Rra.
Отметим, что конструкция типа (1.4) ранее успешно использовалась для построения точных решений нелинейных параболических систем со степенными нелинейностями [3;4] и для нелинейной системы уравнений в частных производных первого порядка [5]. Решения некоторых нелинейных гиперболических уравнений с функциональным разделением переменных получены в [7]. Построение точных решений нелинейных систем в частных производных важно и в ряде других приложений, например в теории кинетических систем [10; 12; 13].
2. Редукция к системе ОДУ
После подстановки функций (1.3) в систему уравнений (1.2), с учетом очевидных равенств
А (¡ЗАи + ¡ли2 + ей) = ¡ЗА(Аи) + 2 ц (иАи + |Уи\2) + еАи,
|УИ^(х)|2 = (А2х,х) + 2(АВ,х) + |В|2, ДИ^(х) =ЬгА — след матрицы А, и несложных преобразований приходим к равенствам
"1
(ф'( - 2^ А)ф2 - а1ф2)
2(Ах,х) + (В,х)+С
= 2цхф\ [(А2х, х) + 2(АВ, х) + |В|2] +2[гфтА)1р1<р1+е1-ф1(Ъ А) + ац<р2,
{фЧ - 2ц2(ъ А)ф2 - а2фг)
^(Ах,х) + (В,х)+С
+ ^ =
= 2ц2ф22 [(А2*, х) + 2(АВ, х) + |В|2] +21л2(%г А)ф2^2 + е2ф2{1г А) + а2^.
Очевидно, что если числовая симметрическая матрица А, постоянный вектор В € Кга и константа С € К удовлетворяют системе алгебраических уравнений (САУ)
А = 2(т А2, В = 2(7 АВ, С = а\В\2, (2.1)
где сг/0— константа разделения, то последние равенства сводятся к следующей системе обыкновенных дифференциальных уравнений (ОДУ) второго порядка:
ф'1 - \\ф\ - а\ф2 = 0, ф2 - \2ф22 - а2ф\ = 0, (2.2)
<р" = Огфг^г + а^2 + ^ = 62ф2(р2 + а2(р 1 + и2ф2. (2.3)
Здесь введены обозначения 9\ = 2/л\ ^г А), 02 = 2/л2(Ьг А), = е\(1т А), г/2 = е2(Ьг А),
Л1 = 2/л (^гА + , Л2 = 2/л2 + ^ . (2.4)
Таким образом, проведенными рассуждениями установлена справедливость следующего утверждения.
Теорема 1. Система нелинейных уравнений типа Буссинеска (1.2) имеет точные решения (1.3), где функция Ш{х) может быть выбрана произвольным полиномом вида (1.4) с коэффициентами, удовлетворяющими САУ (2.1), функции ф¿(¿), г = 1,2 находятся из систем ОДУ (2.2), (2.3).
Отметим, что САУ вида (2.1) ранее была получена и исследована в [3;5]. Поэтому в этой работе мы приведем только пример функций (1.4), коэффициенты которой удовлетворяют САУ (2.1).
Пример 1. В трехмерном случае п = 3 для построения анизотропных по пространственным переменным точных решений системы (1.2) можно использовать, например следующие функции
Wi(x, у, z) = [5х2 + 8у2 + 5z2 - 4ху - 8xz - 4yz] + (2.5)
+ l\X - 2(11 + l2)y + l2z +а^51\ + 8hh +
w2(x, y, z) = \bx2 + 5y2 + 8z2 + 8xy + 4xz - 4yz] + (2.6)
+ l\X + l2y + 2(h — l2)z + а^51\ - 8l\l2 + bl22
Ws(x, y, z) = \8x2 + 5y2 + 5z2 - 4xy + 4xz + 8yz] + (2.7)
+2(l2 - h)x + hy + l2z + a(^,l21 - 8l\l2 + Ы22
Для построения радиально-симметричных решений будем использовать функцию
W0(x, V,z) = j^ [(х + 2al\)2 + (у + 2al2)2 + (z + 2al3)2] . (2.8) Здесь а ф 0, г = 1, 2, 3 — произвольные параметры.
3. Разрешимость нелинейной системы ОДУ (2.2)
В системах ОДУ (2.2), (2.3), полученных в результате редукции системы нелинейных уравнений (1.2), нелинейной является только система уравнений (2.2), построение общего решения которой маловероятно в силу её нелинейности и высокой размерности. Поэтому в этом разделе все внимание уделим построению частных точных решений нелинейной системы ОДУ (2.2).
Элементарными преобразованиями система ОДУ (2.2) сводится к одному нелинейному ОДУ четвертого порядка относительно функции ф\(¿) следующего вида:
(ф'{ - Xгф2)" - ^ (ф'{ - Xгф2)2 - а\а2ф\ = 0. (3.1)
Нетрудно убедиться, что если параметры уравнения (3.1) связаны равенством aiAf — СК2Л2 = 0, то уравнение (3.1) имеет частное решение i¡)\(t) = y(t), которое удовлетворяет ОДУ второго порядка
= (3.2)
В [3; 5] показано, что решением матричного уравнения (2.1) является 1 т
матрица А = —— SEmS , где Ет — диагональная матрица, у которой 2(7
на диагонали произвольным образом расположены m € {1,2,... , гг} единиц и п — т, нулей, S — произвольная ортогональная матрица. При этом
ТП
tr А = —, m<n, neu, п > 2. (3.3)
2а v 7
С учетом формулы (3.3) выражения для параметров (2.4) запишутся как
Ai = MîL±2), Л2 = Mîi±2). (3.4)
а а
Таким образом, соотношение на параметры ai Л2 — СК2Л2 = 0 , в терминах параметров системы (1.2), примет вид
а\ц\ — (12^2 = 0. (3.5)
6(7
Если Ai = 6, Л1СК1 = —4Лг или ш = -, a\Ui + 4u2 = 0, то как
m+ 2
указано в справочнике [2] ОДУ (3.2) имеет частное точное решение
y(t) = —77---, где ¿о — произвольная постоянная. Это частное
sin ¿(t — to)
решение приводит к частному решению системы ОДУ (2.2) следующего вида
да) = • 2(! Ш =--■ t. .у (3.6)
Sin' (t — to) 0>\ Sin (t — to)
Из теоремы 1 для этого частного решения приходим к справедливости утверждения.
Утверждение 1. Пусть параметры системы нелинейных уравнений типа Буссинеска (1.2) связаны соотношением (3.5) и равенствами
6(7
¡11 = ——, oliUi + 4/Х2 = 0, m € {1,2,... ,п}, 771 Zi
где сг/0— произвольная постоянная. Тогда система (1.2) имеет частные точные многомерные решения
^(Ах,х) + (В,х) + С
+ <Pi(t),
±04х,х) + (В,х) + С
где числовая симметрическая матрица А, постоянный вектор В и константа С удовлетворяют САУ (2.1), функции ф*(Ь), г = 1,2 выражаются формулами (3.6), а функции (р определяются из линейной неавтономной системы ОДУ второго порядка (2.3).
В общем случае ОДУ (3.2) сводится к следующей квадратуре
I , = (3.7)
¡л.и + ^ + с,
где С\, Ьо — произвольные постоянные. Вычислив этот интеграл, получим, что зависимость у от переменной £ задается неявно равенством, содержащим эллиптический интеграл первого рода. Однако, если кон-
Л гу3
станта С\ принимает значения С\ = 0 и С\ = —ддз-' т0 легк0 выписать явные частные точные решения уравнения (3.2) в элементарных функциях. Кроме того, если параметры А1, сц и постоянная С\ выбираются определенным образом, то можно получить точные решения уравнения (3.2) в эллиптических функциях Якоби.
Пусть С\ = 0, тогда интеграл (3.7) вычисляется в элементарных функциях и в этом случае система ОДУ (2.2) имеет частные точные решения следующих видов
I \iCXi
«'М^Н^Н"1)- (3'8)
Л ГУ3
Пусть С\ =---3-, тогда интеграл (3.7) также вычисляется в элемен-
ЗА2
тарных функциях и в этом случае система ОДУ (2.2) имеет частные точные решения следующих видов
2 I 1 ¡\lOti
Хюц 3-2со82 I *0) I Ф?а) = —^- \_2—(з.9)
С учетом этих результатов и использованием теоремы 1 приходим к справедливости утверждения
Утверждение 2. Пусть параметры системы нелинейных уравнений типа Буссинеска (1.2) связаны равенством (3.5). Тогда она имеет частные точные многомерные решения
= ф\1 и2{ъ,1)=ф?{1)
(Ах,х) + (В,х) + С
-(Ах,х) + (В,х) + С
-(Лх,х) + (В,х) + С
-(Ах,х) + (В,х)+С
где числовая симметрическая матрица А, постоянный вектор В и константа С удовлетворяют САУ (2.1), функции ф*3^), г,.] = 1,2 выражаются формулами (3.8), (3.9), а функции (р \(1;), (р2(Ь) определяются из линейной неавтономной системы ОДУ второго порядка (2.3).
Пример 2. Система нелинейных уравнений типа Буссинеска (1.2), с соотношением на параметры (3.5) и дополнительным условием е\ = е2 = 0, в трехмерном координатном пространстве имеет частные точные, анизотропные по пространственным переменным, параметрические семейства решений
щ(х,у,г,1) = ^ о) ) +1 ) \¥г(х,у,г),
цгага ( 3 - 2сое2 ~ ¿о) ы(х,у,г,Ь) =-Т—---
(* - ¿о)
где функции \Уг(х,у,г,1), г = 1,2,3 задаются формулами (2.5)—(2.7), а ф 0, ¿о — произвольные постоянные.
Пусть, теперь, параметры А1, Л2, СК1 и постоянная С\ таковы, что А1 > 0 и кубический трехчлен в интеграле (3.7) представим в виде (у — а\){у — а2)(у — аз), где вещественные постоянные а\ < а2 < аз определяются из следующей системы алгебраических уравнений
аг + а2 + аз + ^- = 0, а1а2 + а\аз + а2аз = 0, С\ + а1а2аз = 0. (3.10) 2А2
Тогда вычислив интеграл (3.7), находим
^з1^) = («2 - аг) вп2 (Т, к) + аь
\lCXi
ФГЮ =
Л2
(аг — а^вп (Т, к) + а\
(3.11)
Здесь введено обозначение Т = л/^1^3 а1) 8п(Т, /г) — эллип-
6(7
тический синус Якоби, к = \ —-- — модуль эллиптической функ-
V а3 - а\
ции. С учетом приведенных результатов и использованием теоремы 1 приходим к справедливости утверждения.
Утверждение 3. Пусть параметры системы нелинейных уравнений типа Буссинеска (1.2) связаны равенством (3.5). Тогда она имеет частные точные многомерные решения
и3(х,1)=ф*31(1)
1
(Ах,х) + (В,х) + С
-(Ах,х) + (В,х) + С
где числовая симметрическая матрица А, постоянный вектор В и константа С удовлетворяют САУ (2.1), функции фф^(Ь) выражаются формулами (3.11), в которых вещественные постоянные а\ < а2 < аз удовлетворяют системе алгебраических уравнений (3.10), а функции р> 1^), определяются из линейной неавтономной си-
стемы ОДУ второго порядка (2.3).
Отметим, что в приведенных утверждениях 1-3 для случая п = 3 (х € К3) в качестве \¥(х,у,г) можно взять любую из функций приведенных в примере 1. Так для функций определяемых формулами (2.5)—(2.7) получим анизотропные по пространственным переменным точные решения, а для функции (2.8) можно выписать радиально-сим-метричное по пространственным переменным х, у, г точное решение.
4. Заключение
Таким образом, получены формулы новых точных многомерных решений систем нелинейных уравнений типа Буссинеска. Полученные явные выражения точных многомерных решений, которые выражаются в элементарных и эллиптических функциях Якоби, могут иметь не только теоретическое, но и прикладное значение, поскольку их можно использовать для тестирования, настройки и адаптации численных методов и алгоритмов построения приближенных решений краевых задач
для нелинейных систем уравнений в частных производных четвертого порядка большой размерности.
Список литературы
1. Солитоны и нелинейные волновые уравнения / Р. Додд, Дж. Эйлбек, Дж. Гиббон, X. Моррис. М. : Мир, 1988. 694 с.
2. Камке Э. Справочник по обыкновенным дифференциальным уравнениям. М. : Наука, 1971. 576 с. https://dx.doi.org/10.1007/978-3-663-05925-7
3. Косов А. А., Семенов Э. И. О точных многомерных решениях системы уравнений реакции-диффузии со степенными нелинейностями // Сибирский математический журнал. 2017. Т. 58, № 4. С. 796-812. https:// 10.17377/smzh.2017.58.408
4. Косов А. А., Семенов Э. И. О точных многомерных решениях одной нелинейной системы уравнений реакции-диффузии // Дифференциальные уравнения. 2018. Т. 54, № 1. С. 108-122. https://10.1134/S0374064118010090
5. Kosov A. A., Semenov Е. I., Tirskikh V. V. On Exact Multidimensional Solutions of a Nonlinear System of First Order Partial Differential Equation // Известия Иркутского государственного университета. Серия Математика. 2019. Т. 28. С. 53-68. https://doi.org/10.26516/1997-7670.2019.28.53
6. Павлов М. В. Уравнение Буссинеска и преобразование Миуры // Фундаментальная и прикладная математика. 2004. Т. 10, № 1. С. 175-182.
7. Полянин А. Д., Журов А. И. Решения с функциональным разделением переменных двух классов нелинейных уравнений математической физики // Доклады АН. 2019. Т. 486, № 3. С. 287-291.
8. Полянин А. Д., Зайцев В. Ф. Справочник по нелинейным уравнениям математической физики: Точные решения. М. : Физматлит, 2002. 432 с.
9. Полянин А. Д., Зайцев В. Ф., Журов А. И. Методы решения нелинейных уравнений математической физики и механики. М. : Физматлит, 2005. 256 с.
10. Скубачевский A. J1. Уравнения Власова - Пуассона для двухкомпонентной плазмы в однородном магнитном поле //УМН. 2014. Т. 69, № 2. С. 107-148. https://doi.org/10.4213/rm9579.
11. Galactionov V. A., Svirshchevskii S. R. Subspaces of nonlinear partial differential equations in mechanics and physics. Chapman & Hall/CRC, 2007. 493 p.
12. Steady state solutions of the Vlasov-Maxwell system and their stability / Y. Markov, G. Rudykh, N. Sidorov, A. Sinitsyn, D. Tolstonogov // Acta Appl. Math. 1992. Vol. 28, N 3. P. 253-293.
13. Сидоров H. А., Синицын А. В. Стационарная система Власова - Максвелла в ограниченных областях // Нелинейный анализ и нелинейные дифференциальные уравнения. М. : Физматлит, 2003. Р. 50-88.
Александр Аркадьевич Косов, кандидат физико-математических наук, ведущий научный сотрудник, Институт динамики систем и теории управления им. В. М. Матросова СО РАН, Российская Федерация, 664033, г. Иркутск, ул. Лермонтова, 134, тел.: (3952) 427100, e-mail: [email protected] Эдуард Иванович Семенов, кандидат физико-математических наук, старший научный сотрудник, Институт динамики систем и теории управления им. В. М. Матросова СО РАН, Российская Федерация, 664033, Иркутск, ул. Лермонтова, 134, тел.: (3952) 453099, e-mail: [email protected], ORCID iD https://orcid.org/0000-0002-9768-9945
Владимир Викторович Тирских, кандидат физико-математических наук, доцент, Иркутский государственный университет путей сообщения, Российская Федерация, 664074, г. Иркутск, ул. Чернышевского, 15, тел. (3952)638311, e-mail: [email protected]
Поступила в редакцию 12.07.19
Multidimensional Exact Solutions of a System of Nonlinear Boussinesq Type Equations
A. A. Kosov1, E. I. Semenov1, V. V. Tirskikh2
lMatrosov Institute for System Dynamics and Control Theory of SB of RAS, Irkutsk, Russian Federation
2 Irkutsk State Transport University, Irkutsk, Russian Federation
Abstract. We study the system of two nonlinear partial differential equations of the fourth order. The right parts of the system of equations contain multidimensional analogs of Boussinesq equation, expressed in terms of two-fold Laplace operators and squares of gradients of the required functions, as well as linear functions of the relationship. This kind of equations, similar to Navier-Stokes equations, encountered in problems of hydrodynamics. The paper proposes to search for a solution in the form of anzatz containing quadratic dependence on spatial variables and arbitrary functions on time. The use of the proposed anzatz allows to decompose the process of finding the solution components depending on the spatial variables and time. To find the dependence on spatial variables it is necessary to solve the algebraic system of matrix, vector and scalar equations. The General solution of this system of equations in parametric form is found. To find the time-dependent components of the solution of the initial system, a system of nonlinear ordinary differential equations arises. This system is reduced to one fourth-order equation for which particular solutions are found. A number of examples of the constructed exact solutions of the initial system of Boussinesq equations, including those expressed in terms of Jacobi functions in time and anisotropic in spatial variables, are given.
Keywords: nonlinear system, nonlinear Boussinesq equations, reduction, exact solutions
References
1. Dodd R.K., Eilbeck J.C., Gibbon J.D., Morris H.C. Solitons and nonlinear waves equations. London, Academic Press, 1982, 630 p.
2. Kamke E. Differentialgleichungen: Lösungsmethoden und Lösungen I. Gewöhnliche Differentialgleichungen. B.G. Teubner, Leipzig, 1977. https://dx.doi.org/10.1007/978-3-663-05925-7
3. Kosov A.A., Semenov E.I. Multidimensional exact solutions to the reaction-diffusion system with power-law nonlinear terms. Siberian Math. J., 2017, vol. 58, no. 4, pp. 619-632. https://doi.org/10.1134/S0037446617040085
4. Kosov A.A., Semenov E.I. On Exact Multidimensional Solutions of a Nonlinear System of Reaction-Diffusion Equations. Differential Equations, 2018, vol.54, no. 1, pp. 106-120. https://10.1134/S0012266118010093
5. Kosov A.A., Semenov E.I., Tirskikh V.V. On Exact Multidimensional Solutions of a Nonlinear System of First Order Partial Differential Equation. The Bulletin
of Irkutsk State University. Series Mathematics, 2019, vol. 28, pp. 53-68. https://doi.org/10.26516/1997-7670.2019.28.53
6. Pavlov M.V. Uravnenie Bussineska i preobrazovanie Miury [Boussinesq equation and Miura transform] Fundamental'naja i prikladnaja matematika [Fundamental and Applied Mathematics], 2004, vol. 10, no. 1, pp. 175-182. (in Russian)
7. Polyanin A.D., Zhurov A.I. Functional separable solutions of two classes of nonlinear mathematical physics equations. Doklady Mathematics, 2019, vol. 99, no. 3, pp. 321-324. https://10.31857/S0869-56524863287-291
8. Polyanin A.D., Zaitsev V.F. Handbook of Nonlinear Partial Differential Equations. Chapman & Hall/CRC Press, Boca Raton, London, New York, 2012, 1912 p.
9. Polyanin A.D., Zaitsev V.F., Zhurov A.I. Metody resheniya nelinejnyx uravnenij matematicheskoj fiziki i mexaniki [Methods for solving nonlinear equations of mathematical physics and mechanics]. M.: Fizmatlit Publ., 2005, 256 p. (in Russian)
10. Skubachevskii A.L. Vlasov - Poisson equations for a two-component plasma in homogeneous magnetic field. Russian Math. Surveys, 2014, vol. 69, no. 2, pp. 291335. http://dx.doi.org/10.1070/RM2014v069n02ABEH004889
11. Galactionov V.A., Svirshchevskii S.R. Subspaces of nonlinear partial differential equations in mechanics and physics. Chapman & Hall/CRC, 2007, 493 p.
12. Markov Y., Rudykh G., Sidorov N., Sinitsyn A., Tolstonogov D. Steady state solutions of the Vlasov - Maxwell system and their stability. Acta Appl. Math., 1992, vol. 28, no. 3, pp. 253-293.
13. Sidorov N.A., Sinitsyn A.V. Stacionarnaya sistema vlasova-maksvella v ogranichennyx oblastyax [The stationary Vlasov - Maxwell system in bounded domains] «Nelinejnyj analiz i nelinejnye differencialnye uravneniya» [Nonlinear analysis and nonlinear differential equations]. Moscow, Fizmatlit Publ., 2003, pp. 50-88. (in Russian)
Alexander Kosov, Candidate of Sciences (Physics and Mathematics), Leading Researcher; Matrosov Institute for System Dynamics and Control Theory SB RAS, Post Box 292, 134, Lermontov st., Irkutsk, 664033, Russian Federation; tel.: (3952) 427100, e-mail: [email protected]
Edward Semenov, Candidate of Sciences (Physics and Mathematics), Senior Researcher; Matrosov Institute for System Dynamics and Control Theory SB RAS, Post Box 292, 134, Lermontov st., Irkutsk, 664033, Russian Federation; tel.: (3952) 453099, e-mail: [email protected], ORCID iD https://orcid.org/0000-0002-9768-9945
Vladimir Tirskikh, Candidate of Sciences (Physics and Mathematics), Associate Professor, Irkutsk State Transport University, 15, Chernyshevsky st., Irkutsk, 664074, Russian Federation; tel.: (3952) 638311, e-mail: [email protected]
Received 12.07.19
