We consider the perturbed dual risk model with constant interest and a threshold dividend strategy. Firstly, we investigate the moment-generation function of the present value of total dividends until ruin. Integrodifferential equations with certain boundary conditions are derived for the present value of total dividends. Furthermore, using techniques of sinc numerical methods, we obtain the approximation results to the expected present value of total dividends. Finally, numerical examples are presented to show the impact of interest on the expected present value of total dividends and the absolute ruin probability.

1. Introduction

In insurance mathematics, the classical risk model has been the center of focus for decades. The surplus in the classical model at time can be expressed as where is the initial surplus, is the premium rate, and is the aggregate claims by time , usually modeled by a compound Poisson process. In recent years, quite a few interesting papers have been written on a model which is dual to the classical insurance risk model. See, for example, Avanzi et al. [1], Avanzi and Gerber [2]. The surplus of a dual classical risk model at time is

In this model, the premium rate should be viewed as an expense rate and claims should be viewed as profits or gains. While not very popular in insurance mathematics, this model has appeared in various literature (see Cramér [3, Section  5.13] and Seal [4, pages 116–119] and Takács [5, pages 152–154]). There were many possible interpretations for this model. For example, we can treat the surplus as the amount of the capital of a business engaged in research and development. The company paid expenses for research, and occasional profit of random amounts (such as the award of a patent or a sudden increase in sales) arises according to a Poisson process. A similar model was used by Bayraktar and Egami [6] to model the capital of a venture capital investment. Another model was an annuity business. The company issues payments continuously to annuitants, while the gross reserve of an annuitant was released as emerging profit when he died. Yang and Zhu [7] generalized the dual model into a regime-switching setting and calculated bounds for ruin probabilities. One of the current topics of interest in insurance mathematics is the problem of maximizing the expected total discounted dividends until ruin, which goes back to de Finetti [8], and it had also been studied by Bühlmann [9, Section  6.4], Gerber [10, Sections  7 and 8], and Gerber [11, Section 10.1]. In Avanzi et al. [1], the authors studied the expected total discounted dividends until ruin for the dual model under the barrier strategy by means of integrodifferential equations. They derived explicit formulas when profits or gains followed an exponential or a mixture of exponential distributions and showed that the optimal value of the dividend barrier under the dual model was independent of the initial surplus. Avanzi and Gerber [2] studied a dual model perturbed by diffusion and discussed how the optimal value of the dividend barrier can be determined. Albrecher et al. [12] studied a dual model that also paid taxes when the surplus was at a running maximum and calculated the expected total discounted dividends before ruin for exponentially distributed profits.

Considering the perturbed dual risk model, the surplus of an insurer has the following form where is the initial surplus and is the constant rate. is a sequence of independent income size random variables with a common distribution function with and density function . is the Poisson income-number process with an intensity and is defined as , where the i.i.d. interincome times have a common exponential distribution with a parameter . is a standard Brownian motion with and is a constant, representing the diffusion volatility parameter. In addition, , and , are mutually independent.

Dividend strategy for insurance risk models was first proposed by de Finetti [8] to reflect more realistically the surplus cash flowed in an insurance portfolio. Various threshold strategies have been studied by many authors, including Gerber and Landry [13], Ng [14], and Lin and Pavlova [15]. Among them, Ng [14] showed that the threshold strategy was optimal when the dividend rate was bounded and the individual claim distribution was exponential. Recently, the threshold dividend strategy has been considered in the class of compound Poisson process perturbed by diffusion and its generalizations; readers may refer to Avanzi and Gerber [2], Gao and Liu [16], Wan [17], and the references therein.

Then, we consider the modification of the surplus process by a threshold strategy with a threshold level . When is below , no dividends are paid and the surplus decreases at the original rate . But when is above , the surplus would decrease at a different rate and dividends are paid at rate . Incorporating the threshold strategy into (3) yields the surplus process which can be expressed by where , , and .

In this paper, we consider that the insurance company earns credit interest with a constant force . In the mean time, the insurer will receive interest. Incorporating interest into (4) yields the surplus process which can be expressed by

Let denote the cumulative amount of dividends paid out up to time and the force of interest; then is the present value of all dividends until , where denoted by is the time of ruin. An alternative expression for is with denoting the indicator function. It is obvious that .

In the sequel, we will be interested in the following moment generating function: (for those values of where it exist) and the expected total discounted dividends until ruin is Throughout this paper, we assume that and are sufficiently smooth functions in and , respectively.

2. Integrodifferential Equations

2.1. Integrodifferential Equations for

In this section, we will give the integrodifferential equations satisfied by the moment generating function , respectively. It is easy to see that behaves differently with different initial surplus. Hence, for notation convenience, we set In the case of , we write as , and the other one is similar.

In the following, we firstly derive the integrodifferential equations satisfied by .

Theorem 1. When , and  for , with the following boundary conditions:

Proof. When , consider such that the surplus cannot reach level by time; that is, . By conditioning on the time and amount of first claim, if it occurs by , and whether the claim causes ruin, one gets By Taylor’s expansion, Also, due to Substituting the above expansion into (16), dividing both sides of (16) by , and letting , we can get (11). Using the same way as (11), we can easily prove (12) when , so the condition (12) is correct, by the method used in Gao and Liu [16], and it is easy to check that the boundary conditions hold. This ends the proof of Theorem 1.

