Complex Boundary Value Problems of Nonlinear Differential Equations: Theory, Computational Methods, and Applications
View this Special IssueResearch Article  Open Access
Numerical Method for a MarkovModulated Risk Model with TwoSided Jumps
Abstract
This paper considers a perturbed Markovmodulated risk model with twosided jumps, where both the upward and downward jumps follow arbitrary distribution. We first derive a system of differential equations for the GerberShiu function. Furthermore, a numerical result is given based on Chebyshev polynomial approximation. Finally, an example is provided to illustrate the method.
1. Introduction
The risk model with twosided jumps was first proposed by Boucherie et al. [1] and has been further investigated by many authors during the last few years. For example, Kou and Wang [2] studied the Laplace transform of the first passage time and the overshoot for a perturbed compound Poisson model with double exponential jumps. Xing et al. [3] extended the results of Kou and Wang [2] to the case that the surplus process with phasetype downward and arbitrary upward jumps. Zhang et al. [4] assumed that the downward jumps follow arbitrary distribution and the upward jumps have a rational Laplace transform. They derived the Laplace transform of the GerberShiu function by using the roots of the generalized Lundberg equation. Under the assumption that the upward jumps follow Laplace distribution and arbitrary downward jumps, Chi [5] obtained a closedform expression for the GerberShiu function by applying WienerHopf factorization technique. The applications of the model in finance were also discussed. Jacobsen [6] studied a perturbed renewal risk model with phasetype interclaim times and twosided jumps, where both the jumps have rational Laplace transforms. Based on the roots of the CramérLundberg equation, the joint Laplace transform on the time to ruin, and the undershoot at ruin were given. However, in all the aforemetioned papers, the topic that the jumps in both directions are arbitrary distributions is still not discussed. The Markovmodulated risk model (Markovian regime switching model) was first proposed by Asmussen [7] to extend the classical risk model. Since then, it has received remarkable attention in actuarial mathematics, see, for example, Zhu and Yang [8, 9], Zhang et al. [4], Ng and Yang [10], Li and Lu [11], Lu and Tsai [12], and references therein. Motivated by the papers mentioned above, in this paper, we will study the Markovmodulated risk model with twosided jumps.
Let be a homogenous, irreducible, and recurrent Markov process with finite state space . Denote the intensity matrix of by with and for . Let be a sequence of independent random variables representing the jumps, and be a standard Brownian motion with . Here we assume that the premium rates, claim interarrival times, the distributions of the jumps, and the diffusion parameter are all influenced by the environment process . When , the premium rate is , jumps arrive according to a Poisson process with intensity , the diffusion parameter is , and the size of the jumps which arrives at time follows the distribution with density and finite mean . Then the Markovmodulated diffusion risk model is defined by where is the initial surplus. If we denote the stationary distribution of by , then the positive security loading condition is given by
In this paper, we further assume that the jumps in (1.1) are twosided. The upward jumps can be explained as the random income (premium or investment), while the downward jumps are interpreted as the random loss. In this case, the density function is given by where , , is the indicator function, and are two arbitrary functions on .
Let ( otherwise) be the time to ruin. For , let be the GerberShiu function at ruin given that the initial state is , where is a nonnegative penalty function, is the surplus immediately prior to ruin, and is the deficit at ruin. Without loss of generality, we assume that . Thus for . When , (1.4) reduces to the Laplace transform of the time to ruin when and , (1.4) reduces to the probability of ruin
The purpose of this paper is to present some numerical results on the GerberShiu function for the Markovmodulated diffusion risk model with arbitrary upward and downward jumps. In Section 2 we derive a system of integrodifferential equations and approximate solutions for . Numerical example is given in the last section.
2. Integrodifferential Equations and Approximate Solution
Theorem 2.1. For , satisfies the following integrodifferential equation where with boundary conditions
Proof. Similar to Ng and Yang [10].
Remark 2.2. When , (2.1) is identical to in Zhang et al. [4].
Clearly, (2.1) is a system of second order linear integrodifferential equations of FredholmVolterra type. As is well known, it is very difficult to find analytical solution of this system. Motivated by AkyüzDascioglu [13], we will study an alternative system defined on by Chebyshev collocation method. First, we transform the interval to . Following Diko and Usábel [14], we set , that is, . Furthermore, we assume that is an arbitrary strictly monotone, twice continuously differentiable function throughout the paper.
Theorem 2.3. Let be a monotone increase function and for . Then satisfies the following integrodifferential equation where with boundary conditions
Proof. By the definitions of function and , we have
Substituting (2.7) and into (2.1) and simplifying lead to (2.4). The boundary conditions are direct result of the boundary conditions in Theorem 2.1. This completes the proof.
Remark 2.4. The existence of the solution for the system of integrodifferential equations (2.4) can be found in Fariborzi and Behzadi [15].
According to AkyüzDascioglu [13], and its derivatives have truncated Chebyshev series expression where , are shifted Chebyshev polynomials of the first kind and are the unknown coefficients to be determined.
Let , , . Then (2.8) can be written in the matrix form where for odd , and for even .
Similarly, the kernel functions and can be expanded to univariate Chebyshev series where with and are Chebyshev coefficients determined by Clenshaw and Curtis [16].
Theorem 2.5. For , an approximate expression for is given by where the column vector can be determined by the following systems where matrix with elements matrix with elements and are collocations.
Proof. Using (2.8) and (2.12), one obtains Substituting (2.18) into (2.4), we have which is identical to (2.15) in form. Substituting the collocations into (2.19) leads to (2.15). and can be obtained by (2.6).
Example 2.6. To illustration our method, we use the example of Zhang et al. [4]. Let , , , , , , , the downward jumps are exponentially distributed with parameter , and the upward jump density is given by . We set and the collocation points are .
Figure 1 shows that the approximate solution is very near to the exact solution for any initial surplus . We remark that the horizontal axis in Figure 1 is and .
From Table 1 we can see that the errors between the approximate solutions and the exact solutions decrease when increases. The initial surplus can also influence the approximate solution: the bigger need a bigger to decrease the error.

