CS476/676 Mar 6, Today s Topics. American Option: early exercise curve. PDE overview. Discretizations. Finite difference approximations

Similar documents
Numerical Methods in Option Pricing (Part III)

FINITE DIFFERENCE METHODS

American Equity Option Valuation Practical Guide

Computational Finance Finite Difference Methods

Lecture 4 - Finite differences methods for PDEs

arxiv: v1 [q-fin.cp] 1 Nov 2016

Option Pricing Models for European Options

CS 774 Project: Fall 2009 Version: November 27, 2009

FE610 Stochastic Calculus for Financial Engineers. Stevens Institute of Technology

Definition Pricing Risk management Second generation barrier options. Barrier Options. Arfima Financial Solutions

Binomial model: numerical algorithm

PDE Methods for the Maximum Drawdown

4. Black-Scholes Models and PDEs. Math6911 S08, HM Zhu

Near-Expiry Asymptotics of the Implied Volatility in Local and Stochastic Volatility Models

Queens College, CUNY, Department of Computer Science Computational Finance CSCI 365 / 765 Fall 2017 Instructor: Dr. Sateesh Mane.

Lecture 11: Ito Calculus. Tuesday, October 23, 12

Black-Scholes model: Derivation and solution

Solving the Black-Scholes Equation

Price sensitivity to the exponent in the CEV model

A Note about the Black-Scholes Option Pricing Model under Time-Varying Conditions Yi-rong YING and Meng-meng BAI

Chapter 5 Finite Difference Methods. Math6911 W07, HM Zhu

Project 1: Double Pendulum

32.4. Parabolic PDEs. Introduction. Prerequisites. Learning Outcomes

Exercises for Mathematical Models of Financial Derivatives

MAFS Computational Methods for Pricing Structured Products

NUMERICAL METHODS OF PARTIAL INTEGRO-DIFFERENTIAL EQUATIONS FOR OPTION PRICE

The Black-Scholes Equation

Short-time-to-expiry expansion for a digital European put option under the CEV model. November 1, 2017

A distributed Laplace transform algorithm for European options

Computational Finance

A Worst-Case Approach to Option Pricing in Crash-Threatened Markets

WKB Method for Swaption Smile

Volatility Trading Strategies: Dynamic Hedging via A Simulation

Hedging Credit Derivatives in Intensity Based Models

HIGHER ORDER BINARY OPTIONS AND MULTIPLE-EXPIRY EXOTICS

EFFICIENT MONTE CARLO ALGORITHM FOR PRICING BARRIER OPTIONS

Lecture 4. Finite difference and finite element methods

An Asymptotic Expansion Formula for Up-and-Out Barrier Option Price under Stochastic Volatility Model

2.1 Mathematical Basis: Risk-Neutral Pricing

arxiv: v1 [q-fin.pm] 13 Mar 2014

An Asymptotic Expansion Formula for Up-and-Out Barrier Option Price under Stochastic Volatility Model

Stochastic Differential Equations in Finance and Monte Carlo Simulations

Introduction to Financial Mathematics

Research Article Exponential Time Integration and Second-Order Difference Scheme for a Generalized Black-Scholes Equation

6. Numerical methods for option pricing

ADAPTIVE PARTIAL DIFFERENTIAL EQUATION METHODS FOR OPTION PRICING

MSC FINANCIAL ENGINEERING PRICING I, AUTUMN LECTURE 6: EXTENSIONS OF BLACK AND SCHOLES RAYMOND BRUMMELHUIS DEPARTMENT EMS BIRKBECK

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

Fractional Black - Scholes Equation

Boundary conditions for options

ANALYSIS OF THE BINOMIAL METHOD

Computational Finance Binomial Trees Analysis

Aspects of Financial Mathematics:

CONVERGENCE OF NUMERICAL METHODS FOR VALUING PATH-DEPENDENT OPTIONS USING INTERPOLATION

Numerical Evaluation of American Options Written on Two Underlying Assets using the Fourier Transform Approach

Finite Element Method

Write legibly. Unreadable answers are worthless.

King s College London

