# Valuation of endowment-insurance equity-linked contracts for stocks with exotic dynamics.

1. Introduction

The fair value of these kinds of products involves considering two separate sources of randomness: one stemming from the stochastic nature of the dynamics of the stock markets and a different one due to the uncertainty in mortality. Here we consider a simple model of financial market consisting of two securities: a savings account [B.sub.t] which evolves via d[B.sub.t] = [r.sub.t][B.sub.t]dt, where [r.sub.t] is the instantaneous interest rate of the market. For convenience it is assumed to be deterministic, but not necessarily constant as befits a contract held for a long time. The second instrument in the market, to which the policy is linked, is a given stock whose t-price, [X.sub.t], varies according to definite stochastic dynamics. The prototype model for stocks-price evolution [13, 14] assumes that [X.sub.t] is a geometric Brownian motion (GBM) model; that is, [X.sub.t] satisfies the stochastic differential equation (SDE)

d[X.sub.t] = [mu][X.sub.t]dt + [sigma][X.sub.t]d[W.sub.t]. (1)

Here [mu] is the mean return rate and [sigma] is the volatility, while [W.sub.t] is a Brownian motion under the empirical or real world probability. Despite the fact that this simple model describes well the basic dynamical properties, there are nonetheless several stylized facts that the model fails to capture.

In this paper we modify the dynamics so as to incorporate the possibility that the traded stock has a strong support at some level, say c. Such a feature may stem from a market consensus under which heavy buy orders are triggered when the stock price hits this level. Investing in such stock may well be an appealing feature to the insurance company as it partially hedges the exposure to a possible plummeting of the price. In Section 4 we propose a dynamical model which incorporates such a feature by modifying appropriately the "infinitesimal volatility" term; by using Feller's theory, we discuss the nature of the barrier. We next show that under some election of parameters the SDE that drives the dynamics can be solved in an exact way (see (39) to (42)). We suppose that the benefit is received either at the insurer's demise or at maturity, whichever comes first.

The structure of the paper is the following. In Section 2 we give the basic setup and determine explicit formulas for the fair premium required by the insurance firm to hedge the exposure to the evolution of the risky asset for an endowment insurance contract under quite general conditions. The fair value is given in terms of expectations with respect to the risk-free measure under which the discounted price of any security happens to be a martingale. In Section 3, assuming a simple GBM dynamics, concrete valuation formulas for several contracts of actuarial interest are obtained (see formulae (9)-(10), (17)-(20), and (26)-(28)). Recall that the classical work of Brennan and Schwartz [1] and Boyle and Schwartz [2] corresponds to a pure endowment contract with payoff at maturity [[PSI].sub.T] [equivalent to] max{[X.sub.T], k}. Remarkably, when the guaranteed amount is taken as the interest accrued with rate [delta] for a principal [X.sub.t], then the premium contingent on death can also be obtained in closed form, as we show below. The second case we study corresponds to a floating strike lookback put (FLOP) which entitles the owner to sell the stock at the highest realized price, that is, to a path dependent payoff [[PSI].sup.t.sub.T] [X.sup.t.sub.T] - [X.sub.T], where [X.sup.t.sub.T] [equivalent to] [max.sub.t [less than or equal to] s [less than or equal to] T] [X.sub.s] is the running maxima of the process. Look back options guarantee a nonloss outcome and as such are interesting in dynamic investment fund protection.

In Section 4 we obtain and solve a natural stochastic equation that models evolution in a situation when strong resistances are present (see (39) to (42)), whereupon we obtain the fair insurer's liability for several endowment-insurance contracts. The presence of barriers is found to complicate the prizing problem.

Section 5 is devoted to study the partial differential equation (PDE) that the premium satisfies. We derive under appropriate conditions a modified Kolmogorov Backwards equation (68) and sketch a brief comparison between replicating portfolios techniques and direct martingale approaches.

2. Contract Characteristics and Valuation

As we have already pointed out, unit-linked contracts involve two sources of randomness, one stemming from the associated financial market and another corresponding to mortality expectations. We now pass to discuss some features of the latter. Consider a given individual aged d at time 0 and let [tau] : [bar.[OMEGA]] [right arrow] R be the time at which decease occurs. Here, [tau] is a random variable defined on a complete probability space ([bar.[OMEGA]], [bar.G], [bar.P]). Under natural assumptions the survival function [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] given by

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (2)

where t [??] h(T) the conditional density of [tau] and [??] [equivalent to] T - t. Further the mortality intensity is taken in the Gompertz-Makeham form [mu](t + d, t) = [alpha] + [beta][[lambda].sup.(t+d)] (see [15, 16]). We finally introduce the filtration [[bar.F].sub.t] [equivalent to] [sigma]([1.sub.{[tau] > s}], 0 [less than or equal to] s [less than or equal to] t).

