MELLIN TRANSFORM METHOD FOR THE VALUATION OF SOME VANILLA POWER OPTIONS WITH NON-DIVIDEND YIELD Chuma

This paper presents the Mellin transform method for the valuation of some vanilla power options with non-dividend yield. This method is a powerful tool used in the valuation of options. We extend the Mellin transform method proposed by Panini R. and Srivastav R.P. [15] to derive the price of European and American power put options with non-dividend yield. We also derive the fundamental valuation formula known as the Black-Scholes model using the convolution property of the Mellin transform method. To provide a sufficient numerical analysis, we compare the results generated by the Mellin transform method for the valuation of American power put option for n = 1 which pays no dividend yield to two other numerical methods namely Crank Nicolson finite difference method [2] and binomial model [3] for options valuation against Black-Scholes analytical pricing formula [1]. The numerical experiment shows that the Mellin transform method is efficient, easy to implement, agree with the values of Black-Scholes [1], Crank Nicolson finite difference method [2] and binomial model [3]. Hence the Mellin transform method is a better alternative method compared to the Crank Nicolsion finite difference and binomial model for the valuation of some vanilla power options. Received: June 2, 2014 c © 2014 Academic Publications, Ltd. url: www.acadpubl.eu §Correspondence author 80 C.R. Nwozo, S.E. Fadugba AMS Subject Classification: 44A15, 44A99, 60H30, 91G20, 91G99


Introduction
The derivative market has becomes extremely popular, this popularity even exceeds that of the stock exchange.Many problems in mathematical finance entail the computation of a particular integral.In many cases these integrals can be valued analytically and in some cases they can be valued using numerical integration, or computed using a partial differential equation.Option price estimation as the most interesting of the derivatives has many approaches and diverse qualities.
An option is an instrument whose value derives from that of another asset; hence it is called a "derivative".In other words an option on an underlying asset is an asymmetric contract that is negotiated today with the following conditions in the future.The holder has either the right, but not the obligation to buy, as it is the case with the European call option, or the possibility to sell, as in the case of the European put option, an asset for a certain price at a prescribed date in the future.The American type of option can be exercised at any time up to and including the date of expiry.The distinctive features of American option is its early exercise privilege.The pricing of American options has been the subject of extensive research in the last decades.There is no known closed form solution and many numerical and analytic approximations have been proposed.Black F. and Scholes M. [1] published their seminal work on option pricing in which they described a mathematical frame work for finding the fair price of a European option.They used a no-arbitrage argument to describe a partial differential equation which governs the evolution of the option price with respect to the maturity time and the price of the underlying asset.In the recent years the complexity of numerical computational in financial theory and practice has increased greatly, putting more demands on computation speed and efficiency.The subject of numerical methods in the area of options valuation and hedging is very broad.A wide range of different types of contracts are available and in many cases there are several candidate models for the stochastic evolution of the underlying state variables.
We present an overview of Mellin transform method for the valuation of options in the context of Black F. and Scholes M. [1].The Mellin transform is an integral transform named after the Finnish mathematician Hjalmar Mellin .Mellin transforms in option theory were introduced by Panini R. and Srivastav R.P. [14].They derived the expression for the free boundary and price of an American perpetual put as the limit of a finite-lived option.In this contribution Frontczak F. and Schöbel [7] extended the results obtained in Panini R. and Srivastav R.P. [14] and showed how the Mellin transform approach can be used to derive the valuation formula for perpetual American put options on dividend-paying stocks.Panini R. and Srivastav R.P. [15] derived integral equation representations for the price of European and American basket put options using Mellin transform techniques.Samuelson P.E.[16] derived a closed-form expression for the free boundary and price of a perpetual American put option using Mellin transform techniques, Frontczak F. and Schöbel [8] extended a framework based on the Mellin transforms and showed how to modify the approach to value American call options on dividend paying stocks, Zieneb A.E and Rokiah R.A [21] derived a closed form solution for a continuous arithmetic Asian option by means of partial differential equation.They also provided a new method for solving arithmetic Asian options using Mellin transforms in a stock price.Vasilieva O. [17] introduced a new method of pricing Multi-options using Mellin transforms and integral equations.
For mathematical backgrounds, Mellin transform method in derivatives pricing and some numerical methods for options valuation see [2], [3], [4], [5], [6], [8] [9], [10], [11], [12], [13], [16], [17], [18], [19], [20] just to mention a few.In this paper, we focus on Mellin transform method for the valuation of some vanilla power options with non-dividend yield.This paper is structured as follows: Section 2 presents the power options and a review of classic integral transforms.Section 3 presents some operational fundamental properties of Mellin transform method, derivation of Black-Scholes Formula for the valuation of European power put option using Mellin transform method and Mellin transform method for the valuation of European and American power put options.Section 4 concludes the paper.