Forward Monte-Carlo Scheme for PDEs: Multi-Type Marked Branching Diffusions

MATH4143: Scientific Computations for Finance Applications Final exam Time: 9:00 am - 12:00 noon, April 18, Student Name (print):

PAijpam.eu ANALYTIC SOLUTION OF A NONLINEAR BLACK-SCHOLES EQUATION

From Discrete Time to Continuous Time Modeling

Linear and Nonlinear Models for the Underlying Asset V(P) and the NPV Equation

Computational Finance

Option Valuation with Sinusoidal Heteroskedasticity

Optimal exercise price of American options near expiry

Risk Neutral Pricing Black-Scholes Formula Lecture 19. Dr. Vasily Strela (Morgan Stanley and MIT)

Numerical Methods for Optimal Stochastic Control in Finance

Weak Reflection Principle and Static Hedging of Barrier Options

1 The Hull-White Interest Rate Model

Lecture 15: Exotic Options: Barriers

MASM006 UNIVERSITY OF EXETER SCHOOL OF ENGINEERING, COMPUTER SCIENCE AND MATHEMATICS MATHEMATICAL SCIENCES FINANCIAL MATHEMATICS.

Advanced Topics in Derivative Pricing Models. Topic 4 - Variance products and volatility derivatives

SPDE and portfolio choice (joint work with M. Musiela) Princeton University. Thaleia Zariphopoulou The University of Texas at Austin

MATH3075/3975 FINANCIAL MATHEMATICS TUTORIAL PROBLEMS

1 Geometric Brownian motion

The Evaluation of American Compound Option Prices under Stochastic Volatility. Carl Chiarella and Boda Kang

Math 623 (IOE 623), Winter 2008: Final exam

Numerical Methods for Stochastic Differential Equations with Applications to Finance

Economathematics. Problem Sheet 1. Zbigniew Palmowski. Ws 2 dw s = 1 t

Options. An Undergraduate Introduction to Financial Mathematics. J. Robert Buchanan. J. Robert Buchanan Options

Introduction Credit risk

The Black-Scholes Model

Continuous Time Mean Variance Asset Allocation: A Time-consistent Strategy

The Volatility Smile Dynamics Implied by Smile-Consistent Option Pricing Models and Empirical Data

Multiple Optimal Stopping Problems and Lookback Options

A Study on Numerical Solution of Black-Scholes Model

BACHELIER FINANCE SOCIETY. 4 th World Congress Tokyo, Investments and forward utilities. Thaleia Zariphopoulou The University of Texas at Austin

A Simple Numerical Approach for Solving American Option Problems

Sample Path Large Deviations and Optimal Importance Sampling for Stochastic Volatility Models

Homework Set 6 Solutions

The Black-Scholes Equation using Heat Equation

A Consistent Pricing Model for Index Options and Volatility Derivatives

Monte Carlo Methods. Prof. Mike Giles. Oxford University Mathematical Institute. Lecture 1 p. 1.

Fast Convergence of Regress-later Series Estimators

AMERICAN OPTION PRICING UNDER STOCHASTIC VOLATILITY: A SIMULATION-BASED APPROACH

STOCHASTIC INTEGRALS

1 Explicit Euler Scheme (or Euler Forward Scheme )

Deriving and Solving the Black-Scholes Equation

( ) since this is the benefit of buying the asset at the strike price rather

Transcription:

CS476/676 Mar 6, 2019 1 Today s Topics American Option: early exercise curve PDE overview Discretizations Finite difference approximations

CS476/676 Mar 6, 2019 2 American Option American Option: PDE Complementarity problem V t + 1 2 σ2 S 2 2 V S 2 +(r q)s V S V(S,t) payoff(s) 0 rv 0, (V(S,t) payoff(s))( V t + 1 2 σ2 S 2 2 V S 2 +(r q)s V S rv) = 0 0 < S < +, 0 t < T and V(S,T) = payoff(s).