In a unit-linked contract the premium [v.sub.t] paid at time t by policyholders is invested in a equity fund. Let [X.sub.t] be the value at time t of a unit of the fund and ([OMEGA], [F.sub.[infinity]], P) be a probability space big enough to contain the filtration [F.sub.s] generated by the stock: [F.sub.t] [equivalent to] [sigma]([X.sub.s], 0 [less than or equal to] s [less than or equal to] t). Let [[PSI].sup.1).sub.T] be the reward payable at maturity T if the policyholder is alive. Further, if the insurer's demise happens at a time [tau] before maturity, the policy entitles the beneficiaries to a payment [[PSI].sup.2).sub.T] at the decease time. In other words, we consider an endowment insurance, mixing a pure endowment and a term insurance contract. Here [[PSI].sup.1).sub.T] are supposed to be right continuous processes adapted to the filtration [F.sub.s]. These conditions guarantee that [[PSI].sup.1,2)] are progressively measurable and that [[PSI].sup.1).sub.[tau]] is [F.sub.[tau]] measurable. Note that in our basic setup surrender is not permitted; however, it could be easily incorporated by treating the surrender decision like a premature decease, as several authors do.

We assume that the insurance company is risk-neutral with respect to mortality. This assumption means that it does not request any compensation for assuming mortality risk. The exposure could be hedged to some extent whenever a large number of contracts are written yearly.

We will also assume that our market is efficient; that is, the existence of the Harrison and Pliska [17] risk-neutral probability on ([OMEGA], [F.sub.[infinity]]) under which discounted prices of self-financing strategies [v'.sub.t] [equivalent to] [v.sub.t]/[B.sub.t] are martingales with respect to the history of the process up to time t.

The pricing of equity-linked life insurance policies is a classical problem in the actuarial literature and, particularly for pure endowment contracts, explicit formulas are known. However, far less is known as regards explicit analytical formulas for endowment insurance contracts even for the simplest GBM dynamics. Here we determine the insurer's liability required by the insurance firm to hedge the payoff for given [[PSI].sup.1)], [[PSI].sup.2)], and a general mortality intensity [mu](d + t, t).

Let [v.sub.t] be the price of such a contract written at time t at which [X.sub.t] [equivalent to] x is known and maturing at time T. Let E [equivalent to] [E.sup.Px[bar.P]] denote the expectation with respect to the product measure P x [bar.P]. Let [partial derivative] be the demise time or maturity, whichever comes first, [partial derivative] = T [conjunction] [tau]. Then, we have

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (3)

Here [v.sub.'.sub.t] [equivalent to] [v.sub.t]/[B.sub.t] are the deflated prices and we used the martingale property and the optional stopping theorem (note that the optional stopping theorem can be applied since [bar.P]([partial derivative] < [infinity]) = 1, [E.sup.P]([sup.sub.0 [less than or equal to] s [less than or equal to] T] [[PSI].sup.i).sub.s]) < [infinity]). Further,

[THETA] [equivalent to] [[[PSI].sup.1).sub.T]/[B.sub.T]] [1.sub.{[tau] [greater than or equal to] T}] + [[PSI].sup.2).sub.[tau]]/[B.sub.[tau]]] [1.sub.{[t [less than or equal to] [tau] < T}] (4)

is a function naturally defined on the product probability space ([OMEGA] x [bar.[OMEGA]], [F.sub.[infinity]] [disjunction] [[bar.F].sub.[infinity]], P x [bar.P]), adapted to the augmented filtration: [F.sub.t] [disjunction] [[bar.F].sub.t] [equivalent to] [sigma]([F.sub.t] [union] [[bar.F].sub.t]). Equation 2) involves

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (5)

where we assume that [F.sub.t], [[bar.F].sub.t] are independent filtrations, as it is reasonable to assume that the risk stemming from the market has no influence on the mortality risk.

Similarly we have in terms of the conditional density h of [tau] (see (2))

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (6)

Thus the insurer's liability at time t is made up of two terms, in correspondence with the benefits at maturity or at decease:

[v.sub.t] = [sub.[??][P.sub.d+t][v.sup.1).sub.t] + [v.sup.2).sub.t], (7)

where

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (8)

For obvious reasons, the term [1.sub.{t[less than or equal to][tau]}] will be dropped in the sequel. For guaranteed unit-linked contracts (GULC) the benefit at maturity depends on the value of the associated portfolio but there is a minimum guaranteed amount if the stock price falls below a fixed level; this can be taken to correspond to the capital accrued at a fixed interest rate [delta], the "technical rate." A fairly common example (see, [12]) is given by the choice [[PSI].sup.1, 2).sub.[partial derivative]] = [X.sub.[partial derivative]] + [X.sub.t]([e.sup.[delta]([partial derivative]-t)] - 1), where contingent on an insurance event happening (maturity [partial derivative] = T or death [partial derivative] = [tau]) the insured receives the stock plus the interest accrued with rate [delta]. It is remarkable that in such situation it is possible to derive a closed expression for the fair value with all generality, independently of the dynamics of the stock. Indeed, using that [X.sub.t]/[B.sub.t] is a martingale we find that the maturity and mortality premiums for such a GULC are given by

[[PSI].sup.1).sub.t] = [X.sub.t](1 + [[B.sub.t]/[B.sub.T]]([e.sup.[delta][??]] - 1)),