Power Options
Power option is defined as a contingent claim on the product of powers of several assets.It can also be seen as a class of options in which the payoff at expiry is related to the n th power of the underlying price of the asset, where n ∈ {Z + − 0}.American power option is an option that can be exercised before or at the expiry date with non-linear payoff.European power option is an option that can be exercised only at the expiry date with non-linear payoff.
European power option comes in two forms namely European power call and put options.A European power call option is an option with non-linear payoff given by the difference between price of the underlying asset at maturity raised to a strictly positive power and the exercise price.A European 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 a European power option on the underlying price of the asset S t with exercise price K and time to expiry T , we have the payoff for the European power call option as The payoff for the European power put option is given as Remark 2.0.1: • For n > 0, European power option allows parties to negotiate the underlying price of the asset, exercise price, expiration date and other features.It also gives investors the opportunity to trade on a large scale with expanded or eliminated position limit and is of practical interest since Over-The-Counter (OTC) traded options exhibit such a payoff structure.
• The payoff curves for European power call and put options become concave and thus have a negative time value for n < 1 and n > 1 respectively, i.e. current price of the underlying asset (S n t ) < payoff (E n p (S, T )).
• For n = 1, the payoffs for European power call and put options in (1) and (2) become (3) and ( 4) respectively Equations ( 3) and (4) give the payoffs for regular European call and put options respectively and can also be expressed as • As n increases, the value of European power put option becomes very large.
Theorem 2.0.1:Let S n t denote the price of the underlying asset, σ the volatility, r the risk interest rate, n the power of the option and W t the Wiener process or Brownian motion.If the underlying price of the asset S n t follows a random process in then the explicit formula for the evolution of the underlying price of the asset is given by Differentiating (7) we have Recall from Ito's lemma using (5) for any derivative u(S n t , t) we have From (5), we can write for substituting ( 7), ( 8) and ( 10) into (9) and rearranging the terms, we have Therefore, Thus, ln S n t is a Brownian motion with drift parameter nr − 1 2 nσ 2 and variance parameter nσ 2 .To derive an explicit formula for the evolution of the stock price, we integrate (12) from 0 to T to obtain ln (14) can also be written as where Z ∼ N (0, 1).Therefore the stock dynamic follows a log-normal distribution.For n = 1, (15) becomes (16) shows that plain vanilla option follows a log-normal distribution.The partial differential equation for any derivative f written on S n t is given by Remark 2.0.2: • This result shows that the price of the underlying asset follows a random process.Also an explicit formula for the evolution of the underlying price of the asset can be obtained.
• For n = 1, ( 17) is known as the regular Black-Scholes partial differential equation which is given by The European power put option with non-dividend yield can be obtained by setting f = E n p (S n t , T ) = E n p in (17) as with boundary conditions lim Also we have for the European power call option with non-dividend yield: with boundary conditions lim The solutions to ( 18) and ( 22) can be obtained using the same procedures for solving Black-Scholes partial differential equation.

A Review of Classic Integral Transforms
We present here some integral transforms namely Fourier, Laplace and Mellin transforms and their inversion formulas.Under a suitable constant c ∈ ℜ and appropriate condition for the function f we have the following: • The Fourier transform and its inversion formula are given respectively by • The Laplace transform and its inversion formulas are given respectively by and • The Mellin transform and its inversion formula are given respectively by and The Fourier transform is popular in pure mathematics and finance.The Laplace transform is a widely used integral transform in mathematics with many applications in physics and engineering.The Mellin transform is used often in statistics, in analytic number theory and in the analysis of algorithms.All the three transforms and their inversion formulas are equivalent by elementary substitution.With the change of variable x = e s , it is observed that the Mellin transform is closely related to the Laplace transform and the Fourier transform.In particular, if F(f (x), ω) and L(f (x), ω) denote the two-sided Fourier and Laplace transform respectively, then we have However, there are numerous applications where it has been established that it is more convenient to operate directly with the Mellin transform rather than the Laplace-Fourier version.Next we discuss the Mellin transform for the valuation of some vanilla options.

The Mellin Transform Method for the Valuation of Some Vanilla Power Options with Non-Dividend Yield
This section presents some operational fundamental properties and Mellin transform method for the valuation of some vanilla power options namely European and American power options with non-dividend yield.