Now we consider the case of . For , define . Thus, is an It stochastic integral. Denote by the variance process of . We have for . Let ; then

Set for , . By the time change of Brownian motion, we have that is local standard Brownian motion with running for an amount of time .

Let and . Then, we have the following results.

Theorem 2. When , and when , with the following boundary conditions:

Proof. For , assume that such that . Define and . We have for all . Then, by the strong Markov property, we get By the assumption of independence, one gets where Noting that , a.s. and using It’s formula, we have From (26) and (29), we can get (20) for all . Hence, (20) holds in .
Similarly, using the above method, we get (21). Conditions (22) and (23) are obvious. By (14) and (15) and by the weak convergence method used in [17], it is easy to check that the boundary conditions (24) and (25) hold. This ends the proof of Theorem 2.

2.2. Integrodifferential Equations for

Theorem 3. When , and when , with the following boundary conditions:

3. Numerical Analysis

The second order system of integrodifferential equations such as in Theorem 3 does not have an explicit solution except in some special cases. In this section, we propose the approximate solution of this equation via the use of a collocation based on sinc methods. The sinc method is highly efficient numerical method developed by Frank Stenger. It is widely used in various fields of numerical analysis such as interpolation, quadrature, approximation of transforms, and solution of integral, ordinary differential, and partial differential equations. An introduction to the sinc approximation theory can be found in the Appendix.

3.1. Numerical Solution of the Expected Present Value of Total Dividends

To construct an approximation on the interval , we consider the following conformal map:

The function also provides a one-to-one transformation of onto the real line . The sinc grid points are defined for and , by

In order to adopt the sinc method procedure, we arrange the systems in Theorem 3 into the following integrodifferential equation:

Using the properties of convolution, the above equation can be further written as with the following boundary conditions:

Set ; then and satisfies with the following boundary conditions: where

Let be defined by (33). The sinc grid points are defined for by

Then, by using Theorems A.3 and A.6 in the Appendix, we have where , is the th element of matrix . denotes an approximate value of , and , , , and . and for are defined as Definition A.2 in the Appendix.

Since , by Theorem A.3 in the Appendix, then it is convenient to take , . Moreover, by some simple calculations, we have , with and being defined in (33) and (41), respectively.

Having replaced the integral term on the left-hand side of system (38) with the right-hand side of (42) and having substituted for , where are sinc grid points and also by replacing by as in (43), we obtain the following collocation result: where .

Then, by Theorem A.7, we have

Replacing (45) in (44), we rewrite (44) as

Having multiplied the resulting equations by , we have

Now, since , , , and , we obtain the collocation result as

We set and , where denotes the th element of the matrix . Also, we denote and , are square matrices of order . So, the system (48) can be given in matrix form as where ,  , and

The above linear system contains equations with unknown coefficient , . Solving this linear system, we obtain the approximate solution of the system as follows:

4. Numerical Example

In this section, we consider some numerical samples to illustrate the performance of sinc method and investigate how much the values of are affected by interest . An example is solved under the assumption that the claim size density is given by .

Example 1. Let , , , , , , , , , , , and , the result is shown in Figure 1.
From Figure 1, we see that is an increasing function with respect to .


The sinc method is a highly efficient numerical method developed by Stenger, the pioneer of this field, people in his school, and others [1822]. It is widely used in various fields of numerical analysis such as interpolation, quadrature, approximation of transforms, and solution of integral, ordinary differential, and partial differential equations.

sinc methods are based on the use of the cardinal function, , which is sinc expansion of function , defined by where the step size , and the function is defined on the whole real line, , by

For any , the translated sinc functions with evenly spaced nodes are given as

The sinc functions are cardinal for the interpolating points in the sense that

Definition A.1 (see [19, P73, De. 1.5.2]). Let denote a smooth one-to-one transformation of an arc  with end-point and onto , such that and . Let denote the inverse map, so that
Given , , and a positive number , define the sinc points by and a function by
Observe that increases from to as traverses from to .
Corresponding to positive numbers , , and , let denote the family of all functions defined on for which such that the Fourier transform satisfies the following relation: for all .
Another important family of functions is , with , and . It consists of all those functions defined on such that and where is defined by .
Let and denote two positive integers, such that , or where denotes the greatest integer function. Moreover, we denote .

Definition A.2. Given three positive integers , , and as above, let and denote linear operators acting on functions defined on given by where denote the sinc points. Set and we then define an matrix by , with denoting the th element of . Finally, let denote the uniform norm on that is,

Theorem A.3 (see [19, P85, Th. 1.5.13]). Let ; then, as , If , then it is convenient to take , .

Theorem A.4 (see [19, P87, Th. 1.5.14]). Let ; then,

Theorem A.5 (see [19, P93, Th. 1.5.19]). Let ; then, for all , where with denoting the transpose of .

Theorem A.6 (see [19, P95-96, Th. 1.5.20]). Let ; suppose that can be diagonalized with and with and that is a diagonal matrix. If for some positive independent of , one has for all , then there exists a constant independent of such that where with .

Theorem A.7 (see [18, P106, Th. 4.1]). Let be a conformal one-to-one transformation of an arc . Then,


This work was supported in part by the National Natural Science Foundation of China under Grant nos. 61370096 and 61173012 and the Key Project of Natural Science Foundation of Hunan Province under Grant no. 12JJA005.