Selçuk J. Appl. Math. Selçuk Journal of Vol. 10. No. 1. pp. 135-140, 2009 Applied Mathematics
A Comparative Study of Black-Scholes Equation Refet Polat
Department of Mathematics, Faculty of Science and Letters Ya¸sar University Bornova, Izmir, Turkey
e-mail: refet.p olat@ yasar.edu.tr
Received: March 19, 2009
Abstract. In this paper, we analyzed Options and Black —Scholes Models for the valuing and pricing of commodities. In particular, we examined the numer-ical solution techniques of American Option Problems. For the comparison of the results pertaining to different methods, we used classical methods utilizing chronological order. Then we compared the results of these methods and tried to determine the most efficient method.
Key words: Options, American, Options, Black-Scholes Equation, Finite-difference Method.
2000 Mathematics Subject Classification: 26D15. 1.Introduction
The Black-Scholes Model was developed by Fisher Black and Myran Scholes in 1973 [2]. Firstly it was used to value European options that are not paying any financial profit. In the case of American options, a volatility event which will be an important research area for mathematicians occurs.
BLACK —SCHOLES EQUATION
The Black-Scholes equation in general form can be given as
(1) + 1 2 222 2 + − = 0 where
is current stock price, is option price, is volatility, is option expiration date,
The numerical methods [10], which are developed for solving American Options problem (1), are the following: Finite-difference method, direct discretization of PDE, implicit method, explicit method, Wilmott, Howison and Dewyne Model and Linear Complementarity Theta-weighted method. Instead of implementing the direct discretization of Scholes equation, Wilmott transformed Black-Scholes equations into heat equations and came up with a solution [10,12,13]. In this study, Wilmott’s results are compared with the numerical methods that are mentioned above.
2. Reduction of Black-Scholes Equation to Parabolic Equation In Wilmott’s paper [12], the reduction of Black-Scholes equation given in (1) is reduced to heat equation as follows:
In equation (1) Let us assign K as S’s appropriate basic value i.e.
(2) = Then = () ⇒ = 1 2 2 = − 1 2 (3) = = 1 2 2 = µ ¶ = µ ¶ + µ ¶ = µ 2 2 ¶ + 2 2 (4) = 1 2 2 2 − 1 2 If we insert (3) and (4) in place of (1)
+ ( − ) 1 + 1 2 22 1 2 µ 2 2 − ¶ − = 0 or + ( − ) + 1 2 22 2− 1 2 2 − = 0 or
(5) ⇒ + µ − −122 ¶ + 1 2 22 2 − = 0
We reach the equation (5). If we multiply this equation by 2 2 2 2 + 2 2( − ) − + 2 2− 2 2 = 0 ⇒ 22 + µ 2 2− 2 2 ¶ − + 2 2− 2 2 = 0 1= 2 2 and 2= 2 ( − ) 2 ;
Rearranging the above equation becomes
(6) 2 + (2− 1) − 1 = −2 2 Let us arrange the time variable again:
(7) = 1
2
2
( − )
For convenience, we take independent of x and t. Hence, t from equation (7); = −22 so; = = µ −2 2 ¶ ⇒ = − 2 2
is reached. If we insert this equation in (6)
(8) 2 2 + (2− 1) − 1 =
is reached. In order to solve equation (8), let us apply the separation of variables method.
( ) = ()( ). In this case we reach;
= ( )
0()
= ()
0( )
If we insert these equations in (8)
( )( 00() − 1 () + (2− 1) 0() = ()0( ) ⇒ 00() + (2− 1) ()0() − () = 0( ) ( ) = ⇒ 00() + (2− 1) 0() − 1 () = () and 0( ) ( ) = ⇒ () = + 1 ⇒ () = 1 In order to solve 00() + ( 2− 1) 0() − (1+ ) () = 0 () = Then the characteristic equation is
2+ (2− 1) − (1+ ) = 0 ⇒ 12= 1 − 2∓ q (2− 1)2+ 41 2 = 1 2(1 − 2) + r 1 4(1− 1) 2 + 1 () = − 1 2(2−1)+ 1 4(2−1)2+1 ( ) = ()( ) (9) ( ) = − 1 2(2−1)+ 1 4(2−1)2+1+ ( )
The solution (9) is the semi-analytical solution for the Black-Scholes Partial Dif-ferential Equation. This solution will be compared with the numerical solutions mentioned in introduction.
3. Numerical Experiment and The Result
The analysis of our results shows that among the different methods [10] used– namely, Explicit, Implicit, Brennan & Schwartz First Model [3], Brennan & Schwarz Second Model [4], Courtadon [5], Wilmott Howison Dewyne Theta weighted and Linear Complementarity Theta weighted method [12]–Wilmott’s method is the most efficient one [10]. The table given below verifies that all obtained results are diverse due to the applied methods. These empiric results show that the stability of the equation differs from method to method; but, in our opinion, Wilmott’s method may be taken as the most appropriate one. Therefore the reason that we made our comparison with respect to this solution is that the method is a semi-analytic solution. In the Table 1, it is observed that the best value is obtained with linear complementarily methods, but the complexity of the method is very high which makes it less effective. In the other numerical methods, complexity looks pretty appropriate. Nevertheless the values differ from method to method, which shows the unstability of the equation.
Cpu: Pentium 4, 2.2
References
1. C. Albenese, S. Jaimungal and D. Rubinsov, “The Model of lines for Option Pricing with Jump”, 2001.
2. F. Black and M. Scholes, “The pricing of options and corporate liabilities, Journal of Political Economy, 81 637-59, 1973.
3. M. Brennan and E. Schwartz, “The valuation of American put options, Journal of Finance, 1977
4. M. Brennan and E. Schwartz, “Finite-difference methods and jump processes arising in the Pricing of contingent claims: A synthesis, Journal of Financial and Quantitative Analysis, 1978
5. G. Courtadon, “A more accurate finite difference approximation for the valuation of Options, Journal of Financial and Quantitative Analysis, 1982
6. H. P Mc Kean, Appendix, “A Free Boundary Problem for the Heat Equation Arising from a Problem in Mathematical Economics.” Industr. Manage. Rev., 6: 32-39, 1965. 7. Hull, John, “Options, Future and other Derivative Securities”, Willeys, 1997. 8. J. J McCutcheon-W. F. Scott, “An Introduction to the Mathematics of Finance”, 1989.
9. P. Van Moerbeke, “On Optimal Stopping and Free Boundary Problem”, Arch. Rational Mech. And. 60: 101-148, 1976.
10. S. Sukha, “Finite-Difference Methods for Pricing the American Put Options”, 2001.
11. W. Allegretto, G. Barone-Adesi and R. J. Elliot, “Numerical Evaluation of the Critical Price and American Options”, Europen J. Finance, 1: 69-78, 1995.
12. Wilmott, Paul-Dewynne, Jeff-Howison, Sam, “Option Pricing Mathematical Modals and Computation”, Oxford Financial Press, 1998.
13. Wilmott, Paul-Dewynne, Jeff-Howison, Sam, “The Mathematics of Financial Derivatives” A student Introduction, 1997.