# Questions tagged [numerical-methods]

The tag has no usage guidance.

118 questions
Filter by
Sorted by
Tagged with
242 views

### Anyone has detailed explanation on how to use epstein-zin preferences in asset pricing models

I'd be interested to know how Epstein-Zin preferences are used in, say, consumption-based asset pricing models. I'm looking for specific derivations (how you get the SDF) and possible numerical ...
28 views

### Quick Discretization question for finite difference and finite element methods

Assume we have the discretization in space $x_1, x_2, ... , x_M$ and time $t_1, t_2, ... , t_N$ for a finite difference or finite element method for option pricing and we want to solve for the option ...
21 views

### boundary conditions in finite element method

In the appendix A of this paper, https://citeseerx.ist.psu.edu/viewdoc/download?doi=10.1.1.227.5073&rep=rep1&type=pdf, a finite element method is demonstrated to price a straddle. The same ...
69 views

54 views

### Good ways to select best decision among N decisions, each with a profit/loss distribution? [closed]

I'm working on a problem where an asset owner (e.g., owner of a factory, power plant, etc.) can take a number of possible decisions (say 10). Each of those 10 decisions entails certain actions, but ...
123 views

### CIR model. Is there a closed-form solution or even a good proxy of analytical solution?

Is there a closed-form (analytical) solution for the Cox-Ingersoll-Ross SDE \begin{equation} dr_t=k_r(\theta_r-r_t)dt+\sigma_r\sqrt{r_t}dW_t\tag{1} \end{equation} ? Notice that $\{r_t\}$ is our ...
73 views

### How to compute returns from cumulative returns in Python? [closed]

If X is a $T\times N$ pandas DataFrame of multivariate asset returns, the cumulative returns can be computed in python as (1 + X).cumprod() - 1 How can I reverse this operation so that I go ...
76 views

### How to simulate asset prices/returns that display market regimes?

Are there any techniques that can make a multivariate random number generating process for stock prices/returns, like geometric Brownian motion via Cholesky, also include the simulation of a finite ...
14 views

### Is it easier to estimate the mean of volume/dollar bar returns than time bars?

Financial data, by default, is sampled using constant intervals of time, i.e. we collect the prices of assets every day, and then take the log-differences of prices to compute daily returns. These are ...
151 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 ...
155 views

104 views

### L2 Assumptions of the Longstaff Schwartz method

In page 121 of the original LS Paper they use the fact that the space of functions they are dealing with (payoffs of American options), belong to the $\mathcal L^2$ space. They use this assumption ...
53 views

### Optimized search for yield-to-worst of a callable bond

Suppose that I need to find the yield-to-worst of a callable bond, and that the option is American (call any time). The bond may have step-up coupons and/or non-constant call price (oprion strike). ...
206 views

### QuantLib returns slightly different bondYield when backtested

I am just starting to get familiar with QuantLib (in particular, fixed rate bond pricing functions). I read a number of examples, from which I am able to calculate bond price and bond yield. The ...
29 views

### How to approximate expectation and variance of an integral from a discrete Time series financial dataset?

I have discrete time series financial data, with time($u$), price($S$) and someVariable($q$) which looks something like this. ...
178 views

### In Carr-Madans option pricing method, why do they use FFT?