[[PSI].sup.2).sub.t] = [X.sub.t] (1 - [sub.[??][P.sub.d+t] + [[integral].sup.T.sub.t] dsh(s) [B.sub.t]/[B.sub.s] ([e.sup.[delta](s-t)] - 1)). (9)

When the interest and mortality rates r and [mu] are constants and [X.sub.t] is assumed to have standard Black and Scholes dynamics one recovers the result of Shen and Xu [12] obtained by PDE techniques.

In a generic case, explicit evaluation of [v.sup.1).sub.t] and [v.sup.1).sub.t] can be a difficult matter. Those cases when [[PSI].sup.2).sub.t] is either a Markov process with a time homogeneous transition function or when it can be represented as [[PSI].sup.2).sub.[tau]] = [[PSI].sub.2]([tau], [X.sub.[tau]]) in terms of a given measurable function [[PSI].sub.2] : [R.sup.2] [right arrow] R are of particular interest; in this case [v.sup.2).sub.t] can be simplified as follows. Let [[theta].sub.s] be the standard time-shift map on the path space acting on sample paths [omega] [member of] [OMEGA] via ([[theta].sub.s] [omicron] [omega]) = [omega](t + s). Then, Markov's property yields

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (10)

3. Valuation of GUL in a Black and Scholes Framework

3.1. Endowment Insurance Contracts. In the sequel we study the valuation of two different types of GULC of actuarial interest. The first case generalizes the result of Brennan and Schwartz [1] and Boyle and Schwartz [2] from a pure endowment case to an endowment insurance contract, where the payoff contingent on death occurring before maturity is given by [[PSI].sup.2).sub.[tau]] = max{[X.sub.[tau]], [X.sub.t][e.sup.[delta]([tau]-t)]}. This means that the initial capital, accrued at a fixed interest rate [delta], is guaranteed. We also study the nonarbitrage price for a floating strike look-back put (FLOP) which entitles the owner to sell the stock at the highest realized price before maturity, that is, to a terminal payoff [[PSI].sup.1).sub.T] = [X.sup.t.sub.T], where [X.sup.t.sub.T] [equivalent to] [max.sub.t [less than or equal to] s [less than or equal to] T] [X.sub.s] is the running maxima of the process. Look back options guarantee a nonloss outcome and as such are interesting in dynamic investment fund protection.

In this section we assume classical stochastic dynamics so that under the martingale probability [X.sub.t] evolves via a geometric Brownian motion. Let T [??] [X.sup.t,x.sub.T] be the solution that starts from x at time t, t < T:

[X.sup.t,x.sub.T] = [X.sub.t] exp{[sigma]([W.sub.T] - [W.sub.t]) + [[integral].sup.T.sub.t] [r.sub.s]ds - [[sigma].sup.2](T - t)/2}, (11)

where we allow for a time varying short term interest rate [r.sub.t], since a life insurance policy could typically be expected to be held for a long time. However, to obtain closed formulas for the premium contingent on death [v.sup.2).sub.T] we assume in this section that r is constant.

The premium for a contract with payoff [[PSI].sup.2).sub.T] = max{[X.sub.T], [X.sub.t][e.sup.[delta](T-t)]} is written at time t such that [X.sub.t] = x follows a minor modification of the classical Black-Scholes-Merton formula: the equality in law [X.sup.t,x.sub.T] = [X.sup.0,X.sub.[??]], where [??] = T-t yields that [v.sup.1).sub.t] [equivalent to] [v.sup.1)](T|t, x) = [v.sup.1)](T - t|0, x), where

[v.sup.1)] (T|0, x) = x ([PHI]([m.sub.+] [square root of (T)]) + [e.sup.-r[??]] [THETA] (-[m.sub.-] [square root of (T)])), (12)

[PHI] is the normal distribution function, and

[m.sub.[+ or -]] [equivalent to] [[??][sigma]] [+ or -] [[sigma]/2], [??] [equivalent to] r - [delta]. (13)

The limit behavior for long values of T is interesting; depending on whether [??] > 0, [??] = 0, or [??] < 0, three different possibilities are found for the premium: [v.sup.1)](T|t, x) tends, respectively, to either x [equivalent to] [v.sup.1)](t | t, x), to 2x, or to [infinity]. The result is easy to understand. The higher the guaranteed rate the more interesting the contract becomes. Further, when [delta] is higher than the market rate r, the discounted benefit tends to [infinity] asymptotically in time and so does the premium.

As commented, it turns out that an analytical formula can be derived for the demise component premium when r and [mu] are independent of time. Here [[PSI].sup.2).sub.[tau]] = max[[X.sub.[tau]], [X.sub.t][e.sup.[delta]([tau]-t)]}. In this case, using (10) we find that [v.sup.2)](T|t, x) = [v.sup.2)](T-t|0, x), where

[v.sup.2)] (T-t|0, x)

= x[mu] [[integral].sup.T.sub.0] ([e.sup.-[mu]s][PHI] ([m.sub.+] [square root of (S)]) + [e.sup.-([mu]+[??])S][PHI] (- [m.sub.-] [square root of (S)])) ds (14)

and hence it involves the integral

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (15)

By interchanging integrals we find I to be given by

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (16)

It follows that the demise contribution to the premium is given by

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (17)

where [eta] [equivalent to] [m.sup.2.sub.+] + 2[mu] and [m.sub.[+ or -]] are defined above (cf. (13)).

Thus, the full premium for a contract paying [[PSI].sub.[partial derivative]] = max{[X.sub.[partial derivative]], [X.sub.t][e.sup.[delta]([partial derivative]-t)]} at the policyholder's demise ([partial derivative] = [tau]) or else at expiry of the policy ([partial derivative] = T) is given by

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (18)

