Numerical Evaluation of European Option on a Non Dividend Paying Stock

Similar documents
The Financial System. Instructor: Prof. Menzie Chinn UW Madison

Fugit (options) The terminology of fugit refers to the risk neutral expected time to exercise an

Section 6 Short Sales, Yield Curves, Duration, Immunization, Etc.

American basket and spread options. with a simple binomial tree

SOCIETY OF ACTUARIES FINANCIAL MATHEMATICS. EXAM FM SAMPLE SOLUTIONS Interest Theory

Dynamic Relationship and Volatility Spillover Between the Stock Market and the Foreign Exchange market in Pakistan: Evidence from VAR-EGARCH Modelling

Pricing and Valuation of Forward and Futures

Interest Rate Derivatives: More Advanced Models. Chapter 24. The Two-Factor Hull-White Model (Equation 24.1, page 571) Analytic Results

STOCHASTIC LOCAL VOLATILITY

FITTING EXPONENTIAL MODELS TO DATA Supplement to Unit 9C MATH Q(t) = Q 0 (1 + r) t. Q(t) = Q 0 a t,

A valuation model of credit-rating linked coupon bond based on a structural model

Online appendices from Counterparty Risk and Credit Value Adjustment a continuing challenge for global financial markets by Jon Gregory

Normal Random Variable and its discriminant functions

MORNING SESSION. Date: Wednesday, May 4, 2016 Time: 8:30 a.m. 11:45 a.m. INSTRUCTIONS TO CANDIDATES

Deriving Reservoir Operating Rules via Fuzzy Regression and ANFIS

Chain-linking and seasonal adjustment of the quarterly national accounts

Lecture 11 Partial Differential Equations. Partial Differential Equations (PDEs). What is a PDE? Examples of Important PDEs. Classification of PDEs.

Pricing Model of Credit Default Swap Based on Jump-Diffusion Process and Volatility with Markov Regime Shift

Lab 10 OLS Regressions II

Improving Forecasting Accuracy in the Case of Intermittent Demand Forecasting

The Proposed Mathematical Models for Decision- Making and Forecasting on Euro-Yen in Foreign Exchange Market

Real-World Pricing for a Modified Constant Elasticity of Variance Model

Recall from last time. The Plan for Today. INTEREST RATES JUNE 22 nd, J u n e 2 2, Different Types of Credit Instruments

Agricultural and Rural Finance Markets in Transition

OPTIMAL EXERCISE POLICIES AND SIMULATION-BASED VALUATION FOR AMERICAN-ASIAN OPTIONS

A Solution to the Time-Scale Fractional Puzzle in the Implied Volatility

Gaining From Your Own Default

Baoding, Hebei, China. *Corresponding author

Floating rate securities

Differences in the Price-Earning-Return Relationship between Internet and Traditional Firms

In-Arrears Interest Rate Derivatives under the 3/2 Model

Michał Kolupa, Zbigniew Śleszyński SOME REMARKS ON COINCIDENCE OF AN ECONOMETRIC MODEL

Stochastic Local Volatility

A New N-factor Affine Term Structure Model of Futures Price for CO 2 Emissions Allowances: Empirical Evidence from the EU ETS

Comparing Sharpe and Tint Surplus Optimization to the Capital Budgeting Approach with Multiple Investments in the Froot and Stein Framework.

Accuracy of the intelligent dynamic models of relational fuzzy cognitive maps

Valuing Surrender Options in Korean Interest Indexed Annuities

Noise and Expected Return in Chinese A-share Stock Market. By Chong QIAN Chien-Ting LIN

Correlation of default

Centre for Computational Finance and Economic Agents WP Working Paper Series. Amadeo Alentorn Sheri Markose

Quarterly Accounting Earnings Forecasting: A Grey Group Model Approach

Economics of taxation

Network Security Risk Assessment Based on Node Correlation

Financial Innovation and Asset Price Volatility. Online Technical Appendix

Terms and conditions for the MXN Peso / US Dollar Futures Contract (Physically Delivered)

Pricing under the Real-World Probability Measure for Jump-Diffusion Term Structure Models

Estimating intrinsic currency values

The UAE UNiversity, The American University of Kurdistan

Correlation Smile, Volatility Skew and Systematic Risk Sensitivity of Tranches

STOCK PRICES TEHNICAL ANALYSIS

Albania. A: Identification. B: CPI Coverage. Title of the CPI: Consumer Price Index. Organisation responsible: Institute of Statistics

