CC BY
17
УДК 519.21
DOI: 10.18698/1812-3368-2018-2-12-22
ГРУППОВОЙ метод поиска функции римана
ДЛЯ НЕКОТОРЫХ УРАВНЕНИЙ ЭПИДЕМИИ
А.В. Мастихин М.Н. Шевченко
МГТУ им. Н.Э. Баумана, Москва, Российская Федерация
Аннотация
На примере марковского процесса эпидемии Вейса рассмотрена задача поиска функции Римана для стационарного первого уравнения Колмогорова относительно экспоненциальной (двойной) производящей функции вероятностей перехода. С помощью групповых методов найдены четырехмерная алгебра Ли сим-метрий для этого гиперболического уравнения в частных производных и функция Римана
Ключевые слова
Инфинитезимальный оператор, марковский процесс, экспоненциальная производящая функция, первое уравнение Колмогорова, функция Римана
Поступила в редакцию 16.02.2017 © МГТУ им. Н.Э. Баумана, 2018
Введение. В математической теории эпидемий базовыми примерами являются общая эпидемия, т. е. эпидемия Бартлетта — Мак-Кендрика и эпидемия Вейса, или простая эпидемия. В первой из них пара частиц двух типов (инфицированная частица — переносчик инфекции — и частица, восприимчивая к инфекции), взаимодействуя, переходит в пару инфицированных частиц. Естественно определяемый процесс тем не менее оказывается достаточно сложным для исследования, некоторые результаты получены асимптотическими методами. В эпидемии Вейса пара частиц, взаимодействуя, переходит в одну инфицированную частицу, один переносчик удаляется из популяции. Этот процесс хорошо изучен и допускает обобщения. Так, Дж. Гани предложил процесс заболевания с двумя стадиями и тремя типами частиц. Метод экспоненциальной (двойной) производящей функции вероятностей перехода для всех трех процессов приводит к одному и тому же гиперболическому уравнению.
Цель настоящей работы — проанализировать гиперболическое уравнение с помощью теории симметрий. Задача поиска функции Римана для стационарного первого уравнения Колмогорова относительно экспоненциальной (двойной) производящей функции вероятностей перехода рассмотрена на примере марковского процесса эпидемии Вейса. Четырехмерная алгебра Ли симметрий для этого гиперболического уравнения в частных производных и ее одномерная подалгебра, инвариантная относительно задачи Гурса, найдена с помощью групповых методов. Это позволяет получить функцию Римана, применяемую при выводе предельных теорем.
Определение процесса. Уравнение Колмогорова для двойной производящей функции. Пусть t ) = ),..., (0), I е[0, да) — однородный во вре-
мени марковский процесс со счетным множеством состояний Nn = {а = = (аь... , <Xn ): 0,1,2,—, i — 1,.., nj. Состояние системы, характеризуемое вектором a = (ai,...,an), означает наличие совокупности из а1-частиц типа T1, а2-частиц типа T2,..., an-частиц типа Tn. Возможные переходы системы из одного состояния в другое представляются соответствующей схемой взаимодействий. Обозначим переходные вероятности как [1]
Paß(t) = P{£(t) = aß(0) = ß}, a, ßeNn.
Введем экспоненциальную (двойную) производящую функцию переходных вероятностей (|s| < 1)
z a
F (t, 5, z ) = E— Paß (t )sP
a,ß a !
и линейные дифференциальные операторы с постоянными коэффициентами
hk Ц Иp , k=1, -1,
соответствующие распределению вероятностей
j pk > 0, Ypk = 1, pk = 0 j , k = 1,..., I,
определяемому процессом [2]. Двойная производящая функция переходных вероятностей F(t, s, z) удовлетворяет дифференциальному уравнению
б 1 k ( ( б Л бек ^
— F(t,s,z) = ^XkZs hk I —I--r F(t, s, z), F(0, s, z) = esz,
ot k=1 ^ \ 6z) dzs y
называемому первым уравнением Колмогорова для двойной производящей функции. Здесь sk — число a-частиц; Xk — положительный коэффициент, характеризующий интенсивность взаимодействия частиц из k-комплекса.
Пусть у е Nn — поглощающие состояния процесса. Вероятности попадания из состояния а в поглощающее состояние у называются финальными вероятностями [3].
Для нахождения финальных вероятностей qay используют производящую функцию
Фа( s ) = £ qay st, |s| < 1,
У
и экспоненциальную производящую функцию
za
Ф (z, s ) = Z Ф a(s ).
W !
а. ^ •
Рассуждая так же, как и при доказательстве теоремы 3.1, приведенной в работе [2], можно показать, что
Ф (z, 5 ) = lim Т (t, 5, z)
t ^rn
и Ф(г,5) удовлетворяет стационарному первому уравнению
I J Г дЛ дек л
Z^kZsk к =1
v dz J Qz
у v У (У/, у
ek
Ф (5, Z )= 0.
Примеры. Процесс эпидемий Вейса, Бартлетта — Мак-Кендрика, Гани.
Рассмотрим стационарные уравнения на трех примерах марковских процессов, интерпретируемых как процессы эпидемии [4]. Для процесса эпидемии Вейса, следуя изложенному выше алгоритму, приведем подробную постановку задачи. Начнем с определения процесса. На дискретной четверти плоскости N2 ={а = (аьа2): аьа2 = 0, 1, 2,...} рассмотрим однородный во времени марковский процесс £,(t) = (^1 (t), 2 (t)), [0,да), с переходными вероятностями Р((р^2) = Р £ (t) = (р1, р2) | £ (0) = («1, а2)}. При & ^0 полагаем
раьЙ = 1" (°Ча2 + ^) ^ + 0 (■^);
P,
а1а 2 At+0 (At);
PtÜ )=^At + 0 (At ),
где ц>0. Событие t ) = (а1, а2 ) интерпретируется как наличие совокупности
из а1-частиц типа T1 и а2-частиц типа T2. Частицы рассматривают как больные (типа T2) и здоровые (типа T1) особи. Принятая запись для схемы взаимодействий T + T2 ^ T1, T1 ^ 0. Для приложений важна задача об оценке числа ^(а1,а2) финальных частиц типа T2 (оставшихся здоровых особей). Случайная величина ^(аьа2) имеет распределение {?(оу2)2), У2 = 0,..., а2Ь где Я((°ч,°)2) = lim P((a b"2)(t).
Аналитический подход [2] к решению задачи о финальных вероятностях ^(™у2)2) для поглощающих состояний (0, у2), где у2 = 0,1,2,..., приводит для экспоненциальной производящей функции
Ф (Z1, Z2,5)= Е ^ Е 5»
а1, а2 = 0 • ^2 • у2 = 0
к дифференциальному уравнению в частных производных гиперболического типа
(1)
с граничными условиями
Ф (0, z2, s ) = ez2s, Ф (zb0, s ) = ez1.
(2)
Для марковских процессов с взаимодействием частиц, таких как процессы эпидемии Вейса [5], эпидемии Гани [6] и процесса общей эпидемии Бартлетта — Мак-Кендрика [7], задача о финальных вероятностях приводит к задачам Гурса для гиперболических уравнений.
Приведем стационарные уравнения для эпидемии Гани и процесса общей эпидемии Бартлетта — Мак-Кендрика соответственно:
Здесь ц,, р — коэффициенты, определяемые процессом как интенсивности взаимодействия частиц. Для первых двух процессов соответствующие уравнения были решены методом Римана [8, 9], для третьего была получена функция Ри-мана [10]. Для всех трех процессов уравнения Колмогорова могут быть сведены к одному гиперболическому уравнению, рассматриваемому ниже.
Групповые свойства уравнения. Рассмотрим гиперболическое уравнение
и для краткости обозначим его F = 0. Вычислим группу Ли, допускаемую этим дифференциальным уравнением.
Теорема 1. Группа Ли уравнения (3) обладает четырехмерным базисом из инфинитезимальных операторов.
< Пусть p1 = ux, q1 = uy, ... Рассмотрим переменные я, у, u, p1, q1, ... и ограничение операторов полных производных на (3):
или
u — u = 0,
(3)
Определяющее уравнение [11, 12]
ÔF
If (ф) |f=о = Z~— £>аф = 0,
а ÖW0
составленное для (3), имеет вид
ц
БуБх^--ф = 0.
У
Вычислим последовательно значения операторов на ф: ^ ц
Ахф = фх + р1фи + р2фр1 + фql — и;
У
Ц
БуБх ф = фху + qlфxu +—иф хр1 + q2фxq1 +
У
М* ( Ц
+ — ифи + р1 I фиу + qlфuu + — ифир1 + q2фuql
У I У
ц, I ц
--р1фР1 + p2 I фP1 y + ?1фpiu + — Мфpp + ?2фpiqi
У l y
Ц Ц | Ц ( Ц
+ 1 —7 +—q1 |ф?1 +— ul ф„1 y + ^1ф„1и +— Мф„1 Р1 + ^2ф?1А1
У2 У ) У l У
Теперь запишем новый вид уравнения (4)
Фху + qlф xu + — Мф xp1 + д2фх?1 +
У
ц f ц,
+ — ифи + p1 I фиу + qlфuu + — Мфир1 + q2фuql
У I У
ц f ц,
+ _ plфp1 + p2 I фp1 y + q19p1u + — ифp1 p1 + q2фp1q1
y l y
f ^qл
-—+-q1
У2 У
фщ + ~ и
У
(
фq1y + qlфq1u + — Uфq1p1 + q2фq1q1
У
л
(4)
-Ьф = 0. (5) У
Последовательно рассмотрим коэффициенты при независимых переменных и приравняем их нулю. При произведении переменных р^2 имеем фрж = 0, следовательно, ф распадается на сумму двух функций ф = А(х, у,и, р1 ) + Б(х,у,и, q1). При р2 : БУ (АР1) = 0, таким образом, АР1 зависит от х; при q2 : Вх (Бт ) = 0, следовательно, Б^ зависит от у; ф имеет вид ф = а(х) р1 +Р(у)q1 + С(х, у,и). Подставим в (5) и получим
Cxy ^ q1Cxu H— иаx--uCu + p1 (Cuy + qCuu ) +
У У
Я У
p1a + ß
(
ц ц
—7 и+—q1 У У
+—фУ -—(а p1 +ßq1 + C) = 0. У У
(6)
Получаем, продолжая расщепление, при р^ : Сху = 0; при рх : Сиу = 0; при ql : Схи = 0, следовательно, С = схи + у(х, у). Подставляя в (6), получим
уХу +—иах -—исх —иР+—иру -—(схи + у) = 0. (7)
у у у у у
При u = 1
х Ci_iL р+£р y _£ Ci = о.
y y y2 У У
После приведения подобных слагаемых получаем уравнение
^ х ~ Р у ,
у
дифференцируя которое по х, определяем ахх = 0 и а = с2х + с3; дифференцируя по у, после преобразования находим уравнение Эйлера
у2Руу - уРу +Р = 0,
которое имеет решение Р = с4 у + с5 у 1п у. Подставляя полученные из дифференциальных следствий результаты в исходное уравнение, запишем связь между константами
с5 =-с2. (8)
Отметим, что, собирая в уравнении (7) коэффициенты при 1, получаем
ц
У ху —у = 0, у
т. е. симметрию, задаваемую любым решением исходного уравнения и связанную с ними бесконечномерную алгебру Ли . В итоге (с учетом (8)) был получен следующий вид производящей функции:
ф = а(х)р1 +Р(у+у(х, у) = (с2х + сз )р1 +(с4у -с2у 1пу+ си + у =
= сх_ф1 + с2_ф2 + сз_фз + с4_ф4 +У, (9)
8
где фх = и, соответствующий оператор Хх = и—; ф2 = хрх - у 1п yql, соответ-
ди
д . 5
ствующий оператор Х2 = х--у 1п у —; ф3 = рх, соответствующий оператор
дх ду
8 8 8 Х3 = —; ф4 = yql, соответствующий оператор Х4 = х--у—. Инфинитези-
8х 8х 8у
мальные операторы Хх, Х2, Х3, X 4 составляют базис четырехмерной алгебры Ли Ь4, что и требовалось доказать. ►
8 8 8
Для элементов Х = £,--ь^--ь^— алгебры Ли Ь4 умножение определя-
8х 8у 8и
ется по формуле
[хг, х, ] = (* ,)-х, & +(хг (л,)-х, (л, +(хг (С,)-х, (С, ))£-.
Таблица умножения образующих элементов приведена ниже.
Таблица умножения образующих элементов
Элемент X1 X 2 X 3 X 4
X! 0 0 0 0
X 2 0 0 -X3 X 4
X 3 0 X 3 0 0
X 4 0 -X4 0 0
д
Отметим, что симметрия, соответствующая оператору х1 = и—, имеется у
ди
каждого линейного уравнения второго порядка, и поэтому в работе [12] набор образующих алгебры Ли рассмотрен с точностью до х1; в такой классификации полученная алгебра Ли является трехмерной.
Приложение. Вывод функции Римана. Метод функции Римана состоит в сведении задачи интегрирования гиперболического уравнения
д2и ди ди
I [и ] = + а (х > У + Ь (х > У + с (х' У )и = 0 (10)
дхду ох ду
к построению вспомогательной функции V, которая называется функцией Римана, если она удовлетворяет сопряженному уравнению
II [V1 = -^— а (х, у )— - Ь (х, у) — + с (х, у )v = 0 (11)
1 J дхду v п дх V /;ду v п
и следующим условиям на характеристиках:
^Т|х=хо = а(хо, у)V; ^|у=у0 = Ь(х, уо)v. (12)
ду дх
Если решена задача (11), (12), то решение задачи Коши для уравнения (10) дается явным интегральным представлением [13]. По теореме 4.1, приведенной в работе [12], если уравнение допускает четырехмерную алгебру Ли, то задача Гурса допускает одномерную инвариантную алгебру Ли и потому сводится к обыкновенному дифференциальному уравнению. Будем искать его, следуя алгоритму, изложенному в работе [14]. Из найденных в теореме 1 операторов х^ х2, х3 составим (единственно возможным образом) линейную комбинацию
х = (х -х0 у (1пу - 1пу0 , дх ду
оставляющую инвариантными характеристики, а также условия на характеристиках
Групповой метод поиска функции Римана для некоторых уравнений эпидемии V ( х0, х0; у, у0 ) = 1, V ( х, х0; у0, у0 ) = 1.
Отметим, что для второй характеристики, инвариантной относительно оператора, является ее запись в виде 1п у / у0 = 0. Критерий инвариантности
д д с
функции Р относительно оператора Х = --ь^--\-С,— состоит в удовле-
дх ду ди
творении ею уравнения ХР = 0.
В рассматриваемом случае это уравнение
дР дР
(х " х0 )-— у (1п у "1п у0 ) —=
дх ду
Решая его как уравнение в частных производных
йх йу
x - Xo y (in y - lnyo)'
получаем, что инвариантами оператора X являются первый интеграл 9 = ( x — x0) in y / y0 и v. В уравнении (3), совпадающим со своим сопряженным,
( y \
выполним замену v = V(6) = V (x-x0 )ln— и запишем уравнение типа Бес-
I Уо )
селя
ev "(е)+v '(e)-^v (е) = о.
Замена z = -02/ (4ц) приводит его к каноническому уравнению Бесселя zVZz + VZZ + zV = 0, отсюда V = J0(z), где J0(z) — функция Бесселя нулевого порядка; V (0) = J0 -ц0 ). Окончательно запишем следующее утверждение. Теорема 2 [8-10]. Уравнение эпидемии (3) обладает функцией Римана
v ( x, xo; y, yo ) = J 0 ( 2yJ-^(x - xo )ln y / yo ).
Далее выводят [8, 9] решение стационарного уравнения (3) и рассматривают предельные теоремы, относящиеся для процессов эпидемии к теоремам «порогового» типа, которые применяют для определения пороговой численности инфицированных частиц, превышение которого означает начало эпидемии [15]. Отметим, что граничные условия для уравнения Бартлетта — Мак-Кендрика не являются условиями на характеристиках, поэтому приведенный метод Римана к нему неприменим.
Заключение. Предложенный в работе групповой подход к исследованию уравнений на двойную производящую функцию процесса дает новый взгляд на аналитические методы и результаты теории эпидемий. Применение полученной в теореме 1 алгебры Ли хорошо вписывается в сложившийся набор используемых для исследования процессов эпидемии методов, не только подтверждает уже полученные результаты для процессов Вейса и Гани, но и позволяет вы-
явить причину возникших трудностей исследования процесса Бартлетта — Мак-Кендрика. В следующей работе с групповой точки зрения будет рассмотрен метод разделения переменных для системы первого и второго уравнений Колмогорова процессов эпидемии.
ЛИТЕРАТУРА
1. Гихман И.И., Скороход А.В. Введение в теорию случайных процессов. М.: Наука, 1977. 568 с.
2. Калинкин А.В. Марковские ветвящиеся процессы с взаимодействием // Успехи математических наук. 2002. Т. 57. № 2. С. 23-84.
3. Севастьянов Б.А. Ветвящиеся процессы. М.: Наука, 1971. 436 с.
4. Эпидемии процесс // Математическая энциклопедия. Т. 5. М.: Советская энциклопедия, 1985. 623 с.
5. Weiss G. On the spread of epidemics by carries // Biometrics. 1965. Vol. 21. No. 2. P. 481-490. DOI: 10.2307/2528105 URL: http://www.jstor.org/stable/2528105
6. Gani J. Approaches to the modelling of AIDS // Stochastic processes in epidemic theory. Springer, 1990. Pp. 145-154.
7. Bartlett M.S. Some evolutionary stochastic processes // J. of Royal Statistical Society. Ser. B (Methodological). 1949. Vol. 11. No. 2. P. 211-229.
8. Калинкин А.В. Финальные вероятности ветвящегося процесса с взаимодействием частиц и процесс эпидемии // Теория вероятностей и ее применение. 1998. Т. 43. № 4. С. 773-780.
9. Мастихин А.В. Финальное распределение для марковского процесса эпидемии Гани // Математические заметки. 2007. Т. 82. № 6. С. 873-884.
10. Мастихин А.В. Функция Римана для некоторых уравнений Колмогорова // Инженерный вестник. 2014. № 12. URL: http://engsi.ru/doc/745888.html
11. Виноградов А.М., Красильщик И.С., ред. Симметрии и законы сохранения уравнений математической физики. М.: Факториал Пресс, 2005. 379 c.
12. Овсянников Л.В. Групповой анализ дифференциальных уравнений. М.: Наука, 1978. 339 c.
13. Бицадзе А.В., Калиниченко Д.Ф. Сборник задач по уравнениям математической физики. М.: Наука, 1985. 312 с.
14. Ибрагимов Н.Х. Опыт группового анализа обыкновенных дифференциальных уравнений. М.: Знание, 1991. 48 c.
15. Kalinkin A.V., Mastikhin A.V. A limit theorem for a Weiss epidemic // J. Appl. Probab. 2015. Vol. 52. No. 1. P. 247-257.
Мастихин Антон Вячеславович — канд. физ.-мат. наук, доцент кафедры «Высшая математика» МГТУ им. Н.Э. Баумана (Российская Федерация, 105005, Москва, 2-я Бауманская ул., д. 5, стр. 1).
Шевченко Маргарита Николаевна — старший преподаватель кафедры «Высшая математика» МГТУ им. Н.Э. Баумана (Российская Федерация, 105005, Москва, 2-я Бауманская ул., д. 5, стр. 1).
Просьба ссылаться на эту статью следующим образом:
Мастихин А.В., Шевченко М.Н. Групповой метод поиска функции Римана для некоторых уравнений эпидемии // Вестник МГТУ им. Н.Э. Баумана. Сер. Естественные науки. 2018. № 2. С. 12-22. БО!: 10.18698/1812-3368-2018-2-12-22
GROUP METHOD IN SEARCHING RIEMANN FUNCTION FOR SOME EPIDEMIC EQUATIONS
A.V. Mastikhin M.N. Shevchenko
Bauman Moscow State Technical University, Moscow, Russian Federation
Abstract
In the stochastic theory of epidemics, the general stochastic epidemic (i. e. Bartlett — Mac-Kendric epidemic process) and the simple stochastic epidemic (i. e. Weiss epidemic process) are the basic models; both of them are considered on the set of states N2.In Bartlett — Mac-Kendric epidemic process an ineracting pair of infected and susceptible particles transform into a pair of two infected particles (or two carriers). This process is rather complicated; some results were obtained by asymptotic methods. In Weiss epidemic process an interacting pair of infected and susceptible particles transform into one carrier, i. e. one infected individual is removed from the population. This Markov process is more accessible to study; there are many generalizations of Weiss model on the case of Nn. For instance, J. Gany introduced a carrier borne epidemic process with two stages of infection on N3 which may be of some relevance to the spread of AIDS. We consider the first stationary Kolmogorov equation for exponential (double) generating function of transitions probabilities of Markov epidemic process. For all three processes considered above we have the same stationary Kolmogorov equation. It is possible to solve it by using Riemann function as the solution of a characteristic value problem. In present paper we study this hyperbolic equation by the symmetry method. The 4-dimentional Li algebra was obtained in the theorem 1. It contains a subalgebra, which is invariant due to the characteristic problem. So, we can reduce the hyperbolic equation to Bessel differential equation. Then, it is easy to construct the solution in an integral form and to obtain limit theorems. We also discuss the possibility of using Riemann method in the case of general epidemic process
Keywords
Infinitesimal symmetry generator, Markov process, exponential (double) generating function, first Kolmogorov equation, Riemann function
Received 16.02.2017 © BMSTU, 2018
REFERENCES
[1] Gikhman I.I., Skorokhod A.V. Vvedenie v teoriyu sluchaynykh protsessov [Introduction to the random processes theory]. Moscow, Nauka Publ., 1977. 568 p.
[2] Kalinkin A.V. Markov branching processes with interaction. Russian Mathematical Surveys, 2002, vol. 57, no. 2, pp. 241-304. DOI: 10.1070/RM2002v057n02ABEH000496
[3] Sevast'yanov B.A. Vetvyashchiesya protsessy [Branching processes]. Moscow, Nauka Publ., 1971. 436 p.
[4] Epidemii protsess. Matematicheskaya entsiklopediya. T. 5 [Epidemic process. In: Mathematical encyclopedia. Vol. 5]. Moscow, Sovetskaya entsiklopediya Publ., 1985. 623 p.
[5] Weiss G. On the spread of epidemics by carries. Biometrics, 1965, vol. 21, no. 2, pp. 481-490. DOI: 10.2307/2528105 Available at: http://www.jstor.org/stable/2528105
[6] Gani J. Approaches to the modelling of AIDS. Stochastic processes in epidemic theory. Springer, 1990. Pp. 145-154.
[7] Bartlett M.S. Some evolutionary stochastic processes. J. of Royal Statistical Society. Ser. B (Methodological), 1949, vol. 11, no. 2, pp. 211-229.
[8] Kalinkin A.V. Final probabilities for a branching process with interaction of particles and an epidemic process. Theory of Probability and its Applications, 1999, vol. 43, no. 4, pp. 633-640. DOI: 10.1137/S0040585X97977203
[9] Mastikhin A.V. Final distribution for Gani epidemic Markov processes. Mathematical Notes, 2007, vol. 82, iss. 5-6, pp. 787-797. DOI: 10.1134/S0001434607110223
[10] Mastikhin A.V. Riemann functions for some Kolmogorov equations. Inzhenernyy vestnik [Engineering Bulletin], 2014, no. 12 (in Russ.). Available at: http://engsi.ru/doc/745888.html
[11] Vinogradov A.M., Krasil'shchik I.S., ed. Simmetrii i zakony sokhraneniya uravneniy matematicheskoy fiziki [Symmetry and energy conservation laws of equations of mathematical physics]. Moscow, Faktorial Press Publ., 2005. 379 p.
[12] Ovsyannikov L.V. Gruppovoy analiz differentsial'nykh uravneniy [Group analysis of differential equations]. Moscow, Nauka Publ., 1978. 339 p.
[13] Bitsadze A.V., Kalinichenko D.F. Sbornik zadach po uravneniyam matematicheskoy fiziki [Problem book of mathematical physics equations]. Moscow, Nauka Publ., 1985. 312 p.
[14] Ibragimov N.Kh. Opyt gruppovogo analiza obyknovennykh differentsial'nykh uravneniy [Group analysis experience of ordinary differential equations]. Moscow, Znanie Publ., 1991. 48 p.
[15] Kalinkin A.V., Mastikhin A.V. A limit theorem for a Weiss epidemic. J. Appl. Probab., 2015, vol. 52, no. 1, pp. 247-257.
Mastikhin A.V. — Cand. Sc. (Phys.-Math.), Assoc. Professor, Department of Higher Mathematics, Bauman Moscow State Technical University (2-ya Baumanskaya ul. 5, str. 1, Moscow, 105005 Russian Federation).
Shevchenko M.N. — Assist. Professor, Department of Higher Mathematics, Bauman Moscow State Technical University (2-ya Baumanskaya ul. 5, str. 1, Moscow, 105005 Russian Federation).
Please cite this article in English as:
Mastikhin A.V., Shevchenko M.N. Group Method in Searching Riemann Function for some Epidemic Equations. Vestn. Mosk. Gos. Tekh. Univ. im. N.E. Baumana, Estestv. Nauki [Herald of the Bauman Moscow State Tech. Univ., Nat. Sci.], 2018, no. 2, pp. 12-22 (in Russ.). DOI: 10.18698/1812-3368-2018-2-12-22
