Questions tagged [finite-difference-method]

Filter by
Sorted by
Tagged with
8
votes
2answers
135 views

Improve Finite Difference Scheme

I understand how to derive and implement standard finite difference schemes. I wonder how to improve such a standard FD scheme? For example, when solving the standard Black-Scholes equation, the ...
5
votes
1answer
306 views

Other numerraire choices when applying Feynman Kac

all of the books and notes I have seen on the Feynman Kac formula mostly applied to Risk neutral measure, i.e. different interest rate models, stochastic volatility, etc. I think risk neutral measure ...
5
votes
0answers
115 views

Benchmark a Libor Market Model implementation

Assume I have implemented a solution of the Libor Market model PDE in terms of the Finite Difference method. What is a good strategy for validating and benchmarking the results of this implementation? ...
5
votes
0answers
199 views

Finite Difference with SVI Vol Model

I am attempting to implement a local vol pricing model in finite difference for equity index options. I have followed Gatheral's Lectures and fitted an SVI Model bringing me to the following local ...
4
votes
2answers
424 views

Smoothing of the payoff function as a terminal condition for numerical option pricing

I am interested in using a 4th order finite difference method in (underlying asset) space to price a European call basket option. I have developed the solver and everything works as expected, except ...
4
votes
1answer
290 views

Finite Difference method in Matlab for SABR volatility model fails to provide correct option values

Currently, I'm trying to implement a Finite Difference (FD) method in Matlab for my thesis (Quantitative Finance). I implemented the FD method for Black-Scholes already and got correct results. ...
4
votes
0answers
378 views

binomial trees and finite differences

I was reading Tavella Randall book and their explanation why binomial trees are a particular example of finite differences. I started having additional questions. So, they way they do that is saying ...
3
votes
1answer
203 views

Errors on Finite Differences + Implicit Scheme + Black & Scholes