ESSAYS ON MONETARY POLICY AND INTERNATIONAL TRADE. A Dissertation HUI-CHU CHIANG

Fairing of Polygon Meshes Via Bayesian Discriminant Analysis

DYNAMIC ECONOMETRIC MODELS Vol. 8 Nicolaus Copernicus University Toruń 2008

Convexity Adjustments in Inflation linked Derivatives using a multi-factor version of the Jarrow and Yildirim (2003) Model

Bank of Japan. Research and Statistics Department. March, Outline of the Corporate Goods Price Index (CGPI, 2010 base)

Convertible Bonds and Stock Liquidity. Author. Published. Journal Title DOI. Copyright Statement. Downloaded from. Griffith Research Online

Tax Dispute Resolution and Taxpayer Screening

An Inclusion-Exclusion Algorithm for Network Reliability with Minimal Cutsets

Time-domain Analysis of Linear and Nonlinear Circuits

Jean-Paul Murara, Västeras, 26-April Mälardalen University, Sweden. Pricing EO under 2-dim. B S PDE by. using the Crank-Nicolson Method

Permanent Income and Consumption

Prediction of Oil Demand Based on Time Series Decomposition Method Nan MA * and Yong LIU

Convexity adjustments in inflation-linked derivatives

Credit Value Adjusted Real Options Based Valuation of Multiple-Exercise Government Guarantees for Infrastructure Projects

A Novel Application of the Copula Function to Correlation Analysis of Hushen300 Stock Index Futures and HS300 Stock Index

Asset Pricing and Option Trading for Harvesting Rights to Renewable Resources

Adjusted-Productivity Growth for Resource Rents: Kuwait Oil Industry

Boğaziçi University Department of Economics Money, Banking and Financial Institutions L.Yıldıran

OPTIMIZED CALIBRATION OF CURRENCY MARKET STRATEGIES Mustafa Onur Çağlayan 1, János D. Pintér 2

Assessment of The relation between systematic risk and debt to cash flow ratio

The impact of intellectual capital on returns and stock prices of listed companies in Tehran Stock Exchange

Impact of Stock Markets on Economic Growth: A Cross Country Analysis

INSTITUTE OF ACTUARIES OF INDIA

THE APPLICATION OF REGRESSION ANALYSIS IN TESTING UNCOVERED INTEREST RATE PARITY

Improving Earnings per Share: An Illusory Motive in Stock Repurchases

A Framework for Large Scale Use of Scanner Data in the Dutch CPI

Methodology of the CBOE S&P 500 PutWrite Index (PUT SM ) (with supplemental information regarding the CBOE S&P 500 PutWrite T-W Index (PWT SM ))

Cointegration between Fama-French Factors

Using Fuzzy-Delphi Technique to Determine the Concession Period in BOT Projects

Management of financial and consumer satisfaction risks in supply chain design

A Novel Particle Swarm Optimization Approach for Grid Job Scheduling

A Hybrid Method for Forecasting with an Introduction of a Day of the Week Index to the Daily Shipping Data of Sanitary Materials

Online Technical Appendix: Estimation Details. Following Netzer, Lattin and Srinivasan (2005), the model parameters to be estimated

Lecture 12: Asymptotics

Optimal Combination of Trading Rules Using Neural Networks

Volatility Forecasting Models and Market Co-Integration: A Study on South-East Asian Markets

Conditional Skewness of Aggregate Market Returns

Global regional sources of risk in equity markets: evidence from factor models with time-varying conditional skewness

IJEM International Journal of Economics and Management

Return Calculation Methodology

THE TYRANNY OF THE IDENTITY: GROWTH ACCOUNTING REVISITED

The Macrotheme Review A multidisciplinary journal of global macro trends

A Hybrid Method to Improve Forecasting Accuracy Utilizing Genetic Algorithm An Application to the Data of Operating equipment and supplies

UC San Diego Recent Work

iii) pay F P 0,T = S 0 e δt when stock has dividend yield δ.

Chapter 9 - Level 2 - Course FM Solutions

Mixtures of Normal Distributions: Application to Bursa Malaysia Stock Market Indices

The Empirical Research of Price Fluctuation Rules and Influence Factors with Fresh Produce Sequential Auction Limei Cui

The Comparison among ARMA-GARCH, -EGARCH, -GJR, and -PGARCH models on Thailand Volatility Index

Transcription:

