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{\partial V}{\partial M} + c \frac{\partial V}{\partial \phi} - rV = 0$$

with the terminal condition $V(T, S, M, \phi) = g(S, M, \phi)$

If I try to apply a Crack Nicholson method

$$2 \frac{df}{dx} = \frac{f(t+1, x_{i+1}) - f(t+1, x_i)}{\delta x} + \frac{f(t, x_{i+1}) - f(t, x_{i})}{\delta x}$$

to $S$, $M$ and $\phi$ the equation gets too complicated to solve.

So, should I only apply this to the $S$ variable and approximate the derivatives for $M$ and $\phi$ only usingw backward time values? How would you approach this?

  • $\begingroup$ What do you mean too complicated to solve? To complicated for you to write down, the computer computes this too slowly, memory explodes, unstable, impossible to debug, etc? $\endgroup$
    – oliversm
    Jul 7, 2020 at 12:57

1 Answer 1


You may want to have a look on Alternating Direction Implicit for solving multi-dimension PDE on finite difference method. The linear system will still be tridiagonal matrix.


Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge that you have read and understand our privacy policy and code of conduct.

Not the answer you're looking for? Browse other questions tagged or ask your own question.