Pricing Arithmetic Asian Options under Hybrid Stochastic and Local Volatility

Abstract

Recently, hybrid stochastic and local volatility models have become an industry standard for the pricing of derivatives and other problems in finance. In this study, we use a multiscale stochastic volatility model incorporated by the constant elasticity of variance to understand the price structure of continuous arithmetic average Asian options. The multiscale partial differential equation for the option price is approximated by a couple of single scale partial differential equations. In terms of the elasticity parameter governing the leverage effect, a correction to the stochastic volatility model is made for more efficient pricing and hedging of Asian options.

(ProQuest: ... denotes non-US-ASCII text omitted.)

Min-Ku Lee 1 and Jeong-Hoon Kim 2 and Kyu-Hwan Jang 2

Academic Editor:K. S. Govinder

1, Department of Mathematics, Sungkyunkwan University, Suwon, Gyeonggi-do 440-746, Republic of Korea
2, Department of Mathematics, Yonsei University, Seoul 120-749, Republic of Korea

Received 31 July 2013; Revised 24 November 2013; Accepted 26 November 2013; 8 January 2014

This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

1. Introduction

Since the well-known work of Black and Scholes [1] on the classical vanilla European option, there has been concern about the pricing of more complicated exotic options. An exotic option is a derivative which has a payoff structure more complex than commonly traded vanilla options. They are usually traded in over-the-counter market or embedded in structured products. Also, the pricing of them tends to require more complex methods than the classical Black-Scholes approach.

This paper is concerned with one of the exotic options called an Asian option. This option is a path dependent option whose final payoff depends on the paths of its underlying asset. More precisely, the payoff is determined by the average value of underlying prices over some prescribed period of time. The name of "Asian" options is known to come from the fact that they were first priced in 1987 by David Spaughton and Mark Standish of Bankers Trust when they were working in Tokyo, Japan (cf. [2]). The main motivation of creating these options is that their averaging feature could reduce the risk of market manipulation of the underlying risky asset at maturity. Since Asian options reduce the volatility inherent in the option, the price of these options is usually lower than the price of classical European vanilla options. Note that there are two types of Asian options depending on the style of averaging: one is the arithmetic average Asian option and the other is the geometric average Asian option.

Since there is no general analytical formula for the price of Asian option, a variety of techniques have been developed to approximate the price of this option. Subsequently, there has been quite an amount of literature devoted to studying this option. For instance, Geman and Yor [3] computed the Laplace transform of the price of continuously sampled Asian options. However, there is a problem of slow convergence for low volatility or short time-to-maturity cases as indicated by Fu et al. [4]. Linetsky [5] derived a new integral formula for the price of Asian options but with the same convergence problem. Apart from this approximation technique, there are Monte Carlo simulation approach and partial differential equation (PDE) approach. Monte Carlo simulation methods typically require variance reduction techniques and also have to take into account the discretization error caused by discrete sampling. Refer to Kemna and Vorst [6] for a discussion of the pricing of Asian options with Monte Carlo methods. On the other hand, the PDE methods must deal with an extra state variable representing the running sum of the underlying process, which leads to an issue for reducing the dimension of the PDE. Refer to Ingersoll [7] and Rogers and Shi [8] and Vecer [9] for the PDE methods.

It is well known that the constant volatility assumption, on which a review of literature quoted above is based, for the underlying asset price is severely in contrast with many empirical studies which demonstrate the skew or smile effect of implied volatility, fat-tailed and asymmetric returns distributions, and the mean-reversion of volatility. Thus a number of alternative underlying models have been proposed. The constant elasticity of variance model by Cox [10], a stochastic volatility model by Heston [11] or Fouque et al. [12], and a Levy model by Carr et al. [13] are among those representative ones that can reproduce an empirically reasonable outcome.

So, it is desirable to study Asian options based on these alternative models. In fact, there are a number of recent studies along the lines of this type of extension. For instance, B. Peng and F. Peng [14] for the CEV model, Fouque and Han [15] for a stochastic volatility model, and Lemmens et al. [16] for a Levy type model are among those works extending the price for the Black-Scholes framework with constant volatility. As long as we understand, up to now, however, there has been no work for Asian options based on a hybrid stochastic and local volatility model. Recently, hybrid stochastic and local volatility models have become an industry standard for the pricing of derivatives and several financial institutions have incorporated those models into their systems [17]. Therefore, it is worth studying Asian options on the hybrid models.

