American Journal of Applied Mathematics
Volume 4, Issue 4, August 2016, Pages: 181-185

Efficient Pricing of European Options with Stochastic Interest Rate Using Fourier Transform Method

Michael Chimezie Anyanwu, Roseline Ngozi Okereke

Department of Mathematics, Michael Okpara University of Agriculture, Umudike, Abia State, Nigeria

Email address:

(M. C. Anyanwu)
(R. N. Okereke)

To cite this article:

Michael Chimezie Anyanwu, Roseline Ngozi Okereke. Efficient Pricing of European Options with Stochastic Interest Rate Using Fourier Transform Method. American Journal of Applied Mathematics. Vol. 4, No. 4, 2016, pp. 181-185. doi: 10.11648/j.ajam.20160404.13

Received: May 6, 2016; Accepted: May 20, 2016; Published: July 13, 2016

Abstract: The inclusion of stochastic interest rate is an essential element of any realistic option pricing method. Therefore, the purpose of this paper is to incorporate interest rates in the Fourier transform method for pricing European options in exponential Levy models. With the assumption of stochastic independence between the underlying log asset price and the stochastic interest rate, we obtain a pricing of pure discount bond available in the market. Our method of valuation is to apply eigenfunction expansion to the variable that describes the evolution of the interest rate, and Fourier transform to the variable that describes the log asset price.

Keywords: Fourier Transform, Eigenfunction Expansion, Stochastic Interest Rate, Levy Process, Simple Harmonic Oscillator

1. Introduction

The first acceptable contribution to the field of quantitative modelling of derivative securities is given by Black and Scholes [1], in which the price of a European option on an underlying asset with price process given by a geometric Brownian motion was given in close-form. However, in most exponential Levy models and stochastic volatility models, analytic formula for the option price does not exist. Hence, alternative methods are sought. One of such approaches is the Fourier transform and inversion methods, since it has been observed that when characteristic function of the density function is mapped to the option payoff in Fourier space, option prices are easily calculated much easier for some of the much complex processes. The idea behind Fourier transform is to take an integral of the payoff function over the probability distribution function, where the probability distribution function is obtained by inverting the corresponding Fourier transform.

In mathematical finance, Fourier transform method was first used in Stein and Stein [3] volatility model, to obtain the distribution of the underlying asset price. This was followed by Heston [2], where through characteristic function, obtains an analytical formula for the valuation of European options when the volatility of the underlying asset price is stochastic. Following these two papers, quite a lot have being written on the application of Fourier transform and inversion methods in the valuation of more complex contingent claims. For example, Carr and Madan [4] introduces fast Fourier transform algorithm for the numerical valuation of the Fourier integral encountered in option pricing. Since then, fast Fourier transform algorithm has become an efficient mathematical tool used in mathematical finance in the valuation problems, e.g. see [10, 7, 8, 12, 13]. A more general pricing framework was developed in Lewis [9] by separating the underlying asset price from the option payoff through the use of Plancherel-Parseval theorem. Through his method, a variety of pricing formula for different types of contingent claims can be obtained. For other applications of Fourier transform methods to option pricing, see [5,6,8] and bibliographies therein.

All the existing Fourier transform methods in option pricing known in literature make a simplifying assumption that the interest rate in the market remain constant throughout the life of the option. This assumption is not generally true, especially when pricing options of long maturity, as interest rates are bound to change, and this level of change is not known before time (that is, it is stochastic). Therefore, the contribution of this paper to the existing literature on Fourier transform method for option pricing is to incorporate stochastic interest in pricing European style contingent claim under the assumption that the underlying log asset price  is an exponential Levy process. Therefore, to price the European option, we assume that the interest rate, , is a stochastic function of the Markov process  when state space is D = .

The method of valuation we adopt here is the joint application of the Fourier transform and eigenfunction expansion. To the best of our knowledge, this is the first time the two methods are being applied together. The idea is that we use the Fourier transform for the variable  that models the log-asset of the firm, and apply the eigenfunction expansion with respect to the variable  that describes the evolution of interest rate.

