ВЕСТНИК ТОМСКОГО ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА
2016 Управление, вычислительная техника и информатика № 2 (35)
УДК 519.2
Б01: 10.17223/19988605/35/4
К.И. Лившиц, Л.Ю. Сухотина
ВЕРОЯТНОСТЬ РАЗОРЕНИЯ СТРАХОВОЙ КОМПАНИИ ПРИ ГИПЕРЭКСПОНЕНЦИАЛЬНЫХ РАСПРЕДЕЛЕНИЯХ СТРАХОВЫХ ПРЕМИЙ И СТРАХОВЫХ ВЫПЛАТ ДЛЯ РАЗЛИЧНЫХ МОДЕЛЕЙ СТРАХОВАНИЯ
Находятся вероятности разорения страховой компании и среднее условное время до разорения для модели Крамера-Лундберга со стохастическими премиями, модели Крамера-Лундберга с ММР-потоком страховых выплат при гиперэкспоненциальных распределениях страховых премий и страховых выплат.
Ключевые слова: вероятность разорения; среднее условное время до разорения; модель Крамера-Лундберга; модель Крамера-Лундберга со стохастическими премиями; гиперэкспоненциальное распределение.
Стандартной задачей актуарной математики является задача вычисления вероятности разорения страховой компании, т.е. вероятности ситуации, когда страховая компания не может исполнять свои финансовые обязательства ввиду отсутствия денежных средств при различных предположениях о потоках, поступающих в компанию страховых премий и страховых выплат, производимых страховой компанией [1-3]. Сложность состоит, как правило, в нахождении явных решений соответствующих систем интегро-дифференциальных уравнений. Поэтому представляет интерес определение явного вида для вероятности разорения для какого-то достаточно широкого класса распределений страховых премий и страховых выплат. В настоящей работе находятся вероятности разорения и среднее условное время до разорения для некоторых стандартных моделей деятельности страховой компании в предположении, что страховые премии и страховые выплаты имеют гиперэкспоненциальное распределение. Аналогичная идея используется, например, в задачах управления запасами, где считается, что объемы потребляемых ресурсов имеют гиперэкспоненциальное распределение [4].
1. Модель Крамера-Лундберга
Простейшей моделью деятельности страховой компании является модель Крамера-Лундберга, которая строится при следующих предположениях [1-3]: процесс поступления страховых премий в компанию считается детерминированным, за время t приращение капитала компании равно С, где С - количество средств, поступивших в компанию за единицу времени; страховые выплаты - независимые, одинаково распределенные случайные величины с плотностью распределения ¥(х) и средним значением а ; моменты наступления страховых выплат образуют пуассоновский поток интенсивности X . Поток страховых выплат не зависит от поступления страховых премий. Пусть в момент времени t капитал компании равен S^). Определим Т = шт^ :5^) < 0} и Т = м, если S^) > 0 Vt. Случайная
величина Т - момент разорения. Тогда вероятность предельного разорения страховой компании при условии, что ее капитал в начальный момент равен 5":
Р(5) = Рг {Т <м}. (1)
Как показано, например, в [1-3], вероятность разорения Р(5) удовлетворяет интегро-дифференциальному уравнению
СР (5) = ХР(5) - Р(5 - х)¥( х)с1х - х)с1х (2)
0 X
с очевидным граничным условием Р(+<х>) = 0. Для существования решения необходимо выполнение условия
С = (1 + 9)Ха, (3)
где 9> 0 - нагрузка страховой премии. В принципе, для нахождения решения уравнения (2) можно применить преобразование Лапласа. Сложность состоит в нахождении обратного преобразования.
Будем теперь предполагать, что страховые выплаты имеют гиперэкспоненциальное распределение п -го порядка
У (5) = ¿4 а ^ (4)
к =1
с параметрами ак > 0, Ак > 0, причем
к к=1
Среднее значение величины страховой выплаты в этом случае будет равно
. А
а
к=1 а к
X Ак = 1. (5)
= Х~. (6)
а
Распределение (4) в нашей задаче можно интерпретировать следующим образом. Имеется всего п различных типов страхования. При наступлении страхового случая страховая выплата с вероятностью Ак соответствует к-му типу страхового договора и имеет плотность распределения аке-ак. В силу линейной
независимости функций е-а распределение (4) можно использовать при подходящем выборе параметров для аппроксимации произвольного монотонного распределения У (5). При п = 2 уравнение (2) с плотностью распределения (4) решалось в работе [1].
Учитывая граничное условие Р(+да) = 0, уравнение (2) можно переписать в виде
5 да да да
СР( 5) = х| Р(5 - х)|т( у )СуСх + X | у)с1ус1х. (7)
0 х Б х
Решение уравнения (7) будем искать в виде
Р(5) = 1Р/ ^, (8)
,=1
где параметры у, > 0, Р подлежат определению. Подставляя выражение (8) в уравнение (7) и прирав-
нивая коэффициенты при линейно независимых функциях е акБ и е , получим, что выражение (8) есть решение уравнения (7), если выполняются следующие условия:
С- XI—^ = 0, (9)
к=1 ак -У, п 1 1
х—!—р =—. (10)
,=1 ак- у, ак
Таким образом, величины у, должны быть положительными корнями алгебраического уравнения
г (7) = с- х£ = о. (11)
к=1 а к— 7
Покажем, что уравнение (11) имеет ровно п различных положительных корней. Будем считать, что
п А
а1 < а2 <... < ап. Имеем, во-первых, /(0) = С - Ха = 9Ха > 0. Во-вторых, /'(0) = - ХХ_2 < 0. Таким об-
к=1 ак
разом, в окрестности точки г = 0 функция /(г) монотонно убывает. При г — а1 - 0 функция /(г) — -да. Следовательно, на отрезке (0, а1) уравнение (11) имеет корень у1. Далее, при г — а1 +0 /(г) —^ +да, при г — а2- 0 /(г) — -да. На отрезке (а1, а2) функция /(г) непрерывна. Следовательно, на отрезке (а1,а2) уравнение (11) имеет корень у2. Перебирая остальные отрезки (а,.,а1.+1), убеждаемся, что уравнение (11) имеет ровно п различных положительных корней.
Коэффициенты Р. должны являться решением системы уравнений (11). Определитель матрицы системы |^(ак -у..)-1 ^ есть определитель Коши [5], равный
П (а- а.)(у- у) -. (12)
А = det
ак -Уj
ПП(ак "Уj)
к=1 j=1
Поэтому система уравнений (10) имеет, причем единственное, решение. Покажем, что решение системы (10) положительно. Пусть система решается по формулам Крамера [5]. Тогда
А -
р =_m m = 1 n
1 m а ' ' '
где определитель Am получается из определителя А заменой его m-го столбца на столбец из свободных членов уравнений (10). Рассмотрим вначале определитель А (12). Так как при
7 /Л /Л /-,„4 П(П -1)
к < j ак - а < 0, у -ук > 0, то в числителе выражения (12) стоит произведение из - отрица-
2
тельных сомножителей. Так как уj <ак при j < к, то в знаменателе выражения (12) стоит произведение из n(n_2+1) отрицательных сомножителей. Поэтому А > 0. Определитель Am вычисляется по формуле (12) при уm = 0. Поэтому как в числителе, так и в знаменателе выражения (12) добавится m -1 отрицательный сомножитель и, следовательно, Am > 0.
Второй характеристикой, позволяющей анализировать перспективы страховой компании, является среднее время до разорения t(S) при условии, что разорение произошло и начальный капитал компании равен S. Как показано в [3], среднее условное время t(S) определяется соотношением
t (S) = , (13)
P( S)
где функция T (S) удовлетворяет интегро-дифференциальному уравнению
S
CT (S) = XT (S) - Xj T (S - x)¥( x) dx - P( S) (14)
0
с граничным условием T (+да) = 0. Учитывая граничное условие, уравнение (14) можно переписать в виде
S да да
CT(S) = XjT(S - x)jV(y)dydx + J P(y)dy. (15)
0 x S
В нашем случае функции ^(S) и P(S) определяются соотношениями (4) и (8) соответственно, поэтому решение уравнения (15) будем искать в виде
T(S) = £ (Uj + VjS)e YjS. (16)
j=1
Подставляя выражения (4), (8) и (16) в уравнение (15) и приравнивая коэффициенты при е-а^, e y S и Se 1]S и учитывая (9), получим, что коэффициенты Uj и Vj должны удовлетворять соотношениям
P,
V■ =-j-, (17)
j » A
Ху j Si^
Ё Л- Uj =1^77 Vj. (18)
j=1 ак у j j=1 (ак у j)
2. Модель Крамера-Лундберга со стохастическими премиями
Второй достаточно хорошо известной моделью страхования является модель Крамера-Лундберга со стохастическими премиями [6, 7]. В этом случае предполагается, что поток страховых премий, поступающих в компанию, является пуассоновским с интенсивностью X, премии - независимые, одинаково распределенные случайные величины с плотностью распределения Ф(x) и средним значением a . Страховые выплаты также образуют пуассоновский поток с интенсивностью ц, выплаты - независимые, одинаково распределенные случайные величины с плотностью распределения ^(x) и средним значением b. Поток страховых выплат не зависит от потока страховых премий. Как показано в [6, 7], в этом случае вероятность разорения (1) определяется уравнением
S м
(X + ц) P(S) = xj™ P(S + x )Ф( x)dx + P( S - x)¥( x )dx + x )dx (19)
0 S
с граничным условием P(+ro) = 0. Для существования решения уравнения необходимо выполнение условия
Xa = (1 + 0) цЬ, (20)
где 0 > 0 - нагрузка страховой премии.
Получим решение уравнения (19) в предположении, что распределения страховых премий и выплат являются гиперэкспоненциальными:
m n
Ф^) = £ Ak ake-«kS, S) = £ Bk^-kS. (21)
k=1 k=1
Решение уравнения (19) будем искать в виде
P(S) = £TP.e-1jS. (22)
j =1
Подставляя выражения (22) и (21) в уравнение (19) и приравнивая коэффициенты при линейно независимых функциях e--kS и e 1 ]S , получим, что выражение (22) есть решение уравнения (19), если выполняются условия
x + ц- xjr -^Ь- - ц£ -B-^ = 0, (23)
k=1 « +ij k=1 -k-ij n 1 1
P =-L. (24)
Таким образом, величины 1. должны быть положительными корнями алгебраического уравнения
f (z) = xjr «Л- + ц£ -Bk-k- - X - ц = 0. (25)
k=1 «k + z k=1 -k - z
Покажем, что уравнение (25) имеет ровно n различных положительных и m -1 различных отрицательных корней. Имеем, во-первых, что f (0) = 0, во-вторых, f '(0) = -Xa + цЬ < 0. Значит, в окрестности точки z = 0 f (z) = -(Xa -цЬ) z + o(z). Функция f (z) непрерывна на [0, ) и при z —^ — 0 f (z) —+м . Поэтому на отрезке [0, ) уравнение (25) имеет корень 11. Далее, функция f (z) непрерывна на (Р1, -2). При z — + 0 f (z) — -м, при z — -2 - 0 f (z) — +м. Поэтому на отрезке (Р1,Р2) уравнение (25) имеет корень 12. Аналогично устанавливается наличие остальных n- 2
положительных и m -1 различных отрицательных корней уравнения. Так как уравнение (24) по форме совпадает с уравнением (10), то решение уравнения (24) положительно.
Среднее условное время до разорения t(S) для данной модели определится соотношением (13), где теперь функция T (S) является решением уравнения [6]:
м S
(X + ц)Г (S) = Xj T (x - S )dx + ц| T (x)¥( S - x)dx + P( S) (26)
S
0
с граничным условием Т (+<») = 0. Решение уравнения (26) будем опять искать в виде
Т (5) = £ + ¥^е~3, (27)
3=1
где показатели степеней уудовлетворяют соотношениям (23). Подставляя соотношения (27), (21) и (22) в уравнение (26), получим систему соотношений на коэффициенты и. и ¥3
Р
¥■ =-3-, (28)
3 Ху 4а, + цу_ВР^
¿К +у, )2 ^(Р, - у, )2
я 1 я ¥
^^ и, =£ ТТ"^. (29)
3=1 Рк -уз 3=1 (Рк -уз )
Аналогично вышеизложенному могут быть найдены вероятность разорения и среднее условное время до разорения для модели Крамера-Лундберга со стохастическими премиями и постоянными не страховыми расходами [8].
3. Модель Крамера-Лундберга с ММР потоком страховых выплат
Одним из возможных обобщений модели Крамера-Лундберга является допущение того, что интенсивность потока страховых выплат или скорость поступления страховых премий могут скачкообразно изменяться в случайные моменты времени. Такие модели были рассмотрены, например, в работах [9-12]. В настоящей работе будем считать, что интенсивность потока страховых выплат Х^) является
однородной цепью Маркова с непрерывным временем и двумя состояниями Х1 и Х2. Переход из состояния в состояние задается матрицей инфинитезимальных характеристик
"-Р1 Р _ Р2 -Р2
где Р1, Р2 > 0. Тогда финальные вероятности состояний Х1 и Х2 равны соответственно
В =
(30)
^ =—^ (31)
1 Pl +Р2 2 Pl +Р2 и средняя интенсивность потока страховых выплат в стационарном режиме
Х0 = Х1л1 +Х2 я2. (32)
Страховые выплаты - независимые одинаково распределенные случайные величины с плотностью распределения ¥(x) и средним значением а . Наконец, будем считать, что страховые премии поступают непрерывно во времени с постоянной скоростью C. Для нормального функционирования страховой компании необходимо выполнение условия [11]:
C = (1 + 0) Х0а, (33)
где 0 > 0. При 0< 0 компания разоряется.
Обозначим, как и ранее, через Т момент разорения страховой компании, и пусть
P(S) = Рг{Г <м|S(0) = S, Х(0) = Х,}, i = 1,2 (34)
есть вероятности разорения страховой компании при условии, что в начальный момент ее капитал равен S и значение интенсивности потока выплат равно Хi. Как показано в [11], вероятности р (S) удовлетворяют системе уравнений
S м
CP(S ) = (Х1 +Р1) P(S ) -P1P2(S ) -Х1J р( x)^ (S-x)dx- Х1J^ ( x)dx, (35)
0 S
S м
CP2(S ) = -P2 P1 (S ) + (Х2 + P2) P2(S ) - Х 2 J P2( x)^ (S - x)dx - Х 2 J^ ( x)dx
0S
с граничными условиями Pi (+м) = 0.
Пусть страховые выплаты имеют гиперэкспоненциальное распределение п -го порядка (8). Решения системы уравнений (35) будем искать в виде
2п
Р (?) = 2Р^. (36)
к =1
Подставляя соотношения (36) и (8) в уравнения (35) и приравнивая коэффициенты при линейно независимых функциях е"^ и е"Ук?, получим систему соотношений
^ а 1 —
2-— Рк1 = 1, 1 = 1, п, (37)
к=1 а 1 "У к
2-— Р 2 = 1, 1 = 1, п, (38)
к=1 а 1 "У к
А■ а ,■
(^2 +Р2 + Сук "^22^"^)Рк2 "Р2^к1 = 0, к = 1,2п, (39)
1=1 а 1 "У к
А, а,
(X, +Р1 + С у к " X Рк1 " РР 2 = 0, к = 1,2п. (40)
1=1 а 1 "У к
При фиксированных ук относительно коэффициентов Рк1 и Рк2 имеем систему из 6п уравнений с 4п неизвестными. Для совместности системы необходимо, чтобы ранг получившейся системы уравнений был равен 4п [5]. Для этого необходимо, чтобы коэффициенты соответствующих уравнений систем (39) и (40) были пропорциональны, т.е. чтобы выполнялись условия
А,а . " А,а .
(X +Р1 +СУк " X2^-^)(Х2 +Р2 +СУк " X" Р1Р2 = 0, к = 1,2п. (41)
1=1 а 1 "У к 1=1 а 1 "У к
Таким образом, величины Ук должны являться корнями уравнения
/ (=(А, +Р1+С2" х 2 о^сх+Р2+с^" Х2 2 о^) - Р1Р2=0. (42)
1=1 а 1 2 1=1 а 1 2
Покажем, что уравнение (42) имеет 2п различных положительных корней.
Корни уравнения (42) совпадают с корнями многочлена степени 2п +2. И значит, уравнение (42) имеет всего 2п +2 корня. При 2 = 0 / (2) = 0. Следовательно, 2 = 0 - корень уравнения (42). Далее,
/ '(0) = (Р1 +Р2)(С" Х0а) >0. Поэтому при малых 2 <0 /(2) <0. С другой стороны, при 2»1
/ (2) ~ 22 и, следовательно, / (2) > 0. Поэтому уравнение (42) имеет хотя бы один отрицательный корень.
Перепишем функцию / (2) в виде
/ (2)=2 2(с" х 2—~—)(с" Х2 22 —^М+2 (Р1+Р2)(с" х0)22 (43)
1=1а 1 2 1=1а 1 2 1=1а 1 2
и обозначим
/0( 2) = С" X 22 а1. (44)
1=1 а 1"2
Как было показано ранее, уравнение /0(2) = 0 имеет п различных корней к1,к2,...,кп, удовлетворяющих неравенствам
0 <к <а1 <к2 <а2 <... <к <а .
112 2 п п
Воспользуемся этим. Имеем теперь, во-первых, /(0) = 0, /'(0) > 0; во-вторых,
/ к) = к2(х " Х)(Х " Х)(2-^ )2 < 0,
1=1 а 1 " к
так как или Х1 — Х0 < 0, или Х2 — Х0 < 0 . Поэтому уравнение /(2) = 0 имеет корень у1 е (0,к1) . Далее, п а
при 2 ^а1 /(2) ~ (V--—)2 и, следовательно, /(2) ^+со. Таким образом, уравнение /(2) = 0 име-
-=1 а - — 2
ет корень у2 е (к^ а1). Аналогично
/(к2) = кЦХ, — Х)(Х2 — ——-)2 < 0,
-=1 а - к2
откуда вытекает существование корня уравнения у3 е (а1, к2). Продолжая этот процесс, убеждаемся в существовании 2п различных положительных корней уравнения (43), удовлетворяющих неравенствам
0 <У1 <к1 <У2 <а1 <Уз <к2 <У4 <... <кп <У2п <ап. Определив величины у - , коэффициенты Рк1 и Рк 2 получим теперь как решение системы уравнений
2 п
а-
V-Рк1 = 1, - = 1, п, (45)
к=1 а - — У к
1 2п а п а _
- V-— Ук (С — х V-'—)Рк1 = 0, - = 1, п, (46)
Р1 к=1а -— у к и а — у к
Р2 = (1+7Т (С — X V -^1-))Р1. (47)
Р1 /=1 а,- — у к
Пример. В качестве примера рассмотрим случай гиперэкспоненциального распределения при п = 2 . Параметры распределения: а1 = 1, а2 = 2, Д = 0,5, А2 = 0,5 ; параметры потока страховых выплат: X = 2, X = 4, Р1 = 3, Р2 = 5 ; нагрузка страховой премии 0 = 0,2. Корни уравнения (42) в этом случае равны: у1 = 0,19, у2 = 0,866,у3 = 1,683, у4 = 1,839. Графики вероятностей разорения Р(£) и Р2(£) приведены на рис. 1. Как и следовало ожидать, вероятность разорения для состояния с большей интенсивностью выплат больше, чем вероятность разорения для состояния с меньшей интенсивностью выплат.
Р,я
05
06
0<
\
Х/\ "чХ \\
^Х V я
X
] 5
Рис. 1. Зависимость вероятностей разорения Р{ ) от начального капитала
Заключение
В работе получены расчетные формулы, позволяющие вычислить вероятность разорения и среднее условное время до разорения страховой компании при гиперэкспоненциальных распределениях страховых выплат и страховых премий для модели Крамера-Лундберга, модели Крамера-Лундберга со стохастическими премиями, модели Крамера-Лундберга с ММР потоком страховых выплат. В принци-
пе, тот же подход может быть использован для вычисления вероятностей разорения и для более сложных моделей, например, для случая, когда и моменты поступления страховых премий образуют ММР поток [13].
ЛИТЕРАТУРА
1. Panjer H.Y., Wilmot G.E. Insurance Risk Models. Society of Actuaries, 1992. P. 442.
2. Бауэрс Н., Гербер Х., Джонс Д., Несбит С., Хикман Д. Актуарная математика. М. : Янус К, 2001. 656 с.
3. Глухова Е.В., Змеев О.А., Лившиц К.И. Математические модели страхования. Томск : Изд-во Том. ун-та, 2004. 180 с.
4. Назаров А.А., Бронер В.И. Система управления запасами с гиперэкспоненциальным распределением объемов потребления
ресурсов // Вестник Томского государственного университета. Управление, вычислительная техника и информатика. 2016. № 1(34). С. 43-50.
5. Лившиц К.И. Линейная алгебра и аналитическая геометрия. Томск: Изд-во НТЛ, 2011. Ч. I. 252 с.
6. Livshits K.I. Probability of Ruin of an Insurance Company for the Poisson Model // Russian Physics Journal. 1999. V. 42, No. 4.
P. 394-399.
7. Бойков А.В. Модель Крамера-Лундберга со стохастическими премиями // Теория вероятностей и ее применение. 2002. Т. 47,
вып. 3. С. 549-553.
8. Livshits K.I., Yakimovich K.Yu. Cramer-Lundberg Model with Stochastic Premiums and Continuous Non-insurance Costs // Com-
munications in Computer and Information Science. Springer, 2014. V. 487. P. 251-260.
9. Jasinlewicz H. Probability of Ruin with Variable Premiun Rate in a Markovian Environment // Insurance: Mathematics end Econom-
ics. 2001. V. 29. P. 291-296.
10. Lu Y., Li S. On the Probability of Ruin in a Markov-modulated Risk Model // Insurance: Mathematics end Economics. 2005. V. 37(3). P. 522-532.
11. Лившиц К. И., Бублик Я. С. Вероятность разорения страховой компании при дважды стохастическом потоке страховых выплат // Вестник Томского государственного университета. Управление, вычислительная техника и информатика. 2010. № 1(10). С. 66-77.
12. Лившиц К.И., Бублик Я.С. Распределение условного времени до разорения страховой компании при дважды стохастическом потоке страховых выплат // Вестник Томского государственного университета. Управление, вычислительная техника и информатика. 2010. № 4(13). С. 15-23.
13. Лившиц К. И., Бублик Я. С. Вероятность разорения страховой компании при дважды стохастических потоках страховых премий и страховых выплат // Вестник Томского государственного университета. Управление, вычислительная техника и информатика. 2011. № 4(17). С. 64-73.
Лившиц Климентий Исаакович, д-р техн. наук, профессор. E-mail: [email protected] Томский государственный университет
Сухотина Лариса Юрьевна, канд. физ.-мат. наук, доцент. E-mail: [email protected]
Томский государственный университет Поступила в редакцию 2 февраля 2016 г.
Livshits Klimenty I., Suhotina Larisa Yu. (Tomsk State University, Russian Federation).
Ruin probability of an insurance company with hyperexponential distribution of insurance premiums and insurance payments for different insurance models.
Keywords: probability of ruin; mean conditional time to ruin; Cramer-Lundberg model; Cramer-Lundberg model with stochastic premiums; hyperexponential distribution.
DOI: 10.17223/19988605/35/4
In this paper the calculated formulas which allow to calculate the ruin probability of an insurance company for some insurance models under assumption that insurance premiums and insurance payments have the hyperexponential distribution are obtained.
n
For Cramer-Lundberg model if the insurance payments distribution is T (S) = ^ Ak a ke"a'kS, rate cash flows is C and the insur-
k=1
ance payments flow intensity is X , then the probability of ruin is defined as
p(S ) = £ Pj
-JjS
J
J=1
where y , are the simple positive roots of equation
f ( z) = c—xjr = 0,
k=1 — k—z
n 1 1
and the coefficients P. are the solution of equations ^-P. =-.
,=i ak "Yj a*
For Cramer-Lundberg model with stochastic premiums if the insurance premiums and insurance payments distributions are
m n
®(S) = X Akake-akS, T(S) = X BkPke-PkS , respectively, the insurance premiums and insurance payments flows intensities are
k=1 k=1
X and ^ , then the probability of ruin is determined by the same relation, but y j are the simple positive roots of equation
f (z)=XX A---B^--x- 0,
k=1 ak + Z k =1 Pk _ z
n 1 1
and the coefficients Pj are the solution of equations X-Pj = —.
j=1 Pk - yj Pk
For Cramer-Lundberg model with MMP payments flow with two states X1 and X2 and the final probabilities of states P2 /(P1 +P2), P1 /(P1 +P2) the probabilities of ruin
2n
Pj (S) = X Pje-ykS, j = 1,2,
k=1
where y j are the simple positive roots of equation
n A a n A a
f (z) = (X1 +P1 + Cz - X1 +P2 + Cz - X2 " P1P2 = 0,
j=1 a j z j=1 a j z
and the coefficients Pkj are the solution of equations
2n a 1 2n a n A __y n A
X-— Pk1 = 1; T"X-— yk(C-X1X-i—)Pk1 = 0, j = 1,n; Pk2 = (1 + "PL(C-X1 ^_^))Pk1.
k=1a j- y k P1 k=1a j- y k '=1 a- y k P1 '=1 a- y k
REFERENCES
1. Panjer, H.Y. & Wilmot, G.E. (1992) Insurance Risk Models. Society of Actuaries.
2. Bowers, N.L., Gerber, H.U., Hekman, J.C., Jones, D.A. & Nesbitt, C.J. (2001) Aktuarnaya Matematika [Actuarial Mathematics].
Translated from English by V.K. Malinovskiy. Moscow: Yanus K.
3. Glukhova, E.V., Zmeev, O.A. & Livshits, K.I. (2004) Matematicheskie model' strakhovaniya [Mathematical Models of Insurance].
Tomsk: Tomsk State University.
4. Nazarov, A.A. & Broner, V.I. (2016) Inventory model with hyperexponential distribution of demand's batch size. Vestnik Tomskogo
gosudarstvennogo universiteta. Upravlenie, vychislitel'naya tekhnika i informatika - Tomsk State University Journal of Control and Computer Science. 1(34). pp. 43-50. DOI: 10.17223/19988605/34/5
5. Livshits, K.I. (2011) Lineynaya algebra i analiticheskaya geometriya [Linear Algebra and Analytical Geometry]. Tomsk: NTL.
6. Livshits, K.I. (1999) Probability of Ruin of an Insurance Company for the Poisson Model. Russian Physics Journal. 42(4). pp. 394-
399. DOI: 10.1007/BF02509675
7. Boykov, A.V. (2002) The Cramer-Lundberg model with stochastic premium process. Teoriya veroyatnostey i ee primenenie. 47(3).
pp. 549-553. DOI: 10.4213/tvp3693
8. Livshits, K.I. & Yakimovich, K.Yu. (2014) Cramer-Lundberg Model with Stochastic Premiums and Continuous Non-insurance Costs.
Communications in Computer and Information Science. 487. pp. 251-260. DOI: 10.1007/978-3-319-13671-4_30
9. Jasinlewicz, H. (2001) Probability of Ruin with Variable Premiun Rate in a Markovian Environment. Insurance: Mathematics and
Economics. 29. pp. 291-296. DOI: 10.1016/S0167-6687(01)00090-7
10. Lu, Y. & Li, S. (2005) On the Probability of Ruin in a Markov-modulated Risk Model. Insurance: Mathematics and Economics. 37(3). pp. 522-532. DOI: 10.1007/s11009-007-9044-4
11. Livshits, K.I. & Bublik, Ya.S. (2010) Ruin probability of an insurance company under double stochasticpayment current. Vestnik Tomskogo gosudarstvennogo universiteta. Upravlenie, vychislitel'naya tekhnika i informatika - Tomsk State University Journal of Control and Computer Science. 1(10). pp. 66-77. (In Russian).
12. Livshits, K.I. & Bublik, Ya.S. (2010) Distribution of the conditional time before ruin of an insurancecompany under double stochastic payment current. Vestnik Tomskogo gosudarstvennogo universiteta. Upravlenie, vychislitel'naya tekhnika i informatika - Tomsk State University Journal of Control and Computer Science. 4(13). pp. 15-23. (In Russian).
13. Livshits, K.I. & Bublik, Ya.S. (2011) Ruin probability of an insurance company under double stochastic insurance premium and insurance payment currents. Vestnik Tomskogo gosudarstvennogo universiteta. Upravlenie, vychislitel'naya tekhnika i informatika -Tomsk State University Journal of Control and Computer Science. 4(17). pp. 64-73. (In Russian).