Hot answers tagged


You might want to set $a= \epsilon - d$ and write $\epsilon>0$ as a constraint. I guess $\textbf{lsqnonlin}$ is the suitable fonction for what you intend to do. I personnally like to use and play around with $\textbf{fmincon}$, which allows more flexibility and performs well, if you are willing to provide Jacobian and/or Hessian in algorithms options


Yes it can be done. However, bear in mind that a naive explicit FD scheme is not unconditionally stable (see CFL stability condition). As far as your initial/boundary conditions issue is concerned: [Time domain] Use terminal condition $V(S,T)=h(T)$ where $h(T)$ figures the payoff of the target derivative claim at maturity $T$ (e.g. $h(T)=(S-K)^+$ for a ...


Who gave you that idea? You absolutely can use Finite Differences for other PDEs. They are routinely used to solve hyperbolic PDEs (wave equation, both first and second order) and elliptic PDEs (steady state diffusion/heat equation). You can even mix and match the equation types and create PDEs that have characteristic of both hyperbolic and parabolic ...


Rather than thinking about the steps, think about the piecewise regions where your value is constant. When using the explicit scheme, time zero option value at any stock price for your simple digital option is basically just a function of which antecedent nodes (accounting for backwards timestepping) were above or below the strike. Slight modifications of ...

Only top voted, non community-wiki answers of a minimum length are eligible