The study of eigenfunction expansion of linear operators was motivated by the seminal work of Mckean [14], where the spectral representation for the transition density of one dimensional diffusion process is obtained. Since then, eigenfunction expansion has been shown to be an important methodology for derivative pricing. For example, Gorovoi and Linetsky [20] models the Japanese interest rate using eigenfunction expansion method, to obtain analytical solutions for zero coupon bonds and bonds options under the Vasicek [19] and CIR [18] processes for the shadow rate. Lewis [27] uses eigenfunction expansion to price generalized European style options on stocks that pay constant dividends, and finds an estimate for the yield curve associated with the Merton [26] economy. Davydov and Linetsky [16] develops eigenfunction expansion approach for valuing options on scalar diffusion processes. Boyarchenko and Levendorskii [17] applies eigenfunction expansion to multifactor models of interest rate. For more applications of eigenfunction expansion method to derivative pricing, see in applying eigenfunction expansion, the problem of derivative pricing is reduced to that of solving a single eigenvalue-eigenfunction equation. Once the equation is solved, the approximate value of the derivative can be computed. In this paper, we do not do eigenfunction expansions with the infinitesimal generator of the pricing semigroup, rather we do expansion with infinitesimal generator  is the infinitesimal generator of the process . This ultimately leads to solving an infinite system of disjoint first order ordinary differential equation for the bond price, which is easy to solve. In this case, we obtain an infinite series that gives the bond price.

The remaining of this article is organized as follows: in the remaining part of section one, we discuss the preliminaries, which include Fourier transform and characteristic function, Levy process and its theory and eigenfunction expansion. In section two, our main method for pricing European option with stochastic interest rate is presented. This is the joint application of eigenfunction expansion and Fourier transform. In section three, we process our numerical result, discuss the result and conclude the article.

1.1. Fourier Transform and Characteristic Function

Let  be a piecewise continuous real function which satisfies the condition  The Fourier transform of  defined for any complex variable   is given by


This transform is known to exist and is analytic for all  in the strip  The inverse of the transform through which the function  can be recovered is


1.2. Brief on Levy Processes

Levy processes are becoming very popular in mathematical finance because they have been discovered to describe the observed reality of financial markets more accurately than models with Brownian motion such as geometric Brownian motion of Black and Scholes [1]. In real world, it is known that asset price paths do exhibit jump or spikes, and these have to be taken into account in order to manage risks efficiently. Moreover, the empirical distribution of asset returns exhibits fat tails and skewness. These show that the distribution of asset returns are far from being normal, as predicted in Brownian motion models. Levy processes are known to possess the appropriate properties that accurately describe all these properties of asset price and its returns in the real world and risk neutral world.

By definition, a Levy process  is stochastic processes with stationary and independent increments, and it is continuous in probability. A general 1-dimensional Levy process can be represented as


Where is a 1-dimensional Brownian motion,  is a compound Poisson process which includes the jumps of , and is a compensated compound Poisson process which includes the jumps of  with  The characteristic function of the distribution of  can be represented in the form

The function ψ is called the characteristic exponent for . the general form of characteristic exponent of any Levy process can be deduced from the Levy. Khintchine formula


where  is a measure on \0 satisfying