The premium (17) simplifies when the insurance company is committed to pay a technical interest rate [delta] equal to the short rate r. In this case, [m.sub.+] = -[m.sub.-] = [sigma]/2, [eta] = ([[sigma].sup.2]/4) + 2[mu], and (17) yields

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (19)

In Figure 1 we plot [v.sup.2)](T|t, x) as a function of T - t. As expected, the premium vanishes if [mu] = 0. Otherwise it increases towards x(1 + [sigma]/2 [square root of ([eta])]) as the graph shows.

Further, the corresponding full premium (18) can be written in a quite neat way as

[v.sub.t] [equivalent to] [v.sup.1).sub.t] + [v.sup.2).sub.t] = x[1 + [[sigma]/[square root of ([eta])]([PHI]([square root of ([eta][??])])- 1/2)]. (20)

3.2. Floating Strike Lookback Options

3.2.1. Maturity Component ([[PSI].sup.1).sub.T] = [max.sub.0 [less than or equal to] s [less than or equal to] T][X.sub.s]). European floating strike lookback options are a different interesting kind of GULC. For a put option (flop) they entitle the owner to sell the stock at the highest realized price; thus, they imply a terminal payoff [[PSI].sub.[partial derivative]] = [X.sub.[partial derivative]] - [X.sub.[partial derivative]], where [X.sup.t.sub.T] [equivalent to] [max.sub.t [less than or equal to] s [less than or equal to] T] [X.sub.s] is the running maxima of the process. Lookback options guarantee a nonloss outcome and as such are interesting in dynamic investment fund protection.

Let [W.sup.v.sub.t] be BM with drift v and let [M.sub.[??]] [equivalent to] [max.sub.t[less than or equal to]s[less than or equal to]T][W.sup.v.sub.s] be the running maxima. The distribution of the latter can be obtained in terms of the joint density of BM without drift and its running maximum as follows.

For standard BM, [W.sub.t] is obviously [W.sup.t.sub.T] [equivalent to] [max.sub.t[less than or equal to]s[less than or equal to]T][W.sub.s] = [W.sup.0.sub.T-t] [equivalent to] [W.sub.T-t]. If ([W.sup.v.sub.T], [W.sup.v.sub.T]) are BM with drift rate v and its running maxima, their joint density may be retrieved by the known [18] joint density of [W.sub.t] (BM without drift) and [W.sub.T], the Radon-Nikodym transformation dP = exp(v[W.sup.v.sub.[??]] - ([v.sup.2]/2)[??])dQ along with the Cameron-Martin-Girsanov theorem. One finds, upon evaluation of some integrals, that

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (21)

where

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (22)

The distribution of the running maxima of GBM follows from that of the running maxima of BM with drift [M.sub.[??]] by noting that if [X.sup.t,x.sub.s], t < s is the solution to (1) with drift r passing through x at time t then

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (23)

where [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] is increasing and [W.sup.v] is a BM with drift v [equivalent to] r/[sigma] - [sigma]/2. It follows that

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (24)

Therefore,

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (25)

By using (22) and upon tedious integration we obtain that the the premium for a flop contract written at t = 0 and maturity at T is given by

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (26)

where we define now [m.sub.[+ or -]] = r/[sigma] [+ or -] [sigma]/2 (note that they coincide with (13) when [delta] = 0). In Figure 2 we compare this premium with that of the Brennan-Schwartz-Boyle formula (12) corresponding to the case [delta] = 0. Notice how (26) is significantly larger than the premium (12). The humped form is in this case characteristic of the latter but not necessarily of the max-premium.

3.2.2. Floating Strike Options: Demise Component ([[PSI].sup.2).sub.T] = [max.sub.0[less than or equal to]t[less than or equal to][tau]] [X.sub.s]). We now evaluate the demise component to the premium for a Flop contract written at t = 0 when the hazard function p is time independent. Here [[PSI].sup.2).sub.T] = [max.sub.0[less than or equal to]t[less than or equal to][tau]][X.sub.s] and it follows from (10) that

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (27)

where we recall that [m.sub.[+ or -]] [equivalent to] r/[sigma] [+ or -] [sigma]/2, [eta] = [m.sup.2.sub.+] + 2[mu]. By recalling (16) we obtain the fair premium as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (28)

4. Price Dynamics When Supports Are Present

Here we describe a model of price dynamics that incorporate the possibility of the existence of a strong lower support at some constant level c, where 0 < c < [x.sub.0] [equivalent to] [X.sub.0]. To meet this requirement one must modify appropriately the stock dynamics. We assume the existence of a risk neutral probability [P.sup.*] and a process [W.sup.*.sub.t] which is a BM with respect to the latter such that the risk neutral evolution of [X.sub.t] satisfies the SDE:

d[X.sub.t] = [r.sub.t][X.sub.t]dt + b (t, [X.sub.t]) d[W.sup.*.sub.t], [X.sub.0] = [x.sub.0] > c, (29)

where the infinitesimal variance coefficient b(t, x) is to be defined appropriately.