This paper studies the pricing of an arithmetic Asian option under a hybrid stochastic and local volatility model which was introduced by Choi et al. [18], where the volatility is given by the product of a multiscale stochastic process and a power (the elasticity of variance) of the underlying's price. The hybrid nature of this volatility enables us to capture the leverage effect produced by the constant elasticity of variance (CEV) model as well as the smile effect of implied volatility, fat-tailed and asymmetric returns distributions, the tendency of the volatility process to revert towards a long-term mean at a certain rate, and a degree of correlation between the randomness of volatility and the randomness of underlying's price produced by a "pure" stochastic volatility (SV) model. This hybrid model is called the SVCEV model. So, this paper will generalize both [14, 15] into an approximation problem for the price of arithmetic average Asian options based on the SVCEV model.

This paper is structured as follows. In Section 2, we review the SVCEV model introduced by [18]. A multiscale partial differential equation for the price of an arithmetic average Asian option is obtained in Section 3. Section 4 is devoted to obtain an approximated option price under the condition of fast mean-reverting volatility. In Section 5, the effect of hybrid structure on the pricing of the Asian option is illustrated by a numerical study.

2. The SVCEV Model

In this section, we establish a partial differential equation for the price of Asian floating strike call option based on the SVCEV model.

As introduced by [18], the SVCEV model for the underlying's price is given by the stochastic differential equations (SDEs): [figure omitted; refer to PDF] where μ and m are some constants, γ is greater than - 1 , α and β are positive constants, and _{W t} and _{B t} are correlated Brownian motions such that d Y9; W , B _{YA; t} = ρ d t for some ρ . If ρ > 0 , then _{S t} may fail to be a true martingale since _{Y t} may go infinite. Refer to [19]. Therefore, ρ < 0 is assumed here. It is observed in most financial markets that there is a minus correlation or leverage effect between stock price and volatility impact. However, some commodity markets show the opposite effect (the inverse leverage effect). So, the model under the negative correlation condition can apply to many financial markets but has limits for some commodity markets. Also, generally speaking, the correlation ρ may rely on time but it is assumed to be a constant for simplicity. In fact, in most real situations, it is taken to be such. We do not specify the function f but it has to satisfy a growth condition to avoid the nonexistence of moments of _{S t} . It is assumed in this paper that 0 < _{c 1} ...4; f ...4; _{c 2} < ∞ for some constants _{c 1} and _{c 2} .

From the I to formula, the solution Y of the second equation in (1) is an ergodic process given by the Ornstein-Uhlenbeck process: [figure omitted; refer to PDF] and so _{Y t} ~ N ( m + ( _{Y 0} - m ) ^{e - α t} , ( ^{β 2} / 2 α ) ( 1 - ^{e - 2 α t} ) ) , which has an invariant distribution given by N ( m , ^{β 2} / 2 α ) . Later, we will use notation Y9; · YA; for the average with respect to this invariant distribution; that is, [figure omitted; refer to PDF] for arbitrary function g .

Now, we take the process _{Y t} to be a fast mean-reverting process as in Fouque et al. [20]. This means that we take into account a fast time scale volatility factor as a major driving force for the volatility movement of the underlying asset. This assumption can be justified by an empirical analysis of high-frequency S&P 500 index data. The analysis confirms that volatility is fast mean-reverting when looked at over the time scale of a derivative contract although it reverts slowly to its mean in comparison to the tick-by-tick fluctuations of the index value. See [20] for details. Also, this assumption will provide us with analytic tractability for an approximation for the price of the Asian option as also shown in [15]. So, in terms of a parameter, say ... , which is taken as positive and small, we assume that [figure omitted; refer to PDF] where ν ~ ...AA; ( 1 ) which implies that the long run magnitude of volatility fluctuations remains fixed as a constant. Then, under a risk-neutral probability measure ^{P *} , we have [figure omitted; refer to PDF] where ^{W t *} and ^{B t *} are two Brownian motions under the measure ^{P *} whose correlation is given by d Y9; ^{W t *} , ^{B t *} YA; = ρ d t and Λ ( y ) denotes the combined market price of risk given by [figure omitted; refer to PDF] Here, λ is the market price of volatility risk, which is assumed to be a bounded function depending on y .

3. A PDE for Option Price

In this section, we utilize the generalization of Vecer's dimension reduction technique given by Fouque and Han [15] to derive a two-space dimensional PDE representation for the option price.

