CC BY
12
Том 26, № 136
2021
© Ланеев Е.Б., Анисимов В.А., Лесик П.А., Ремезова В.И., Романов A.A., Хегай А.Г., 2021 DOI 10.20310/2686-9667-2021-26-136-394-403 ШШ
УДК 519.6
■8
OPEN fil ACCESS
Об одной некорректно поставленной краевой задаче для метагармонического уравнения в круговом цилиндре
Евгений Борисович ЛАНЕЕВ, Виктор Александрович АНИСИМОВ, Полина Александровна ЛЕСИК, Виктория Ивановна РЕМЕЗОВА, Андрей Андреевич РОМАНОВ, Анна Георгиевна ХЕГАЙ
ФГАОУ ВО «Российский университет дружбы народов» 117198, Российская Федерация, г. Москва, ул. Миклухо-Маклая, 6
Аннотация. Рассматривается смешанная по краевым условиям задача для метагармонического уравнения в области, представляющей собой часть кругового цилиндра. Эту цилиндрическую область с одной стороны ограничивает поверхность общего вида, на которой заданы условия Коши, т. е. заданы функция и ее нормальная производная. Другая граница цилиндрической области свободна. На боковой поверхности цилиндрической области заданы однородные краевые условия первого рода. Задача некорректно поставлена и ее приближенное решение, устойчивое к погрешности в данных Коши, построено с применением методов регуляризации. Рассматриваемая задача сведена к интегральному уравнению Фредгольма первого рода. На основе решения интегрального уравнения, полученного в виде ряда Фурье по собственным функциям первой краевой задачи для уравнения Лапласа в круге, построено явное представление точного решения поставленной задачи. Устойчивое решение интегрального уравнения построено методом регуляризации Тихонова. В качестве приближенного решения интегрального уравнения рассматривается экстремаль функционала Тихонова. На основе этого решения строится приближенное решение задачи в целом. Приведена теорема сходимости приближенного решения поставленной задачи к точному при стремлении к нулю погрешности в данных Коши и при согласовании параметра регуляризации с погрешностью в данных. Результаты работы могут быть использованы для математической обработки данных тепловидения в ранней диагностике в медицине.
Ключевые слова: некорректно поставленная задача, метагармоническое уравнение, функции Бесселя, интегральное уравнение первого рода, метод регуляризации Тихонова
Для цитирования: Ланеев Е.Б., Анисимов В.А., Лесик П.А., Ремезова В.И., Романов А.А., Хегай А.Г. Об одной некорректно поставленной краевой задаче для метагармонического уравнения в круговом цилиндре // Вестник российских университетов. Математика. 2021. Т. 26. № 136. С. 394-403. DOI 10.20310/2686-9667-2021-26-136-394-403.
© Е.В. Laneev, V.A. Anisimov, P. A. Lesik, V.I. Remezova, A. A. Romanov, A. G. Ivhegai, 2021 DOI 10.20310/2686-9667-2021-26-136-394-403 Ш
On an ill-posed boundary value problem for a metaharmonic equation in a circular cylinder
Evgeniy B. LANEEV, Viktor A. ANISIMOV, Polina A. LESIK, Viktoriya I. REMEZOVA, Andrey A. ROMANOV, Anna G. KHEGAI
Peoples' Friendship University of Russia (RUDN University) 6 Miklukho-Maklay St., Moscow 117198, Russian Federation
Abstract. In this paper, we consider a mixed problem for a metaharmonic equation in a domain in a circular cylinder. The cylindrical area is bounded on one side by an arbitrary surface on which the Cauchy conditions are set, i. e. the function and its normal derivative are set. The other border of the cylindrical area is free. On the lateral surface of the cylindrical domain, homogeneous boundary conditions of the first kind are given. The problem is ill-posed and its approximate solution, stable to errors in the Cauchy data, is constructed using regularization methods. The problem is reduced to a first kind Fredholm integral equation. Based on the solution of the integral equation obtained in the form of a Fourier series by the eigenfunctions of the first boundary value problem for the Laplace equation in a circle, an explicit representation of the exact solution of the problem is constructed. A stable solution of the integral equation is obtained by the method of Tikhonov regularization. The extremal of the Tikhonov functional is considered as an approximate solution. Based on this solution, an approximate solution of the problem as a whole is constructed. A theorem on convergence of the approximate solution of the problem to the exact one as the error in the Cauchy data tends to zero and the regularization parameter is matched with the error in the data, is given. The results can be used for mathematical processing of thermal imaging data in early diagnostics in medicine.
Keywords: ill-posed problem, metaharmonic equation, Bessel function, integral equation of the first kind, method of Tikhonov
Mathematics Subject Classification: 35R25, 35R30.
For citation: Laneev E.B., Anisimov V.A., Lesik P.A., Remezova V.I., Romanov A.A., Khegai A.G. Ob odnoy nekorrektno postavlennoy krayevoy zadache dlya metagarmonicheskogo uravne-niya v krugovom tsilindre [On an ill-posed boundary value problem for a metaharmonic equation in a circular cylinder]. Vestnik rossiyskikh universitetov. Matematika - Russian Universities Reports. Mathematics, 2021, vol. 26, no. 136, pp. 394-403. DOI 10.20310/2686-9667-2021-26136-394-403. (In Russian, Abstr. in Engl.)
Введение
В работе рассматривается некорректно поставленная смешанная краевая задача для метагармонического уравнения в круговом цилиндре с условиями Коши на поверхности общего вида. Такая задача возникает в медицинской диагностике как задача обработки термографических данных [1] с целью выявления патологий у пациента, которые могут быть сопоставлены аномальным стационарным источникам тепла в теплопроводящей среде со стационарным распределением температуры, удовлетворяющим стационарному уравнению теплопроводности, т. е. уравнению Лапласа. Учет влияния кровотока приводит к метагармоническому уравнению [2]. Функция температуры на поверхности тела, регистрируемая как термограмма, дает представление о плотности распределения источников тепла внутри тела. Как правило, такое представление сильно искажено за счет удаленности источников тепла от поверхности исследуемого тела. Уточненную информацию об источниках можно получить, анализируя распределение температуры вблизи источников, решая задачу продолжения метагармонической функции с границы в сторону источников по известной функции и нормальной производной на границе. Следуя работам [3, 4, 5], в которых решена соответствующая задача для уравнения Лапласа в цилиндрической области прямоугольного сечения [3, 4] и в области кругового цилиндра [5], краевая задача для метагармонического уравнения приведена к линейному интегральному уравнению первого рода, устойчивое решение которого строится на основе метода регуляризации Тихонова [6]. Рассматриваемая здесь задача аналогична задаче продолжения поля ньютоновского потенциала с плоской поверхности [7].
1. Постановка задачи
Имеется цилиндрическая область кругового сечения
D(F,H ) = {(r,p,z):0 < r < a, 0 2n, F (r,<p) <z<H},
ограниченная поверхностью общего вида
S = {(r,ip,z) : 0 < r < a, 0 < ^ < 2n, z = F(r, <p) < H}. В этой области рассмотрим краевую задачу для метагармонического уравнения
Au(M) - k2u(M) = 0, M е D(F, H), k = const, u\s = f,
du = (1.1)
dn s g,
u\r=a = 0,
включающую условия Коши на поверхности S и однородные условия первого рода на боковой поверхности цилиндра. Граница z = H открыта.
Будем считать, что функции f и g непрерывны на S и обеспечивают существование решения u е C2(D(F, H)) П C 1(D(F,H)).
Смешанная задача (1.1) с условиями Коши некорректно поставлена. Решение задачи неустойчиво по отношению к погрешности в данных f и д. Получим явное выражение точного решения задачи в зависимости от точных данных Коши f и g.
2. Точное решение краевой задачи
Пусть p(M, P) — функция источника первой краевой задачи для метагармонического уравнения
Av(P) - k2v(P) = -p(P), P e D
v\ r=a = 0,
v ^ 0 при \z\ ^ то в бесконечном круговом цилиндре
= {(r,p,z) : 0 <r < a, 0 < p < 2п, -то < z < то}.
Функция источника есть сумма фундаментального решения метагармонического уравнения и метагармонической по P функции W(M, P)
p(M, P) = eXp{-krMP} + W(M, P), (2.1)
4пгмр
где Гмр — расстояние между точками M и P, и удовлетворяет граничным условиям. Функция источника (2.1) может быть представлена в виде разложения по собственным функциям оператора Лапласа в круге радиуса a
1 ^^ ^^ Jn(—n ~~м) Jn(—n ~~р)
P(M, P) = -tn eXP{ knm \ zM - zP|}-k aV ( m)]2 a COS П(РМ - PP),
Па n=0 m=\ knm Jn (-n )]
e0 = 0.5, tn = 1 при n = 0.
Здесь Jn — функция Бесселя, -f, m =1, 2, 3,... — нули функции Jn, а также введено обозначение _
knm =\j (-f )2 + k2. (2.2)
Пусть М Е Б (Г, И). Применяя формулы Грина в области Б (Б, И) к функции и(Р) — решению задачи (1.1) и функции источника р(М,Р), получим
и(М)= У [дп(РММ,Р) - и(Р)(М,Р)]<Ър, М Е Б(Б,И). (2.3)
дВ(Е,И)
Учитывая однородные граничные условия для р и и на боковой поверхности цилиндрической области Б (Г, И), получим
u(M)= I g(P)p(M,P) - f (P)(M,P)
da
p
S
+ J
n(H)
где второй интеграл берется по кругу
(P )P(M-P)- u(P) dp [M.F-)
daP, (2.4)
П(Н ) = {(r, p,z) : 0 < r < a, 0 <p< 2n, z = H}.
(2.5)
Обозначим
f г дт
daP, zM < H, (2.6)
ф(м)=/ g(P)т(м,р) - f (p)дт(M,P)
5
v(M)= ^ [Щ(Р)<р(М,Р) - и(Р)(М,Р)](ар, гы < Н, (2.7)
П(Н)
тогда решение задачи (1.1) получим в виде
и(М) = ^(М) + Ф(М), М е Б(Р,Н), (2.8)
где функция Ф вычисляется по известным функциям f и д. Функцию V можно рассматривать как решение задачи
Дv(M) - k2v(M) = 0, М е Б(—то,Н),
v\z=н = vн, v\r=a = 0,
V ^ 0 при г ^ —то.
Отсюда следует, что если решение задачи (1.1) существует, то функция v может быть представлена в виде ряда Фурье
оо оо
v(M) = Y^ J2eXP{knm(ZM - H[(vHYnm COS nm + (vH)Snm Sin , (2.9)
n=0 m=1
2n a
(vh)cnm = —2r jn 2 dm rdrvH(r,m)Jn(^m-)cos nт, na2[J'n(tfm)\2 J J a
0 0
eo = 0.5, tn =1 при n = 0, (2.10)
2n a
(vh)Sm = —2r ví 2 dm rdrvH(r,m)Jn(»m-)sinnт,
na2[Jn(№)\ J J a
00
причем ряд (2.9) сходится равномерно области d(—to, H — е) при любом е > 0, так как
i r i i i
|(vh)nmexp{knm(zM — H)]Jn{^m~) cos nm1 ^ |(vh)nm| exp{—knme}.
Аналогичная оценка справедлива и для (г)н)Пт.
Таким образом, из представления (2.8) решения задачи (1.1) и (2.9) следует, что для получения явного выражения для точного решения задачи (1.1) достаточно выразить функцию Vн в (2.9) через заданные функции f и д.
Из (2.9), (2.10) функция v может быть выражена через Vн в виде интеграла
v(M) = У 0(М,Р^н(Р)гр(грскрр, М е Б(—то,Н), (2.11)
п(н)
где
2 ~ ~ ^(кт—к(кт -)
G(M, р) = exp{knm(zM — H)}-J ( m)\2 a cosп(тм — mp),
П n=0 m=l a [ n(^n )]
t0 = 0.5, tn =1 при n = 0. (2.12)
Покажем, что функция Ун удовлетворяет интегральному уравнению Фредгольма первого рода. Пусть М Е Б(-то, Г), где
Б(-то,Г) = {(г, р,г):0 < г < а, 0 <р< 2п, -то (т,р)}.
Применяя формулу Грина в области Б(Г,И) к функции и(Р) — решению задачи (1.1) и функции р(М,Р) вида (2.1), аналогично (2.3) получим
^(P)p(M, P) - u(P)-dp(M, P) dap, M Е D(F, H). ion dnP
(2.15)
0 = J
dD(F,H)
Отсюда с учетом однородных граничных условий для p и u и обозначений (2.6) и (2.7) получим
v(M) = -Ф(М), M Е D(-to, F). (2.13)
Пусть b < min F(r, p) и M Е n(b), где n(b) — область вида (2.5) при z = b, тогда из
(r,v)
(2.13) и (2.11) получим интегральное уравнение первого рода
J G(M,P)vh(P)dxpdyp = -Ф(M), M Е П(Ь). (2.14)
П(Н)
Из уравнения (2.14) с учетом разложения (2.12) при zm = b получаем соотношение между коэффициентами Фурье единственного решения Vh и коэффициентами Фурье правой части
-(VH)Cnm exP{-knm(H - b)} = ^m^
-(VH)Cm exP{ knm(H - b)} = Km(b),
где фCm(b), фCm(b) — коэффициенты Фурье функции Ф^)|men(b) :
2n a
(^)Cnm(b)= na2JJ^rn)]2 J dP J rd^(r,P,b)Jn{vma) COS np, n C 0 0
e0 = 0.5, cn = 1 при n = 0,
2n a
($)Cm(b) = na2[J2(^m)]^ dP / rdr^rpp^Jn^ a ) Sin np.
0 0
Отметим, что формулы (2.15) характеризуют убывание коэффициентов Фурье фnm(b) с ростом n и к, если функции f и д таковы, что обеспечивают существование решения задачи (1.1) и, следовательно, функции vh вида (2.9). Подставляя коэффициенты Фурье (Фн)nm из (2.15) в ряд (2.9), получим функцию v в области D(-to,H)
те те
v(M) = - ^ ^ exp{knm(H - b)} JnOCa)[(ф)Пто(Ь) cos np + (<l)nm(b) sin np]. (2.16)
n=0 m=1
Ряд (2.16), как и ряд (2.9), сходится равномерно в D(-to, H — е) при любом е > 0, если решение задачи (1.1) существует при данных f и д.
Формула (2.8), где функции v и Ф вида (2.16) и (2.6) соответственно, дает явное выражение для решения задачи 1.1.
3. Построение приближенного задачи в случае приближенных данных Коши
Пусть функции f и д в задаче (1.1) заданы с погрешностью, то есть вместо f и д заданы функции f5 и д*, такие что
— f 1112(50 ^ 6, \\д* — д\\ь2[я) ^ 6.
Построим приближенное решение задачи (1.1), сходящееся к точному решению при 6 ^ 0. Функция Ф вида (2.6) в этом случае может быть получена приближенно:
Фй (M )=/ gs (P )f(M, P) - f (P) ^ (M,P)
dap. (3.1)
S
Применяя неравенство Коши-Буняковского к разности функций (3.1) и (2.6) при M Е n(b), b < min F(r,f), получим оценку погрешности правой части интегрального
уравнения (2.14)
|Ф'(M) - Ф(М)|мод ^ max ( / f2(M,P)dap)1/2\\g' - g\\b2{s)
м en(b) J
S
+ MS) (/ ^M P)Td*p)1,2\\fl - f \U(S) < a. SP
В качестве приближенного решения интегрального уравнения (2.14) с приближенной правой частью (3.1) будем рассматривать экстремаль функционала Тихонова [6] со стабилизатором нулевого порядка
M>] = \\Gw + Ф'||L2(n(b)) + а\Н\|2(п(Я)), а > 0. (3.2)
Экстремаль функционала может быть получена как решение уравнения Эйлера для функционала (3.2)
G*Gw + aw = -G* Ф',
которое, в свою очередь, с использованием разложения (2.12) ядра интегрального оператора G может быть представлено в виде алгебраических уравнений относительно коэффициентов Фурье функции w
exp{-2knm(H - b)}wnm + awnm = - exp{-knm(H - Ь)}Ф6пт(Ь), (3.3)
где Фnm(b) — в одном случае косинус-коэффициенты, в другом — синус-коэффициенты Фурье функции Ф'(M)|Men(b)
2п a
^nm(b) = (Ф' )Cnm(b) = df\ ^тФв (т^Ь)^ ^) COS Wf,
n n 0 0
e0 = 0.5, tn =1 при n = 0,
2n a
^nm (b) = (Ф S )nm(b) = na2[J2(^m)]^ dff ЫТФ (T, b) Jn ^ ) Sin Wf.
n n 0 0
Решая алгебраические уравнения относительно коэффициентов Фурье экстремали и подставляя их вместо коэффициентов Фурье функции Ун в ряд (2.9), найдем приближение Уёа к функции у в области Б (-то, И) :
S (M) = _ ^^ exP{ knm (zM - b)}Jn(-ma ЖФ * )Cnm(b) COS nP + (Ф * )Пп (b) sin nP] 'Ja( ) 1 + a exp{2knm(H - b)}
n=0 m=1
где knm определено формулой (2.2). Функция (3.4) отличается от точной функции v вида (2.16) регуляризирующим множителем (1 + a exp{2knm(H - b)}) , обеспечивающим сходимость ряда.
В соответствие с (2.8) приближенное решение задачи (1.1) построим в виде
ui(M)= va(M) + Ф*(M), M e D(F,H), (3.5)
где vOa и Ф* — функции вида (3.4) и (3.1) соответственно.
Для построенного приближенного решения (3.5) задачи (1.1) имеет место следующая теорема сходимости к точному решению.
Теорема 3.1. Пусть решение задачи (1.1) существует. Тогда для любого a = a(8) такого, что a(8) ^ 0 и 8/^/a(8) ^ 0 при 8 ^ 0 функция ua(S) вида (3.5) равномерно
сходится при 8 ^ 0 к точному решению в D(F+е, H - е), 0 < е < 0,5(H - max F (r, p)).
(r,^)
Доказательство теоремы 3.1 повторяет доказательство соответствующей теоремы в [3].
Формулы (3.5), (3.4) и (3.1) дают, таким образом, приближенное решение поставленной задачи (1.1).
4. Приложение результатов к обратной задаче термографии
Построенное решение задачи (1.1) может быть использовано для решения обратной задачи термографии [8] в приложении к задачам математической обработки термограмм в тепловизионных исследованиях в медицине [1]. При моделировании участка тела пациента к задаче (1.1) приводит модель теплопроводящего тела цилиндрической формы, содержащего стационарные источники тепла. С учетом кровотока [2] стационарное распределение температуры внутри тела описывается метагармоническим уравнением. На боковой поверхности цилиндра поддерживается постоянная температура, что соответствует первому краевому условию в (1.1), а на поверхности S имеет место конвективный теплообмен с внешней средой, описываемый законом Ньютона. В этом случае поток тепла через поверхность S, т. е. нормальная производная, прямо пропорционален разности температур внутри тела и снаружи, что соответствует третьему краевому условию. Если распределение температуры на поверхности S (термограмма) может быть измерено как функция f, например, тепловизионными методами, то в рамках описанной модели оказывается известной и нормальная производная, что приводит к рассмотренной здесь задаче (1.1) восстановления пространственного распределения температуры, аномалии которого могут быть связаны с патологиями внутренних органов пациента. Кроме того, след функции распределения температуры на плоскости z = H, более близкой к источникам, чем поверхность S, может рассматриваться как скорректированная термограмма, более точно передающая изображение тепловыделяющей структуры, связанной с патологиями.
References
[1] Г. Р. Иваницкий, "Тепловидение в медицине", Вестник РАН, 76:1 (2006), 48-58. [G.R. Ivanitskii, "Thermovision in medicine", Herald of the Russian Academy of Sciences, 76:1 (2006), 48-58 (In Russian)].
[2] J. P. Agnelli, A. A. Barrea, C.V. Turner, "Tumor location and parameter estimation by thermography", Mathematical and Computer Modelling, 53:7-8 (2011), 1527-1534.
[3] Е. Б. Ланеев, Б. Васудеван, "Об устойчивом решении одной смешанной задачи для уравнения Лапласа", Вестник РУДН. Серия Прикладная математика и информатика, 1999, №1, 128133. [Е. B. Laneev, B. Vasudevan, "On a stable solution of a mixed problem for the Laplace equation", PFUR Reports. Series: Applied mathematics and computer science, 1999, № 1, 128133 (In Russian)].
[4] Е. Б. Ланеев, "О построении функции Карлемана на основе метода регуляризации Тихонова в некорректно поставленной задаче для уравнения Лапласа", Дифференциальные уравнения, 54:4 (2018), 483-491; англ. пер.:Е. B. Laneev, "Construction of a Carleman Function Based on the Tikhonov Regularization Method in an Ill-Posed Problem for the Laplace Equation", Differential Equations, 54:4 (2018), 476-485.
[5] Е. Б. Ланеев, Д. Ю. Быков, А. В. Зубаренко, О.Н. Куликова, Д. А. Морозова, Е.В. Шу-нин, "Об одной некорректно поставленной краевой задаче для уравнения Лапласа в круговом цилиндре", Вестник российских университетов. Математика, 26:133 (2021), 35-43. [Е.B. Laneev, D. Yu. Bykov, A. V. Zubarenko, O.N. Kulikova, D. A. Morozova, E. V. Shunin, "On an ill-posed boundary value problem for the Laplace equation in a circular cylinder", Russian Universities Reports. Mathematics, 26:133 (2021), 35-43 (In Russian)].
[6] А.Н. Тихонов, В. Я. Арсенин, Методы решения некорректных задач, Наука, М., 1979. [A. N. Tikhonov, V. Ya. Arsenin, Methods for Solving Ill-posed Problems, Nauka Publ., Moscow, 1979 (In Russian)].
[7] А. Н. Тихонов, В. Б. Гласко, О. К. Литвиненко, В. Р. Мелихов, "О продолжении потенциала в сторону возмущающих масс на основе метода регуляризации", Изв. АН СССР. Физика Земли, 1968, №1, 30-48. [A.N. Tikhonov A.N., V.B. Glasko, O.K. Litvinenko, V.R. Melikhov, "On the continuation of the potential towards the perturbing masses based on the regularization method", Izv. AN SSSR. Fizika Zemli, 1968, № 1, 30-48 (In Russian)].
[8] Е. Б. Ланеев, М. Н. Муратов, "Об одной обратной задаче к краевой задаче для уравнения Лапласа с условием третьего рода на неточно заданной границе", Вестник РУДН. Серия Математика, 10:1 (2003), 100-110. [Е. B. Laneev, M.N. Muratov,, "An inverse problem to a boundary value problem for the Laplace equation with a condition of the third kind on an inexactly specified boundary", PFUR Reports. Series: Mathematics, 10:1 (2003), 100-110 (In Russian)].
Информация об авторах
Ланеев Евгений Борисович, доктор физико-математических наук, профессор Математического института им. С.М. Никольского. Российский университет дружбы народов, г. Москва, Российская Федерация. E-mail: elaneev@yandex.ru ORCID: https://orcid.org/0000-0002-4255-9393
Анисимов Виктор Александрович, студент магистратуры, Математический институт им. С.М. Никольского. Российский университет дружбы народов, г. Москва, Российская Федерация. E-mail: dm.yurievich@mail.ru
Information about the authors
Evgeniy B. Laneev, Doctor of Physics and Mathematics, Professor of S.M. Nikol'skii Mathematical Institute. Peoples' Friendship University of Russia (RUDN University), Moscow, Russian Federation. E-mail: elaneev@yandex.ru ORCID: https://orcid.org/0000-0002-4255-9393
Viktor A. Anisimov, Master's Degree. S.M. Nikol'skii Mathematical Institute. Peoples' Friendship University of Russia (RUDN University), Moscow, Russian Federation. E-mail: dm.yurievich@mail.ru
Лесик Полина Александровна, аспирант, Математический институт им. С.М. Никольского. Российский университет дружбы народов, г. Москва, Российская Федерация. E-mail: polinalesik@yandex.ru
Ремезова Виктория Ивановна, студент, Математический институт им. С.М. Никольского. Российский университет дружбы народов, г. Москва, Российская Федерация. E-mail: remezova.98@mail.ru
Романов Андрей Андреевич, аспирант, Математический институт им. С.М. Никольского. Российский университет дружбы народов, г. Москва, Российская Федерация. E-mail: an1romanov@gmail.com
Хегай Анна Георгиевна, студент магистратуры, Математический институт им. С. М. Никольского. Российский университет дружбы народов, г. Москва, Российская Федерация. E-mail: annhegay98@gmail.com
Конфликт интересов отсутствует.
Для контактов:
Ланеев Евгений Борисович E-mail: elaneev@yandex.ru
Поступила в редакцию 30.08.2021 г. Поступила после рецензирования 15.11.2021 г. Принята к публикации 27.11.2021 г.
Polina A. Lesik, Post-Graduate Student. S.M. Nikol'skii Mathematical Institute. Peoples' Friendship University of Russia (RUDN University), Moscow, Russian Federation. E-mail: polinalesik@yandex.ru
Viktoriya I. Remezova, Student. S.M. Nikol'skii Mathematical Institute. Peoples' Friendship University of Russia (RUDN University), Moscow, Russian Federation. E-mail: remezova.98@mail.ru
Andrey A. Romanov, Post-Graduate Student. S.M. Nikol'skii Mathematical Institute. Peoples' Friendship University of Russia (RUDN University), Moscow, Russian Federation. E-mail: an1romanov@gmail.com
Anna G. Khegai, Master's Degree. S.M. Nikol'skii Mathematical Institute. Peoples' Friendship University of Russia (RUDN University), Moscow, Russian Federation. E-mail: annhegay98@gmail.com
There is no conflict of interests.
Corresponding author:
Evgeniy B. Laneev E-mail: elaneev@yandex.ru
Received 30.08.2021 Reviewed 15.11.2021 Accepted for press 27.11.2021