Remark 1. Valuation under exotic dynamics where the stock is driven by a SDE whose variance coefficient b(t, x) is not a linear function of x has often been the subject of financial literature, even as early as 1976. A familiar case is the CEV model of Cox [19], where b(x) = [sigma][x.sup.[rho]] and a and [rho] < 1 are constants. A more complete account of prizing under these dynamics is given in Delbaen and Shirakawa [20]. The extension of CEV models to b(x) having also jump singularities is discussed by Decamps et al. [21]. Our stochastic evolution equation (see (29)) is reminiscent, to some extent, of models describing stochastic interest rate dynamics, such as the classical CIR model [22] where the underlying is related to a Bessel process. For other models in this regard see Schroder [23], Geman and Yor [24], and Goovaerts and de Schepper [25]. Pricing of equity linked products with the stock following some exotic dynamics driven by Levy processes appears in Jaimungal and Young [26].

We find it reasonable to assume that there exists positive probability to attain the boundary; we suppose that this event "triggers" bid orders and hence that [X.sub.t] ricochets upon hitting the boundary. This requirement yields that the point x = c must be what, in the terminology of Feller's boundaries classification, is termed a regular boundary. We also require that for large x, b(t, x) displays a linear dependence on x, as happens with the GBM (1). The obvious choice b(t,x) = [sigma](x - c) is not acceptable since under such dynamics c is a natural, nonattainable barrier. The simplest choice that meets all these requirements is given by taking b(t, x) = [sigma](t)[square root of ([x.sup.2] - [c.sup.2])], where [sigma](t) is an arbitrary function; that is,

d[X.sub.t] = [r.sub.t][X.sub.t]dt + [sigma](t) [square root of ([X.sup.2.sub.t] - [c.sup.2])]d[W.sup.*.sub.t], [X.sub.0] = [x.sub.0] > c. (30)

Note first that under the risk neutral probability [P.sup.*] the deflated process [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] is a martingale. For ease of notation in the sequel we drop the symbol *. Note also that the square root branch point x = c prevents [X.sub.t] to reach the region [0, c) but nevertheless there exists positive probability to attain the barrier x = c.

When r and a are constants these ideas are substantiated by appealing to Feller's theory. Let [[tau].sub.c] be the first hitting time of the barrier. The behaviour of the Feller functions depends on the parameter [??][equivalent to] r/[[sigma].sup.2]. The classical scale and speed functions s(x), m(x) are given by

s(x)= [[integral].sup.x] [([y.sup.2] - [c.sup.2]).sup.-[??]] dy,

[[[sigma].sup.2]/2] m(x) = [[integral].sup.x] [([y.sup.2] - [c.sup.2]).sup.[??]-1] dy. (31)

We also consider the Feller functions

[SIGMA](c, x) [equivalent to] [[integral].sup.x.sub.c] dm(y)(s(y) - s(c)),

[OMEGA](c, x) [equivalent to] [[integral].sup.x.sub.c] dm(y)(s(y) - s(c)). (32)