In this paper, a payoff function for arithmetic average Asian options is given by [figure omitted; refer to PDF] where h is a function satisfying [figure omitted; refer to PDF] for any real number α . Note that, when _{K 1} = 0 , it becomes a payoff for fixed strike Asian options, whereas, when _{K 2} = 0 , it becomes a payoff for floating strike Asian options. So, the option price P ( t , x , y ; T , _{K 1} , _{K 2} ) is defined by [figure omitted; refer to PDF] under a risk-neutral measure ^{P *} .

First, we would like to replicate the averaged process ( 1 / t ) ^{∫ 0 t} ... _{S u} d u with a portfolio [figure omitted; refer to PDF] Here, _{α t} is assumed to be a nonrandom function. By the self-financing strategy, we note that [figure omitted; refer to PDF] The assumption of the nonrandomness of _{α t} yields [figure omitted; refer to PDF] Thus, using (11) and (12), one can obtain [figure omitted; refer to PDF] By integration of (13), we obtain [figure omitted; refer to PDF] which can be written as [figure omitted; refer to PDF] For technical reason, if we choose the trading strategy _{α t} as [figure omitted; refer to PDF] and the initial portfolio price _{X 0} as [figure omitted; refer to PDF] then the final portfolio price _{X T} becomes ( 1 / T ) ^{∫ 0 T} ... _{S t} d t - _{K 2} and so the general payoff function of arithmetic average Asian options given by (8) becomes h ( _{X T} - _{K 1 S T} ) . Refer to [9, 21] or [15] for details. In terms of the portfolio process _{X t} , the price of an arithmetic average Asian option at time t = 0 is given by [figure omitted; refer to PDF]

Next, we change the probability measure ^{P *} into the probability measure ^{P ~ *} defined by [figure omitted; refer to PDF]

Then we are ready to obtain the option price as the solution of a partial differential equation.

Theorem 1.

