УДК 519.213
DOI: 10.14529/mmph170303
ИССЛЕДОВАНИЕ ВЕРОЯТНОСТНОГО РАСПРЕДЕЛЕНИЯ ТИПА ГИПЕРБОЛИЧЕСКОГО КОСИНУСА
М.С. Токмачев
Новгородский государственный университет им. Ярослава Мудрого, г. Великий Новгород, Российская Федерация E-mail: [email protected]
Исследовано полученное в результате характеризации трехпараметри-ческое вероятностное распределение, являющееся обобщением известных распределений: однопараметрического распределения гиперболического косинуса (секанса) и двухпараметрического распределения Майкснера. Приведено доказательство его безграничной делимости. По характеристической функции в общем виде восстановлена плотность распределения вероятностей, выраженная через бета-функцию комплексно-сопряженных аргументов. Наряду с единой формулой, при целых значениях параметра т для функции плотности выведены соотношения в элементарных функциях.
Ключевые слова: распределение типа гиперболического косинуса; характеристическая функция; безгранично делимое распределение; бета-функция.
Введение
Вероятностным распределением типа гиперболического косинуса называют трехпараметри-ческое распределение с характеристической функцией
f(t) = ich ßt - i msh ßt\ , где m ß, m e R, m > 0, ßФ 0. (1)
v m ß m )
Впервые функция (1) получена автором при решении задачи характеризации распределений свойством постоянства регрессии квадратичной статистики Q на линейную статистику Л [1]:
Е(01 Л) = Е(0, где 0 = XIа]кХ]Хк; Л = £Х} . (2)
]=1 к=1 ]=1
Х1,Х2,...,Хп - независимые, одинаково распределенные случайные величины. В зависимости от соотношения коэффициентов статистики б условие постоянства регрессии (2) является характе-ризационным (характеристическим) для ряда известных распределений [2-4] и, в частности, распределения типа гиперболического косинуса.
Указанное трехпараметрическое распределение является обобщением двухпараметрического распределения Майкснера (Ж. Ме1хпег) [5, 6] с характеристической функцией
/(г) = ( Л г - iвsh г )~т , т > 0, ве R.
Распределение Майкснера получается, если в (1) положить //= т и переобозначить параметр т1Р=т т = в. Наличие третьего параметра позволяет получить большее разнообразие распределений: каждое двухпараметрическое распределение Майкснера с конкретными параметрами (в, т) является семейством распределений с параметром // . Например, при в = 2, т = 1 получаем тройку параметров (т = 2/, т = 1, /), где /е R, /ф 0. В частности, эти параметры (т, т, /3) могут быть (2; 1; 1) - распределение Майкснера, а также (1; 1; 0,5), (-1; 1; -0,5), (4; 1; 2) и т. д. Заметим, что в трехпараметрическом распределении математическое ожидание случайной величины равно параметру т. При этом параметр т является параметром не только сдвига, но и участвует в формировании дисперсии и других моментов.
В частности, при т = 0, т = 1 из (1) получаем характеристическую функцию однопараметрического распределения гиперболического косинуса (секанса). При сдвиге случайной величины, распределенной по закону гиперболического косинуса, на величину Л получим двухпараметрическое распределение, известное как распределение Чампернауна [7].
Таким образом, распределение типа гиперболического косинуса оказывается обобщением двухпараметрического распределения Майкснера, которое, в свою очередь, является обобщением однопараметрического распределения гиперболического косинуса на случай /Ф 0.
В данной работе восполняются пробелы в доказательной части теории: представлено обоснование найденного ранее распределения и выведен ряд соотношений.
Исследование характеристической функции
Докажем, что функция / ^) вида (1), полученная как результат характеризации распределений свойством постоянства регрессии, действительно является характеристической функцией. Теорема 1. Функция
1 1 (3)
q(t)
, ß — , ß '
ch — t -i — sh — t
т Ь т
где л, Ь, т - действительные постоянные, Ь, т Ф 0 , является характеристической функцией некоторого распределения.
Доказательство теоремы 1. Для доказательства теоремы потребуются две леммы. Лемма 1. Функция
q(t) = chßt - i— shßt, m ß m
(4)
где /, Ь, т - действительные константы, р, т Ф 0 , является целой функцией порядка 1.
Доказательство леммы 1. Для доказательства леммы 1 достаточно проверить условия Ко-ши-Римана применительно к функции комплексной переменной д(2) вида (4), которые справедливы на всей комплексной плоскости. Порядок целой функции д(2) определяют функции еЬ Ь 2
т
и Ь 2, которые являются целыми функциями первого порядка.
т
Лемма 2. Задана функция
q(z) = ch ß z - i— sh ß z, m ß m
где — , ß, m - действительные константы, ß, m Ф 0 . Тогда а) q( z) имеет лишь чисто мнимые нули; б) все нули функции q(z) простые; в) z = 0 не является нулем функции q(z).
Доказательство леммы 2. Положим, q(z) = 0 . Решая указанное уравнение с использованием
формулы cth z = i ctg iz, получаем ctg iß z = — . Отсюда находим нули функции q(z):
m ß
im
f
zь. =--
ß
m
arcctg —+ pk g ß .
\
k e Z.
(5)
Из (5) следует справедливость утверждений а) и б) леммы 2. Утверждение в) проверяется непосредственно. Лемма 2 доказана.
Переходим к доказательству теоремы 1. Для целой функции д(2) вида (4) применим теорему Вейерштрасса о факторизации целых функций [8]:
q( z) = zS0 e
Soa z)
П
к=-
1 --L
, zk;
z!zk
mk
Упростим данное выражение, используя утверждения леммы 1: так как 2 = 0 не является нулем функции д(2), то = 0; так как все корни 2к простые, то все тк равны 1, к е Ъ . Следовательно,
+¥ (
q(z) = eh(z> П
к=-¥ V
1 --L
zk;
z!zk
e
e
Исходя из (5), обозначим zk =--, k е Z, где Бк - действительные и различные при раз-
B
ных k . Тогда
q(z) = eh(z) Д(1 -iBkZ).
,iBkz
(6)
k =-¥
Поскольку, согласно лемме 1, д(2) - целая функция первого порядка, то Т](2) - многочлен степени не выше первой, а именно: Т](2) = аг + Ь . Из условия д(0) = 1 и из (6) следует, что Ь = 0 . Таким образом, (6) принимает вид
д(2) = еа П(1 - ТБк2) етБк2. к=-¥
Перейдем к действительной переменной (, тогда
Ч($) = еа? П(1 - Бк()
к=-¥
Для определения а продифференцируем функцию ), заданную соотношением (7):
( +¥
q'(t) = aeat Д (1 - Bkt)
¿Bkt
V k =-¥
( +¥
+ e
П(! - iBkt)
JBkt
V k=
(7)
(8)
Из (8) следует д'(0) = а + ТБ, где Б - действительное число. Однако, из (4) получаем
д'(0) = —тЛ. Два последних соотношения приводят к выводу, что а - чисто мнимое (возможно, т
а = 0). Положим, а = ТА , где А - действительная постоянная. Подставив это значение в (7), приходим к выражению
q(t) = еш П (1 - iBkt)
Тогда
1
k=-¥ +¥ (
■ = е
—iAt
П
1
\
V 1 - iBkt J
iBkt
-iBkt
q(t) k=-¥
Функции-сомножители в правой части (9) являются характеристическими функциями:
e~At, е" 1
(9)
iAt " iBkt - характеристические функции вырожденного распределения,
1 - iBkt
- характеристическая функция гамма-распределения при Bk > 0 и характеристиче-
ская функция распределения, сопряженного с гамма-распределением при Бк < 0 . Введем функцию Нп ^) следующим образом
К (t) = e-At П
k=-n
1
v 1 - iBkt J
,~iBkt
(10)
Согласно известному свойству произведения характеристических функций, функция Нп (^) при любом конечном п является характеристической. Из (9), (10) следует, что
1
q(t)
= lim hn (t).
(11)
По следствию из теоремы непрерывности характеристических функций ([9], следствие 2 тео-
ремы 3.6.1) предельная функция
1
q(t)
также является характеристической функцией. Таким об-
разом, теорема 1 доказана. Отметим, что соотношения (10), (11) представляют структуру характеристических функций вида (3).
e
Исследуем найденное распределение на безграничную делимость. Теорема 2. Функция
1 1
^ оЬ / - г /
т / т
является безгранично делимой характеристической функцией. Доказательство теоремы 2. Согласно теореме 1, функция
1
q(t)
является характеристиче-
ской функцией. Обратимся к доказательству теоремы 1.
Из (10) следует, что характеристическая функция Нп (^) является произведением конечного числа безгранично делимых характеристических функций гамма-распределения и вырожденного распределения. Следовательно ([9], теорема 5.3.2), Нп ^) также безгранично делима.
Характеристическая функция
1
q(t)
безгранично делима, поскольку, согласно (11), является
пределом последовательности безгранично делимых характеристических функций Нп (^) ([9], теорема 5.3.3).
Теорема 2 доказана.
Теорема 3. Функция f (t) =
f
оЬ — I - г — sh — I I при любом т > 0 является безгранично ч т / т )
делимой характеристической функцией.
Доказательство теоремы 3. Поскольку любая положительная степень безгранично делимой характеристической функции сама является характеристической функцией ([9], следствие
1 Г 1 Г
теоремы 5.3.3), а функция- безгранично делима, согласно теореме 2, то /(^) = - , т > 0
д(0 у ))
также является характеристической функцией, причем безгранично делимой. Теорема 3 доказана.
Вывод соотношений для плотности распределения вероятностей
Плотность распределения рт (х) как обратное преобразование Фурье характеристической функции /(^) вида (1) при /> 0 имеет вид [10]:
1 +Т 1
(Х) = 2p ^ e"tXf (t) dt = i
(ch ßt - i ßt)m
imx
2m-2 mßm-1 ( ß- iß\2p ( m imx m imx^
B
т
где г =
, В(р; д) - бета-функция, а сомножитель с мнимой единицей имеет вид:
ß + iß,
-; — + 2 2ß 2 2ß
A °
ß- iß ß + iß,
e
2ßß arctg 2-2
ß2 -ß2
2ßß
arctg—---+psign ß
, ß -ß
-sign ß
при ß2-ß2 > 0, при ß2-ß2 < 0, при ß2-ß2 = 0.
При целых т функцию плотности в (12) можно выразить в элементарных функциях. Теорема 4. При т = 1 функция плотности рт (х) имеет вид:
(12)
(13)
m
—itx
e
m
e
л( x) =
A 2b
2J b2+m2 ch— УИ И 2Ь
(14)
Справедливость (14) следует из (12) при использовании известных соотношений:
B|1 -z;2 + z I = -
2 2 J cos—z
cos iz = ch z .
Теорема 5. При m = 2 функция плотности pm (x) имеет вид:
x
2 xAb
p2(x)="
(b2 +m2)sh
—x b
(15)
Справедливость (15) следует из (12) при использовании известных соотношений
B (1 + z; 1 - z )=-
—z
sin iz = ish z.
sin —z
Также отметим, что при целых т случайная величина X с характеристической функцией вида (1) представляет собой сумму т независимых случайных величин X1 + X2 +... + Xm , каждая из которых обладает характеристической функцией
f (t) =
f
v
ch —t - i msh —t m b m J
\—1
Следовательно, функция p2(x) вида (15) как плотность вероятностей распределения суммы независимых одинаково распределенных случайных величин X1,X2 является сверткой:
P2 (x) = | Pl (у) Pl (x — У^у , где д (у) = pl у Теорема 6. При т = 3, 5, 7,... функция плотности pm (x) имеет вид:
mx
A 23
m b 2 , 2
pm (x) ="
m
m
# . ..2V2",
m-1 2
. —mx
ch- n=!
П [(2n - 1)2b
2 o2 , 2 2
" + m x
(16)
2(Ь +m2)2(m-1)! - 2b
Доказательство теоремы 6. Итак, полагаем m = 2k +1 (k = 1, 2,...). Тогда
B
f m imx m imx ^ ; — +
v
2 2b 2 2b
f
= B
1 imx
1 imx
k +---; k + — +
2 2b 2 2b
Г
( (
k +
V v
1 imx
2 - 2b.
Г
(( k +
1 imx
- + -
Л
V
2 2b JJ r(k + z) Г( k + z)
r(2k +1)
(2k)!
-¡-у \ , 1 imx _
где Г (...) - гамма-функция, z = , z и z - комплексно-сопряженные величины.
По свойствам гамма-функции
r(k + z) = (k -1 + z)(k - 2 + z) ... zr(z), Г | 1 -©] Г | 1 + w| = ——
V 2 J V 2 J cos— w
получаем
r(k + z)r(k + z) =
(1 imx V 3 imx ^ ( 2k -1 imx
2 - 2b Jl 2 - 2b J... V ~2 2b.
x
x
X
( 1 imx Л ( 3 imx Л (2k -1 imx Л
— +- - +- ... -+-
2 2b А 2 2b J ( 2 2b
Г
(1 imx 1
imx
Л
, Г - + ■
2 2b J {2 2b
2 2 Л
1 m x — +
4 4b2 J! 4 4b
9 m x — +
2~2 Л ((2k -1)2 m2x2 Л
+
4b2
p
J cos
ipmx
Ж
В последнем выражении перейдем к функции гиперболического косинуса, cos iw = chw, за-
меним значение k на
m -1 2
и вынесем общие множители слагаемых в скобках. Тогда
B
m imx ; m imx Л _ p(fi + mx)(9b + mx)...((m - 2y b + mx)
2 2b '2 2b J ^m-xom-u ,s. . Pmx
(17)
2m-1b (m - 1)!ch
2b
Следовательно, исходя из выражения (12) и соотношений (13), (17), для функции плотности вероятностей при т = 3,5, 7,... приходим к формуле (16). Теорема 6 доказана.
Теорема 7. При т = 4, 6, 8,... функция плотности рт (х) имеет вид:
m 1
mx m
Pm (x) _-
m2 x
A 2b
m
?2 . "2\2
. pmx sh- "_!
П(4n2b2 + m2x2).
(18)
2(b2 + m )2 (m-1)^- 2b Доказательство теоремы 7. Полагаем m _ 2k (k _ 2, 4,...). Тогда
-, ikx ikx ^ k--; k + —
. b b
в
f m imx m imx Л ^
---; — +-I _ B
2 2b 2 2b
Г
k -
ikx b
v
Л ( Г
k +
ikx
\
m
b J Г(k - z)Г(k + z) ikx , ...
-_-, где z _ —, k _ — .
(2k -1)! b 2
H2k)
Вновь используем свойства гамма-функции: Щ + z) _ (k - 1 + z)(k - 2 + z) ... (1 + z)Г(1 + z), Щ - z) _ (k -1 - z)(k - 2 - z) ... (1 - z)H1 - z)
Г(1 - z )Г(1 + z )_-
pz
sin pz
и получаем Г(
(k - z )Г( k + z)_((k -1)2 - z2 )((k - 2 )2 - z2) ... (12 - z2 )г(1 - z )Г(1 + z)_
(,-z2)(4-z2)...((k-1)2 -z2)-pz
sin pz
2 2 Л
1 +
k 2 x b2
22
4 +
k 2 x b2
22
k2 x
(k -1)2 + 42 b2
p-
ikx b
J sin p-
ikx b
В последнем выражении перейдем к функции гиперболического синуса, sin iw = i sh w, и вынесем общие множители слагаемых в скобках. Тогда
B
( m imx m + imx ) _ (02 + k2x2)(4p2 + k2x2)...((k - 1)2b2 + k2x2)pkx
2 2b 2 2b
b2k-1(2k - 1)!sh
pkx b
(19)
Таким образом, согласно соотношениям (13), (19), для функции плотности вероятностей, заданной (12), при т = 4, 6, 8,... приходим к формуле (18). Теорема 7 доказана. Можно доказать, что в общем случае
г^ч 2 b2 +m2 2m 2
e(X)_m; _-——; g _—; g _— +
m ms m
(2m Л2 2
v ms J
_ - + J?,
m
4
1
где s2 = V(X), g = As(X), g2 = Ex(X). Исходя из указанных соотношений, легко найти зависимость параметров ß, ß , m от первых моментов:
m = E(X); m = -2ß = ^-; ß2 = ms2 -ß2. gS g -g
В статистическом анализе данных найденные соотношения позволяют использовать метод моментов. При обработке реальных данных из области медицины и здравоохранения (различные показатели физического состояния и показатели заболеваемости населения) с помощью разработанного программного обеспечения установлено согласие многих из них с распределением типа гиперболического косинуса [11].
Заключение
Можно найти примеры применения указанного распределения кроме вероятностных. В частности, в [12, 13] представлено множество нетривиальных интегралов, вычисляемых на основе моментов распределения при различных параметрах. В [14] также из найденных взаимосвязи и структуры моментов распределения типа гиперболического косинуса сформировано структурированное множество в виде бесконечной числовой призмы, исследован ряд её сечений, в частности, в связи с числами Стирлинга первого рода и коэффициентами полиномов Бесселя. В [15] в качестве сечений числовой призмы представлены и систематизированы как широко известные классические, так и числовые треугольники, и числовые последовательности, ранее в литературе не встречавшиеся. Для них найдены многие интересные свойства и соотношения.
Работа выполнена при финансовой поддержке проектной части государственного задания в сфере научной активности Министерства образования и науки Российской Федерации, проект № 1.949.2014/К.
Литература
1. Токмачев, М.С. Характеризация распределения типа гиперболического косинуса свойством постоянства регрессии / М.С. Токмачев // Деп. в ВИНИТИ 21.06.94. - № 1542-В94. - 11 с.
2. Токмачев, М.С. Постоянство регрессии квадратичной статистики на линейную статистику / М.С. Токмачев // Вестник НовГУ. - 1995. - № 1. - С. 139-141.
3. Каган, А.М. Характеризационные задачи математической статистики / А.М. Каган, Ю.В. Линник, С Р. Рао. - М.: Наука, 1972. - 656 с.
4. Клебанов, Л.Б. Когда квадратичная статистика имеет постоянную регрессию на выборочное среднее? / Л.Б. Клебанов // Теория вероятн. и её примен. - 1979. - ХХ1У, 3 - C. 646-648.
5. Lai, C.D. Meixner classes and Meixner hypergeometric distributions / C.D. Lai // Australian & New Zealand Journal of Statistics. - 1982. - Vol. 24. - P. 221-233.
6. Lai, C.D. A characterization of gamma, Meixner hypergeometric and negative binomial distributions based on canonical measures / C.D. Lai // Annals of the Institute of Statistical Mathematics. -1982. - Vol. 34. - Issue 1. - P. 359-363.
7. Вадзинский, Р.Н. Справочник по вероятностным распределениям / Р.Н. Вадзинский. -СПб.: Наука, 2001. - 295 с.
8. Смирнов, В.И. Курс высшей математики / В.И. Смирнов - M.: Наука, 1974. - Т. 3, Ч. 2. -672 с.
9. Лукач, Е. Характеристические функции / Е. Лукач. - М.: Наука, 1979. - 424 с.
10. Токмачев, М.С. Распределение типа гиперболического косинуса / М.С. Токмачев, А.М. Токмачев // Вестник НовГУ. - 2001. - № 17. - С. 85-88.
11. Токмачев, М.С. Прикладной аспект обобщенного распределения гиперболического косинуса / М.С. Токмачев // Вестник НовГУ. Сер.: Техн. науки. - 2005. - № 34. - С. 96-99.
12. Токмачев, М.С. Некоторые интегралы, связанные с распределением типа гиперболического косинуса / М.С. Токмачев // Математика в вузе и в школе: Труды XXIV Международной научно-методической конференции. - СПб.: Петербургский гос. ун-т путей сообщения, 2012. -С. 185-186.
13. Токмачев, М.С. Вычисление интегралов от функций некоторого класса с вероятностной интерпретацией / М.С. Токмачев // Вестник НовГУ. Сер.: Физико-математические науки. - 2014. - № 80. - С. 42-46.
14. Токмачев, М.С. О числовых множествах и последовательностях в связи с распределением типа гиперболического косинуса / М.С. Токмачев // Вестник НовГУ. Сер.: Физико-математические науки. - 2015. - № 3-2(86). - С. 35-39.
15. Токмачев, М.С. Множество подмножеств в структуре некоторой числовой призмы: монография / М.С. Токмачев // Деп. В ВИНИТИ 09.06.2016. - № 91-В2016. - 90 с.
Поступила в редакцию 30 ноября 2016 г.
Bulletin of the South Ural State University Series "Mathematics. Mechanics. Physics" _2017, vol. 9, no. 3, pp. 18-26
DOI: 10.14529/mmph170303
THE STUDY OF PROBABILITY DISTRIBUTIONS OF THE HYPERBOLIC COSINE TYPE
M.S. Tokmachev
Yaroslav-the-Wise Novgorod State University, Veliky Novgorod, Russian Federation E-mail: [email protected]
Probability distribution, obtained earlier by the author in a result of the characterization of distributions by the property of the constant regression of quadratic statistics on a linear form and called the hyperbolic cosine type distribution, is considered. Along with the classic normal distribution, gamma distribution, negative binomial and some other distributions, the given three-parameter distribution is related to the class of probability distributions, united by a common characterization condition. The obtained distribution is a generalization of the famous two-parameter Meixner distribution.
In the article, a proof that a function appeared in the result of given characterization is indeed characteristic is given. The structure of this function, in connection with gamma distribution and the corresponding distribution conjugated with the gamma distribution, along with the constant distribution, is presented. Infinite decomposability of distribution is proved.
Based on characteristic function, the probability density function, expressed through the beta function of complex conjugate arguments is recovered in general terms. Along with a unified formula, correlations in elementary functions are deduced for the distribution density functions at integer values of parameter m. Density formulas at odd and even values of the parameter are similar on the structure of cofactors: exponent function, hyperbolic secant or cosecant correspondingly, and polynomial factors.
The distribution under study has multiple applications, not only in probability problems. Moments of distribution at specified parameters change as polynomials of some class with corresponding coefficients. These coefficients can be constructed as number sets (number triangles and number sequences), both known and new with setting of a row of functional relationships.
Keywords: distribution of the hyperbolic cosine type; characteristic function; infinitely divisible distribution; beta function.
References
1. Tokmachev M.S. Dep. in VINITI21.06.94, no. 1542^94, 11 p. (in Russ.).
2. Tokmachev M.S. Postoyanstvo regressii kvadratichnoy statistiki na lineynuyu statistiku (Constancy of regression of quadratic statistics with linear statistics). Bulletin of the Novgorod state University, 1995, no. 1, pp. 139-141 (in Russ.).
3. Kagan A.M., Linnick Y.V., Rao S.R. Kharakterizatsionnye zadachi matematicheskoj statistiki (Characterization problems of mathematical statistics), Moscow, Nauka Publ., 1972, 656 p. (in Russ.).
4. Klebanov L.B. Kogda kvadratichnaya statistika imeetpostoyannuyu regressiyu na vyborochnoe srednee? (When the quadratic statistic has constant regression on the sample mean?). Teoriya veroyat-nostey i eye primenenie (Theory of probability and its applications), 1979, XXIV, 3. pp. 646-648. (in Russ).
5. Lai C.D. Meixner classes and Meixner hypergeometric distributions. Aust. J. Stat., 1982, Vol. 24, pp. 221-233. DOI: 10.1111/j.1467-842X.1982.tb00828.x
6. Lai C.D. A characterization of gamma, Meixner hypergeometric and negative binomial distributions based on canonical measures. Ann. Inst. Stat. Math., 1982, Vol. 34, Issue 1, pp. 359-363. DOI: 10.1007/BF02481035
7. Vadzinskiy R.N. Spravochnik po veroyatnostnym raspredeleniyam (Handbook of Probabilistic Distributions), St. Petersburg, Nauka Publ., 2001, 295 p. (in Russ.).
8. Smirnov V.I. Kurs vysshey matematiki [Course of Higher Mathematics], Moscow, Nauka Publ., 1974, Vol. 3, part 2, 672 p. (in Russ.).
9. Lukacs E. Characteristic functions. Hafner Publishing Co., New York, 1970, 350 p.
10. Tokmachev M.S., Tokmachev A.M. Raspredelenie tipa giperbolicheskogo kosinusa (Distribution type hyperbolic cosine). Bulletin of the Novgorod state University, 2001, no. 17, pp. 8588. (in Russ.).
11. Tokmachev M.S. Prikladnoj aspekt obobshhennogo raspredelenija giperbolicheskogo kosinusa (Applied aspect of the generalized hyperbolic cosine distribution). Bulletin of the Novgorod state University, 2005, no. 34, p. 96-99 (in Russ.).
12. Tokmachev M.S. Nekotorye integraly, svyazannye s raspredeleniem tipa giperbolicheskogo kosinusa (Some integrals associated with the distribution type hyperbolic cosine). Matematika v vuze i v shkole: Trudy XXIVMezhdunarodnoy nauchno-metodicheskoy konferentsii (Mathematics in high school and in school: Proceedings of the XXIV International Scientific and Methodical Conference), St. Petersburg, Peterburgskiy gos. universitet putey soobshcheniya Publ., 2012, pp. 185-186. (in Russ.).
13. Tokmachev M.S. Vychislenie integralov ot funkcij nekotorogo klassa s verojatnostnoj in-terpretaciej (Evaluation of integrals for some class of functions with probabilistic interpretation). Vestnik NovGU. Ser.: Fiziko-matematicheskie nauki (Bulletin of the Novgorod state University. Ser.: Physical and mathematical Sciences), 2014, no. 80, pp. 42-46 (in Russ.).
14. Tokmachev M.S. O chislovyh mnozhestvah i posledovatel'nostyah v svyazi s raspredeleniem tipa giperbolicheskogo kosinusa (About number sets and sequences in connection with a distribution of the hyperbolic cosine type). Vestnik NovGU. Ser.: Fiziko-matematicheskie nauki (Bulletin of the Novgorod state University. Ser.: Physical and mathematical Sciences), 2015, no. 3-2 (86), pp. 35-39. (in Russ).
15. Tokmachev M.S. Dep. in VINITI 09. 06. 2016, no. 91^2016, 90 p. (in Russ.).
Received November 30, 2016