CS476/676 Mar 6, 2019 3 The early exercise boundary {S t} is implicitly determined. Assume that V(S, t) solves the complementarity problem. For a put, S t = max{s : V(S,t) = payoff(s)}. For a call (with dividend yield), S t = min{s : V(S,t) = payoff(s)}.

CS476/676 Mar 6, 2019 4 Pricing European Options: PDE Approach Assume ds t = µs t dt+σ(s t,t)s t dz t (assume no dividend for simplicity) Any European option value V(S, t) satisfies V t + 1 2 σ2 S 2 2 V S +(r q)s V rv = 0, 2 S 0 < S < +, 0 t < T A PDE has an infinite number of solutions. It needs proper combinations of initial/final conditions to uniquely determine a solution. Numerically, boundary conditions are also required.

CS476/676 Mar 6, 2019 5 Assume that σ,r,q are constants. Consider transformations: S = Ke x, t = T 2τ σ 2, ˆδ = 2r 2(r q) σ2, δ = σ 2, ( V(S,t) = V Ke x,t 2τ ) def = v(x,τ) σ 2 v(x,τ) = Ke 1 2 (δ 1)x (1 4 (δ 1)2 +ˆδ)τ u(x,τ)

CS476/676 Mar 6, 2019 6 It can be shown that the BS PDE is equivalent to which is called a heat equation. u τ = u xx It models the diffusion of heat on an (uniform) infinite bar from an initial temperature distribution. The solution to a heat equation is infinitely smooth, even if initially the distribution is discontinuous. European option value function V(S, t) is infinitely smooth

CS476/676 Mar 6, 2019 7 A partial differential equation with two independent variables can be written as Au xx +Bu xy +Cu yy +F(u x,u y,u,x,y) = 0 Note. Subsequently subscripts denote partial derivatives with respect to the independent variables. Parabolic: B 2 4AC = 0 (BS Eqns) Hyperbolic: B 2 4AC > 0 Elliptic: B 2 4AC < 0 Different types of PDE have different properties.

CS476/676 Mar 6, 2019 8 The heat equations (thus BS equation) has an explicit solution. But for a LVF model and American options, numerical methods based on discretizations are required to compute their values.

CS476/676 Mar 6, 2019 9 Numerical methods for option pricing: Lattice methods Actually simple explicit finite difference methods Not very efficient, very complex for certain cases Monte-Carlo Method of choice for > 3 underlying stochastic factors Generally difficult to handle American early exercise (ideas: Longstaff-Schwartz) 3 stochastic factors, not very efficient

CS476/676 Mar 6, 2019 10 Numerical PDE Method of choice for < 3 stochastic factors Complex features easy to handle (i.e. complex American type contracts) Get value, hedging parameters for whole range of asset prices, not just a single point Useful for simulating hedging strategies Can handle more complex nonlinear PDE problems arising in finance Can build up general library so that many different contracts can be quickly priced

CS476/676 Mar 6, 2019 11 In folloing discussion, we consider the transformation τ = T t (time to expiry) and BS PDE has the form: V τ = σ2 2 S2 V SS +rsv S rv V = value of option T = expiry time of option τ = T t σ = volatility S = asset price r = risk free interest rate Note: V(S,τ) now represents the option value when time to expiry is τ and underlying price is S.

CS476/676 Mar 6, 2019 12 Numerical Methods for PDE: an overview Initial condition: V(S,0) = payoff(s). Use a finite computational domain [0,S max ] [0,T] to approximate [0,+ ] [0,T] Apply some boundary conditions at S = 0 and S = S max ; the effect of the boundary conditions is small if S 0 is sufficiently away from the boundary

CS476/676 Mar 6, 2019 13 Discretize the continuous region by a finite set of grids {(S i,τ n ) : n = 0, N,i = 0,,M} Discretize PDE at each grid point to obtain an algebraic equation relates values at n with values at n+1 As N,M +, computed values approximate the solution to the PDE

CS476/676 Mar 6, 2019 14 Discretization PDE Methods Finite difference method Replace derivatives by finite differences, Taylor series error analysis Standard FD analysis requires smooth payoffs Alternative FE analysis allows discontinuous payoffs FD,FE give the same discrete equations