Acknowledgments
This work is supported by the Natural Science Foundation of Shandong (no. ZR2010AQ015), the Tianyuan fund for Mathematics (no. 11226251), and the Natural Science Foundation of Qufu Normal University (no. 2012ZRB01473).
References
 R. J. Boucherie, O. J. Boxma, and K. Sigman, “A note on negative customers, $GI/G/1$ workload, and risk processes,” Probability in the Engineering and Informational Sciences, vol. 11, no. 3, pp. 305–311, 1997. View at: Publisher Site  Google Scholar
 S. G. Kou and H. Wang, “First passage times of a jump diffusion process,” Advances in Applied Probability, vol. 35, no. 2, pp. 504–531, 2003. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 X. Xing, W. Zhang, and Y. Jiang, “On the time to ruin and the deficit at ruin in a risk model with doublesided jumps,” Statistics & Probability Letters, vol. 78, no. 16, pp. 2692–2699, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 Z. Zhang, H. Yang, and S. Li, “The perturbed compound Poisson risk model with twosided jumps,” Journal of Computational and Applied Mathematics, vol. 233, no. 8, pp. 1773–1784, 2010. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 Y. Chi, “Analysis of the expected discounted penalty function for a general jumpdiffusion risk model and applications in finance,” Insurance: Mathematics & Economics, vol. 46, no. 2, pp. 385–396, 2010. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 M. Jacobsen, “The time to ruin for a class of Markov additive risk process with twosided jumps,” Advances in Applied Probability, vol. 37, no. 4, pp. 963–992, 2005. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 S. Asmussen, “Risk theory in a Markovian environment,” Scandinavian Actuarial Journal, no. 2, pp. 69–100, 1989. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 J. Zhu and H. Yang, “Ruin theory for a Markov regimeswitching model under a threshold dividend strategy,” Insurance: Mathematics & Economics, vol. 42, no. 1, pp. 311–318, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 J. Zhu and H. Yang, “On differentiability of ruin functions under Markovmodulated models,” Stochastic Processes and Their Applications, vol. 119, no. 5, pp. 1673–1695, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 A. C. Y. Ng and H. Yang, “On the joint distribution of surplus before and after ruin under a Markovian regime switching model,” Stochastic Processes and Their Applications, vol. 116, no. 2, pp. 244–266, 2006. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 S. Li and Y. Lu, “The decompositions of the discounted penalty functions and dividendspenalty identity in a Markovmodulated risk model,” ASTIN Bulletin, vol. 38, no. 1, pp. 53–71, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 Y. Lu and C. C. L. Tsai, “The expected discounted penalty at ruin for a Markovmodulated risk process perturbed by diffusion,” North American Actuarial Journal, vol. 11, no. 2, pp. 136–149, 2007. View at: Google Scholar
 A. AkyüzDascioglu, “A Chebyshev polynomial approach for linear FredholmVolterra integrodifferential equations in the most general form,” Applied Mathematics and Computation, vol. 181, no. 1, pp. 103–112, 2007. View at: Publisher Site  Google Scholar
 P. Diko and M. Usábel, “A numerical method for the expected penaltyreward function in a Markovmodulated jumpdiffusion process,” Insurance: Mathematics & Economics, vol. 49, no. 1, pp. 126–131, 2011. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 M. A. Fariborzi Araghi and S. S. Behzadi, “Numerical solution of nonlinear VolterraFredholm integrodifferential equations using homotopy analysis method,” Journal of Applied Mathematics and Computing, vol. 37, no. 12, pp. 1–12, 2011. View at: Publisher Site  Google Scholar
 C. W. Clenshaw and A. R. Curtis, “A method for numerical integration on an automatic computer,” Numerische Mathematik, vol. 2, pp. 197–205, 1960. View at: Publisher Site  Google Scholar  Zentralblatt MATH
Copyright
Copyright © 2012 Hua Dong and Xianghua Zhao. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.