MAFS Computational Methods for Pricing Structured Products
|
|
- Clarissa Reed
- 5 years ago
- Views:
Transcription
1 MAFS550 - Computational Methods for Pricing Structured Products Solution to Homework Two Course instructor: Prof YK Kwok 1 Expand f(x 0 ) and f(x 0 x) at x 0 into Taylor series, where f(x 0 ) = f(x 0 ) f (x 0 ) + 4 x f (x 0 ) + O( x 3 ) f(x 0 x) = f(x 0 ) xf (x 0 ) + x f (x 0 ) + O( x 3 ) We determine α, α 1 and α 0 such that α f(x 0 ) + α 1 f(x 0 x) + α 0 f(x 0 ) = f (x 0 ) + O( x ) Collecting like terms, we obtain (α + α 1 + α 0 )f(x 0 ) + ( α α 1 ) xf (x 0 ) + (4α + α 1 ) x f (x 0 ) = f (x 0 ) + O( x ) The corresponding linear system of equations for α, α 1 and α 0 is α α 1 = 1/ x α 0 0 The solution of the system gives α = 1, α 1 = x and α 0 = 3 The corresponding finite difference formula is the one-sided backward difference formula for the first order derivative, namely f (x 0 ) f(x 0 ) 4f(x 0 x) + 3f(x 0 ) By changing x to x, we deduce that the one-sided forward difference formula for the first order derivative is given by f (x 0 ) f(x 0 + ) + 4f(x 0 + x) 3f(x 0 ) Local truncation error of the Crank-Nicolson scheme V ( x, (n + 1)) V ( x, n) = σ V (( + 1) x, (n + 1)) V ( x, (n + 1)) + V (( 1) x, (n + 1)) 4 x V (( + 1) x, n) V ( x, n) + V (( 1) x, n)) + x 1 ) (r σ V (( + 1) x, (n + 1)) V (( 1) x, (n + 1)) V (( + 1) x, n) V (( 1) x, n) + + r V ( x, (n + 1)) + V ( x, n) 1
2 Expanding ( ( each term using the Taylor series expansion at the intermediate time level x, n + 1 ) ), we obtain = V ( x, (n + 1)τ) V ( x, n) { V + V τ + 1 ( ) V + 1 τ 6 V V τ + 1 ( V τ = V τ V 4 τ 3 () + O( 4 ) 3 V τ 3 ) 1 3 V 6 τ 3 ( ) 3 + ( ) } 3 / + Here, we adopt the convention that any derivative ( with no ( argument specified would implicitly implies that the derivative is evaluated at x, n + 1 ) ) Note that and V (( + 1) x, n) V ( x, n) + V (( 1) x, n)/( x) = V x 4 V ( x, n) + x 1 x ( x, n) + 4 O( x4 ) = V x 3 V x τ + 1 ( ) 4 V x τ + + x 4 V 1 x 5 V 4 x 4 τ + 1 ( ) 6 V x 4 τ + + O( x 4 ), V (( + 1) x, (n + 1)) V ( x, (n + 1)) + V (( 1) x, (n + 1))/( x) = V x + 3 V x τ + 1 ( ) 4 V x τ + + x V 1 x + 5 V 4 x 4 τ + 1 ( ) 6 V x 4 τ + + O( x 4 ) Combining the results, we have σ V (( + 1) x, (n + 1)) V ( x, (n + 1)) + V (( 1) x, (n + 1)) 4 x V (( + 1) x, n) V ( x, n) + V (( 1) x, n) + x = σ V x + O( ) + O( x )
3 Similarly, we obtain 1 = 1 (r σ ) V (( + 1) x, (n + 1)) V (( 1) x, (n + 1)) V (( + 1) x, n) V (( 1) x, n) +, ) (r σ V x + O( ) + O( x ) and r V ( x, (n + 1)) + V ( x, n) = rv + O() ( ( By noting that V x, n + 1 ) ) satisfies the Black-Scholes equation, we obtain the local truncation error of the Crank-Nicholson scheme = O( ) + O( x ) 3 Instead of imposing artificial boundary conditions for the bond price at r = r min and r = r max, we enforce the condition that the option values along the boundary nodes of the computational domain remain to be governed by the bond price equation We apply backward difference operators to approximate the continuous differential operators: B r rmax=(n+1) r B r rmax=(n+1) r The corresponding explicit FTCS becomes B N+1 Bn N+1 = σ r max B N+1 5B N + 4B N 1 B N r, 3B N+1 4B N + B N 1 r B n N+1 5Bn N + 4Bn N 1 Bn N r + α(β r max ) 3Bn N+1 4Bn N + Bn N 1 r r max B n N+1 In a similar manner, we now apply the forward difference operators to approximate the differential operators: B r B 0 5B 1 + 4B B 3, rmin r B r 3B 0 + 4B 1 B rmin r We obtain the following explicit FTCS at = 0: B 0 B n 0 = σ r min B n 0 5B n 1 + 4B n B n 3 r + α(β r min ) 3Bn 0 + 4B n 1 B n r 3 r min B n 0
4 4 For the given pricing formulation of the floating strike lookback put, the binomial parameters are determined by equating the mean and variance of the discrete random walk and the continuous price process up to O( t): ) (1 α) x α x = (q r σ t, (1 α) x + α x = σ t We obtain x = σ t and This gives (1 α) x = ) (q r σ t = α = 1 + x ( r q + 1 ) σ n 1 V 1 ) (q r σ x σ 1 x n V x Here, α denotes the probability of down-move in the binomial tree The binomial scheme takes the form: V n = αv 1 + (1 α)v+1 qv n t, V n 1 = αv 1 + (1 α)v+1, q t The term qv in the governing equation for V is similar to the discount term rv in the 1 usual Black-Scholes equation This gives rise to the discount factor in the above 1 + q t binomial pricing formula When = 0, the binomial formula involves the grid point at = 1, which is outside the computational domain We then approximate the Neumann V boundary condition: (0, t) = 0 using the one-sided finite difference formula, that is, V 1 = V n 1 V 1 x 0 The numerical boundary value is given by V n 0 = 1 αv 0 + (1 α)v1 1 + q t This binomial scheme is similar to the Cheuk-Vorst scheme To complete the construction of the binomial scheme, we adopt the terminal payoff condition: V N = e x 1, 0 4
5 5 We write the FTCS scheme in the form of the standard -level-4-point scheme: ( ) σ V = S S + rs V+1 n S ( ) + 1 σ S S r V n ( ) σ + S S rs V n S 1 In order to avoid spurious oscillations, it suffices to have all the coefficients to be positive That is, σ S S + rs S > 0 1 σ S S r > 0 σ S S rs S > 0 The first inequality is always satified The last two inequalities are satisfied provided that S < σ S r and < 1 r + σ S S 6 Let the first barrier be an up-stream barrier B H while the second barrier be the downstream barrier B L The sequential barrier option is equivalent to a one-sided up-and-out barrier with a rebate paid upon knock-out The rebate is a one-sided down-and-out barrier option We assume a call payoff of the sequential barrier option Let B n and R n denote the numerical option value of the sequential barrier option and the rebate barrier option at the (, n) th node, respectively Design of the computational domain For the sequential barrier option, we set the right boundary to coincide with the upstream barrier B H The left boundary must lie sufficiently far to the left end For the down-and-out barrier option (treated as the rebate upon breaching the up-barrier in the sequential barrier option), we set the left boundary to coincide with the downstream barrier B L The right boundary must lie sufficiently far to the right end Here, we choose M such that M >> N 5
6 Boundary conditions (i) At = 0, the sequential barrier option is deep out-of-the-money so that the option value is close to zero We set B n 0 = 0, for all n (ii) At = N, which corresponds to the up-barrier B H, the sequential barrier option is apparently knocked out, receiving the down-and-out barrier option as the rebate Therefore, we have B n N = R n N, for all n (iii) At = M, the (rebate) down-and-out barrier option is deep in-the-money so that it is almost like a forward contract Accordingly, the second order derivative of the price function with respect to the stock price is close to zero That is, RM n = 5Rn M 1 4Rn M + Rn M 3, for all n (iv) At = L, the rebate barrier option is knocked out at x = B L with zero value That is R n L = 0, for all n Both price functions of the sequential barrier option and the rebate barrier option satisfy the Black-Scholes equation We start with computation of the rebate barrier option, then continue with the sequential barrier option The explicit finite difference schemes in both option calculations take an identical form: R = B = ( where µ = σ and c = x µ + c Rn +1 + (1 µ)r n + µ c Rn 1 e r, = L + 1,, M 1, n = 0, 1,, µ + c Bn +1 + (1 µ)b n + µ c Bn 1 e r, = 1,, N 1, n = 0, 1,,, ) r q σ x 6
7 7 In the continuation region where V > h, the penalty term vanishes However, if the option remains unexercised when it should be optimally exercised (in the region of optimal stopping), we have V < h As a result, the penalty term becomes highly dominant since the penalty parameter ρ is typically chosen to be a very large positive parameter Intuitively, we observe that V should assume the exercise payoff h in the limit ρ Suppose we use S as the state variable and let V n denote the numerical option value at the node ( S, n), where S = S + S 0 and S 0 is the initial stock price The Crank-Nicolson scheme takes the form: V V +1 V n = σ S V + V 1 + V n S V +1 + (r q)s V 1 + V n S r V + V n +1 V n S +1 V 1 n S + ρ max{h(s ) V + V n + V n 1 Due to the non-linearity exhibited in the penalty term, we cannot apply the Thomas algorithm Indeed, we need to use an iteration method to solve the resulting non-linear algebraic system of equations 8 The exponential distribution is the distribution of the random time between two successive umps of a Poisson process with rate 1/θ Inverting the exponential distribution gives This can be implemented as X = θ ln(1 U) X = θ ln(u) since U and 1 U have the same distribution As a summary, with U U(0, 1), we have X = θ ln U exp(θ) 9 Write ϵ 1 α x 1 ϵ = α 1 α 0 x ϵ 3 α 31 α 3 α 33 x 3 where α M = α 1 α 0 α 31 α 3 α 33 The entries in M are determined by the relation: MM T = Σ This yields α MM T 11 α 11 α 1 α 31 α = α 11 α 1 α1 + α α 1 α 31 + α α 3 = , α 31 α 11 α 1 α 31 + α α 3 α31 + α3 + α , 0}
8 and accordingly, α 11 = 1, α 11 α 1 = 06, α 31 α 11 = 05 α 1 + α = 1, α 1 α 31 + α α 3 = 07, α 31 + α 3 + α 33 = 1 We take α 11 = 1 (we choose the positive root without loss of generality), then Also, α 31 = 05, then α 1 = 06, α = 1 06 = 08 (choosing the positive root) (06)(05) + 08α 3 = 07 so that α 3 = = Lastly, α 33 = 1 α 31 α 3 = 1/ (choosing the positive root) Hence 10 Since c AV = c + c so that ( ) c + c var(c AV ) = var ϵ x 1 ϵ = ϵ / x x 3 = 1 4 var(c) var( c) + 1 cov(c, c) As c is generated using ϵ (i) while c is generated using ϵ (i), we expect to have var(c) = var( c) so that var(c AV ) = 1 var(c) + 1 cov(c, c) (A) Now, we apply the following criterion of determining the trade-off between computational work units and variances: σ1/σ < W /W 1 Since the amount of computational work to compute c AV is about twice that of c, the control variate is preferred in terms of computational efficiency provided that Based on Eq (A), Ineq (B) is equivalent to var(c AV ) < var(c) (B) cov(c, c) < 0 Since we have chosen ϵ (i) for computing c i, the chances are high that c i and c i are negatively correlated Hence, the antithetic variates method improves computational efficiency 8
9 11 E Ŷ = E t ( Z Z ) = 0; ( ) E Ŷ = te Z Z = t 4 E Ŷ Ŵ = t E Z 1 (Z Z 3 ) = t Correlation coefficient between Ŷ t E Ŷ Ŵ E Ŷ E Ŵ = = var( Ŷ ) var( Ŵ ) 1 We consider t 3 var(x t ) = E = 3 t ( ) = t ; 3 and Ŵ ( W t t ) T W T E W t t T W T = EW t t T EW tw T + t T EW T = t t T EW t(w T W t ) + W t + t T T Since the Brownian increments over non-overlapping time intervals are independent, so Hence Let Y t = EW t (W T W t ) = 0 var(x t ) = t t ( T t + t T T = t t T = t 1 t ) T t ( ) 1 t T Z, Z N(0, 1), we have ( EY t = EX t = 0 and var(y t ) = t 1 t ) = var(x t ) T Also, Y 0 = X 0 = 0 and Y T X T = 0; so t ( ) 1 t T Z, with Z N(0, 1), is a realization of X t 13 (a) Note that var( X Y ) = 1 N var(x Y ) = 1 N var(x) + var(y ) cov(x, Y ) We achieve a reduction of the variance for the control variate estimator provided that var(x) > var(x Y ), which is equivalent to cov(x, Y ) > var(y ) If Y is very close to X, then we have the satisfaction of the inequality This would imply the significant reduction of the variance of the crude Monte Carlo estimator by using its control variate variable (b) The efficiency of the control variate hinges on the fact that we can directly simulate the difference X i Y i as one random variable It is supposed that we know its exact distribution, so one needs to generate one random variable 9
10 To obtain the confidence interval for the control variate Monte Carlo estimator, one starts with that of the crude Monte Carlo estimator for EX Y and add EY to the interval For example, we obtain an approximate 95%-confidence interval by X Y 196 ˆσ X Y, XY ˆσ X Y, N N where ˆσ X Y = 1 N 1 N X i Y i 1 N i=1 N (X Y ) =1 10
6. Numerical methods for option pricing
6. Numerical methods for option pricing Binomial model revisited Under the risk neutral measure, ln S t+ t ( ) S t becomes normally distributed with mean r σ2 t and variance σ 2 t, where r is 2 the riskless
More informationAdvanced Numerical Methods
Advanced Numerical Methods Solution to Homework One Course instructor: Prof. Y.K. Kwok. When the asset pays continuous dividend yield at the rate q the expected rate of return of the asset is r q under
More informationMATH6911: Numerical Methods in Finance. Final exam Time: 2:00pm - 5:00pm, April 11, Student Name (print): Student Signature: Student ID:
MATH6911 Page 1 of 16 Winter 2007 MATH6911: Numerical Methods in Finance Final exam Time: 2:00pm - 5:00pm, April 11, 2007 Student Name (print): Student Signature: Student ID: Question Full Mark Mark 1
More informationFE610 Stochastic Calculus for Financial Engineers. Stevens Institute of Technology
FE610 Stochastic Calculus for Financial Engineers Lecture 13. The Black-Scholes PDE Steve Yang Stevens Institute of Technology 04/25/2013 Outline 1 The Black-Scholes PDE 2 PDEs in Asset Pricing 3 Exotic
More informationUtility Indifference Pricing and Dynamic Programming Algorithm
Chapter 8 Utility Indifference ricing and Dynamic rogramming Algorithm In the Black-Scholes framework, we can perfectly replicate an option s payoff. However, it may not be true beyond the Black-Scholes
More informationNumerical Methods in Option Pricing (Part III)
Numerical Methods in Option Pricing (Part III) E. Explicit Finite Differences. Use of the Forward, Central, and Symmetric Central a. In order to obtain an explicit solution for the price of the derivative,
More informationMASM006 UNIVERSITY OF EXETER SCHOOL OF ENGINEERING, COMPUTER SCIENCE AND MATHEMATICS MATHEMATICAL SCIENCES FINANCIAL MATHEMATICS.
MASM006 UNIVERSITY OF EXETER SCHOOL OF ENGINEERING, COMPUTER SCIENCE AND MATHEMATICS MATHEMATICAL SCIENCES FINANCIAL MATHEMATICS May/June 2006 Time allowed: 2 HOURS. Examiner: Dr N.P. Byott This is a CLOSED
More informationMAFS5250 Computational Methods for Pricing Structured Products Topic 5 - Monte Carlo simulation
MAFS5250 Computational Methods for Pricing Structured Products Topic 5 - Monte Carlo simulation 5.1 General formulation of the Monte Carlo procedure Expected value and variance of the estimate Multistate
More informationCh 5. Several Numerical Methods
Ch 5 Several Numerical Methods I Monte Carlo Simulation for Multiple Variables II Confidence Interval and Variance Reduction III Solving Systems of Linear Equations IV Finite Difference Method ( 有限差分法
More informationMATH4143: Scientific Computations for Finance Applications Final exam Time: 9:00 am - 12:00 noon, April 18, Student Name (print):
MATH4143 Page 1 of 17 Winter 2007 MATH4143: Scientific Computations for Finance Applications Final exam Time: 9:00 am - 12:00 noon, April 18, 2007 Student Name (print): Student Signature: Student ID: Question
More information2.1 Mathematical Basis: Risk-Neutral Pricing
Chapter Monte-Carlo Simulation.1 Mathematical Basis: Risk-Neutral Pricing Suppose that F T is the payoff at T for a European-type derivative f. Then the price at times t before T is given by f t = e r(t
More informationFrom Discrete Time to Continuous Time Modeling
From Discrete Time to Continuous Time Modeling Prof. S. Jaimungal, Department of Statistics, University of Toronto 2004 Arrow-Debreu Securities 2004 Prof. S. Jaimungal 2 Consider a simple one-period economy
More informationThe Binomial Lattice Model for Stocks: Introduction to Option Pricing
1/33 The Binomial Lattice Model for Stocks: Introduction to Option Pricing Professor Karl Sigman Columbia University Dept. IEOR New York City USA 2/33 Outline The Binomial Lattice Model (BLM) as a Model
More informationComputational Finance Finite Difference Methods
Explicit finite difference method Computational Finance Finite Difference Methods School of Mathematics 2018 Today s Lecture We now introduce the final numerical scheme which is related to the PDE solution.
More informationFINITE DIFFERENCE METHODS
FINITE DIFFERENCE METHODS School of Mathematics 2013 OUTLINE Review 1 REVIEW Last time Today s Lecture OUTLINE Review 1 REVIEW Last time Today s Lecture 2 DISCRETISING THE PROBLEM Finite-difference approximations
More informationIEOR E4703: Monte-Carlo Simulation
IEOR E4703: Monte-Carlo Simulation Simulation Efficiency and an Introduction to Variance Reduction Methods Martin Haugh Department of Industrial Engineering and Operations Research Columbia University
More informationChapter 5 Finite Difference Methods. Math6911 W07, HM Zhu
Chapter 5 Finite Difference Methods Math69 W07, HM Zhu References. Chapters 5 and 9, Brandimarte. Section 7.8, Hull 3. Chapter 7, Numerical analysis, Burden and Faires Outline Finite difference (FD) approximation
More informationQueens College, CUNY, Department of Computer Science Computational Finance CSCI 365 / 765 Spring 2018 Instructor: Dr. Sateesh Mane.
Queens College, CUNY, Department of Computer Science Computational Finance CSCI 365 / 765 Spring 218 Instructor: Dr. Sateesh Mane c Sateesh R. Mane 218 19 Lecture 19 May 12, 218 Exotic options The term
More informationKing s College London
King s College London University Of London This paper is part of an examination of the College counting towards the award of a degree. Examinations are governed by the College Regulations under the authority
More informationAmerican Equity Option Valuation Practical Guide
Valuation Practical Guide John Smith FinPricing Summary American Equity Option Introduction The Use of American Equity Options Valuation Practical Guide A Real World Example American Option Introduction
More informationKing s College London
King s College London University Of London This paper is part of an examination of the College counting towards the award of a degree. Examinations are governed by the College Regulations under the authority
More informationResults for option pricing
Results for option pricing [o,v,b]=optimal(rand(1,100000 Estimators = 0.4619 0.4617 0.4618 0.4613 0.4619 o = 0.46151 % best linear combination (true value=0.46150 v = 1.1183e-005 %variance per uniform
More informationAdvanced Topics in Derivative Pricing Models. Topic 4 - Variance products and volatility derivatives
Advanced Topics in Derivative Pricing Models Topic 4 - Variance products and volatility derivatives 4.1 Volatility trading and replication of variance swaps 4.2 Volatility swaps 4.3 Pricing of discrete
More informationGamma. The finite-difference formula for gamma is
Gamma The finite-difference formula for gamma is [ P (S + ɛ) 2 P (S) + P (S ɛ) e rτ E ɛ 2 ]. For a correlation option with multiple underlying assets, the finite-difference formula for the cross gammas
More informationVersion A. Problem 1. Let X be the continuous random variable defined by the following pdf: 1 x/2 when 0 x 2, f(x) = 0 otherwise.
Math 224 Q Exam 3A Fall 217 Tues Dec 12 Version A Problem 1. Let X be the continuous random variable defined by the following pdf: { 1 x/2 when x 2, f(x) otherwise. (a) Compute the mean µ E[X]. E[X] x
More information1.1 Basic Financial Derivatives: Forward Contracts and Options
Chapter 1 Preliminaries 1.1 Basic Financial Derivatives: Forward Contracts and Options A derivative is a financial instrument whose value depends on the values of other, more basic underlying variables
More informationOption Pricing Models for European Options
Chapter 2 Option Pricing Models for European Options 2.1 Continuous-time Model: Black-Scholes Model 2.1.1 Black-Scholes Assumptions We list the assumptions that we make for most of this notes. 1. The underlying
More informationBinomial model: numerical algorithm
Binomial model: numerical algorithm S / 0 C \ 0 S0 u / C \ 1,1 S0 d / S u 0 /, S u 3 0 / 3,3 C \ S0 u d /,1 S u 5 0 4 0 / C 5 5,5 max X S0 u,0 S u C \ 4 4,4 C \ 3 S u d / 0 3, C \ S u d 0 S u d 0 / C 4
More informationLecture 4. Finite difference and finite element methods
Finite difference and finite element methods Lecture 4 Outline Black-Scholes equation From expectation to PDE Goal: compute the value of European option with payoff g which is the conditional expectation
More informationAMH4 - ADVANCED OPTION PRICING. Contents
AMH4 - ADVANCED OPTION PRICING ANDREW TULLOCH Contents 1. Theory of Option Pricing 2 2. Black-Scholes PDE Method 4 3. Martingale method 4 4. Monte Carlo methods 5 4.1. Method of antithetic variances 5
More informationThe Multistep Binomial Model
Lecture 10 The Multistep Binomial Model Reminder: Mid Term Test Friday 9th March - 12pm Examples Sheet 1 4 (not qu 3 or qu 5 on sheet 4) Lectures 1-9 10.1 A Discrete Model for Stock Price Reminder: The
More information1. In this exercise, we can easily employ the equations (13.66) (13.70), (13.79) (13.80) and
CHAPTER 13 Solutions Exercise 1 1. In this exercise, we can easily employ the equations (13.66) (13.70), (13.79) (13.80) and (13.82) (13.86). Also, remember that BDT model will yield a recombining binomial
More informationOption Pricing Models. c 2013 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 205
Option Pricing Models c 2013 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 205 If the world of sense does not fit mathematics, so much the worse for the world of sense. Bertrand Russell (1872 1970)
More informationMonte Carlo Methods in Structuring and Derivatives Pricing
Monte Carlo Methods in Structuring and Derivatives Pricing Prof. Manuela Pedio (guest) 20263 Advanced Tools for Risk Management and Pricing Spring 2017 Outline and objectives The basic Monte Carlo algorithm
More informationTrinomial Tree. Set up a trinomial approximation to the geometric Brownian motion ds/s = r dt + σ dw. a
Trinomial Tree Set up a trinomial approximation to the geometric Brownian motion ds/s = r dt + σ dw. a The three stock prices at time t are S, Su, and Sd, where ud = 1. Impose the matching of mean and
More informationReal Options and Game Theory in Incomplete Markets
Real Options and Game Theory in Incomplete Markets M. Grasselli Mathematics and Statistics McMaster University IMPA - June 28, 2006 Strategic Decision Making Suppose we want to assign monetary values to
More informationFinancial Risk Management
Financial Risk Management Professor: Thierry Roncalli Evry University Assistant: Enareta Kurtbegu Evry University Tutorial exercices #4 1 Correlation and copulas 1. The bivariate Gaussian copula is given
More informationFINANCIAL OPTION ANALYSIS HANDOUTS
FINANCIAL OPTION ANALYSIS HANDOUTS 1 2 FAIR PRICING There is a market for an object called S. The prevailing price today is S 0 = 100. At this price the object S can be bought or sold by anyone for any
More informationAmerican options and early exercise
Chapter 3 American options and early exercise American options are contracts that may be exercised early, prior to expiry. These options are contrasted with European options for which exercise is only
More informationUniversity of California, Los Angeles Department of Statistics. Final exam 07 June 2013
University of California, Los Angeles Department of Statistics Statistics C183/C283 Instructor: Nicolas Christou Final exam 07 June 2013 Name: Problem 1 (20 points) a. Suppose the variable X follows the
More informationCHAPTER 6 Numerical Schemes for Pricing Options
CHAPTER 6 Numerical Schemes for Pricing Options In previous chapters, closed form price formulas for a variety of option models have been obtained. However, option models which lend themselves to a closed
More informationDefinition 9.1 A point estimate is any function T (X 1,..., X n ) of a random sample. We often write an estimator of the parameter θ as ˆθ.
9 Point estimation 9.1 Rationale behind point estimation When sampling from a population described by a pdf f(x θ) or probability function P [X = x θ] knowledge of θ gives knowledge of the entire population.
More informationStochastic Differential Equations in Finance and Monte Carlo Simulations
Stochastic Differential Equations in Finance and Department of Statistics and Modelling Science University of Strathclyde Glasgow, G1 1XH China 2009 Outline Stochastic Modelling in Asset Prices 1 Stochastic
More informationThe Binomial Lattice Model for Stocks: Introduction to Option Pricing
1/27 The Binomial Lattice Model for Stocks: Introduction to Option Pricing Professor Karl Sigman Columbia University Dept. IEOR New York City USA 2/27 Outline The Binomial Lattice Model (BLM) as a Model
More informationEquity correlations implied by index options: estimation and model uncertainty analysis
1/18 : estimation and model analysis, EDHEC Business School (joint work with Rama COT) Modeling and managing financial risks Paris, 10 13 January 2011 2/18 Outline 1 2 of multi-asset models Solution to
More informationOptions Pricing Using Combinatoric Methods Postnikov Final Paper
Options Pricing Using Combinatoric Methods 18.04 Postnikov Final Paper Annika Kim May 7, 018 Contents 1 Introduction The Lattice Model.1 Overview................................ Limitations of the Lattice
More informationQueens College, CUNY, Department of Computer Science Computational Finance CSCI 365 / 765 Fall 2017 Instructor: Dr. Sateesh Mane.
Queens College, CUNY, Department of Computer Science Computational Finance CSCI 365 / 765 Fall 217 Instructor: Dr. Sateesh Mane c Sateesh R. Mane 217 13 Lecture 13 November 15, 217 Derivation of the Black-Scholes-Merton
More informationIEOR E4703: Monte-Carlo Simulation
IEOR E4703: Monte-Carlo Simulation Simulating Stochastic Differential Equations Martin Haugh Department of Industrial Engineering and Operations Research Columbia University Email: martin.b.haugh@gmail.com
More informationInfinite Reload Options: Pricing and Analysis
Infinite Reload Options: Pricing and Analysis A. C. Bélanger P. A. Forsyth April 27, 2006 Abstract Infinite reload options allow the user to exercise his reload right as often as he chooses during the
More informationLecture Quantitative Finance Spring Term 2015
and Lecture Quantitative Finance Spring Term 2015 Prof. Dr. Erich Walter Farkas Lecture 06: March 26, 2015 1 / 47 Remember and Previous chapters: introduction to the theory of options put-call parity fundamentals
More informationIntroduction Dickey-Fuller Test Option Pricing Bootstrapping. Simulation Methods. Chapter 13 of Chris Brook s Book.
Simulation Methods Chapter 13 of Chris Brook s Book Christopher Ting http://www.mysmu.edu/faculty/christophert/ Christopher Ting : christopherting@smu.edu.sg : 6828 0364 : LKCSB 5036 April 26, 2017 Christopher
More informationStochastic Modelling in Finance
in Finance Department of Mathematics and Statistics University of Strathclyde Glasgow, G1 1XH April 2010 Outline and Probability 1 and Probability 2 Linear modelling Nonlinear modelling 3 The Black Scholes
More informationQueens College, CUNY, Department of Computer Science Computational Finance CSCI 365 / 765 Fall 2017 Instructor: Dr. Sateesh Mane.
Queens College, CUNY, Department of Computer Science Computational Finance CSCI 365 / 765 Fall 2017 Instructor: Dr. Sateesh Mane c Sateesh R. Mane 2017 14 Lecture 14 November 15, 2017 Derivation of the
More informationChapter 5. Statistical inference for Parametric Models
Chapter 5. Statistical inference for Parametric Models Outline Overview Parameter estimation Method of moments How good are method of moments estimates? Interval estimation Statistical Inference for Parametric
More informationComputational Finance. Computational Finance p. 1
Computational Finance Computational Finance p. 1 Outline Binomial model: option pricing and optimal investment Monte Carlo techniques for pricing of options pricing of non-standard options improving accuracy
More informationAN OPERATOR SPLITTING METHOD FOR PRICING THE ELS OPTION
J. KSIAM Vol.14, No.3, 175 187, 21 AN OPERATOR SPLITTING METHOD FOR PRICING THE ELS OPTION DARAE JEONG, IN-SUK WEE, AND JUNSEOK KIM DEPARTMENT OF MATHEMATICS, KOREA UNIVERSITY, SEOUL 136-71, KOREA E-mail
More informationRandom Tree Method. Monte Carlo Methods in Financial Engineering
Random Tree Method Monte Carlo Methods in Financial Engineering What is it for? solve full optimal stopping problem & estimate value of the American option simulate paths of underlying Markov chain produces
More informationOptimal Search for Parameters in Monte Carlo Simulation for Derivative Pricing
Optimal Search for Parameters in Monte Carlo Simulation for Derivative Pricing Prof. Chuan-Ju Wang Department of Computer Science University of Taipei Joint work with Prof. Ming-Yang Kao March 28, 2014
More informationCONVERGENCE OF NUMERICAL METHODS FOR VALUING PATH-DEPENDENT OPTIONS USING INTERPOLATION
CONVERGENCE OF NUMERICAL METHODS FOR VALUING PATH-DEPENDENT OPTIONS USING INTERPOLATION P.A. Forsyth Department of Computer Science University of Waterloo Waterloo, ON Canada N2L 3G1 E-mail: paforsyt@elora.math.uwaterloo.ca
More informationHomework Assignments
Homework Assignments Week 1 (p 57) #4.1, 4., 4.3 Week (pp 58-6) #4.5, 4.6, 4.8(a), 4.13, 4.0, 4.6(b), 4.8, 4.31, 4.34 Week 3 (pp 15-19) #1.9, 1.1, 1.13, 1.15, 1.18 (pp 9-31) #.,.6,.9 Week 4 (pp 36-37)
More informationEFFICIENT MONTE CARLO ALGORITHM FOR PRICING BARRIER OPTIONS
Commun. Korean Math. Soc. 23 (2008), No. 2, pp. 285 294 EFFICIENT MONTE CARLO ALGORITHM FOR PRICING BARRIER OPTIONS Kyoung-Sook Moon Reprinted from the Communications of the Korean Mathematical Society
More information( ) since this is the benefit of buying the asset at the strike price rather
Review of some financial models for MAT 483 Parity and Other Option Relationships The basic parity relationship for European options with the same strike price and the same time to expiration is: C( KT
More informationMATH3075/3975 FINANCIAL MATHEMATICS TUTORIAL PROBLEMS
MATH307/37 FINANCIAL MATHEMATICS TUTORIAL PROBLEMS School of Mathematics and Statistics Semester, 04 Tutorial problems should be used to test your mathematical skills and understanding of the lecture material.
More informationProject 1: Double Pendulum
Final Projects Introduction to Numerical Analysis II http://www.math.ucsb.edu/ atzberg/winter2009numericalanalysis/index.html Professor: Paul J. Atzberger Due: Friday, March 20th Turn in to TA s Mailbox:
More informationA distributed Laplace transform algorithm for European options
A distributed Laplace transform algorithm for European options 1 1 A. J. Davies, M. E. Honnor, C.-H. Lai, A. K. Parrott & S. Rout 1 Department of Physics, Astronomy and Mathematics, University of Hertfordshire,
More informationContents Critique 26. portfolio optimization 32
Contents Preface vii 1 Financial problems and numerical methods 3 1.1 MATLAB environment 4 1.1.1 Why MATLAB? 5 1.2 Fixed-income securities: analysis and portfolio immunization 6 1.2.1 Basic valuation of
More informationWeek 7 Quantitative Analysis of Financial Markets Simulation Methods
Week 7 Quantitative Analysis of Financial Markets Simulation Methods Christopher Ting http://www.mysmu.edu/faculty/christophert/ Christopher Ting : christopherting@smu.edu.sg : 6828 0364 : LKCSB 5036 November
More informationMATH 476/567 ACTUARIAL RISK THEORY FALL 2016 PROFESSOR WANG
MATH 476/567 ACTUARIAL RISK THEORY FALL 206 PROFESSOR WANG Homework 5 (max. points = 00) Due at the beginning of class on Tuesday, November 8, 206 You are encouraged to work on these problems in groups
More informationStrategies for Improving the Efficiency of Monte-Carlo Methods
Strategies for Improving the Efficiency of Monte-Carlo Methods Paul J. Atzberger General comments or corrections should be sent to: paulatz@cims.nyu.edu Introduction The Monte-Carlo method is a useful
More informationFinal Exam Key, JDEP 384H, Spring 2006
Final Exam Key, JDEP 384H, Spring 2006 Due Date for Exam: Thursday, May 4, 12:00 noon. Instructions: Show your work and give reasons for your answers. Write out your solutions neatly and completely. There
More informationNEWCASTLE UNIVERSITY SCHOOL OF MATHEMATICS, STATISTICS & PHYSICS SEMESTER 1 SPECIMEN 2 MAS3904. Stochastic Financial Modelling. Time allowed: 2 hours
NEWCASTLE UNIVERSITY SCHOOL OF MATHEMATICS, STATISTICS & PHYSICS SEMESTER 1 SPECIMEN 2 Stochastic Financial Modelling Time allowed: 2 hours Candidates should attempt all questions. Marks for each question
More informationExotic Derivatives & Structured Products. Zénó Farkas (MSCI)
Exotic Derivatives & Structured Products Zénó Farkas (MSCI) Part 1: Exotic Derivatives Over the counter products Generally more profitable (and more risky) than vanilla derivatives Why do they exist? Possible
More informationThe Pennsylvania State University. The Graduate School. Department of Industrial Engineering AMERICAN-ASIAN OPTION PRICING BASED ON MONTE CARLO
The Pennsylvania State University The Graduate School Department of Industrial Engineering AMERICAN-ASIAN OPTION PRICING BASED ON MONTE CARLO SIMULATION METHOD A Thesis in Industrial Engineering and Operations
More informationOptimal Stopping. Nick Hay (presentation follows Thomas Ferguson s Optimal Stopping and Applications) November 6, 2008
(presentation follows Thomas Ferguson s and Applications) November 6, 2008 1 / 35 Contents: Introduction Problems Markov Models Monotone Stopping Problems Summary 2 / 35 The Secretary problem You have
More informationFinancial derivatives exam Winter term 2014/2015
Financial derivatives exam Winter term 2014/2015 Problem 1: [max. 13 points] Determine whether the following assertions are true or false. Write your answers, without explanations. Grading: correct answer
More informationLecture 4 - Finite differences methods for PDEs
Finite diff. Lecture 4 - Finite differences methods for PDEs Lina von Sydow Finite differences, Lina von Sydow, (1 : 18) Finite difference methods Finite diff. Black-Scholes equation @v @t + 1 2 2 s 2
More information1 Geometric Brownian motion
Copyright c 05 by Karl Sigman Geometric Brownian motion Note that since BM can take on negative values, using it directly for modeling stock prices is questionable. There are other reasons too why BM is
More informationMonte Carlo Methods. Prof. Mike Giles. Oxford University Mathematical Institute. Lecture 1 p. 1.
Monte Carlo Methods Prof. Mike Giles mike.giles@maths.ox.ac.uk Oxford University Mathematical Institute Lecture 1 p. 1 Geometric Brownian Motion In the case of Geometric Brownian Motion ds t = rs t dt+σs
More information- 1 - **** d(lns) = (µ (1/2)σ 2 )dt + σdw t
- 1 - **** These answers indicate the solutions to the 2014 exam questions. Obviously you should plot graphs where I have simply described the key features. It is important when plotting graphs to label
More informationCRANK-NICOLSON SCHEME FOR ASIAN OPTION
CRANK-NICOLSON SCHEME FOR ASIAN OPTION By LEE TSE YUENG A thesis submitted to the Department of Mathematical and Actuarial Sciences, Faculty of Engineering and Science, Universiti Tunku Abdul Rahman, in
More informationSingular Stochastic Control Models for Optimal Dynamic Withdrawal Policies in Variable Annuities
1/ 46 Singular Stochastic Control Models for Optimal Dynamic Withdrawal Policies in Variable Annuities Yue Kuen KWOK Department of Mathematics Hong Kong University of Science and Technology * Joint work
More information1 The continuous time limit
Derivative Securities, Courant Institute, Fall 2008 http://www.math.nyu.edu/faculty/goodman/teaching/derivsec08/index.html Jonathan Goodman and Keith Lewis Supplementary notes and comments, Section 3 1
More informationShort-time-to-expiry expansion for a digital European put option under the CEV model. November 1, 2017
Short-time-to-expiry expansion for a digital European put option under the CEV model November 1, 2017 Abstract In this paper I present a short-time-to-expiry asymptotic series expansion for a digital European
More informationTrinomial Tree. Set up a trinomial approximation to the geometric Brownian motion ds/s = r dt + σ dw. a
Trinomial Tree Set up a trinomial approximation to the geometric Brownian motion ds/s = r dt + σ dw. a The three stock prices at time t are S, Su, and Sd, where ud = 1. Impose the matching of mean and
More informationMulti-Asset Options. A Numerical Study VILHELM NIKLASSON FRIDA TIVEDAL. Master s thesis in Engineering Mathematics and Computational Science
Multi-Asset Options A Numerical Study Master s thesis in Engineering Mathematics and Computational Science VILHELM NIKLASSON FRIDA TIVEDAL Department of Mathematical Sciences Chalmers University of Technology
More informationPRMIA Exam 8002 PRM Certification - Exam II: Mathematical Foundations of Risk Measurement Version: 6.0 [ Total Questions: 132 ]
s@lm@n PRMIA Exam 8002 PRM Certification - Exam II: Mathematical Foundations of Risk Measurement Version: 6.0 [ Total Questions: 132 ] Question No : 1 A 2-step binomial tree is used to value an American
More informationMartingale Pricing Theory in Discrete-Time and Discrete-Space Models
IEOR E4707: Foundations of Financial Engineering c 206 by Martin Haugh Martingale Pricing Theory in Discrete-Time and Discrete-Space Models These notes develop the theory of martingale pricing in a discrete-time,
More informationHomework Assignments
Homework Assignments Week 1 (p. 57) #4.1, 4., 4.3 Week (pp 58 6) #4.5, 4.6, 4.8(a), 4.13, 4.0, 4.6(b), 4.8, 4.31, 4.34 Week 3 (pp 15 19) #1.9, 1.1, 1.13, 1.15, 1.18 (pp 9 31) #.,.6,.9 Week 4 (pp 36 37)
More informationHomework Set 6 Solutions
MATH 667-010 Introduction to Mathematical Finance Prof. D. A. Edwards Due: Apr. 11, 018 P Homework Set 6 Solutions K z K + z S 1. The payoff diagram shown is for a strangle. Denote its option value by
More informationMultiple Optimal Stopping Problems and Lookback Options
Multiple Optimal Stopping Problems and Lookback Options Yue Kuen KWOK Department of Mathematics Hong Kong University of Science & Technology Hong Kong, China web page: http://www.math.ust.hk/ maykwok/
More informationCHOICE THEORY, UTILITY FUNCTIONS AND RISK AVERSION
CHOICE THEORY, UTILITY FUNCTIONS AND RISK AVERSION Szabolcs Sebestyén szabolcs.sebestyen@iscte.pt Master in Finance INVESTMENTS Sebestyén (ISCTE-IUL) Choice Theory Investments 1 / 65 Outline 1 An Introduction
More informationHedging with Life and General Insurance Products
Hedging with Life and General Insurance Products June 2016 2 Hedging with Life and General Insurance Products Jungmin Choi Department of Mathematics East Carolina University Abstract In this study, a hybrid
More informationCONVERGENCE OF OPTION REWARDS FOR MARKOV TYPE PRICE PROCESSES MODULATED BY STOCHASTIC INDICES
CONVERGENCE OF OPTION REWARDS FOR MARKOV TYPE PRICE PROCESSES MODULATED BY STOCHASTIC INDICES D. S. SILVESTROV, H. JÖNSSON, AND F. STENBERG Abstract. A general price process represented by a two-component
More information4. Black-Scholes Models and PDEs. Math6911 S08, HM Zhu
4. Black-Scholes Models and PDEs Math6911 S08, HM Zhu References 1. Chapter 13, J. Hull. Section.6, P. Brandimarte Outline Derivation of Black-Scholes equation Black-Scholes models for options Implied
More informationLecture 3: Review of mathematical finance and derivative pricing models
Lecture 3: Review of mathematical finance and derivative pricing models Xiaoguang Wang STAT 598W January 21th, 2014 (STAT 598W) Lecture 3 1 / 51 Outline 1 Some model independent definitions and principals
More informationPricing Implied Volatility
Pricing Implied Volatility Expected future volatility plays a central role in finance theory. Consequently, accurate estimation of this parameter is crucial to meaningful financial decision-making. Researchers
More informationExtend the ideas of Kan and Zhou paper on Optimal Portfolio Construction under parameter uncertainty
Extend the ideas of Kan and Zhou paper on Optimal Portfolio Construction under parameter uncertainty George Photiou Lincoln College University of Oxford A dissertation submitted in partial fulfilment for
More informationCS476/676 Mar 6, Today s Topics. American Option: early exercise curve. PDE overview. Discretizations. Finite difference approximations
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
More informationMath 623 (IOE 623), Winter 2008: Final exam
Math 623 (IOE 623), Winter 2008: Final exam Name: Student ID: This is a closed book exam. You may bring up to ten one sided A4 pages of notes to the exam. You may also use a calculator but not its memory
More informationTopic 2 Implied binomial trees and calibration of interest rate trees. 2.1 Implied binomial trees of fitting market data of option prices
MAFS5250 Computational Methods for Pricing Structured Products Topic 2 Implied binomial trees and calibration of interest rate trees 2.1 Implied binomial trees of fitting market data of option prices Arrow-Debreu
More information