The triple  is called the Levy triple for Xt. The proof of this Levy – Khintchine formula can be found in Cont and Tankov [21]. The known exponential Levy models in financial modeling literature correspond to assigning different values to . For example, the case where F( is a finite measure and  gives rise to jump-diffusion Levy processes, e.g. Merton [22] and Kou [23] models. The case where  an infinite measure, and gives rise to pure jump model of infinite activity, e.g. Variance Gamma (VG) model of Madan, et al [24], Normal Inverse Gaussian model of Barndorff-Nielson [25], etc. As an example in this article, we use a VG process whose density of jumps is given by

Here,  is an indicator function of the interval (a, b). The coefficients c+ (respectively, c-) gives the intensity of upwards jumps (respectively, downward jumps). The characteristic exponent  of the VG process can be obtained using (4):

1.3. Eigenfunction Expansion for the Bond Price

Consider a pure discount bond that pays the owner a certain unit amount of currency at maturity date T. The price  of this bond at time  is given by


where  is the interior of , and  is the risk-neutral measure chosen for pricing. In the CIR (1985) model, the variable  satisfies the stochastic differential equation


Where  is one dimensional Brownian motion. In (6), we assume that the feller condition,  is satisfied to ensure that the point 0 is not attained by the process . Applying Feymann-Kac theorem to  we obtain the backward Kolmogorov equation


Subject to  is the infinitesimal generator of the process . We assume that  is diagonalizable so that we can expand  into the series


where  are the orthogonal eigenfunctions of  are the coefficients of expansion given by


 is the measure on  chosen so that the infinitesimal generator  becomes unbounded in the Hilbert space  with discrete spectrum, and

Is the inner product in . We can now expand the bond price  in the series


where . If we substitute (10) into (7), with , we get , n= 0,1,2… Hence


Subject to . Then, (11) is an infinite system of disjoint scalar first order ordinary differential equations for . This is solved to get  for each . therefore, the eigenfunction expansion (10) for the bond price becomes


Hence, if the constant of expansion  converges fast, and if the eigenvalues and eigenfunctions are known, then with moderately few terms, the series (12) can be truncated appropriately to find the approximation to the bond price. Suppose the risk-free interest rate is quadratic function of the stochastic variable  can be transformed to the operator of the simple Harmonic oscillator whose eigenvalues and eigenfunctions are well known. The transformed infinitesimal generator becomes , where  and  are constants to be determined in the transformation,  is the operator of the simple harmonic oscillator.

The normalized eigenfunctions of  can be generally written as

With, and the corresponding eigenvalues are

Where  are the Hermite polynomials, and  means the largest integer which is not greater than . it is verifiable that  is the initial eigenfunction that corresponds to the initial eigenvalue

We use the following theorem from [17] which we state without proof.


1.  The eigenfunctions of the transformed infinitesimal generator  and that of  are the same.

2.  The eigenvalues of  where  are the eigenvalues of  for j = 0,1,2,…

With these, the eigenfunction expansion for the bond price becomes


The detail of this transformation can be found in [17]

The graph below shows how the bond prices from the eigenfunction expansion method approximates closely, the bond prices from the analytical solution.

Figure 1. Approximate and Exact Bond Price for T = .5 years.

2. Pricing European Option with Stochastic Interest Rate

We begin by choosing a risk neutral measure Q, suitable for pricing options on the underlying stock or index  with stochastic interest rate , where  is a Levy process of exponential type . In the presence of stochastic interest rate, the time 0 < T price of a European option with payoff and expiry date  can be written as


Where  are independent, we can write (14) to be


Where  is the risk neutral value of a pure discount bond that pays one unit of currency at maturity date .

Define . Then, the value of the option can be written as


Suppose there exists some  such that  is integrable, we can decompose H into Fourier integral to get




Putting (16) into the pricing formula (14), we get


We can apply Fubinis theorem to change the order of taking expectation and integration to obtain


Using the characteristic function for the distribution of X, we can write (20) as


Where  is the so-called characteristic exponent for the Levy process X. the payoff for European call and put options with strike price K are  respectively. Hence, using (18), we have  for call option and  for put option. The major task in applying Fourier transform to pricing European option is the evaluation of the integral in (21). To do this, we truncate the integral and discretize the integrand. This involves choosing a convenient grid  to obtain


where  and  is the grid size

3. Numerical Example

In this section, we implement the method discussed in We consider this paper. We experiment with the following parameters: For the Levy process, we use  (since we are considering pure jump process), For the CIR interest rate process, we use . We experiment with different stock prices with fixed strike price, K = 100 for both European call and put options. The result of our calculations are presented in Table 1 and Table 2 for call and put options respectively. In this tables, N represent the number of simulations used in (22).

Table 1. Call option values for K = 100, T = 1.

S/N 20 30 40 50 60 70 80 90
120 1.3781 1.3745 1.3753 1.3750 1.3751 1.3751 1.3751 1.3751
150 6.2955 6.3002 6. 2996 6. 2996 6. 2996 6. 2996 6. 2996 6. 2996
180 11.3756 11.3751 11.3758 11.3759 11.3758 11.3758 11.3758 11.3758
200 14.7766 14.7728 14.7727 14.7728 14.7729 14.7729 14.7729 14.7729

Table 2. Put option values for K = 100, T = 1.

S/N 20 30 40 50 60 70 80 90
20 2.9252 2.9255 2.9255 2.9255 2.9255 2.9255 2.9255 2.9255
40 5.2309 5.2310 5.2312 5.2311 5.2312 5.2312 5.2312 5.2312
60 5.6639 5.6667 5.6672 5.6672 5.6672 5.6672 5.6672 5.6672
80 4.2085 4.2084 4.2042 4.2042 4.2043 4.2043 4.2043 4.2043

4. Summary

Pricing European options with short maturity date on a given stock, may not require change in interest rate. But when the maturity date is long, it is necessary to consider the possibility of stochastic change in the level of interest rate in the market. In this article, we have succeeded in incorporating stochastic interest rate in the Fourier transform method for valuing European-style contingent claim on a stock driven by an exponential Levy process. Moreover, the tables above reveals that the method of simplified trapezoid rule adopted here for evaluating the Fourier integral converges fast, thereby allowing the option value to be computed fast.


  1. Black, F. and M. S. Scholes (1973). The Pricing of Options and Corporate Liabilities, Journal of Political Economy 81, 637-654.
  2. Hoston, S. L. (1993) A closed-form Solution for Options with Stochastic Volatility with Application to Bond and Currency Options. Review of Financial Studies 6 (2), 327-343.
  3. Stein, E. and Stein, J. (1991) Stock Price Distributions with Stochastic Volatility: An Analytic Approach. Review of Financial Studies 4, 727-752.
  4. Carr, P. P. and Madan, D. B. (1999) Option valuation using the fast Fourier transform, Journal of Computational Finance 2 (4), 61-73.
  5. Bailey, D. and Swarztrauber, P. (1991). The fractional Fourier transform and applications, SIAM Review 33, 389-404.
  6. Attari, M. (2004) Option Pricing using Fourier Transforms: A Numerically Efficient Simplification. Working Paper, Charles River Associates. Available at SSRN:
  7. Benhamou, E. (2000). Fast Fourier Transform for Discrete Asian Options, EFMA 2001 Lugano Meetings. Available at SSRN:
  8. Hurd, T. R. and Zhou, Z. (2009). A Fourier transform method for spread option pricing, Preprint, arXiv: 0902.3643v1.
  9. Lewis, A. (2001). A simple option Formula for General Jump-Diffusion and other Exponential Levy Processes, Envision Financial Systems and Option City. net, California. Available at
  10. Borak, S., Detlefsen, K. and Hardle, W. (2005) FFT Based Option Pricing. Discussion Paper, Sonderforschungsbereich (SFB) 649, Humboldt Universität Berlin.
  11. Dempster, M. A. and Hong, S. S. (2002) Spread Option Valuation and the Fast Fourier Transform. Technical Report WP26/2000, University of Cambridge.
  12. Itkin, A. (2005) Pricing options with VG model using FFT. arXiv: physics/0503137v1
  13. O’Sullivan, C. (2005) Path Dependent Option Pricing under Levy Processes. EFA 2005 Moscow Meetings Paper. Available at SSRN:
  14. Mckean, Henry P. J. (1956). Elementary solution for certain parabolic differential equations. Transaction of the American mathematical society, 82 (2), 519-548.
  15. Linetsky, V. (2004). The spectral decomposition of the option value. International Journal of Theoretical and Applied Finance, 7 (3), 337-384.
  16. Davydov, V. and Linetsky, V. (2000). Pricing Options on scalar diffusions: An Eigenfunction Expansion Approach, Operations Research 512, 185-209.
  17. Boyarchenko, N. and Levendorskii, S. Z. (2004). Eigenfunction Expansion method in multi-factor quadratic term-structure models. Working paper. Available at SSRN:
  18. Cox, J. C., Ingersol, J. E and Ross, S. A. (1985). A theory of the term structure of interest rates. Econometrica 55, 349-383.
  19. Vasicek, O. A. (1977). An equilibrium characterization of the term structure. Journal of Financial Economics 5, 177-188.
  20. Gorovoi, V. and Linetsky, V. (2004) Black’s model of interest rates as options, eigenfunction expansion and Japanese interest rates. Mathematical Finance, 14 (1), 48-78.
  21. Cont, R. and Tankov, P. (2004). Financial Modeling with Jump Processes. (Financial Mathematics Series), Chapman and Hall/CRC, New York.
  22. Merton, R. C. (1976). Option pricing when the underlying stock returns are discontinuous. Journal of Financial Economics 3, 125-144.
  23. Kou, S. G. (2002). A jump – diffusion model for option pricing. Management Science 48, 1086-1101.
  24. Madan, D. B., Carr, P. P. and Chang, E. C. (1998). The Variance Gamma Process and Option Pricing. European Finance Review 2, 79-105.
  25. Barndorff-Nielson, O. Processes of Normal Inverse Gaussian Type. Finance and Stochastics, 2, (1998), 41-68.

Article Tools
Follow on us
Science Publishing Group
NEW YORK, NY 10018
Tel: (001)347-688-8931