The Mellin Transform Method
The Mellin transform is a complex valued function defined on a vertical strip in the ω-plane whose boundaries are determined by the asymptotic behaviour of f (x) as x → 0 + and x → ∞.Mellin transform is denoted by M(f (x), ω) and defined as where f (x) is a locally Lebesgue integrable function.The largest strip (a, b) in which the integral converges is called fundamental strip.The conditions are when u > v, guarantee the existence of M(f (x), ω) in strip (−u, −v).Thus, the existence is granted for locally integrable functions, whose exponent in the order at 0 is strongly larger than the exponent of the order at infinity.Conversely, if f (x) is an integrable function with fundamental strip (a, b), then if c is such that a < c < b and f (c + it) is integrable, the equality in (32) holds almost everywhere on (0, ∞), where f (x) and f (ω) are known as Mellin transform pair.The inversion formula for the Mellin transform is given as Using (33) and (34) above we have the following conditions for European power put options: The Mellin transform of E n p (S n t , t) is given by where ω is a complex variable with 0 < Re(ω) < ∞.Using (35), the inversion formula of the Mellin transform is defined by

Some Operational Fundamental Properties of the Mellin Transform
If M(f (x), ω) = f (ω), then we have the following properties of the Mellin transform.
• Scaling Property e.g if a = 2, using (40) we have e.g if a = 4, using (41) we have where e.g if f (x) = x 2 , using (43) we have • Raising the independent variable to a Real Power • Inverse of Independent Variable • Multiplication by ln x • Multiplication by the Power of ln x • Convolution Property • Multiplicative Convolution where just to mention a few.

Some Integral Representations for European Power Call and Put Options Payoffs
Consider the power call payoff with strike price K > 0 given by ( 1) as where S n T is used as transform variable.An elementary calculations shows , (ℜω > 1) (54) and the inversion formula gives with c > 0 an arbitrary fixed number.Similarly the power put payoff given by ( 2) as has the same Mellin transform but the formula is valid in a different domain , (ℜω > 1) (56) and the inversion formula gives with d > 0 an arbitrary fixed number.For the case 0 < ℜω < 1, the Mellin transform and its inversion formula for power call payoff are obtained as with 0 < c < 1 an arbitrary fixed number.

Mellin Transform Method for the Valuation of European Power Put Option with Non-Dividend Yield
Consider the European power put option given by ( 18) as Taking the Mellin transform of ( 18) we have where and Substituting the values of (61), ( 62), ( 63), ( 64) and ( 65) into (60) we have Solving (68) integrating from 0 to t using separation of variables, we have Let Ẽn p (ω, 0) = c(ω), where c(ω) is a constant that depends on the boundary conditions, then Now let us consider the terminal condition of European power put option given by ( 24) which is of the form Then the Mellin transform of (24) is obtained as where (70) is independent of n.Therefore the constant c(ω) becomes Substituting ( 71) into (69) yields: Using the inversion formula of Mellin transform, the price of European power put option is obtained as Substituting (72) into (73), we have where φ(ω, t) is given by ( 70), (S n t , t) ∈ {(0, ∞) × [0, T ]}, c ∈ (0, ∞) and {ω ∈ C|0 < Re(ω) < ∞}

Derivation of Black-Scholes Formula for the Valuation of
European Power Put Option using Mellin Transform Method Theorem 3.3.1:Let S n t be the price of the underlying asset, K be the strike price, r be the risk interest rate and T be the time to maturity.Using the convolution property of the Mellin transform we have the Black-Scholes formula for the valuation of European power put option as Proof: Using the convolution property of Mellin Transform method, the price of European power put option is given by Therefore, where φ(v) is to be determined.Let us consider Setting ρ 1 = n 2 σ 2 2 (T − t) on the right side of (78), therefore (78) becomes Thus (77) becomes where G = −ρ 1 (ρ 2 2 + α 2 ).Now φ is the Mellin transform of Using the transformation given by The European power put option can therefore be expressed as ) denote the first and second terms of (86), then we write (86) as where and denote the first and second integrals in (86) respectively.We use the following transformations to evaluate the first and second integrals in (86) respectively: Solving further, the first term of ( 86) is obtained as where Similarly, the second term is evaluated as where We denote the distribution function for a normal variable by N (x) Substituting ( 92) and ( 94) into (87), we have (97) is called the analytical pricing formula for European power put option with non-dividend yield is given by with and as given above.Remark 3.3.1: • This result shows the Black-Scholes formula for the European power put option is obtained using the convolution property of the Mellin transform.
• For n = 1, (97) becomes analytical pricing formula "the so called Black-Scholes model" for European put option given by and