Inernaonal Journal of Compuaonal cence and Mahemacs. IN 0974-389 olume Number 3 (00) pp. 6--66 Inernaonal Research Publcaon House hp://www.rphouse.com Numercal Evaluaon of European Opon on a Non Dvdend Payng ock bha Mshra and K.R. Pardasan Deparmen of Mahemacs M.A.N.I.T Bhopal (M.P.) Inda E-mal: vbha_67@redffmal.com kamalrp@homal.com Absrac In hs paper an aemp has been made o numercally evaluae he value of European opon employng Black-scholes formula. The mporan parameers lke me volaly srke prce neres rae and spo prce have been ncorporaed n he model. Approprae boundary and nal condons have been framed usng fnancal condons of he sock marke. The numercal mehods especally he mplc and explc schemes have been employed o solve Black-choles PDE. The resuls obaned are used o sudy he relaonshps among varous parameers. Keywords: European opon non-dvdend srke prce volaly. Inroducon Black and choles [9] proposed an explc formula for evaluang European call opons whou dvdends whch s sll exensvely used wh underlyng. The numercal soluon of hs equaon has been of paramoun neres due o he governng paral dfferenal equaon whch s very dffcul o generae sable and accurae soluons [94]. A European call opon gves s holder he rgh (bu no he oblgaon) o purchase from he wrer a prescrbed asse for a prescrbed prce a a prescrbed me n he fuure. The prescrbed purchase prce s known as he exercse prce or srke prce (E) and he prescrbed me n he fuure s known as he expry dae (T). The drec oppose of a European call opon s a European pu opon [3 4 5 ]. A European pu opon gves s holder he rgh (bu no he oblgaon) o sell o he wrer a prescrbed asse for a prescrbed prce a a prescrbed me n he fuure. Therefore he value of he opon a expry known as he pay-off funcon s ( T) (-E). Numercal mehod of long sandng for general parabolc free boundary problem s he me dscree. If dfferenal equaon s used o model he value of he opon

6 bha Mshra and K.R. Pardasan hen he resulng problem no longer s formulaed for a consan coeffcen scalar (hea) equaon bu for a one dmensonal non lnear parabolc equaon[ 0]. Fne dfference heory has a long hsory and has been appled for more han 00 years o approxmae he soluons of paral dfferenal equaons n he physcal scences and engneerng. Wha s he relaonshp beween FDM and fnancal mahemacs? To answer hs queson we noe ha he behavor of a sock (or some oher underlyng) can be descrbed by a sochasc dfferenal equaon [6 7]. Then a conngen clam ha depends on he underlyng s modeled by a paral dfferenal equaon n combnaon wh some nal and boundary condons. olvng hs problem means ha we have found he value for he conngen clam [3]. Problem In hs paper we consder how mgh be used o value an European opon on a nondvdend payng sock. The Black-choles paral dfferenal equaon for European opon on non-dvdend payng sock s. r σ r () uppose ha he lfe of he opon s T. we dvde hs no N equally spaced T nervals of lengh. A oal of N mes are consdered..e. N 0 3... T.mlarly max s a sock prce of he opon. We max defne and consder a oal of M equally spaced sock prce of he M opons. 0 3... max. The ( ) pons on he grd s he pon ha correspondng o and sock prce. We wll use he varable o denoe he value of he opon a he ( ) pon. The followng wo mehods are employed for solvng he fundamenal PDE numercally. The mplc Fne Dfference mehod (The backward Euler mehod) The explc Fne Dfference mehod (The Euler mehod) Implc Fne Dfference Mehod: - For he neror pon ( ) on he grd can be approxmaed as () OR (3) Equaon () s known as he forward dfference approxmaon and equaon (3) s known as he backward dfference approxmaon. We use more symmercal

Numercal Evaluaon of European Opon 63 approxmaon by averagng he wo. (4) For we wll use a forward dfference approxmaon so ha he value a me s relaed o he value a me (). (5) The backward dfference approxmaon for a he ( ) pon s gven by equaon () he backward dfference a he () pon s Hence a fne dfference approxmaon for a he ( ) pon s (6) ubsung equaons (4) (5) and (6) no he dfferenal equaon () and gves r r σ For 3.. M- and 03 N-rearrangng erms we oban c b a (7) Where r a σ r b σ r c σ A program has been developed n MATLAB 7.5 for hs mplc fne dfference scheme o evaluae he fnancal opons.

