The Difference between European Option Prices in the Black-Scholes and NIG Models Computed with the DFT

This Demonstration shows the difference between the price of a European call option on stock (no dividend is paid and the interest rate is 0) in the Black–Scholes model and the model (due to Barndorff–Nielson) based on a centered exponential Normal Inverse Gaussian (NIG) Lévy process, as a function of the strike price. There are three model parameters that control the NIG process: steepness, asymmetry, and scale (the fourth parameter, location, is set to 0). The other control parameter is the time to expiry of the option.
Given the values of the NIG parameters one can compute instantaneous volatility, which together with the time to expiry and the initial stock price (set to 60) enables us to compute the price of the call option in the Black–Scholes model. The corresponding call option price for the NIG model is computed by means of the Esscher transfer and the fast Fourier transform.


  • [Snapshot]
  • [Snapshot]
  • [Snapshot]


The NIG process has been shown by several empirical studies to model option prices more accurately than the Black–Scholes model. Unlike the Black–Scholes model, there is no unique non-arbitrage price in the NIG model (and in other incomplete models, e.g. those involving jumps or stochastic volatility), so the price has to be chosen using economic considerations. We use the price given by replacing the real world measure by the equivalent martingale measure defined by means of the Esscher transform (see [2]). The prices for both models are computed for a range of strike prices using the discrete Fourier transform (DFT) following the method described in [1] and generalized in [3]. The method of computing NIG option prices used in this Demonstration is due to Gerber and Shiu [1]. A detailed study of option pricing with Lévy processes, of which the NIG process is an example, can be found in [2].
[1] P. Carr and D. B. Madan, "Option Valuation Using the Fast Fourier Transform," The Journal of Computational Finance, 2(4), 1999.
[2] H. U. Gerber and E. S. W. Shiu, "Option Pricing by Esscher Transforms," Transactions of the Society of Actuaries, 46, 1999 pp. 99–191.
[3] S. Raible, Lévy Processes in Finance: Theory, Numerics, and Empirical Facts, PhD thesis, University of Freiburg, 2000.
    • Share:

Embed Interactive Demonstration New!

Just copy and paste this snippet of JavaScript code into your website or blog to put the live Demonstration on your site. More details »

Files require Wolfram CDF Player or Mathematica.

Mathematica »
The #1 tool for creating Demonstrations
and anything technical.
Wolfram|Alpha »
Explore anything with the first
computational knowledge engine.
MathWorld »
The web's most extensive
mathematics resource.
Course Assistant Apps »
An app for every course—
right in the palm of your hand.
Wolfram Blog »
Read our views on math,
science, and technology.
Computable Document Format »
The format that makes Demonstrations
(and any information) easy to share and
interact with.
STEM Initiative »
Programs & resources for
educators, schools & students.
Computerbasedmath.org »
Join the initiative for modernizing
math education.
Step-by-Step Solutions »
Walk through homework problems one step at a time, with hints to help along the way.
Wolfram Problem Generator »
Unlimited random practice problems and answers with built-in step-by-step solutions. Practice online or make a printable study sheet.
Wolfram Language »
Knowledge-based programming for everyone.
Powered by Wolfram Mathematica © 2017 Wolfram Demonstrations Project & Contributors  |  Terms of Use  |  Privacy Policy  |  RSS Give us your feedback
Note: To run this Demonstration you need Mathematica 7+ or the free Mathematica Player 7EX
Download or upgrade to Mathematica Player 7EX
I already have Mathematica Player or Mathematica 7+