23 Stochastic Ordinary Differential Equations with Examples from Finance
|
|
- Eustace Hall
- 5 years ago
- Views:
Transcription
1 23 Stochastic Ordinary Differential Equations with Examples from Finance Scraping Financial Data from the Web The MATLAB/Octave yahoo function below returns daily open, high, low, close, and adjusted close price data, along with trading volume data for a specified financial instrument and date range. function A = yahoo(symbol, startdate, enddate) % YAHOO: Retrieve daily HLOC financial data from Yahoo! Finance. % A = yahoo(symbol, startdate, enddate) % symbol A string-valued financial instrument name. % startdate The starting date string in the form YYYY-MM-DD % enddate The ending date string in the form YYYY-MM-DD % The return value A is an nx9 numeric matrix, where n is the number of % days data returned, and the columns are: % year, month, day, open, high, low, close, volume, adjusted close a=sscanf(startdate, %f-%f-%f ); b=sscanf(enddate, %f-%f-%f ); request=strcat( &a=,num2str(a(2)), &b=,num2str(a(3)), &c=,num2str(a(1)),... &d=,num2str(b(2)), &e=,num2str(b(3)), &f=,num2str(b(1)),... &g=d&ignore=.csv ); x=urlread(request); y=x(43:end); n=length(regexp(y, \n )); A=sscanf(y, %f-%f-%f,%f,%f,%f,%f,%f,%f,[9,n]); A=flipud(A ); We can, for example, plot the daily adjusted closing price for Google s stock over all of 2009 with: goog=yahoo( GOOG, , )(:,9); plot(goog) Consider the simple initial value problem from Example 22.1: 0 Version April 24, 2014 d dt u(t) = λu(t), λ 0, T t 0, u(0) = u 0, 1
2 or, equivalently in differential form, du(t) = λu(t)dt. (1) This simple ODE models the value of an asset value u(t) as it appreciates over time with a continually compounded interest rate λ. There is no uncertainty in the initial value problem (1); the value of the asset can be computed precisely at any time t. The model is appropriate for financial instruments such as interest-bearing savings accounts and fixed-interest loans. But equation (1) is not good at modeling instruments with risk such as stock prices, as illustrated in the following example. Example 1 Model Google s stock price over 2009 with equation (1). goog=yahoo( GOOG, , )(:,9); n=length(goog); mu = log(goog(n)/goog(1)); t=linspace(0,1,n) ; u=goog(1)*exp(mu*t); plot(1:n,goog, -b,1:n,u, -r ) Although the elementary model derived from equation (1) captures the overall growth in the price of Google s stock over the period, it obviously lacks the many small up and down movements along the way. We will explore modifications to equation (1) that can be used to model a much wider class of problems that admit some level of uncertainty (risk). Brownian motion The botanist Robert Brown observed the erratic motion of pollen particles floating in water in He conducted experiments to rule out self-locomotion and described his results. Although similar observations had been made earlier in history, the random motion of particles is generally called Brownian motion. Some of the greatest thinkers of the 19th and 20th century worked out mathematical models to describe Brownian motion, including Thiele, Einstein, Wiener, and many others. One mathematical formulation of Brownian motion is given by the Wiener process. A Wiener process is a random variable W(t) that depends continuously on t 0 and satisfies the conditions: 1. W(0) = 0 2. W(t) is a continuous function of t. 3. For any 0 s < t, the Brownian increment W(t) W(s) is a normally-distributed random variable with zero mean and variance t s. 2
3 Figure 1: Google s stock price over 2009 (blue) and a simple ODE model (red). 4. Over non-overlapping intervals, the Brownian increments are mutually independent random variables. Computationally, we approximate realizations of W(t) over discrete time intervals by taking a step in a direction sampled from a normal distribution with mean 0 and variance one, scaled by the square root of the step length. The following MATLAB code illustrates this procedure: n=1000; % Number of time steps t=linspace(0,1,n) ; % Time t from 0 to 1. z=randn(n,1); % z contains n pseudorandom variables sampled % from N(0,1) W=(1/sqrt(n)) * cumsum(z); % W is an approximate realized sample path plot (h,b); % Plot the result 3
4 Figure 2: Ten sample paths of the standard Brownian motion process. Figure 2 displays ten sample paths of the standard Brownian motion process computed using the above method. The sample paths are extremely jagged, visually and mathematically. In fact, at no point is W(t) differentiable with respect to t. Louis Bachelier noticed around 1900 that sample paths of Brownian motion looked somewhat like stock prices, and he was the first to use formal ideas from mathematics to model financial markets. Stochastic Differential Equations and Financial Models Recall from equation (22.4) that the solution to the initial value problem (1) is equivalent to evaluating the integral T u(t) = u(0) + λu(t)dt. (2) 0 4
5 Define δt = T/N, for T > 0 and positive integer N, and let t j = jδt. The integral term in (2) can be defined by taking the limit as δt 0 of the Riemann sum N 1 j=0 λu(t j )(t j+1 t j ). Now consider a sum that includes a stochastic component involving the Wiener process W: N 1 j=0 σu(t j )(W(t j+1 ) W(t j )). Similar to the definition of the Riemann integral, the stochastic Itô integral may be formally defined by: T N 1 σu(t)dw(t) = lim σu(t j )(W(t j+1 ) W(t j )). (3) δt 0 0 j=0 The Itô integral (3) represents integration with respect to Brownian motion. We can use it to formulate a stochastic analog of the initial value problem (1): u(t) = u(0) + t 0 λu(s)ds + Expressing this relation in differential form yields t 0 σu(s)dw(s), 0 t T. du(t) = λu(t)dt + σu(t)dw(t), u(0) = u 0, W(0) = 0. (4) Equation (4) is a historically important example from finance. This stochastic differential equation (SDE) represents a modification of the initial value problem (1) that includes a diffusion component driven by Brownian motion. The parameter λ is often referred to as the drift coefficient and σ the diffusion coefficient. The celebrated Black Scholes partial differential equation for pricing options can be derived from (4). The analytic solution of (4) can be shown to be the geometric Brownian motion equation ( u(t) = u 0 exp (λ 1 ) 2 σ2 )t + σw(t). (5) Equation (5) can provide a more realistic model of stock price movements than equation (1). We revisit example (1) to illustrate this, using the sample variance of Google s stock price to derive an an estimate of σ: Example 2 Model Google s stock price over 2009 with equation (4). 5
6 goog = yahoo( GOOG, , )(:,9); n = length(goog); sigma = sqrt(var(log(goog))); lambda = log(goog(n)/goog(1)) + 0.5*sigma^2; t = linspace(0,1,n) ; z = randn(n,1); W = (1/sqrt(n)) * cumsum(z); u = goog(1)*exp((lambda-0.5*sigma^2)*t + sigma*w); plot(1:n,goog, -b,1:n,u, -r ) Figure 3: Google s stock price over 2009 (blue) and a model based on equation (4) (red). We see from figure (3) that the new model presents a more realistic-looking one than that of example (1). The new model is stochastic, thus repeated runs of example (2) will yield different 6
7 approximate solutions. Numerical Methods for SDEs Numerical solutions to equation (4) can be computed using a variation of Euler s method called the Euler Maruyama (EM) method. Define δt = T/N, for T > 0 and positive integer N, and let t j = jδt. The EM method proceeds as: u 0 = u(0) u j+1 = u j + λu j δt + σu j W j, where W j = W(t j ) W(t j 1 ), for j = 1, 2,...,N. Each random Brownian increment W j is computed using W j = t j z j, where z j is sampled from a Gaussian distribution of zero mean and unit variance. Strong and Weak Convergence Computed EM solutions of (4) match the true solution more closely as δt decreases. Note that a solution of (4) is a random variable, as is the difference between it and the exact solution. In order to investigate the convergence properties of numerical methods for SDEs, we need notions of convergence that can handle random variables. One such approach is strong convergence. The numeric method converges strongly of order γ if there exists a constant C so that its solution {u j } N 1 j=0 satisfies E u j u(t) Cδt γ, (6) for any fixed t = jδt [0, T]. The notation EX indicates expected value of the random variable X. It was shown by Gikhman and Skorokhod that the strong order convergence of the EM method is γ = 1/2. We can experimentally investigate strong convergence properties of the EM method by considering the average error at a single fixed point (for example, at the endpoint) for many repeated runs of the EM algorithm. A variation of the standard Runge-Kutta method for ODEs can be applied to problem (4) that has strong order convergence γ = 1. The method proceeds as: u 0 = u(0) u j+1 = u j + λu j δt + σu j W j + 1 ) (σ(u j + σu j δt) σuj ( Wj 2 δt). 2 7
8 Many applications do not require accuracy of specific solutions, but instead focus on solution statistics. The notion of weak convergence order is a useful one for such applications. A method converges weakly with order γ if there exists a constant C such that Ef(u j ) Ef(u(t)) Cδt γ, for all polynomial functions f and for any fixed t = jδt [0, T]. We consider the case in which f is the identity function. The EM method has weak order convergence γ = 1. Monte-Carlo Simulation and Option Pricing Assume that a stock price u(t) evolves according to (4). Consider the European call option with value at expiration time T defined by max{u(t) K, 0}, where K the strike price. Assuming no arbitrage 1 and no short-selling, it can be shown that the expected present value of the option is given by exp( λt)e(max{u(t) K, 0}). We can estimate this quantity by simply averaging max{u(t) K, 0} for many solutions u(t) computed by repeated runs of the EM method. The exact solution of the value European call option, assuming that the present time t = 0, can be shown to be given by the solution of the Black Scholes equation: where, C(u(0)) = u(0)n(d 1 ) K exp( λt)n(d 2 ), (7) d 1 = log(u(t)/k) + (λ + 0.5σ2 )T) σ, T d 2 = d 1 σ T, N(x) = 1 x exp( z 2 /2)dz. 2π inf In light of the exact Black Scholes solution, the so-called Monte-Carlo approach of computing many EM solutions to estimate the value of an option seems overly computationally expensive. However, although an exact solution is known for European-style options contracts, closed-form solutions like (7) are generally not available in options pricing. Such cases are very common and arise, for example, in pricing American options. Methods based on the Monte-Carlo approach similar to the above simple example are widely used to price options in these cases. 1 An arbitrage opportunity exists when there are two or more distinct prices for the same financial instrument available simultaneously. 8
9 Exercises 1. Code the Euler Maruyama method and graphically compare its solution to the analytic solution given by geometric Brownian motion for a stock price series. 2. Using the same stock series chosen for the previous exercise, approximate the strong convergence order of the EM method by computing many solutions and averaging the error against the true solution at the endpoint for several values of δt. Your experiments should compute a value of γ of approximately 1/2. 3. Code the example Runge-Kutta method for SDEs and similarly to the last exercise, investigate its strong order convergence properties. 4. Use the definition of weak order convergence and experimentally estimate the weak order of the EM method as in Exercise 2. Compute many runs of the EM method and compute the solution mean, comparing with the true solution mean for each δt. 5. Let a hypothetical option contract be defined by parameters u(0) = 10, K = 12, λ = 0.05, σ = 0.5, and T = 0.5. Compare the computed option value using the Monte Carlo approach against the exact option value computed using formula (7) for 10, 100, and 1000 runs of the EM method. References [1] Chuck Gartland and Kazim Khan, Lectures on Computational Finance, unpublished notes, Kent State University, [2] Desmond Higham, An Algorithmic Introduction to Numerical Simulation of Stochastic Differential Equations, SIAM Rev. Volume 43, Issue 3, pp (2001). [3] Timothy Sauer, Numerical solution of stochastic differential equations in finance. To appear, Handbook of Computational Finance, Springer. 9
Numerical schemes for SDEs
Lecture 5 Numerical schemes for SDEs Lecture Notes by Jan Palczewski Computational Finance p. 1 A Stochastic Differential Equation (SDE) is an object of the following type dx t = a(t,x t )dt + b(t,x t
More informationContinuous Processes. Brownian motion Stochastic calculus Ito calculus
Continuous Processes Brownian motion Stochastic calculus Ito calculus Continuous Processes The binomial models are the building block for our realistic models. Three small-scale principles in continuous
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 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 informationModeling via Stochastic Processes in Finance
Modeling via Stochastic Processes in Finance Dimbinirina Ramarimbahoaka Department of Mathematics and Statistics University of Calgary AMAT 621 - Fall 2012 October 15, 2012 Question: What are appropriate
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 informationLecture Note 8 of Bus 41202, Spring 2017: Stochastic Diffusion Equation & Option Pricing
Lecture Note 8 of Bus 41202, Spring 2017: Stochastic Diffusion Equation & Option Pricing We shall go over this note quickly due to time constraints. Key concept: Ito s lemma Stock Options: A contract giving
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 informationMonte Carlo Methods for Uncertainty Quantification
Monte Carlo Methods for Uncertainty Quantification Mike Giles Mathematical Institute, University of Oxford Contemporary Numerical Techniques Mike Giles (Oxford) Monte Carlo methods 2 1 / 24 Lecture outline
More informationBROWNIAN MOTION Antonella Basso, Martina Nardon
BROWNIAN MOTION Antonella Basso, Martina Nardon basso@unive.it, mnardon@unive.it Department of Applied Mathematics University Ca Foscari Venice Brownian motion p. 1 Brownian motion Brownian motion plays
More informationModule 4: Monte Carlo path simulation
Module 4: Monte Carlo path simulation Prof. Mike Giles mike.giles@maths.ox.ac.uk Oxford University Mathematical Institute Module 4: Monte Carlo p. 1 SDE Path Simulation In Module 2, looked at the case
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 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 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 informationLecture 8: The Black-Scholes theory
Lecture 8: The Black-Scholes theory Dr. Roman V Belavkin MSO4112 Contents 1 Geometric Brownian motion 1 2 The Black-Scholes pricing 2 3 The Black-Scholes equation 3 References 5 1 Geometric Brownian motion
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 informationRandomness and Fractals
Randomness and Fractals Why do so many physicists become traders? Gregory F. Lawler Department of Mathematics Department of Statistics University of Chicago September 25, 2011 1 / 24 Mathematics and the
More informationNumerical Simulation of Stochastic Differential Equations: Lecture 1, Part 2. Integration For deterministic h : R R,
Numerical Simulation of Stochastic Differential Equations: Lecture, Part Des Higham Department of Mathematics University of Strathclyde Lecture, part : SDEs Ito stochastic integrals Ito SDEs Examples of
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 informationThe University of Chicago, Booth School of Business Business 41202, Spring Quarter 2011, Mr. Ruey S. Tsay. Solutions to Final Exam.
The University of Chicago, Booth School of Business Business 41202, Spring Quarter 2011, Mr. Ruey S. Tsay Solutions to Final Exam Problem A: (32 pts) Answer briefly the following questions. 1. Suppose
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 informationAsian Option Pricing: Monte Carlo Control Variate. A discrete arithmetic Asian call option has the payoff. S T i N N + 1
Asian Option Pricing: Monte Carlo Control Variate A discrete arithmetic Asian call option has the payoff ( 1 N N + 1 i=0 S T i N K ) + A discrete geometric Asian call option has the payoff [ N i=0 S T
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 informationContinuous Time Finance. Tomas Björk
Continuous Time Finance Tomas Björk 1 II Stochastic Calculus Tomas Björk 2 Typical Setup Take as given the market price process, S(t), of some underlying asset. S(t) = price, at t, per unit of underlying
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 informationAn Efficient Numerical Scheme for Simulation of Mean-reverting Square-root Diffusions
Journal of Numerical Mathematics and Stochastics,1 (1) : 45-55, 2009 http://www.jnmas.org/jnmas1-5.pdf JNM@S Euclidean Press, LLC Online: ISSN 2151-2302 An Efficient Numerical Scheme for Simulation of
More informationMath 416/516: Stochastic Simulation
Math 416/516: Stochastic Simulation Haijun Li lih@math.wsu.edu Department of Mathematics Washington State University Week 13 Haijun Li Math 416/516: Stochastic Simulation Week 13 1 / 28 Outline 1 Simulation
More informationPractical example of an Economic Scenario Generator
Practical example of an Economic Scenario Generator Martin Schenk Actuarial & Insurance Solutions SAV 7 March 2014 Agenda Introduction Deterministic vs. stochastic approach Mathematical model Application
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 information1 Mathematics in a Pill 1.1 PROBABILITY SPACE AND RANDOM VARIABLES. A probability triple P consists of the following components:
1 Mathematics in a Pill The purpose of this chapter is to give a brief outline of the probability theory underlying the mathematics inside the book, and to introduce necessary notation and conventions
More informationThe Black-Scholes Model
The Black-Scholes Model Liuren Wu Options Markets Liuren Wu ( c ) The Black-Merton-Scholes Model colorhmoptions Markets 1 / 18 The Black-Merton-Scholes-Merton (BMS) model Black and Scholes (1973) and Merton
More informationComputer labs. May 10, A list of matlab tutorials can be found under
Computer labs May 10, 2018 A list of matlab tutorials can be found under http://snovit.math.umu.se/personal/cohen_david/teachlinks.html Task 1: The following MATLAB code generates (pseudo) uniform random
More informationThe Black-Scholes Model
The Black-Scholes Model Liuren Wu Options Markets (Hull chapter: 12, 13, 14) Liuren Wu ( c ) The Black-Scholes Model colorhmoptions Markets 1 / 17 The Black-Scholes-Merton (BSM) model Black and Scholes
More informationAn Introduction to Stochastic Calculus
An Introduction to Stochastic Calculus Haijun Li lih@math.wsu.edu Department of Mathematics Washington State University Week 2-3 Haijun Li An Introduction to Stochastic Calculus Week 2-3 1 / 24 Outline
More informationSimulating Stochastic Differential Equations
IEOR E4603: Monte-Carlo Simulation c 2017 by Martin Haugh Columbia University Simulating Stochastic Differential Equations In these lecture notes we discuss the simulation of stochastic differential equations
More informationApplications of Stochastic Processes in Asset Price Modeling
Applications of Stochastic Processes in Asset Price Modeling TJHSST Computer Systems Lab Senior Research Project 2008-2009 Preetam D Souza May 26, 2009 Abstract Stock market forecasting and asset price
More informationELEMENTS OF MONTE CARLO SIMULATION
APPENDIX B ELEMENTS OF MONTE CARLO SIMULATION B. GENERAL CONCEPT The basic idea of Monte Carlo simulation is to create a series of experimental samples using a random number sequence. According to the
More informationContinous time models and realized variance: Simulations
Continous time models and realized variance: Simulations Asger Lunde Professor Department of Economics and Business Aarhus University September 26, 2016 Continuous-time Stochastic Process: SDEs Building
More informationComputational Finance Improving Monte Carlo
Computational Finance Improving Monte Carlo School of Mathematics 2018 Monte Carlo so far... Simple to program and to understand Convergence is slow, extrapolation impossible. Forward looking method ideal
More information3.1 Itô s Lemma for Continuous Stochastic Variables
Lecture 3 Log Normal Distribution 3.1 Itô s Lemma for Continuous Stochastic Variables Mathematical Finance is about pricing (or valuing) financial contracts, and in particular those contracts which depend
More informationMonte Carlo Simulations
Monte Carlo Simulations Lecture 1 December 7, 2014 Outline Monte Carlo Methods Monte Carlo methods simulate the random behavior underlying the financial models Remember: When pricing you must simulate
More informationStochastic Runge Kutta Methods with the Constant Elasticity of Variance (CEV) Diffusion Model for Pricing Option
Int. Journal of Math. Analysis, Vol. 8, 2014, no. 18, 849-856 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ijma.2014.4381 Stochastic Runge Kutta Methods with the Constant Elasticity of Variance
More informationLecture 3. Sergei Fedotov Introduction to Financial Mathematics. Sergei Fedotov (University of Manchester) / 6
Lecture 3 Sergei Fedotov 091 - Introduction to Financial Mathematics Sergei Fedotov (University of Manchester) 091 010 1 / 6 Lecture 3 1 Distribution for lns(t) Solution to Stochastic Differential Equation
More informationSTOCHASTIC CALCULUS AND BLACK-SCHOLES MODEL
STOCHASTIC CALCULUS AND BLACK-SCHOLES MODEL YOUNGGEUN YOO Abstract. Ito s lemma is often used in Ito calculus to find the differentials of a stochastic process that depends on time. This paper will introduce
More informationNumerical Simulation of Stochastic Differential Equations: Lecture 1, Part 1. Overview of Lecture 1, Part 1: Background Mater.
Numerical Simulation of Stochastic Differential Equations: Lecture, Part Des Higham Department of Mathematics University of Strathclyde Course Aim: Give an accessible intro. to SDEs and their numerical
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 informationStochastic Processes and Advanced Mathematical Finance. Stochastic Differential Equations and the Euler-Maruyama Method
Steven R. Dunbar Department of Mathematics 203 Avery Hall University of Nebraska-Lincoln Lincoln, NE 68588-0130 http://www.math.unl.edu Voice: 402-472-3731 Fax: 402-472-8466 Stochastic Processes and Advanced
More informationNumerical Methods for Stochastic Differential Equations with Applications to Finance
Numerical Methods for Stochastic Differential Equations with Applications to Finance Matilde Lopes Rosa Instituto Superior Técnico University of Lisbon, Portugal May 2016 Abstract The pricing of financial
More informationRichardson Extrapolation Techniques for the Pricing of American-style Options
Richardson Extrapolation Techniques for the Pricing of American-style Options June 1, 2005 Abstract Richardson Extrapolation Techniques for the Pricing of American-style Options In this paper we re-examine
More informationFinancial Engineering with FRONT ARENA
Introduction The course A typical lecture Concluding remarks Problems and solutions Dmitrii Silvestrov Anatoliy Malyarenko Department of Mathematics and Physics Mälardalen University December 10, 2004/Front
More informationMFE/3F Questions Answer Key
MFE/3F Questions Download free full solutions from www.actuarialbrew.com, or purchase a hard copy from www.actexmadriver.com, or www.actuarialbookstore.com. Chapter 1 Put-Call Parity and Replication 1.01
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 informationErrata, Mahler Study Aids for Exam 3/M, Spring 2010 HCM, 1/26/13 Page 1
Errata, Mahler Study Aids for Exam 3/M, Spring 2010 HCM, 1/26/13 Page 1 1B, p. 72: (60%)(0.39) + (40%)(0.75) = 0.534. 1D, page 131, solution to the first Exercise: 2.5 2.5 λ(t) dt = 3t 2 dt 2 2 = t 3 ]
More informationLévy models in finance
Lévy models in finance Ernesto Mordecki Universidad de la República, Montevideo, Uruguay PASI - Guanajuato - June 2010 Summary General aim: describe jummp modelling in finace through some relevant issues.
More informationThe stochastic calculus
Gdansk A schedule of the lecture Stochastic differential equations Ito calculus, Ito process Ornstein - Uhlenbeck (OU) process Heston model Stopping time for OU process Stochastic differential equations
More informationEnergy Price Processes
Energy Processes Used for Derivatives Pricing & Risk Management In this first of three articles, we will describe the most commonly used process, Geometric Brownian Motion, and in the second and third
More informationAssignment - Exotic options
Computational Finance, Fall 2014 1 (6) Institutionen för informationsteknologi Besöksadress: MIC, Polacksbacken Lägerhyddvägen 2 Postadress: Box 337 751 05 Uppsala Telefon: 018 471 0000 (växel) Telefax:
More informationDr. Maddah ENMG 625 Financial Eng g II 10/16/06
Dr. Maddah ENMG 65 Financial Eng g II 10/16/06 Chapter 11 Models of Asset Dynamics () Random Walk A random process, z, is an additive process defined over times t 0, t 1,, t k, t k+1,, such that z( t )
More informationMonte Carlo Simulation of Stochastic Processes
Monte Carlo Simulation of Stochastic Processes Last update: January 10th, 2004. In this section is presented the steps to perform the simulation of the main stochastic processes used in real options applications,
More informationDeriving and Solving the Black-Scholes Equation
Introduction Deriving and Solving the Black-Scholes Equation Shane Moore April 27, 2014 The Black-Scholes equation, named after Fischer Black and Myron Scholes, is a partial differential equation, which
More informationSYSM 6304: Risk and Decision Analysis Lecture 6: Pricing and Hedging Financial Derivatives
SYSM 6304: Risk and Decision Analysis Lecture 6: Pricing and Hedging Financial Derivatives M. Vidyasagar Cecil & Ida Green Chair The University of Texas at Dallas Email: M.Vidyasagar@utdallas.edu October
More informationOption Pricing under Delay Geometric Brownian Motion with Regime Switching
Science Journal of Applied Mathematics and Statistics 2016; 4(6): 263-268 http://www.sciencepublishinggroup.com/j/sjams doi: 10.11648/j.sjams.20160406.13 ISSN: 2376-9491 (Print); ISSN: 2376-9513 (Online)
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 informationSlides for DN2281, KTH 1
Slides for DN2281, KTH 1 January 28, 2014 1 Based on the lecture notes Stochastic and Partial Differential Equations with Adapted Numerics, by J. Carlsson, K.-S. Moon, A. Szepessy, R. Tempone, G. Zouraris.
More informationStochastic Differential equations as applied to pricing of options
Stochastic Differential equations as applied to pricing of options By Yasin LUT Supevisor:Prof. Tuomo Kauranne December 2010 Introduction Pricing an European call option Conclusion INTRODUCTION A stochastic
More informationAmerican Option Pricing Formula for Uncertain Financial Market
American Option Pricing Formula for Uncertain Financial Market Xiaowei Chen Uncertainty Theory Laboratory, Department of Mathematical Sciences Tsinghua University, Beijing 184, China chenxw7@mailstsinghuaeducn
More informationSTEX s valuation analysis, version 0.0
SMART TOKEN EXCHANGE STEX s valuation analysis, version. Paulo Finardi, Olivia Saa, Serguei Popov November, 7 ABSTRACT In this paper we evaluate an investment consisting of paying an given amount (the
More informationMath 239 Homework 1 solutions
Math 239 Homework 1 solutions Question 1. Delta hedging simulation. (a) Means, standard deviations and histograms are found using HW1Q1a.m with 100,000 paths. In the case of weekly rebalancing: mean =
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 informationAssicurazioni Generali: An Option Pricing Case with NAGARCH
Assicurazioni Generali: An Option Pricing Case with NAGARCH Assicurazioni Generali: Business Snapshot Find our latest analyses and trade ideas on bsic.it Assicurazioni Generali SpA is an Italy-based insurance
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 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 informationProbability in Options Pricing
Probability in Options Pricing Mark Cohen and Luke Skon Kenyon College cohenmj@kenyon.edu December 14, 2012 Mark Cohen and Luke Skon (Kenyon college) Probability Presentation December 14, 2012 1 / 16 What
More informationOption Pricing Formula for Fuzzy Financial Market
Journal of Uncertain Systems Vol.2, No., pp.7-2, 28 Online at: www.jus.org.uk Option Pricing Formula for Fuzzy Financial Market Zhongfeng Qin, Xiang Li Department of Mathematical Sciences Tsinghua University,
More informationJDEP 384H: Numerical Methods in Business
Chapter 4: Numerical Integration: Deterministic and Monte Carlo Methods Chapter 8: Option Pricing by Monte Carlo Methods JDEP 384H: Numerical Methods in Business Instructor: Thomas Shores Department of
More informationModule 10:Application of stochastic processes in areas like finance Lecture 36:Black-Scholes Model. Stochastic Differential Equation.
Stochastic Differential Equation Consider. Moreover partition the interval into and define, where. Now by Rieman Integral we know that, where. Moreover. Using the fundamentals mentioned above we can easily
More informationStochastic Dynamical Systems and SDE s. An Informal Introduction
Stochastic Dynamical Systems and SDE s An Informal Introduction Olav Kallenberg Graduate Student Seminar, April 18, 2012 1 / 33 2 / 33 Simple recursion: Deterministic system, discrete time x n+1 = f (x
More informationComputational Finance
Path Dependent Options Computational Finance School of Mathematics 2018 The Random Walk One of the main assumption of the Black-Scholes framework is that the underlying stock price follows a random walk
More informationMarket Volatility and Risk Proxies
Market Volatility and Risk Proxies... an introduction to the concepts 019 Gary R. Evans. This slide set by Gary R. Evans is licensed under a Creative Commons Attribution-NonCommercial-ShareAlike 4.0 International
More informationDrunken Birds, Brownian Motion, and Other Random Fun
Drunken Birds, Brownian Motion, and Other Random Fun Michael Perlmutter Department of Mathematics Purdue University 1 M. Perlmutter(Purdue) Brownian Motion and Martingales Outline Review of Basic Probability
More information"Pricing Exotic Options using Strong Convergence Properties
Fourth Oxford / Princeton Workshop on Financial Mathematics "Pricing Exotic Options using Strong Convergence Properties Klaus E. Schmitz Abe schmitz@maths.ox.ac.uk www.maths.ox.ac.uk/~schmitz Prof. Mike
More informationNumerical Simulation of Stochastic Differential Equations: Lecture 2, Part 2
Numerical Simulation of Stochastic Differential Equations: Lecture 2, Part 2 Des Higham Department of Mathematics University of Strathclyde Montreal, Feb. 2006 p.1/17 Lecture 2, Part 2: Mean Exit Times
More informationFinancial Risk Modeling on Low-power Accelerators: Experimental Performance Evaluation of TK1 with FPGA
Financial Risk Modeling on Low-power Accelerators: Experimental Performance Evaluation of TK1 with FPGA Rajesh Bordawekar and Daniel Beece IBM T. J. Watson Research Center 3/17/2015 2014 IBM Corporation
More informationMAS3904/MAS8904 Stochastic Financial Modelling
MAS3904/MAS8904 Stochastic Financial Modelling Dr Andrew (Andy) Golightly a.golightly@ncl.ac.uk Semester 1, 2018/19 Administrative Arrangements Lectures on Tuesdays at 14:00 (PERCY G13) and Thursdays at
More informationGeometric Brownian Motions
Chapter 6 Geometric Brownian Motions 1 Normal Distributions We begin by recalling the normal distribution briefly. Let Z be a random variable distributed as standard normal, i.e., Z N(0, 1). The probability
More informationA new approach for scenario generation in risk management
A new approach for scenario generation in risk management Josef Teichmann TU Wien Vienna, March 2009 Scenario generators Scenarios of risk factors are needed for the daily risk analysis (1D and 10D ahead)
More informationEconomathematics. Problem Sheet 1. Zbigniew Palmowski. Ws 2 dw s = 1 t
Economathematics Problem Sheet 1 Zbigniew Palmowski 1. Calculate Ee X where X is a gaussian random variable with mean µ and volatility σ >.. Verify that where W is a Wiener process. Ws dw s = 1 3 W t 3
More informationComputer Exercise 2 Simulation
Lund University with Lund Institute of Technology Valuation of Derivative Assets Centre for Mathematical Sciences, Mathematical Statistics Spring 2010 Computer Exercise 2 Simulation This lab deals with
More informationDepartment of Mathematics. Mathematics of Financial Derivatives
Department of Mathematics MA408 Mathematics of Financial Derivatives Thursday 15th January, 2009 2pm 4pm Duration: 2 hours Attempt THREE questions MA408 Page 1 of 5 1. (a) Suppose 0 < E 1 < E 3 and E 2
More informationCE 513: Statistical Methods in CE
CE 513: Statistical Methods in CE 2017 (Aug-Nov) Budhaditya Hazra Room N-307 Department of Civil Engineering 1 CE 513: Statistical Methods in CE ITO CALCULUS 2 General Ito Integrals If X (t) is a process
More informationMarket Risk: FROM VALUE AT RISK TO STRESS TESTING. Agenda. Agenda (Cont.) Traditional Measures of Market Risk
Market Risk: FROM VALUE AT RISK TO STRESS TESTING Agenda The Notional Amount Approach Price Sensitivity Measure for Derivatives Weakness of the Greek Measure Define Value at Risk 1 Day to VaR to 10 Day
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 informationAnalysing multi-level Monte Carlo for options with non-globally Lipschitz payoff
Finance Stoch 2009 13: 403 413 DOI 10.1007/s00780-009-0092-1 Analysing multi-level Monte Carlo for options with non-globally Lipschitz payoff Michael B. Giles Desmond J. Higham Xuerong Mao Received: 1
More informationBluff Your Way Through Black-Scholes
Bluff our Way Through Black-Scholes Saurav Sen December 000 Contents What is Black-Scholes?.............................. 1 The Classical Black-Scholes Model....................... 1 Some Useful Background
More informationIntroduction to Stochastic Calculus With Applications
Introduction to Stochastic Calculus With Applications Fima C Klebaner University of Melbourne \ Imperial College Press Contents Preliminaries From Calculus 1 1.1 Continuous and Differentiable Functions.
More information[AN INTRODUCTION TO THE BLACK-SCHOLES PDE MODEL]
2013 University of New Mexico Scott Guernsey [AN INTRODUCTION TO THE BLACK-SCHOLES PDE MODEL] This paper will serve as background and proposal for an upcoming thesis paper on nonlinear Black- Scholes PDE
More informationMonte Carlo Methods in Financial Engineering
Paul Glassennan Monte Carlo Methods in Financial Engineering With 99 Figures
More informationComputer Exercise 2 Simulation
Lund University with Lund Institute of Technology Valuation of Derivative Assets Centre for Mathematical Sciences, Mathematical Statistics Fall 2017 Computer Exercise 2 Simulation This lab deals with pricing
More informationby Kian Guan Lim Professor of Finance Head, Quantitative Finance Unit Singapore Management University
by Kian Guan Lim Professor of Finance Head, Quantitative Finance Unit Singapore Management University Presentation at Hitotsubashi University, August 8, 2009 There are 14 compulsory semester courses out
More informationStochastic Calculus - An Introduction
Stochastic Calculus - An Introduction M. Kazim Khan Kent State University. UET, Taxila August 15-16, 17 Outline 1 From R.W. to B.M. B.M. 3 Stochastic Integration 4 Ito s Formula 5 Recap Random Walk Consider
More information