I'm solving the classical Black & Scholes (BS) PDE for a European option using finite difference and the implicit scheme. In other words, I'm trying to solve $\displaystyle\frac{\partial V}{\...
3
votes
2answers
183 views

Is it possible to model path-dependent clauses using finite difference methods?

I'm trying to build a convertible bond pricer. In my case a convertible bond is a complex derivative with call, put and conversion price reset clauses, and all of the clauses are triggered in a path-...
3
votes
1answer
242 views

Binomial Trees vs FDM

Binomial trees as the number of time steps is increased (or equivalently as the time step tends to 0), converge to the exact value for an option. So why do people use FDM for pricing options (for ...
3
votes
1answer
1k views

Local Volatility implementation

The Dupire equation is well-known and mentioned in thousands of articles. Although I could not find a lot of documentation about a consistent and proper way of implementing the formula (The difficulty ...
3
votes
2answers
356 views

von Neumann boundary in the transformed PDE

I have transformed the BSM PDE $$\frac{\partial V}{\partial t} + \frac{\sigma^2}{2}S^2 \frac{\partial^2 V}{\partial S^2} + rS \frac{\partial V}{\partial S} - rV = 0 $$ to $u(\tau,x) = V(T-\tau,S_{0}...
3
votes
0answers
72 views

Operator splitting method on three assets black scholes equation

Currently I am studying finite difference method on derivatives with three (or more) underlyings and little bit confused on operator splitting method because two papers have different result. For the ...
2
votes
1answer
593 views

Brennan-Schwartz algorithm for pricing American options

I'm reading Pricing American Options using LU decomposition by Ikonen and Toivanen (IT). They reference The valuation of American put options by Brennan and Schwartz, and cast it as method that uses ...
2
votes
2answers
5k views

calculate gamma value using finite difference method

I try to use the finite difference method to get the approximately gamma value, but there is an issue I can't solve. First, I set $h$ to 1 basis point of underlying asset value, but the result is not ...
2
votes
1answer
102 views

Greeks: Estimate gamma by Monte Carlo finite difference

When I was using Monte Carlo to calculate the gamma of a vanilla call option by finite difference method, I stuck in this weird situation as below. Consider this, $$ Gamma = \frac{CallPrice(S^{up}_{T})...
2
votes
2answers
179 views

Implicit finite difference method always guarantees positive and stable price of derivative?

For the following black scholes pde $$ f_t + rSf_S+\frac{1}{2}\sigma^2S^2f_{SS} = rf $$ By denoting $f_{i}^{n} = $ Price of derivative at price node $i$ and time node $n$ and assume uniform grid, the ...
2
votes
1answer
228 views

SABR PDE spot/forward upper boundary condition implementation

When running my Finite Difference code, I observe something odd. Although implementing a classical (non-reverting) SABR model, I initialized the variables such that it should be equal to Black-...
2
votes
1answer
142 views

Unable to obtain correct Finite Difference Results

A rather general question regarding a specific problem I am facing with my Matlab implementation of the implicit FD method for this PDE: \begin{equation} \frac{\sigma_s^2}{2}\frac{\partial^2 V}{\...
2
votes
1answer
69 views

Numerical Solution to 3 Dimensional Backward BS PDE

I have a three dimensional backward BS PDE. $$ \frac{\partial V}{\partial t} + a(t) S \frac{\partial V}{\partial S} + \frac{1}{2} \sigma(t, S)^2 \frac{\partial^2 V}{\partial S^2} + b(t, M) \frac{\...
2
votes
2answers
2k views

Local Volatility calculation in Python

I am trying to price Local Volatility in Python using Dupire (Finite Difference Method). I have following set of information Spot: 770.05, Strike: 850, Type: 'C', rfr: 0.0066, time to maturity = ...
2
votes
1answer
102 views

Finite difference methods for (continuously) strike-resettable American options

For simplicity, let us consider an American call/put with a continuously resettable strike price. Current time is $t=0$, maturity is at $t=T$, and the initial strike is $K_0$. We consider a "...
2
votes
3answers
123 views

For using finite difference on PDE, what should the grid be?

If I wish to use finite difference methods to approximate the pricing function $F(t, s)$ for an option (say, a call), what size grid should I use? I mean, it seems to make sense to start the grid at ...
2
votes
1answer
304 views

pricing american put option with fdm

Assume I use some finite difference solver to solve for American type of exercise in BS framework where stock pays dividend discretely. Then at every time iteration, for call option, I firstly adjust ...
2
votes
1answer
113 views

Solving a Non-Linear PDE using a Finite Difference Scheme

I have the following non-linear PDE and I have no idea how to go about solving it using a finite difference scheme in Python. Can someone get me started and/or point me to an algorithm for doing this? ...
2
votes
1answer
64 views

What different techniques exist for modeling exotics near payoff discontinuities in Finite Difference method?

If you are modeling an exotic, like a binary or a barrier, and hedging it with vanillas that have strikes quite close to the exotic's strike, then a large asset step size, for example, $\delta S = \...
2
votes
0answers
42 views

Stability of Finite Difference method for Breeden-Litzenberger

I am trying to derive a risk-neutral density from European call option prices using a second order finite difference scheme. Let $C(K,T)$ be the price of a European call with strike $K$ and expiry $T$ ...
2
votes
0answers
81 views

American-Bermudan-Asian option fixed strike using finite differences

I'm trying to price the same American-Bermudan-Asian option described in Longstaff Schwartz (2001). Specifically, using finite difference methods with an explicit scheme to solve $\begin{aligned} \...
2
votes
0answers
227 views

How do you numerically solve the Dupire Local Volatility PDE in log moneyness-time space?

I am trying to implement a numerical solution to price vanilla calls. I am using the Dupire equation in log moneyness-time (k = ln(F/T)) space as per below PDE I have tried solving it using a fully ...
2
votes
0answers
55 views

How does the diameter of the spatial grid affects the solution of a Crank Nicholson algorithm?

this is my first question so I hope I express myself clearly. I'm trying to implement an Implicit and a Crank Nicolson algorithm for the generic PDE $\partial_\tau u(\tau,x)+a \partial_x^2 u(\tau,x) +...
2
votes
0answers
86 views

Multivariate interpolation for estimating FDM in-between grid points

After implementing some FDM to price some option, there are gaps between our grid points that may be of interest. From reading around, it appears common to use bilinear interpolation to estimate ...
1
vote
1answer
84 views

Finite difference: move forwards or backwards?

In finite differences for the black scholes method, you move backwards in time, since of course you know the prices at time $t = T$, and then you iterate until you get to time $t = 0$. However, why ...
1
vote
1answer
64 views

Explicit Euler stability for the Heat Equation (FDM)

Why the Explicit Euler scheme for the Heat Equation is stable only if $k \leq h^2/2$ ? Here is the difference equation: \begin{equation} \frac{U_j^{n+1}-U_{j}^n}{k} = \frac{1}{h^2}(U_{j+1}^n-2U_j^n+...
1
vote
1answer
130 views

Why can't we use Finite Differences with non-parabolic PDEs?

The title of the question says it all. Why can we only apply the method to parabolic PDEs like the heat equation, and not to ordinary PDEs?
1
vote
1answer
45 views

Canonical text on numerical PDEs in finance

I am looking for a text similar to Glasserman's Monte Carlo Methods in Financial Engineering, but with a focus on numerical methods for PDEs. Glasserman's book seems to cover a lot for what is ...
1
vote
1answer
221 views

Pricing an American derivative with finite differences

I have a basic fundamental question on pricing an American option in the Black-Scholes (BS) framework: I seem to confuse two different approaches to price any early exercise, Write down a linear ...
1
vote
1answer
44 views

How many decimals of accuracy can I expect from FDM and MC (both valuation and risk)

I have implemented some Monte Carlo and FDM code. I can then get greeks by bumping. I am comparing to to exact formulas of price + greeks, and am wondering how many decimals of accuracy I can expect ...
1
vote
1answer
42 views

How are FDE's implemented when one wants one particular price?

Say I want to price a particular call option in the Black Scholes model using finite difference methods. The value process of this option $V(s, t)$ satisfies a PDE. I can use finite difference ...
1
vote
1answer
505 views

Black Scholes Theta Finite difference

I am trying to obtain the Theta from Closed Formula by using Finite Difference methods and I observe some discrepancies. For instance, here with the following parameters: Spot:50, Strike:50, Rate: 0....
1
vote
1answer
154 views

MATLAB exercise on an European call option with time-varying volatility

I have to solve the following exercise: compute and plot the value $V = V(S, t),\ t<T$, ($T=$ maturity) of an European CALL option (with arbitrary $t$, $T$, $K$ (strike price), $r$ (risk-free ...
1
vote
2answers
551 views

Full value function of an American option with QuantLib FD

I am looking at the Equity Option example of QuantLib: http://quantlib.org/reference/_equity_option_8cpp-example.html and more particularly the FDAmericanEngine. However, I am not interested in the ...
1
vote
0answers
249 views

Pricing Knock Out Barrier Options by solving Black Scholes PDE (MATLAB)

This question is based on MATLAB functions. Suppose there is a stock S following the process $dS_t=(r-q)S_tdt+\sigma(S_t,t)dW_t$ r - risk-free rate, q - dividend yield, W - Weiner process The ...
1
vote
0answers
63 views

Optimal allocation problem by finite differences

I am attempting to apply implicit finite difference to solve Merton's problem of optimal portfolio allocation for constant parameters. The equation to solve is the Hamilton-Jacobi-Bellman equation: $$...
1
vote
0answers
42 views

Oscillating errors in finite difference Black Scholes

I am writing an implementation of the explicit finite difference method to price a standard european call option, and comparing the results to the corresponding analytical value to gauge the error ...
1
vote
1answer
143 views

Finite Difference implicit scheme

I'm trying to solve the following PDE numerically using an implicit FD scheme: \begin{equation} \frac{\sigma_s^2}{2}\frac{\partial^2 V}{\partial S^2} + \rho \sigma_S \sigma_\alpha\frac{\partial^2 V}{\...
1
vote
0answers
302 views

Trinomial Tree and finite difference methods

I want to understand the connection between the trinomial tree and the finite difference methods. As far as I understood so far is, if we transform the Black-Scholes-PDE to heat equation, the explicit ...
1
vote
0answers
169 views

methodology confirmation for computing implied risk-neutral CDF from option prices

In this question, the risk-neutral probability distribution $q(S_T=s)$ for the underlying at time $t = T$ is given by the Breeden-Litzenberger identity as: $$ \frac{1}{P(0,T)} \frac{ \partial^2 C }{\...
1
vote
0answers
44 views

Methods Available for Derivative Pricing in Mathematica? [closed]

I am using Mathematica to price options (built in functions, no need to reinvent the wheel, right?). In the documentation, the Binomial method is used as an example of specifying a non-standard method....
0
votes
1answer
1k views

how to price barrier option under local vol model using QuantLib

I use QuantLib in Python. Now I have implied volatility surface data. How can I get the local vol surface than using finite difference method to price a barrier option in QuantLib?
0
votes
2answers
1k views

Which options do not have a closed form pricing formula like BS?

Q1 Are there options that doesn't have a closed form pricing formula like BS? Is this the situation when we have to use Finite Difference Method? Can someone give an example? (I hope this option ...