We see that m(c) and [OMEGA](c, x) are always finite; by contrast s(c) is finite [??] [SIGMA](c, c) is finite [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. Thus if [??] < 1 the point c is a regular boundary. By Feller's test

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (33)

Actually, exploiting further Feller's test we see that

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (34)

We get a sharper result by noting that if [??] [less than or equal to] 1/2 then s([infinity]) = [infinity] which implies [[tau].sub.[infinity]] = [infinity]; thus, Feller's test yields [[tau].sub.c] < [infinity] almost surely [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] whenever [??] [less than or equal to] 1/2.

To determine whether [[tau].sub.c] < [[tau].sub.[infinity]] occurs note that for any [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. Thus letting y [up arrow] [infinity] we have

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (35)

and we conclude that if 0 < [??] < 1 there is positive probability that the support is eventually reached.

We have just seen that if [??] [less than or equal to] 1/2 then [[tau].sub.c] < [infinity] almost surely [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. However even in this case the mean time maybe infinite. We can gain some additional information by noting that if s([infinity]) = [infinity] then

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (36)

We can use this to conclude that [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] if [??] [greater than or equal to] 1/2 and

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (37)

The case [[sigma].sup.2] = 2r is of particular interest since then all Feller functions can be evaluated explicitly as

s(x) = [[[sigma].sup.2]/2] m(x) = log (x + [square root of ([x.sup.2] - [c.sup.2])]/c), [SIGMA](c, x) = [OMEGA](c, x) = [s.sup.2](x)/[[sigma].sup.2]. (38)

Here s and m are the scale and speed functions. Thus, inasmuch as c > 0 we have [SIGMA](c, x) = [OMEGA](c, x) < [infinity] corresponding to a regular boundary.

We now consider the solution of the SDE (29). For general election of the functions r(t) and [sigma](t) the solution remains unknown. However, in the particular case when the time dependent volatility and interest rate satisfy [[sigma].sup.2](t) = 2r(t), then it turns out that the equation admits an analytical solution; we restrict to this situation in the sequel. Indeed, by using Ito's rule, one can prove that a strong solution to (29) satisfying [X.sub.0] = [x.sub.0] is given explicitly by

[X.sub.t] = c cosh ([[infinity].sup.t.sub.0] [square root of (2[r.sub.s])]d[W.sub.s] + v), (39)

where we set

v [equivalent to] [cosh.sup.-1] ([x.sub.0]/c), [phi](t) [equivalent to] 2 log [B.sub.t] = 2 [[integral].sup.t.sub.0] [r.sub.l]dl, [phi] (t, T) [equivalent to] [phi](T) - [phi](t) = 2 [[integral].sup.T.sub.t] [r.sub.l]dl. (40)

Further in terms of a new Brownian motion [[??].sub.t] we have that a weak solution is given by

[X.sub.t] = c cosh [Y.sub.t] where [Y.sub.t] = v + [[??].sub.[phi](t)]. (41)

This follows by noting that [[integral].sup.t.sub.0] [square root of (2[r.sub.s])] d[W.sub.s] is a local martingale and Levy's representation theorem.

In particular these expressions clarify the behavior of the process upon hitting the boundary: [X.sub.t] attains the barrier c whenever the process [Y.sub.t] reaches 0 upon which [X.sub.t] is reflected. In Figure 3 a realization of the stock evolution (39) corresponding to a constant interest rate has been plotted. Notice how the stock eventually hits the support level c = 4 several times.

Explicit valuation formulae under dynamics (see (39)-(42)) could prove quite awkward to obtain. Here we determine the fair premium under these dynamics for several pure endowment contracts, with general deterministic interest rate [r.sub.t].

4.1. GULC of Type (i) ([[PSI].sup.1).sub.T] = max{[X.sub.T], [X.sub.t][e.sup.[delta][??]]}). Let [X.sup.t,x.sub.T], t < T be the solution to (29) at time T that starts from x at time t. We have

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (42)

where we recall that

[zeta] [equivalent to] [cosh.sup.-1] (x/c), [phi](t, s) [equivalent to] 2 [[integral].sup.S.sub.t] [r.sub.l]dl. (43)

Notice that in the sequel we simply write [phi](t, T) [equivalent to] [phi]. It then follows from (9) that whenever the benefit can be written as [[PSI].sup.1).sub.T] = [[PSI].sub.1]([X.sub.T]) for some [[PSI].sub.1]: R [right arrow] R, the fair price is given by

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (44)

In our case [[PSI].sub.1]([X.sub.T]) = max{[X.sub.T], [X.sub.t][e.sup.[delta][??]]} and evaluating the integral (45) we obtain that the fair price for a Brennan-Schwartz-Boyle contract under dynamics (29) is given by

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (45)

where we introduce

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (46)

This expression generalizes the Brennan-Schwartz-Boyle formula (12) to the case when barriers are present. There are interesting differences. Note first that, unlike what happens with (12), [v.sup.1).sub.t] [equivalent to] [v.sup.1)](T|t, x) does not grow linearly with x. Further, for long values of [??], [v.sup.1)] tends either to [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], or to [infinity] depending on whether [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. However convergence takes place with a slower rate than what happens when no resistance is present.

In Figure 4 we plot the maturity premium [v.sup.1.sub.t] in terms of the initial stock price x if the barrier is located at c = 4 and compare with the result obtained when no barrier is present. Notice how (12) (dotted line) always overprices the premium compared with (45) given by the thick solid line. This reflects the fact under the actual dynamics that the exposure of the company to a possible plummeting of the stock is partially hedged by the existence of the barrier. For long x the difference between both expressions grows dimmer. This is easy to understand, since the protection from the barrier diminishes with the distance to the starting point.

The dependence in time of (45) is displayed in Figure 5. Notice how differences between (12) and (45) can be quite marked for moderate maturity times.

Setting c = 0 amounts to having no barrier. In this case one has [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], and M = -[infinity]; most of the terms in (45) drop out and, as expected, we recover (12).

4.2. GULC of Type (ii) ([[PSI].sup.1).sub.T] = [max.sub.t[less than or equal to]s[less than or equal to]T][X.sub.s]). We next consider a contract linked to the evolution of the maximum value of the stock, that is, wherein [[PSI].sup.1).sub.T] = [max.sub.t[less than or equal to]s[less than or equal to]T][X.sub.s]. Accordingly, the fair price of such contract involves the distribution of [max.sub.t[less than or equal to]s[less than or equal to]T] cosh([zeta] + [[integra].sup.S.sub.t] [square root of 2[r.sub.l])]d[W.sub.l]. This entails important difficulties to derive the fair prize accountable to the fact that under dynamics (29), (24) does not hold since [??](z) [equivalent to] c cosh(z + [zeta]) is not increasing. The max distribution can be given in terms of the survival probability for Brownian motion as we now show. Let [M.sub.t] and [m.sub.t] denote respectively running maximum and minimum of BM; then, from (41) one can prove the following equality in law:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (47)

Therefore we have in terms of the joint density f(z, y; t) of ([m.sub.t], [M.sub.t]) that (see (41)-(43))

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (48)

where

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (49)

Recall that the survival function G(z, y; t), or probability that the BM [W.sub.t] has remained within the interval [z, y] up to time t, is given by (see [27] or 28]):

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (50)

where

[a.sub.n1] [equivalent to] (2n + 1) y - 2nz, [a.sub.n2] [equivalent to] (2n + 1)y-2(n+ 1) z, [a.sub.n3] [equivalent to] 2ny - (2n- 1) z, [a.sub.n4] [equivalent to] 2ny - (2n+ 1) z. (51)

By differentiation the joint density f of ([m.sub.t], [M.sub.t]) follows. One could proceed by substituting this f into (49). However this yields a quite messy expression. A more convenient approach is to write (49) in terms of G as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (52)

