In this paper we present the Mellin transform method for the valuation of the American power put option with non-dividend and dividend yields, respectively. We use the Mellin transform method to derive the integral representations for the price and the free boundary of the American power put option. We also extend our results to derive the free boundary and the fundamental analytic valuation formula for perpetual American power put option which has no expiry date. Numerical experiments have shown that the Mellin transform method is a better alternative technique compared to the binomial model (BSM), recursive method (RM) and finite difference method (FDM) for the valuation of the American power put option. In general, the Mellin transform method is accurate, flexible and produces accurate prices for the optimal exercise boundary of the American power put option for a wide range of parameters. Hence the Mellin transform method is mutually consistent and agrees with the values of the analytic option valuation formula called the “Black-Scholes model”.

American Power Option Dividend Yield Mellin Transform Method Non-Dividend Yield Perpetual Power Put Option
1. Introduction

Option valuation has been studied extensively in the last three decades. Many problems in financial mathematics entail the computation of a particular integral. In many cases these integrals can be valued analytically and in some cases they can be computed using a partial differential equation, or valued using numerical integration.

Power option is defined as a contingent claim on the product of powers of several underlying assets. The holder has either the right, but not the obligation to buy, as in the case of the power call option, or the possibility to sell, as in the case of the power put option, an asset for a certain price at a prescribed date in the future. The difference between the American and the European power options is that the European power option can only be exercised at the maturity or expiry date while the American power option can be exercised by its holder at any time before the expiry date. This early exercise feature makes the valuation of the American power option mathematically challenging and therefore, creats a great field of research.

A perpetual American power option is an option that has no expiry date. In other words, this type of power option never expires. In a special case of a plain vanilla perpetual option, a closed form solution for the free boundary and price of the American put was derived by  .

Mellin transforms in option theory were introduced by  ,  extended the results obtained in  and showed how the Mellin transform approach could be used to derive the valuation formula for the perpetual American put options on dividend-paying stocks.  considered the Mellin transform method for the valuation of some vanilla power options with non-dividend yield. They derived the fundamental valuation formula known as the Black-Scholes model using the convolution property of the Mellin transform method. The analytical valuation of the American options was considered by  . An alternative approach to the valuation of American options and applications was considered by  .

For the mathematical background of the Mellin transform method in derivatives valuation see  - , just to mention few. In this paper, we focus on the Mellin transform method for the valuation of the American power put option with non-dividend and dividend yields, respectively, and its extension to power option which has no expiry date, i.e. “perpetual American power put option”. The rest of the paper is structured as follows: in Section 2, we present American power options and the payoffs for power call and put options. Section 3 presents the Mellin transform method for the valuation of the American power put option. Section 4 considers the extension of the Mellin transform method to the valuation of the perpetual American power put option. In Section 5, we present some numerical experiments. Section 6 concludes the paper.

2. American Power Options

The power options can be seen as a class of options in which the payoff at expiry is related to the power of the underlying price of the asset. American power options are options that can be exercised before or at the expiry date with non-linear payoff. The American power option comes in two forms, namely, the American power call option and the American power put option. The American power call option is an option with non- linear payoff given by the difference between the price of the underlying asset at maturity raised to a strictly positive power and the exercise price. The American power put option is an option with non-linear payoff given by the difference between the exercise price and price of the underlying asset at maturity raised to a strictly positive power. For an American power option on the underlying price of the asset with exercise price K and time to expiry T, we have the payoffs for the American power call and put options as

and

respectively.

Remark 1

・ For, the payoffs for American power call and put options in (1) and (2) become the payoffs for plain American call and put options, i.e.

and

respectively.

3. The Mellin Transform Method for the Valuation of the American Power Put Option

There are many methods for the valuation of the American power option leading to different but equivalent mathematical formulations. We consider the derivation of the integral representation for the price of the American power put option and the integral equation to determine the free boundary of the American power put option via the Mellin transform method for the case of both non-dividend and dividend yields.

3.1. American Power Put Option with Non-Dividend Yield

Consider the non-homogeneous Black-Scholes partial differential equation for the American power put option with non-dividend yield given by