Let _{ψ t} : = _{X t} / _{S t} . Then the Asian option price at time t = 0 can be expressed by [figure omitted; refer to PDF] under the measure ^{P ~ *} , where [figure omitted; refer to PDF] If u ( t , ψ , y ; T , _{K 1} , _{K 2} ) : = ^{E ~ *} [ h ( _{ψ T} - _{K 1} ) |" _{ψ t} = ψ , _{Y t} = y ] , then u is given by the solution of the PDE: [figure omitted; refer to PDF] with the terminal condition [figure omitted; refer to PDF]

Proof.

By the Itô formula, we obtain the following two results: [figure omitted; refer to PDF]

Then we have [figure omitted; refer to PDF] where [figure omitted; refer to PDF] which is a Brownian motion under the measure ^{P ~ *} by the Girsanov theorem (cf. [22]).

From (8) and (18), we have [figure omitted; refer to PDF] where [figure omitted; refer to PDF]

If u ( t , ψ , y ; T , _{K 1} , _{K 2} ) : = ^{E ~ *} [ h ( _{ψ T} - _{K 1} ) |" _{ψ t} = ψ , _{Y t} = y ] , the Asian option price at time t = 0 is given by [figure omitted; refer to PDF] Then, by the Feynman-Kac formula (cf. [22]) and (25), one can obtain (22).

Once the solution P ( 0 , s , y ; T , _{K 1} , _{K 2} ) is determined, from the result of Fouque and Han [15], the price P ( t , s , y ; T , _{K 1} , _{K 2} ) is given by [figure omitted; refer to PDF] where u is the solution of the PDE (22) at time t = 0 .

4. Multiscale Analysis

In this section, we are interested in the solution of the multiscale PDE (22) in the asymptotic form u = _{u 0} + ... _{u 1} + ... + ^{... i / 2} _{u i} + ... . Then one can come up with single scale PDEs much easier to solve than the PDE (22) itself as follows.

Substituting the asymptotic form of u into the PDE (22) yields [figure omitted; refer to PDF] with the terminal condition [figure omitted; refer to PDF] where [figure omitted; refer to PDF]

Note that _{[Lagrangian (script capital L)] 0} is the infinitesimal generator of the Ornstein-Uhlenbeck process _{Y t} . Then we have a hierarchy of PDEs as follows: [figure omitted; refer to PDF] with the following terminal condition: [figure omitted; refer to PDF]

In the following two theorems, we obtain PDEs for the leading order term _{u 0} as well as the first correction term _{u 1} .

Theorem 2.

Assume that the term _{u i} does not grow as much as _{u i} ~ ^{e y 2 / 2} with respect to the variable y for i = 0,1 . Then the leading order term _{u 0} is independent of y and satisfies the PDE problem: [figure omitted; refer to PDF] where [figure omitted; refer to PDF]

Proof.

From the ...AA; ( 1 / ... ) term in (34) and the fact that _{[Lagrangian (script capital L)] 0} is the generator of _{Y t} , the growth condition of _{u 0} with respect to the variable y leads to the fact that _{u 0} does not depend on y . So, it is represented as [figure omitted; refer to PDF] without dependence on the variable y . From this fact, equation _{[Lagrangian (script capital L)] 0 u 1} + _{[Lagrangian (script capital L)] 1 u 0} = 0 becomes _{[Lagrangian (script capital L)] 0 u 1} = 0 and so, by the same reason as for _{u 0} , _{u 1} also does not depend on the variable y : [figure omitted; refer to PDF] Since _{[Lagrangian (script capital L)] 1} is an operator in which every term contains the derivative with respect to the variable y , _{[Lagrangian (script capital L)] 0 u 2} + _{[Lagrangian (script capital L)] 1 u 1} + _{[Lagrangian (script capital L)] 2 u 0} = 0 becomes _{[Lagrangian (script capital L)] 0 u 2} + _{[Lagrangian (script capital L)] 2 u 0} = 0 . By applying the Fredholm alternative (cf. Ramm [23]) to this Poisson equation with respect to the operator _{[Lagrangian (script capital L)] 0} , we have Y9; _{[Lagrangian (script capital L)] 2 u 0} YA; = 0 ; that is, [figure omitted; refer to PDF] By adding the boundary condition (35) for _{u 0} to this PDE, we obtain the PDE problem (36) as desired.

Next, we obtain a PDE for the first correction term _{u 1} .

Theorem 3.

Assume that the term _{u i} does not grow as much as _{u i} ~ ^{e y 2 / 2} with respect to the variable y for i = 0,1 . Then the first order correction term _{u 1} is independent of y and satisfies the PDE problem: [figure omitted; refer to PDF] where _{V 2} and _{V 3} are constants, respectively, given by [figure omitted; refer to PDF] Here, [varphi] is the solution of [figure omitted; refer to PDF]

Proof.

Since Y9; _{[Lagrangian (script capital L)] 2 u 0} YA; = 0 and _{[Lagrangian (script capital L)] 0 u 2} + _{[Lagrangian (script capital L)] 2 u 0} = 0 , we obtain [figure omitted; refer to PDF] Then, from the definition of [varphi] , we obtain [figure omitted; refer to PDF] for some function c ( t , ψ ) independent of y .

By applying the Fredholm alternative to the Poisson equation _{[Lagrangian (script capital L)] 0 u 3} + _{[Lagrangian (script capital L)] 1 u 2} + _{[Lagrangian (script capital L)] 2 u 1} = 0 , we have equation Y9; _{[Lagrangian (script capital L)] 1 u 2} + _{[Lagrangian (script capital L)] 2 u 1} YA; = 0 ; that is, [figure omitted; refer to PDF]

Combining (45) and (46), we obtain (41).

5. Approximate Option Price

From the results in Section 4, we formally obtain the first order approximation: [figure omitted; refer to PDF] where _{u 0} and _{u 1} are independent of the unobservable variable y and they are given by the PDEs in Theorems 2 and 3, respectively. Note that _{u 0} solves a homogeneous equation with the nonzero final condition while _{u 1} solves a nonhomogeneous equation with the zero final condition. In terms of accuracy of the approximation (47), if the payoff function h is smooth enough, it follows straightforwardly from [20] that the approximation is of order ... in the pointwise convergent sense. Otherwise, it requires a regularization of the payoff as done for European vanilla options in [24]. Since the relevant argument (with an exact form of the density function for _{α T} ) is omitted in this paper, we limit ourselves to the case of regularized payoffs here.

If we define _{P 0} and _{P 1} by [figure omitted; refer to PDF] respectively, then, from (30), the option price P ( t , s , y ; T , _{K 1} , _{K 2} ) has an approximation given by [figure omitted; refer to PDF]

In this section, we compute numerically the leading order price _{P 0} and the first correction term _{P ~ 1} : = ... _{P 1} by using the finite difference method (the Crank-Nicolson method). The solution has the truncation error ...AA; ( ( Δ t ^{) 2} ) + ...AA; ( ( Δ ψ ^{) 2} ) , where Δ t = 0.005 and Δ ψ = 0.0104 .

Figure 1 shows _{P 0} and _{P ~ 1} at time t = 0 for three different values of γ (the elasticity parameter). The parameter values used in this figure are r = 0.06 , σ ¯ = 0.5 , T = 1 , _{K 1} = 0 , _{K 2} = 2 , _{V ~ 2} ... = ... _{V 2} = - 0.01 , and _{V ~ 3} ... = ... _{V 3} = 0.004 . Note that the parameter ... is absorbed by _{V ~ 2} and _{V ~ 3} . We call the price _{P 0} + _{P ~ 1} the (approximate) SVCEV price for the arithmetic average Asian option. It can reduce to the price for the well-known models. In particular, the price _{P 0} with γ = 0 is the Black-Scholes price. The price _{P 0} + _{P ~ 1} with γ = 0 corresponds to the stochastic volatility model studied by Fouque and Han [15]. Let us call it the SV price. When the leverage effect ( γ < 0 ) takes place, the SVCEV price becomes lower than the SV price. When the inverse leverage effect ( γ > 0 ) takes place, the SVCEV price becomes higher than the SV price. Regardless of the value of γ , the first order stochastic volatility correction term _{P ~ 1} is positive and has the maximum value near s = _{K 2} .

The leading order term _{P 0} and the first order correction term _{P ~ 1} of the price of an arithmetic average Asian option are drawn for three different values of the elasticity of variance parameter γ , where r = 0.06 , σ ¯ = 0.5 , T = 1 , _{K 1} = 0 , _{K 2} = 2 , _{V ~ 2} = - 0.01 , and _{V ~ 3} = 0.004 .

(a) _{P 0}

[figure omitted; refer to PDF]

(b) _{P ~ 1}

[figure omitted; refer to PDF]

Figure 2 shows the correction term _{P ~ 1} with respect to the group parameters _{V ~ 2} and _{V ~ 3} . This figure shows that the correction term _{P ~ 1} decreases as the parameter _{V ~ 2} increases whereas it increases as the _{V ~ 3} increases. Figure 2 also shows that as the parameter γ becomes larger, the correction term _{P ~ 1} becomes larger.

The correction term _{P ~ 1} with respect to the group parameters _{V ~ 2} and _{V ~ 3} ; r = 0.06 , σ ¯ = 0.5 , T = 1 , S = 2 , _{K 1} = 0 , and _{K 2} = 2 .

(a) _{P ~ 1}

[figure omitted; refer to PDF]

(b) _{P ~ 1}

[figure omitted; refer to PDF]

6. Conclusion

A frequent criticism of the stochastic volatility or local volatility models for path dependent options is that they do not produce deltas precise enough for hedging purposes. So, relevant industry experts recommend using a hybrid stochastic local volatility model of their own development for best pricing option products. See, for instance, [25]. In this paper, by transforming the path dependent problem for Asian options in the hybrid SVCEV model into a European vanilla style problem, we approximate the volatility dependent price by the nonvolatility dependent price. The approximate price provides not only a correction to the price of [15] for different values of the elasticity of variance parameter but also a correction to the CEV price based on the assumption of fast mean-reverting stochastic volatility. The elasticity of variance plays an important role in characterizing volatile markets as well as differentiating commodity markets from financial markets. Refer to [26]. The numerically solved Asian option prices show that both the leading order option price (the CEV price) and the corrected price (the approximate SVCEV price) increase as the elasticity of variance goes up. This result can provide more efficient risk hedging by choosing an appropriate elasticity parameter based on the observed market volatility or the chosen commodity.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors would like to thank the anonymous reviewers for their valuable comments and suggestions to improve the quality of the paper. The research was supported by the National Research Foundation of Korea NRF-2013R1A1A2A10006693.

References

[1] F. Black, M. Scholes, "The pricing of options and corporate liabilities," Journal of Political Economy , vol. 81, pp. 637-654, 1973.

[2] P. Wilmott Paul Wilmott on Quantitative Finance , John Wiley & Sons, New York, NY, USA, 2006.

[3] H. Geman, M. Yor, "Bessel processes, Asian option, and perpetuities," Mathematical Finance , vol. 3, no. 4, pp. 349-375, 1993.

[4] M. Fu, D. Madan, T. Wang, "Pricing continuous Asian options: a comparison of Monte Carlo and Laplace transform inversion methods," The Journal of Computational Finance , vol. 2, no. 2, pp. 49-74, 1998.

[5] V. Linetsky, "Spectral expansions for Asian (average price) options," Operations Research , vol. 52, no. 6, pp. 856-867, 2004.

[6] A. G. Z. Kemna, A. C. F. Vorst, "A pricing method for options based on average asset values," Journal of Banking and Finance , vol. 14, no. 1, pp. 113-129, 1990.

[7] J. E. Ingersoll Theory of Financial Decision Making , Rowman & Littlfield, Savage, Md, USA, 1987.

[8] L. C. G. Rogers, Z. Shi, "The value of an Asian option," Journal of Applied Probability , vol. 32, no. 4, pp. 1077-1088, 1995.

[9] J. Vecer, "Unified pricing of Asian options," Risk , vol. 15, no. 6, pp. 113-116, 2002.

[10] J. Cox, "Notes on option pricing I: constant elasticity of variance diffusions," (Reprinted in The Journal of Portfolio Management , vol. 22, pp. 15-17, 1996) Working Paper , Stanford University, Stanford, Calif, USA, 1975.

[11] S. L. Heston, "Closed-form solution for options with stochastic volatility with applications to bond and currency options," The Review of Financial Studies , vol. 6, pp. 327-343, 1993.

[12] J. P. Fouque, G. Papanicolaou, R. Sircar, "Asymptotics of a two-scale stochastic volatility model," Equations aux Derivees Partielles et Applications , of Articles Dedies a Jacques-Louis Lions, pp. 517-526, Gauthier-Villars, Paris, France, 1998.

[13] P. Carr, H. Geman, D. B. Madan, M. Yor, "Stochastic volatility for Lévy processes," Mathematical Finance , vol. 13, no. 3, pp. 345-382, 2003.

[14] B. Peng, F. Peng, "Pricing arithmetic Asian options under the CEV process," Journal of Economics, Finance and Administrative Science , vol. 15, no. 29, pp. 7-13, 2010.

[15] J.-P. Fouque, C.-H. Han, "Pricing Asian options with stochastic volatility," Quantitative Finance , vol. 3, no. 5, pp. 353-362, 2003.

[16] D. Lemmens, L. Z. Liang, J. Tempere, A. de Schepper, "Pricing bounds for discrete arithmetic Asian options under Lévy models," Physica A , vol. 389, no. 22, pp. 5193-5207, 2010.

[17] D. F. DeRosa Options on Foreign Exchange , John Wiley & Sons, Hoboken, NJ, USA, 2011., 3rd.

[18] S.-Y. Choi, J.-P. Fouque, J.-H. Kim, "Option pricing under hybrid stochastic and local volatility," Quantitative Finance , vol. 13, no. 8, pp. 1157-1165, 2013.

[19] L. B. G. Andersen, V. V. Piterbarg, "Moment explosions in stochastic volatility models," Finance and Stochastics , vol. 11, no. 1, pp. 29-50, 2007.

[20] J.-P. Fouque, G. Papanicolaou, K. R. Sircar, K. Solna Multiscale Stochastic Volatility for Equity, Interest Rate and Credit Derivatives , Cambridge University Press, Cambridge, UK, 2011.

[21] S. E. Shreve Stochastic Calculus for Finance II: Continuous-Time Models , Springer, New York, NY, USA, 2008.

[22] B. Øksendal Stochastic Differential Equations , Springer, New York, NY, USA, 2003.

[23] A. G. Ramm, "A simple proof of the Fredholm alternative and a characterization of the Fredholm operators," The American Mathematical Monthly , vol. 108, no. 9, pp. 855-860, 2001.

[24] J.-P. Fouque, G. Papanicolaou, R. Sircar, K. Solna, "Singular perturbations in option pricing," SIAM Journal on Applied Mathematics , vol. 63, no. 5, pp. 1648-1665, 2003.

[25] G. Tataru, F. Travis, "Stochastic local volatility," Quantitative Development Group, Bloomberg Version 1, 2010

[26] J.-H. Kim, J. Lee, S.-P. Zhu, S.-H. Yu, "A multiscale correction to the Black-Scholes formula," Applied Stochastic Models in Business and Industry . In press