In the famous fourier option pricing method by Carr-Madan, (http://citeseerx.ist.psu.edu/viewdoc/download?doi=10.1.1.348.4044&rep=rep1&type=pdf), the crucial formula is They evaluate this by ...
31 views

### How important is it that numerical methods can price for various strikes simultaneously?

I am reading a paper which presents a numerical method to price call options. Call this Method 1. The method can also price several call-options for a range of strikes simultaneously if you want it to,...
1k views

### Simulation of Geometric Brownian Motion in R

Using R, I would like to simulate a sample path of a geometric Brownian motion using \begin{equation*} S(t) = S(0) \exp\left(\left(\mu - \frac{\sigma^{2}}{2}\right)t + \sigma B_{t}\right), \end{...
129 views

### Produce the random variable for an asset from a uniformly distributed random varible

I'm working on a quant interview question from the book called Quant Job Interview Questions And Answers (by Mark Joshi and other authors). I cannot understand the following question(not the answer, ...
61 views

### Calibrate a model parameter with an error function

Suppose I want to find the implied volatility using an option model from market prices. Surely I can find the implied volatility for each strike price ($k$ different strike prices) for a given ...
82 views

### When is a numerical solution the only way to obtain a solution to BS?

I am only now reading into Mathematical Finance, I understand the derivation of the BS equation with vanilla European options. On the next page of my book it starts to delve into obtaining exact ...
105 views

### Fastest way to calculate YTM from bond price

I would like to calculate YTM for every top of the book update on the 10-year note traded on Brokertec. There is no closed form solution so have to use a root finding method like Newton-Rhapson. It ...
82 views

### Numerical Solutions to PDEs with Financial Applications

I am reading a paper by Richard White, Opengamma named Numerical Solutions to PDEs with Financial Applications. There is an implementation codes as stated in paper hosted at https://opengamma.com/...
59 views

### Sensitivity Approximation - Crank Nicolson

I am looking into a new method of calculating sensitivities starting off with a proof of concept with Black Scholes PDE. Suppose I want to calculate Rho and take the derivative of the PDE (heresy!!) ...
87 views

### What are the industry standards and rules of thumb when it comes to numerical methods?

So, as far as I know, we have 3 main numerical methods. Monte Carlo, PDE-methods (FDM), and numerical integration methods (Fourier transforms and so on). How do these methods generally compare to ...
73 views

### How good is a “good accuracy” in pricing?

Say you want to test various numerical algorithms for purposes of pricing. How close do you need to be to some benchmark value (the "actual" price) for your accuracy to be good? Say I am trying to ...
169 views

### Simulation scheme for SABR beside the standard Euler discretization

QUESTION: Beside Euler Scheme, is there another more robust (and preferably easy to implement) way to simulate asset path with SABR dynamics? Simulation that will withstand even for high volatilities....
503 views

### Least Squares Monte Carlo

Could you explain to me in words (no formulas) the concept of the Least Squares Monte Carlo method to price an American style option?
271 views

### How to determine the order of convergence of the Euler-Maruyama method?

To make this simple let us consider the Geometric Brownian Motions. My questions: 1. How can I show that the Euler-Maruyama Method is convergent using GBM? 2. How can I determine the order of ...
631 views

### Implied volatility in Monte Carlo models

Suppose I want to get the implied volatility for a given option, whose process does not generate a closed-form formula. In that framework, how is the IV calculated, given the fact that bisection ...
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 "...
251 views

### Architecture of a global pricing library with immutable payoffs

By global pricing library I mean a library handling equity, rate etc, hybrid products having several models (BS, LV, SV, LSV) having several numerical methods (analytic formula, MC, PDE FD/FE) I ...
159 views

### What is the best source to get 10 milliseconds time-series data for numerical computation?

I am working with 4th order Runge-Kutta method to compute a second order differential equation. For the best accuracy, I need a 10 milliseconds ohlcv time-series data. I know that I can build it ...
249 views

### Is there a more efficient data structure to implement binomial trees than 2d array?

I'm just curious what is the "industry standard" for implementing a binomial tree (if "standards" exist in this case). For simplicity, let's just talk about the simplest trees with recombining nodes. ...
69 views

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 ...
230 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-...
145 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}{\...
144 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}{\...
298 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. ...
Hi I want to simulate in Matlab the following stochastic integral: $x(t) = x(s) e^{-a(t-s)} + \sigma \int_s^t e^{-a(t-u)} dW_1(u)$ with $E[x(t) \vert F_s] = x(s) e^{-a(t-s)}$ \$Var[x(t) \vert ...