where the early exercise function f defined on is given by

The final time condition given by on is called the high contact condition. The other boundary conditions are given by

The free boundary is determined by the smooth pasting conditions given by

and

Applying the Mellin transform to (5), we have that

Setting and. Then (11) yields

The Mellin transform of the early exercise function in (12) is obtained as

Solving further, we have the particular solution of (12) as

The complementary solution to the left hand side of (12) is obtained as

where is the integration constant obtained as

is the Mellin transform of the final time condition and is given by

Using (16) and (17) in (15) we have that

Hence the general solution to (12) is given by

The Mellin inversion of (19) is obtained as

where and

Remark 2

・ The first term in (20) is the integral representation for the price of the European power put option (stems from the minimum guaranteed payoff of the American power put) which pays no dividend yield (see  ). The second term in (20) is called the early exercise premium (the value attributable to the right of exercising the option early) for the American power put option with non dividend yield denoted by. Therefore (20) becomes

where

・ Setting in (21) and using the smooth pasting conditions given by (9) and (10), we have the integral representation for the free boundary of the American power put option with non-dividend yield as

where

We formalized the properties highlighted in Remark 2 in the following results.

Theorem 1 The American power put option which pays no dividend yield satisfies the decom- position

where and, and

Theorem 2 Using the smooth pasting conditions given by and, the free

boundary formulation of the American power put option with non-dividend yield is given by

3.2. American Power Put Option with Dividend Yield

The derivation of the integral representation for the price of the American power put option which pays dividend yield using the Mellin transform method is given in the following result.

Theorem 3 Let be the price of the underlying asset, K be the strike price, r be the risk interest rate, q be the dividend yield and T be the time to maturity. Assume yields dividend, then the integral representation for the price of the American power put option is given by

Proof. Consider the non-homogeneous Black-Scholes partial differential equation for the American power put option with dividend yield given by

where

on and the free boundary of the American power put option with dividend yield. The high contact condition is given by on. The other con- ditions are given by (7) and (8). With the smooth pasting conditions given by

and

The Mellin transform of (24) gives

Putting and, (28) yields

where

Following the same procedures for the case of non-dividend yield, the general solution to (30) is obtained as

The Mellin inversion of (31) is given by

Equation (32) is the integral representation for the price of American power put option with dividend yield, where and

Remark 3

・ The first term in (32) is the integral representation for the price of the European power put option (stems from the minimum guaranteed payoff of the American power put) with dividend yield and the last two terms denote the early exercise premium (the value attributable to the right of exercising the option early) for the American power put option with dividend yield denoted by. Therefore (32) becomes

where

・ Setting in (33) and using the smooth pasting conditions given by (26) and (27), we have the integral representation for the free boundary of the American power put option with dividend yield as

where

From Remark 3, we have the following results.

Theorem 4 The American power put option which pays dividend yield satisfies the decomposition

where and, and

Theorem 5 Using the smooth pasting conditions given by and. Then

the free boundary formulation of the American power put option with dividend yield is given by

The following results present some special cases of (20) and (32).

Theorem 6 If and, then

(i) The integral representation for the American power put option which pays no dividend yield (20) reduces to the integral equation derived by Kim  for the price of the plain American put option given by

where

(ii) The free boundary for the American power put option which pays no dividend yield (22) reduces to the integral equation derived by Kim  for the price of the plain American put option given by

where

Proof. Setting and in (20) yields

where and. Equation (39) can be be written as

where and denote the price of the European put option with no dividend yield and free boundary for the American put option with no dividend yield respectively. Let

where

The early exercise function is given by

and

Using the convolution property of the Mellin transform, (42) becomes

Substituting the value of the early exercise function from (43) and

into (45), we have that

Using the transformation given by

(47) becomes

Substituting (49) into (41) we have the early exercise premium for the American put option with non-dividend yield as

where

Setting, then (50) becomes

where

Substituting (52) into (40) we get the integral equation (35) obtained by Kim  as

Hence (i) is established. For the second reduction, setting in the last integral equation above and

using the smooth pasting conditions given by and, we obtain the free boun-