Note that (50) implies that [[partial derivative].sub.y] G[|.sub.0,y] = [[partial derivative].sub.z]G[|.sub.z,0] = 0. Hence upon appropriate manipulations we find that both terms [I.sub.1] and [I.sub.2] satisfy

[I.sub.1] = [I.sub.2] = [[integral].sup.[infinity].sub.[zeta]] dy cosh y[[partial derivative].sub.y]G(-[zeta] -[zeta]; [phi]). (53)

This observation allows us to simplify the conditional expectation to

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (54)

where [p.sub.t](x) is BM density and we introduce

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (55)

A further simplification is obtained noting that [y.sup.+.sub.n] = -[[??].sup.-.sub.-n] whereupon we have

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (56)

The latter integral can be evaluated in terms of the error function [PHI]; then, upon substitution into (9), we obtain that if [[sigma].sup.2] = 2[r.sub.t] and the short rate [r.sub.t] depends on time in an arbitrary way, then the value of the fair premium is given by

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (57)

where

[s.sub.n] [equivalent to] sign [[alpha].sub.n], [[lambda].sup.[epsilon],[kappa].sub.n] [equivalent to] [epsilon][square root of ([phi])/[[alpha].sub.n] - [([[alpha].sub.n] + [kappa])/[square root of ([phi])]] [zeta]. (58)

The dependence of the premium on time and initial value enters through the variables (49). For a time to maturity ranging from moderate to large only those terms corresponding to n = 0 will contribute significantly. It is then natural to write

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (59)

This representation makes it possible to recover (26) when c [right arrow] 0 and r is constant. Note that in this case [phi] = 2r[??]; further, if n [not member of] 0 then [s.sub.n][[lambda].sub.n] [right arrow] -[infinity] and [PHI]([s.sub.n][[lambda].sup.[epsilon],[kappa].sub.n]) [right arrow] 0. Likewise we have

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (60)

as c [right arrow] 0. Thus, the only nonvanishing term corresponds to [epsilon] = -[kappa] = 1, n = 0, and

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (61)

which is (26) in the case [sigma] = [square root of (2r)].

5. Kolmogorov Backwards Equation for the Premium

In this section we will suppose that the insurer's benefit can be represented as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (62)

Further, we assume that [[PSI].sub.i] : [R.sup.2] [right arrow] R is a piecewise continuous function that does not have a parametric dependence on the "initial value" x [equivalent to] [X.sub.t]. Under these assumptions and conditional on {[tau] > t} having happened, we next derive a PDE that the premium [v.sub.t] [equivalent to] v (T|t, x) satisfies. Notice that, from (9),

[[partial derivative].sub.t] ([B.sub.T][v.sup.1).sub.t]/[B.sub.t]) = [sub.[??] [P.sub.d+t] [[mu] (t + d, t) [PSI] (t, x) + [[partial derivative].sub.t] [psi] (t, x)], (63)

where [psi](t, x) [equivalent to] [E.sup.P]([[PSI].sub.1](T, [X.sub.T]) | [X.sub.t] = x). Recall that if [X.sub.t] is a diffusion with drift coefficient a(t, x) [equivalent to] [r.sub.t]x and variance coefficient b(t, x) then the conditional expectation solves Kolmogorov equation:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (64)

Operating with the infinitesimal generator A on (9) and using (63) we find

(A -[r.sub.t] - [mu](t + d, t)) [v.sup.1).sub.t] = 0. (65)

Likewise we find that

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (66)

Operating again with operator A we get

(A - [r.sub.t] - [mu](t + d, t)) [v.sup.2).sub.t] = -[mu](t + d, t) [[PSI].sub.2](t, x). (67)

It follows that v [equivalent to] [v.sup.1)] + [v.sup.2)] satisfies the backward PDE with final condition

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (68)

where a(t, x) [equivalent to] [r.sub.t]x and b(t, x) = [sigma](t) [square root of ([x.sup.2] - [c.sup.2].

Remark 2. (1) If r and [mu] are constants, there is no resistance: c = 0 and, in addition, [[PSI].sub.2](t, x) = x, then 68) reduces to the equation considered by Shen and Xu [12].

(2) If representation (62) holds but [[PSI].sup.i)] depends parametrically on the "initial values and times": [[PSI].sup.i)] = [[PSI].sub.I](T, [X.sub.T]; t, x), then it does not exist such a clear-cut equationas (68). Notice that this is precisely what happens for the type i GULC considered in this paper. Similarly, the case of path-dependent payoff functions (type ii GULC) may not be covered either with the PDE approach.

http://dx.doi.org/10.1155/2014/314286

Conflict of Interests

The author declares that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The author acknowledges support from the direcion general de investigation cientifica y tecnica under Contracts MTM2009-09676 and MTM2012-38445.

References

[1] M. J. Brennan and E. S. Schwartz, "The pricing of equity-linked life insurance policies with an asset value guarantee," Journal of Financial Economics, vol. 3, no. 3, pp. 195-213, 1976.

[2] P. P. Boyle and E. S. Schwartz, "Equilibrium prices of guarantees under equity-linked contracts," The Journal of Risk and Insurance, vol. 44, no. 4, pp. 639-660, 1977.

[3] A. R. Bacinello and F. Ortu, "Pricing equity-linked life insurance with endogenous minimum guarantees," Insurance: Mathematics and Economics, vol. 12, no. 3, pp. 245-257, 1993.