The Mellin Transform Method for the Valuation of American Power Put Option with Non-Dividend Yield
The key to determining the value of the American power option is to find the critical underlying price of the asset, which specifies the conditions under which the option should be exercised prior to maturity.The difference between European and American power options is that an American power option can be exercised by its holder at any time on or before the expiry date.This early exercise feature makes the pricing of American power option mathematically challenging and created great field of research throughout the last three decades.The valuation of American power option can be seen under several mathematical aspects, leading to different but equivalent mathematical formulations of the problem.In this paper we will make use of free boundary formulation for the valuation of American power put option with non-dividend yield.
We extend the domain of the Black-Scholes partial differential equation by setting the final time condition A n p (S n t , t) = A n p = (K −S n t ) + for S n t < Ŝn t , where the free boundary Ŝn t depends on the expiry time T .The non-homogeneous Black-Scholes equation for the price of American power put option with nondividend yield is given by on (0, ∞) × [0, T ).The boundary conditions are given by lim A n p (S n t , t) = 0 on [0, T ] (106) According to Wilmott et al [18], arbitrage arguments show that the option's price must satisfy smooth pasting conditions at the free boundary Ŝn t .We have the payoff for the American power put option as Differentiating (109) with respect to S n t = Ŝn , we have The Mellin transform of the American power put option in (103) is obtained as From the definition of Mellin transform, the hand side of (111) is given by Solving further, we have the particular solution of (112) as Theorem 3.4.1 [7]: If the underlying price of the asset S n t pays no dividends, the critical underlying price of the asset Ŝn t of an American power put option can be determined approximately by solving the implicit equation below where ) is given by (95).

Numerical Experiment
This section presents numerical experiment and discussion of results for the valuation of American power put option for n = 1 with non-dividend yield using Mellin transform method.For the European power put option the Mellin transform approach is simple and straight forward.

Example
We consider the valuation of American power put option for n = 1 which pays no dividend yield with the following parameters: The values of the Mellin transform columns in the Tables below are obtained from [15].The results generated are shown in the Tables 1, 2, 3, 4, 5 and 6.

Discussion of Results
We compare the performance of the Mellin transform method for the valuation of American put option which pays no dividend yield with other numerical methods namely binomial model [3] with 500 time steps and Crank Nicolson finite difference method [2] with S max = 100, 500 spatial and time grids against Black-Scholes model.Tables 1, 2, 3, 4, 5 and 6 show the variation of the price of the American put option via Mellin transform method, binomial model, Crank Nicolson finite difference method and Black-Scholes analytical pricing formula with the volatility σ = {0.2,0.4}, maturity time T = {0.0833,0.3333, 0.5833} and exercise price K = {35, 40, 45} and the fixed parameters S = 40, r = 0.0488 and c = 2.The results generated in the Tables below demonstrate that the Mellin transform method is mutually consistent, perform very well, accurate and agree with the values of Black-Scholes model (BSM) [1], Crank Nicolson finite difference method [2] and binomial model.It can be seen from the Tables below that the higher the volatility, the higher the price of American put option.Also the longer the time to expiration, the longer the time the asset price has to move until it is in the money.Hence the Mellin transform method [15] is accurate, easy to implement and is a better alternative technique to other approaches like the Crank Nicolson finite difference method [2] and binomial model [3] for the valuation of American put option on a non-dividend yield.

Conclusion
In this paper, we have considered the Mellin transform method for the valuation of some vanilla power options.We have established a formula for the valuation of European power put option using the Mellin transform method on a non-dividend yield consisting of single integral.We also generalize the results of Panini R. and Srivastav R.P. [14] to derive the valuation formula for the price and free boundary of American power put option with non-dividend yield.We also show how to derive the Black-Scholes valuation formula for European power put option.To provide a sufficient numerical analysis, we have compared the results generated by the Mellin transform method for the valuation of American power put option for n = 1 to two other numerical methods the Crank Nicolson finite difference method [2] and binomial model [3] in the context of Black-Scholes model [1].The numerical experiment has shown that the Mellin transform method is accurate, flexible, efficient and produces accurate prices for the critical stock price for a wide range of parameter combinations.The Mellin transform method is also a powerful technique for the valuation of more complex vanilla options and some path dependent options and can be extended to jump diffusion processes, stochastic volatility and other stochastic processes but modification is required on the Mellin transform method for the valuation of call options [8].

Appendix: Table of Results
The comparative results analysis are shown in the Tables below: for the call option and E p (S, T ) = max(K − S T , 0) = (K − S T ) + , for the put option respectively.