dary of the American put option which pay no dividend yield (37) derived by Kim  as

where

Theorem 7 If, then the optimal exercise boundary of the American power put option with with dividend yield is given by

Proof. Let and, (34) becomes

where and. Factorizing and rearranging, (54) becomes

where

and

Notice first that critical stock price is bounded from above i.e.. Taking the limits of (56) and (57) as, we have that

and

respectively. If

We have

Using (63), the limit of (55) is obtained as

Since

and

Then (64) becomes

If

We have that

The first integral can also be written as

Applying the residue theorem of complex number given by

Then the inner integral in (67) becomes

Substituting (69) into (67) yields

Similarly,

Substituting (70) and (71) into (66) for, we have that

Using the l’Hospital rule, for, (64) becomes

Combining (72) and (73)

Hence (53) is established.

Remark 4

The above results confirm the formula of Kim and Yu 

Theorem 8 If the underlying asset price follows a lognormal diffusion process and the interest rate is a positive constant, then the optimal exercise boundary of the American power put option with at maturity is given by

Proof. Let. In order to investigate the behaviour of the optimal exercise boundary of the American power put option with near maturity, we consider (55) which is of the form

If, the limit of the right hand side of (55) as can be evaluated using the l’Hospital’s rule we have that

If, the limit of the right hand side of (55) as is obtained directly as

Combining (75) and (76), we have the optimal exercise boundary of the American power put option with at maturity given by

Hence (74) is established.

Remark 5

From (75), we notice that when the American put can have a positive value at expiration given that it has not been exercised earlier. This indicates that large dividend payouts reduce the incentives of early exercise.

From (76), we deduce that when the American put will have a zero payoff at expiration even if it has not been exercised earlier. This is because it is not possible for the underlying asset price at expiration to fall below K without crossing the exercise boundary at an earlier time.

Theorem 9 The integral representation for the price of the American power put option which pays dividend yield given by (32) can be reduced to integral representation derived by Kim  .

where

Proof. Setting, then (32) becomes

where and, and

Using the procedures of  , (78) can be written as

with the Mellin transforms of and given by

respectively. Using the convolution theorem of the Mellin transform we have that

The price of the American power put option which pays dividend yield can be expressed as

The integral is evaluated as follows

where and. Using the following variables

transformation given by

and

For the first and second integrals in (83) respectively, we have that

Substituting (84) into (82) yields

By changing, (85) becomes

Hence by setting, this proves (77).

4. Application of the Results to Perpetual American Power Put Option Valuation

Now, we apply the results generated for the integral equations in (20) and (32) to power options which have no expiry date. The following results shows the derivation of the expression for the free boundary of perpetual the American power put option and its closed form solution for both non-dividend and dividend yields, using the Mellin transform method.

Theorem 10 (Non-Dividend Yield) If and, then the free boundary of the perpetual American power put option is given by

and the price of the perpetual American power option becomes

where

Proof. The integral representation for the price of the American power put option which pays no dividend yield given by (20) can be expressed as

where and are given by

with

and

respectively. For (90) to hold as, it is necessary that i.e., where is given by (89). The second smooth pasting condition (10) for a perpetual power put can be written as

Differentiating (91) at we have that

where

As and therefore

Also differentiating (92) yields

Taking the limit of (97) as, we have

Therefore,

where with. The limiting cases and are the roots of. Hence (99) becomes

Since, application of the residue theorem given by (68) leads to

Substituting (96) and (101) into (93) yields

Equation (102) is the expression for the free boundary of a perpetual American power put option. Next, we use (102) to derive an expression for the price of perpetual American power put option. Note that the price of a perpetual European power put option is zero, since it can never be exercised. Therefore, taking the limit as in (90), the price of perpetual American put option for is given by

where. Integrate the inner integral, (102) becomes

Once again we apply the residue theorem (68) to get

Equation (105) is the price of a perpetual American power put option obtained as a limit of the price of a finite-lived American power put option.

Theorem 11 (Dividend Yield) If and, then the free boundary of the perpetual American power put option is given by