[4] K. K. Aase and S. A. Persson, "Pricing of unit-linked life insurance policies," Scandinavian Actuarial Journal, vol. 1, pp. 26-52, 1994.

[5] M. J. Brennan and E. S. Schwartz, "Alternative investment strategies for the issuers of equity linked life insurance policies with an asset value guarantee," Journal of Business, vol. 52, pp. 63-93, 1979.

[6] S. Ekern and S.-A. Persson, "Exotic unit-linked life insurance contracts," The GENEVA Papers on Risk and Insurance Theory, vol. 21, no. 1, pp. 35-63, 1996.

[7] P P Boyle and M. R. Hardy, "Reserving for maturity guarantees: two approaches," Insurance: Mathematics and Economics, vol. 21, no. 2, pp. 113-127, 1997

[8] A. Grosen and P L. Jorgensen, "Fair valuation of life insurance liabilities: the impact of interest rate guarantees, surrender options, and bonus policies," Insurance: Mathematics and Economics, vol. 26, no. 1, pp. 37-57, 2000.

[9] T. Moeller, "Hedging equity-linked life insurance contracts," North American Actuarial Journal, vol. 5, no. 2, pp. 79-95, 2001.

[10] C. Bernard, O. Le Courtois, and F. Quittard-Pinon, "Market value of life insurance contracts under stochastic interest rates and default risk," Insurance: Mathematics and Economics, vol. 36, no. 3, pp. 499-516, 2005.

[11] A. R. Bacinello, "Endogenous model of surrender conditions in equity-linked life insurance," Insurance: Mathematics and Economics, vol. 37, no. 2, pp. 270-296, 2005.

[12] W. Shen and H. Xu, "The valuation of unit-linked policies with or without surrender options," Insurance: Mathematics and Economics, vol. 36, no. 1, pp. 79-92, 2005.

[13] F. Black and M. Scholes, "The pricing of options and corporate liabilities," Journal of Political Economy, vol. 81, no. 3, pp. 637-659, 1973.

[14] R. C. Merton, "Theory of rational option pricing," The Bell Journal of Economics and Management Science, vol. 4, no. 1, pp. 141-183, 1973.

[15] P. Marocco and E. Pitacco, "Longevity risk and life annuity reinsurance," in Proceedings of the Transactions of the 26th International Congrese of Actuaries, vol. 6, pp. 453-479, Birmingham, England, 1998.

[16] A. Melnikov and Y. Romaniuk, "Evaluating the performance of Gompertz, Makeham and Lee-Carter mortality models for risk management with unit-linked contracts," Insurance: Mathematics and Economics, vol. 39, no. 3, pp. 310-329, 2006.

[17] J. M. Harrison and S. R. Pliska, "Martingales and stochastic integrals in the theory of continuous trading," Stochastic Processes and their Applications, vol. 11, no. 3, pp. 215-260, 1981.

[18] P. Levy, Processus Stochastiques et Movement Brownien, Gauthier Villars, Paris, France, 1948.

[19] J. C. Cox, "The constant elasticity of variance option pricing model," Journal of Portfolio Management, vol. 23, pp. 15-17, 1996.

[20] F. Delbaen and H. Shirakawa, "A note on option pricing for the constant elasticity of variance model," Asia-Pacific Financial Markets, vol. 9, no. 2, pp. 85-99, 2002.

[21] M. Decamps, A. de Schepper, and M. Goovaerts, "Applications of d-function perturbation to the pricing of derivative securities," Physica A, vol. 342, no. 3-4, pp. 677-692, 2004.

[22] J. Cox, J. Ingersoll, and S. Ross, "A theory of the term structure of interest rates," Econometrica, vol. 53, no. 2, pp. 385-408, 1985.

[23] M. Schroder, "Computing the constant elasticity of variance option pricing," Journal of Finance, vol. 44, pp. 211-219, 1989.

[24] H. Geman and M. Yor, "Bessel processes, Asian options and perpetuities," Mathematical Finance, vol. 3, no. 4, pp. 349-375, 1993.

[25] M. Goovaerts and A. de Schepper, "IBNR reserves under stochastic interest rates," Insurance: Mathematics and Economics, vol. 21, no. 3, pp. 225-244, 1997.

[26] S. Jaimungal and V. R. Young, "Pricing equity-linked pure endowments with risky assets that follow Levy processes," Insurance: Mathematics and Economics, vol. 36, no. 3, pp. 329-346, 2005.

[27] D. Freedman, Brownian Motion and Difussions, Springer, Berlin, Germany, 1971.

[28] D. Revuz and M. Yor, Continuous Martingales and Brownian Motion, Springer, Berlin, Germany, 1991.

Javier Villarroel

Correspondence should be addressed to Javier Villarroel; javier@usal.es

Received 30 August 2013; Accepted 11 November 2013; Published 11 February 2014

Academic Editors: T. Prieto-Rumeau and O. Valero
No portion of this article can be reproduced without the express written permission from the copyright holder.

Title Annotation: Printer friendly Cite/link Email Feedback Research Article Villarroel, Javier The Scientific World Journal Report Jan 1, 2014 7538 An observing system simulation experiment for the western North Pacific region. The effect analysis of strain rate on power transmission tower-line system under seismic excitation. Contracts Endowments Insurance Stock exchanges Stock-exchange