64 bha Mshra and K.R. Pardasan Explc mehod Here we use a backward dfference approxmaon of he me dervave. Ths s he dea behnd he explc mehod whch approxmaes If hs approxmaon s appled n conuncon wh (4) and (6) n he PDE () hen r σ r For 3.. M- and 03 N- and rearrangng erm we oban Here a b c! (8) r σ r a r b ( σ ) r σ r c The program has been developed for explc fne dfference mehod n Malab7.5 o evaluae he fnancal opons. Resuls and Dscusson The value of varous parameers lke me volaly srke prce sock prce and neres rae whch have been used for numercal compuaons are gven n able (). The numercal resuls have been obaned for European call and pu opons by he wo numercal mehods namely Implc and Explc fne dfference mehod. The fgure () & () represens he resuls obaned by mplc mehod whle fgure 3 & 4 represens he resuls obaned by explc mehod. In Fg and 3 we observe ha when me s. yrs (73days) opon prce sar ncreasng when sock prce s 50 and opon prce goes an ncreasng when me s yr (365 days). we see ha hs relaon among opon prce sock prce and me s lnear. Furher he relaonshp beween opon prce and sock prce for pu opon n fgure () and (4) are us oppose of he same relaonshp for call opon n fgure () and (3). There s no sgnfcan dfference n he resuls obaned by he mplc and explc mehod. uch models can be developed o sudy relaonshp among he fnancal parameers for predcon of sock markes.

Numercal Evaluaon of European Opon 65 Table : alues of parameers for European opons [6789]. Parameers ymbol Numercal value ock prce 00 rke /Exercse Prce E 50 Ineres rae R 0% Opon Maury T monhs yr olaly σ.06 Fgure : Graph among Opon Prce sock Prce and me for call Opon by Implc mehod. Fgure : Graph among Opon Prce sock Prce and me for pu opon by mplc mehod. Fgure 3: Graph among Opon Prce sock Prce and me for call opon by Explc mehod. Fgure 4: Graph among Opon Prce sock Prce and me for pu opon by explc mehod. References [] Jula Ankudnova & Mahas Ehrhard(008) On he numercal soluon of nonlnear Black choles equaons Compuers and Mahemacs wh Applcaons vol. 56 799 8.

66 bha Mshra and K.R. Pardasan [] Hull John C. Opons (006) Fuures and Oher Dervaves 6h ed. Upper addle Rver NJ: Prence-Hall chaps. 9. [3] Angelos Dassos & Jayalaxshm Nagaradasarma 004 Prcng Of Asan Opons on Ineres Raes n he CIR Model Proceedng (437) Fnancal Engneerng and Applcaons. [4] Andrew M. Chsholm Dervaves Demysfed (004). A ep-by-ep Gude o Forwards Fuures waps and Opons John Wley & ons Ld. [5] Desmond J. Hgham 004 An Inroducon To Fnancal Opon aluaon Mahemacs ochasc And Compuaon Cambrdge unversy Press. [6] Jaksa Cvanc and Fernando Zapaero (004) Inroducon To The Economcs And Mahemacs Of Fnancal Marke The MIT PRE Cambrdge. [7] Karhan Akcoglu k Mng-Yang Kao & huba Raghavan feb (004) Fas Prcng of European Asan Opons wh Provable Accuracy: ngle-sock and Baske Opons. [8] Leser Ingber & Jennffer K. Wlson (004) ascal Mechancs of Fnancal Markes: Exponenal Modfcaons o Black-choles. [9] N. ukhomln (004) ymmery and new soluons of he Black choles equaon Bolen de la Asocacon Maemaca enezolana () 75 89. [0] Cho. and Marcozz M.D. (00) A Numercal Approach o Amercan Curreny Opon aluaon Journal of Dervaves vol. 9. [] M. Yor (00) Exponenal funconals of Brownan moon and relaed processes prnger Fnance. prnger-erlag Berln. [] G. Alobad (000) Amercan opons and her sraeges Ph.D. hess Unversy of Wesern Onaro London Canada. [3] Rober Kolb John Wley & sons (995) Inc. New York Chcheser undersandng opons [4]. Heson (993) A closed-form soluon for opons wh sochasc volaly wh applcaons o bond and currency opons The Revew of Fnancal udes 6:37 343. [5] J.C. Cox.A. Ross M. Rubnsen 7 (979) Opon prcng:a smplfed approach Journal of Fnancal Economcs 9 63. [6] hp://www.quckmba.com/fnance/black-scholes [7] hp://www.sreeauhory.com/erms/b/black-scholes-model.asp [8] hp://www.oponradngpeda.com/free_black_scholes_model.hm [9] hp://www.hoadley.ne/opons Hoadley Tradng & Invesmen Tool