and the price of perpetual American power put option equals

where

and

Proof. The integral representation for the price of the American power put option which pays dividend yield given by (32) can be expressed as

where, and are given by

with

respectively. The roots of are

and

Thus we write that

with. For (110) to hold as, it is necessary that i.e., where and are given by (109). The second smooth pasting condition (10) for a per-

petual power put which pays dividend yield can be written as

Differentiating (111) at we have that

where

As and therefore

Now differentiating (112) w.r.t and taking the limit we have that

Therefore, by setting we have that

Let and

In the same manner, setting and differentiating (113) w.r.t, we have that

Therefore,

Setting and solving (121) further we have that

Once again by the application of residue theorem (68), then (120) and (122) yield

and

respectively. Substituting (117), (123) and (124) into (114), we obtain

Equation (125) is called the free boundary of the perpetual American power put option which pays dividend yield.

The price for the perpetual American power put option is given by

Using the residue theorem (68), then (126) becomes

5. Numerical Experiments

In this section we present some numerical experiments and discussion of results.

Experiment 1

We consider the valuation of the American power put option for which pays dividend yield with the following parameters:

The result generated is shown in Table 1 below.

Experiment 2

We consider the valuation of the American power put option for which pays non-dividend yield with the following parameters:

The results generated for the price of the American power put option via Black-Scholes model (BSM), binomial model (BM) and the Mellin transform method (MTM) are shown in Tables 2-4 below. Also the results generated for the free boundary of the American power put option are shown in Tables 5-7 below.

Experiment 3

We consider the valuation of the American Power put option with the following parameters:

The comparative results analysis of the Mellin transform method (MTM) in the context of Black-Scholes model (BSM), binomial model (BM), recursive method (RM) and Finite difference method (FDM) are shown in Table 8.

American power put values
1.901.952.02.052.10
0.118.27410.2894.3541.3090.275
0.1518.99712.1476.8093.3161.403
0.2020.16014.1029.1755.5483.125
0.2521.53516.05811.4537.8325.129
0.3023.00817.98113.65310.0777.251
KBlack-Scholes Model (BSM)Binomial Model (BM)Mellin Transform Method (MTM)
350.0060.0060.007
400.8400.8510.852
454.8405.0005.031
KBlack-Scholes Model (BSM)Binomial Model (BM)Mellin Transform Method (MTM)
350.0760.0760.078
401.2951.3101.310
454.9755.0515.058
KBlack-Scholes Model (BSM)Binomial Model (BM)Mellin Transform Method (MTM)
350.2440.2450.247
401.7531.7661.768
455.2315.2855.300
Strike Price, KUnderlying Asset Price, Free Boundary
354031.740
404036.273
454040.810
Strike Price, KUnderlying Asset Price, Free Boundary
354029.783
404034.040
454038.291
Strike Price, KUnderlying Asset Price, Free Boundary
354027.850
404031.830
454035.804
The comparative results analysis of some numerical methods for the valuation of American power put option
KcnrTBSMBMMTMRMFDM
35400.4210.050.0832.1042.1442.1572.1602.168
40400.4210.050.0834.2324.3304.3544.3704.357
45400.4210.050.0837.1447.3647.3847.3907.380
Discussion of Results

From Figure 1 below, we observe that the higher the volatility, the higher the values of the American power put option. Also the higher the power of the American put option, the lower the values of the option. Figures 2-4 below show that the Mellin transform method is mutually consistent, performs very well, accurate and agrees with the values of Black-Scholes model (BSM). In Figure 5 below, we plot the free boundary as a function of the strike price K for different values of volatility. We observe that the higher the volatility, the lower the optimal exercise boundary of the American power put option. Figure 6 below demonstrates that the Mellin transform method is a better alternative technique compared to the Black-Scholes model (BSM), binomial model (BM), recursive method (RM) and finite difference method (BM) for the valuation of the American power put option. Hence the Mellin transform method is a good technique for the valuation of the American power put option.

6. Conclusion

In this paper, we have derived the integral representations for the price and the free boundary of the American power put option for non-dividend and dividend yields using the Mellin transform method. We also extended the integral equation for the price of the American power put option to derive the expression for the free boundary and the price of the perpetual American power put option which pays both non-dividend and dividend yields as the limit of a finite-lived option by means of smooth pasting condition. In general, numerical experiments have shown that the Mellin transform method is accurate, flexible, efficient and produces accurate prices for the optimal exercise boundary for a wide range of parameters.

American power put option values using <xref ref-type="table" rid="table1">Table 1</xref> The comparative results analysis using <xref ref-type="table" rid="table2">Table 2</xref> The comparative results analysis using <xref ref-type="table" rid="table3">Table 3</xref> The comparative results analysis using <xref ref-type="table" rid="table4">Table 4</xref> The free boundaries of American power put option with n = 1 The comparative results analysis using <xref ref-type="table" rid="table8">Table 8</xref>
Cite this paper

Sunday EmmanuelFadugba,Chuma RaphaelNwozo, (2015) Mellin Transform Method for the Valuation of the American Power Put Option with Non-Dividend and Dividend Yields. Journal of Mathematical Finance,05,249-272. doi: 10.4236/jmf.2015.53023

ReferencesSamuelson P.A. ,et al. (1965)Rational Theory of Warrant Pricing Industrial Management Review 6, 13-31.Panini, R. and Srivastav, R.P. (2005) Pricing Perpetual Options Using Mellin Transforms. Applied Mathematics Letters, 18, 471-474. http://dx.doi.org/10.1016/j.aml.2004.03.012Frontczak, R. and Sch&ouml;bel, R. (2008) Pricing American Options with Mellin Transforms. Working Paper.Nwozo, C.R. and Fadugba, S.E. (2014) Mellin Transform Method for the Valuation of Some Vanilla Power Options with Non-Dividend Yield. International Journal of Pure and Applied Mathematics, 96, 79-104. http://dx.doi.org/10.12732/ijpam.v96i1.7Kim, I. (1990) The Analytic Valuation of American Options. The Review of Financial Studies, 3, 547-572. http://dx.doi.org/10.1093/rfs/3.4.547Kim, I. and Yu, G.G. (1996) An Alternative Approach to the Valuation of American Options and Applications. Review of Derivatives Research, 1, 61-85. http://dx.doi.org/10.1007/BF01536395Erdelyi, A., Magnus, W., Oberhettinger, F. and Tricomi, F. (1954) Tables of Integral Transforms, Vol. 1-2. McGraw-Hill, New York.Fadugba, S.E. and Nwozo, C.R. (2015) Integral Representations for the Price of Vanilla Put Options on a Basket of Two-Dividend Paying Stocks. Applied Mathematics, 6, 783-792. http://dx.doi.org/10.4236/am.2015.65074Frontczak, R. and Sch&ouml;bel, R. (2009) On Modified Mellin Transforms, Gauss-Laguerre Quadrature and the Valuation of American Call Options. Tübinger Diskussionsbeitr&uuml;ge, No. 320.Gradshteyn, I. and Ryshik, I. (2007) Table of Integrals Series and Products. 7th Edition, Academic Press, Waltham.Panini, R. and Srivastav, R.P. (2004) Option Pricing with Mellin Transforms. Mathematical and Computer Modelling, 40, 43-56. http://dx.doi.org/10.1016/j.mcm.2004.07.008Vasilieva, O. (2009) A New Method of Pricing Multi-Options Using Mellin Transforms and Integral Equations. Master’s Thesis in Financial Mathematics, School of Information Science, Computer and Electrical Engineering, Halmstad University, Halmstad.Nguyen, T.H. and Yakubovich, S.B. (1991) The Double Mellin-Barnes Type Integrals and Their Applications to Convolution Theory. World Scientific, Series on Soviet Mathematics, Vol. 6, Singapore, 295p.Zemanian, A.H. (1987) Generalized Integral Transformations. Dover Publications, New York.Zieneb, A.E. and Rokiah, R.A. (2011) Analytical Solution for an Arithmetic Asian Option Using Mellin Transforms. International Journal of Mathematical Analysis, 5, 1259-1265.