[ProQuest: [...] denotes non US-ASCII text; see PDF]
Academic Editor:Turgut Özis
Department of Mathematics, Ryerson University, 350 Victoria St., Toronto, ON, M5B 2K3, Canada
Received 2 May 2016; Revised 21 July 2016; Accepted 26 July 2016
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
The objective of the paper is the study of the pricing of basket options under a multivariate Black-Scholes model, using two different polynomial approximations of a related conditional contract.
Our main contribution consists in the proposal of an efficient methodology that combines two approaches previously considered. First, we compute the conditional expected value of the payoff with respect to some of the underlying assets in order to reduce the problem to one-dimensional pricing; see [1]. Then, we use Taylor and Chebyshev polynomials to approximate the resulting conditional price. Numerical findings show that, within an equivalent accuracy, both approximations considerably reduce the computational time necessary to obtain the price, when compared with a standard Monte Carlo method.
While a Taylor expansion of the second order for spread pricing has been considered in [2], an approach based on Chebyshev expansion remains less explored. Recently, it has been considered in [3] for a single-asset option contract and in a multiasset setting in [4]. The latter combines with a Fourier series development, offering an interesting analysis of the error in the approximation. Our method conceptually differs from the works cited above in the way the expansion is carried on.
In addition, our expansions are presented in a general setting with the purpose of illustrating applications in higher dimensions, to other contracts and models, at the expense of somehow more complicated notations.
Basket options are multivariate extensions of univariate European calls or puts. A call basket option takes the weighted average of a group of d stocks (the basket) as the underlying asset and produces a payoff equal to the maximum of zero and the difference between the weighted average and the strike (or the opposite difference for the case of a put). Financial indices such as the S&P 500 and real options based on the difference between gas and oil prices are examples of such contracts.
In the case of European spread options, whose payoffs are given in terms of the spread of both prices at maturity, several approximations have been considered in the works of [2, 5-10], among others, where different ad hoc approaches are studied.
The approach to pricing by Taylor expansions can be traced back to [11], where a method to the price of one-dimensional derivatives is proposed. See also [12] for the use of Taylor in valuation of basket and Asian options, where the expansion is made about the characteristic function. Furthermore, in [10] the Taylor expansion is compared with other pricing techniques, proving to be effective and accurate for most values in the parametric space. In our case, the expansion is done on the function resulting from the conditional price, as opposed to a development based on the conditional strike price, as previously considered. Moreover, our methodology hinges on the calculation of mixed exponential-power moments of a Gaussian distribution and is extended to expansions about any point and order. Our point of view may allow for a better control on the approximation.
Taylor expansions produce reasonable approximations in terms of a simple closed-formula. Nonetheless, they are very sensible to the point around which the development is done. Moreover, as this expansion is local, it may introduce significant errors, albeit infrequent, at values far from the point where the expansion is considered, typically the conditional average price under the risk neutral measure. Significant errors for out-of-the-money contracts or extreme values in the parametric set are reported.
In order to overcome this potential problem, we study developments in terms of Chebyshev polynomials, which offer a uniform convergence of the conditional price on a predetermined closed interval. Interesting alternative approximations in terms of Fourier series can be found in [13, 14].
The organization of the paper is the following: in Section 2, we introduce the model and main notations and derive a Taylor approximation for d-dimensional basket options; in Section 3, we implement Taylor method to spread contracts while in Section 4 we study the case of a Chebyshev approximation and the sensitivity with respect to the spot prices; in Section 4, we discuss the numerical implementation and results; finally, in Section 5 we conclude.
2. Basket Derivatives and Taylor Expansions
We introduce some notations. Let (Ω,F,(Ft)t≥0 ,P) be a filtered probability space. We define the filtration FXt [: =]σ(Xs , 0<=s<=t)⊂Ft as the σ-algebra generated by the random variables {Xs , 0<=s<=t} completed in the usual way. Denote by Q an equivalent martingale risk neutral measure (EMM) and by EQ the expectation under Q. Quantities μa,b (k), MX (u) and MX (u,a,b) represent, respectively, the truncated kth moment and the moment generating function (m.g.f.) truncated on [a,b], while the function N(·) is the cumulative distribution function (c.d.f.) of a standard normal distribution.
The matrix A[variant prime] represents the transpose of matrix A=(aij ), while diag[...](A) is a column vector with components from the diagonal of the matrix A. On the other hand, A1/2 denotes a matrix such that A1/2A1/2 =A. For a d-dimensional vector Y, denote by Y~ the same vector, excluding the first component.
By r, we denote the (constant) interest rate or a vector with components equal to r. The matrix Id is the d×d identity matrix.
For l-times differentiable function f on Rd and a vector L=(l1 ,l2 ,...,ld ) with lk ∈N such that ∑k=1d [...]lk =l, DL f represents its mixed partial derivative of order l differentiated lk times w.r.t. yk .
The process of spot prices is denoted by S=(St(1) ,St(2) ,...,St(d) )t≥0 for 0<=t<=T, while Y=(Yt(1) ,Yt(2) ,...,Yt(d))t≥0 is the asset log-return process; they are related by [figure omitted; refer to PDF] The log-return process follows the dynamic [figure omitted; refer to PDF] where (Bt)t≥0 is a multivariate standard Brownian motion and Σ is a positive definite symmetric matrix with components (σij)i,j=1,2,...,d , with σii =σi2 .
We analyze European Basket options whose payoff at maturity T, for a strike price K, is given by [figure omitted; refer to PDF] where (wj)1<=j<=d are some deterministic weights and x+ =max[...](x,0).
The main examples are spread options, defined for d=2 by a payoff: [figure omitted; refer to PDF] Other related contracts are the so-called 3 : 2 : 1 crack spreads with payoff: [figure omitted; refer to PDF] where ST(1) , ST(2) , and ST(3) are, respectively, the spot prices of gasoline, heating oil, and crude oil.
Exchange options are spread options with K=0; see [15].
We start writing the price of the basket option, denoted by CB , in terms of a conditional price via the following elementary proposition.
Proposition 1.
Let CB be the price of a European Basket contract with maturity at T, strike price K, and payoff h(ST ) under the model given by (1) and (2). Then, [figure omitted; refer to PDF] where, for y∈Rd-1 , [figure omitted; refer to PDF] is the Black-Scholes price of a call option with maturity at T>0, volatility σYT(1) /Y~T =y , spot price S0(1) , and strike price [figure omitted; refer to PDF] with [figure omitted; refer to PDF] where [figure omitted; refer to PDF] and ΣY~ is the covariance matrix of the vector Y~T .
Proof.
From (2), [figure omitted; refer to PDF] in law, where Zd is a random variable with a multivariate standard normal probability distribution in Rd . Moreover, independently of Y~T , the random variable YT(1) has a univariate normal distribution. Thus, we can write [figure omitted; refer to PDF] in law, where Z1 is independent of YT and it has, conditionally on Y~T , a standard univariate normal distribution. Moreover, it is well known (see, e.g., [16]) that μYT(1) /Y~T and σYT(1) /Y~T are, respectively, given by (9).
Next, by the iterated property of the conditional expected value, we have [figure omitted; refer to PDF] where K1 (y)=K/w1 -∑j=2d [...](wj /w1 )S0(j)ey(j) .
Moreover, substituting (12) into (13), we have [figure omitted; refer to PDF] where [figure omitted; refer to PDF] Hence, we can recognize in C(Y~T ) the expression of a Black-Scholes price with the parameters mentioned above.
Remark 2.
Notice that K(y) may be negative. In this case, C(y) is not the Black-Scholes price, but the formula remains valid. Condition wj <=0, j=2,...,d, guarantees its positiveness, as it is the case of spreads and 3 : 2 : 1 crack spread contracts defined above.
Notice that C(y) is smooth enough to apply a kth-order Taylor development around any point y[low *] ∈Rd-1 . After conditioning on the last d-1 underlying assets, it leads to the approximated price: [figure omitted; refer to PDF] where L=(l1 ,l2 ,...,ld-1 ) and Rl ={L∈Nd-1 /l1 +l2 +[...]+ld-1 =l, 0<=lk <=l}.
Taking into account Proposition 1, a natural kth-order Taylor approximation around y[low *] =(y2[low *] ,y3[low *] ,...,yd-1[low *] ) of the price CB of a basket option with payoff h(ST ), defined as CB(k) [: =]e-rTEQ [C(k) (Y~T )], is obtained by [figure omitted; refer to PDF] after replacing (16) into the expression for CB in the proposition above.
Remark 3.
Notice that the approximation CB(k) depends only on the derivatives of the function C(y) with respect to y, which turns out to be the Black-Scholes price composed with the function K(y), and the mixed exponential-power moments of a Gaussian multivariate distribution.
Remark 4.
Sensitivities to the parameters can be computed by a similar approximation, as Greeks for a Black-Scholes option model are known. For example, the Delta with respect to the jth asset can be approximated by [figure omitted; refer to PDF]
3. Pricing Spreads Options by Taylor Approximations
In order to illustrate the method studied in the previous section, we consider the case of a bidimensional spread option under the model given by (1) and (2) with covariance matrix: [figure omitted; refer to PDF] Notice that, for convenience, we have slightly changed some notations.
From (11), the conditional distribution of YT(1) given YT(2) is normal with mean and variance given by [figure omitted; refer to PDF] Thus, we can write [figure omitted; refer to PDF] where Z~N(0,1) independent of YT .
The kth-order Taylor approximation in this case simplifies to [figure omitted; refer to PDF] Moreover, [figure omitted; refer to PDF] where [figure omitted; refer to PDF] Now, marginals of a multivariate normal distribution are normal too. Hence, we have that YT(2) ~N((r-(1/2)σ22 )T,Tσ22 ). The exponential-power moments can be calculated as [figure omitted; refer to PDF] with [figure omitted; refer to PDF] Next, integrating by parts, [figure omitted; refer to PDF] where n!! is the double factorial defined as the product of all odd numbers between 1 and n including both. When the set is empty, by convention, the product is equal to one.
As a consequence, for y[low *] =EQ (YT(2) ), we have [figure omitted; refer to PDF] After gathering all pieces and substituting in (22), we have the following result.
Proposition 5.
The kth-order Taylor approximation of a spread contract with maturity at T and strike price K, under the model described by (1) and (2), is given by [figure omitted; refer to PDF] with [figure omitted; refer to PDF] for m=1,2,...,k and E(0)=1, where EQZν =(ν-1)!! if ν is even or zero if it is odd, and [figure omitted; refer to PDF]
We now need to compute the derivatives of the function C(y) with respect to y. By a straightforward calculation from Black-Scholes formula, [figure omitted; refer to PDF] Also, [figure omitted; refer to PDF] where fZ is the density function of a standard normal random variable and [figure omitted; refer to PDF] Similarly, the second derivative is obtained as [figure omitted; refer to PDF] with [figure omitted; refer to PDF] In particular, when we expand the price around the point ymean =EQ (YT(2) )=(r-(1/2)σ22 )T we have the first and second approximations given, respectively, by [figure omitted; refer to PDF] More generally, expanding around y[low *] ∈R, the first two approximations denoted by p^1 (y[low *] ) and p^2 (y[low *] ), respectively, are given by [figure omitted; refer to PDF]
4. Approximation by Chebyshev Polynomials
We study an alternative approximation of the price via Chebyshev polynomials. For definition and their basic properties see, for example, [17].
Denoting by (Tk (y))k∈N the sequence of Chebyshev polynomials of first type on [-1,1], we consider the nth approximation of the function C(y) on the interval [a,b] described by equation in terms of Chebyshev polynomials, the one given by [figure omitted; refer to PDF] where (Tka,b (x))k∈N is the sequence of Chebyshev polynomials of first type on [a,b] defined by [figure omitted; refer to PDF] and the values (c^k)0<=k<=N are estimators of the corresponding coefficients in the Chebyshev expansion.
Chebyshev polynomials on [a,b] are orthogonal with respect to the scalar product defined as [figure omitted; refer to PDF] with weight function wa,b (x)=(1-(2(x-a)/(b-a)-1)2 )-1/2 .
Notice that [figure omitted; refer to PDF] Then, for k≠0 the coefficients in the expansion can be calculated as [figure omitted; refer to PDF] after changes of variables x=-1+(2/(b-a))(y-a) and x=cos[...]θ.
Also, [figure omitted; refer to PDF] From the trapezoidal rule to approximate Riemann integrals, the coefficients (ck)0<=k<=N can be estimated by an equidistant partition of N points on [0,π].
Chebyshev polynomials of first type can be written in terms of powers of the variable. From [18], [figure omitted; refer to PDF] where [figure omitted; refer to PDF] In a similar way to the case of Taylor polynomials, define the nth order Chebyshev approximation for the basket option price as [figure omitted; refer to PDF] The next theorem provides the Chebyshev approximation for the price of a basket option.
Theorem 6.
The nth order Chebyshev's approximation of the price CB of a European Basket option under the model given by (1) and (2) is given by [figure omitted; refer to PDF] where [figure omitted; refer to PDF] for k=0,1,2,...,n, [figure omitted; refer to PDF]
Proof.
Notice that [figure omitted; refer to PDF] Hence, [figure omitted; refer to PDF] From Proposition 1, we take into account (39) and (45) to get [figure omitted; refer to PDF] Now, [figure omitted; refer to PDF] Substituting the last expression into (53), we get (48).
5. Numerical Results
As a benchmark setting, we consider initial prices of S0(1) =100 and S0(2) =96 dollars, strike price of K=1 dollar, maturity at T=1 year, and an annual interest rate of r=3% under a bivariate Black-Scholes model, with a negative correlation ρ=-0.3 and respective volatilities equal to σ1 =0.3 and σ2 =0.1.
In Figure 1, we show the conditional price function on [-1,1] for the benchmark setting (blue line), together with the first and second Taylor expansions around the mean value of YT(2) ; that is, EQ (YT(2) )=(r-(1/2)σ22 )T (green and red lines, resp.). While the second-order approximation offers a reasonable local fit, significant errors may be found for values far from the mean. These errors may impact the price given by its conditional expected value under the risk neutral probability when, for example, volatilities are high.
Figure 1: Approximations of the conditional price based on Taylor polynomials using the benchmark set of parameters. The blue line describes the conditional price on [-1,1] while the green and red lines show first- and second-order expansions, respectively.
[figure omitted; refer to PDF]
Notice that, due to the concavity of function C, the first approximation underestimates the price. Not surprisingly, the second approximation estimates the price fairly well for values close to the point mean value while it is less accurate for values far from the mean. Although it seems as a drawback of the method, it does not constitute a serious problem as values far from the mean are infrequent, and thus the error in calculating the expected value from a Taylor approximation is fairly small.
A more promising result is obtained when an approximation of the conditional price is done via Chebyshev polynomials. Figure 2 represents, clockwise from top left, approximations of 4th, 10th, and 15th orders, respectively. Expansions of the 10th and 15th orders are practically indistinguishable from the original function. Coefficients in the expansion are calculated following a trapezoidal rule with 100 points on the interval [-1.5,1.5].
Figure 2: Approximations of the conditional price based on Chebyshev polynomials using the benchmark set of parameters. Clockwise, from top left are 4th, 10th, and 15th orders, respectively. The truncation interval is [-1.5,1.5].
[figure omitted; refer to PDF]
In Table 1, we show the impact of correlation on the prices of spread contracts obtained under the benchmark parameter set. Prices are computed by a Monte Carlo approach based on 10 million simulations of the asset prices following a geometric bivariate Black-Scholes model with correlated Brownian motions (column 2). In addition, we implement a second-order Taylor expansion and an approximation by Chebyshev polynomials of order n=15. We consider positive, negative, large, moderated, and weak correlations. In all cases, the Chebyshev approximation shows a notable agreement with Monte Carlo prices at a lesser computational cost.
Table 1: Spread prices for the benchmark parameters and several values of ρ, using Monte Carlo, a Taylor second-order approximation around y[low *] =0, and a Chebyshev approximation with 15 terms, a=-4 and b=0.25.
Correlation | Monte Carlo | Taylor second approx. | Chebyshev, n=15 |
ρ = - 0.1 | 14.2921 | 13.8709 | 14.2906 |
ρ = 0.1 | 13.56278 | 14.78882 | 13.5649 |
ρ = - 0.3 | 14.9734 | 15.0065 | 14.9629 |
ρ = 0.3 | 12.8085 | 12.7901 | 12.7903 |
ρ = - 0.5 | 15.6273 | 15.9238 | 15.6316 |
ρ = 0.5 | 11.9525 | 11.9646 | 11.9566 |
ρ = - 0.7 | 16.2421 | 17.5217 | 16.2521 |
ρ = 0.7 | 11.0315 | 11.1947 | 11.0529 |
Notice that Chebyshev expansion requires extra quadrature to compute the polynomial coefficients as opposed to Taylor where we simply need to evaluate the corresponding derivatives. Hence, the computational complexity in the former multiplies by the number of points used in the trapezoid rule.
In Table 2, we show average routine times for the price calculation following Monte Carlo, Taylor, and Chebyshev expansions. The average is taken upon runs for different parameters in the parametric set. Routines were implemented on a standard PC using MATLAB language. Chebyshev approximation works about 200 times faster than the standard Monte Carlo approach. Taylor approach has an even lesser computational time, but for large asset correlations it is not as accurate as the former.
Table 2: Average routine times for Monte Carlo, Taylor, and Chebyshev approaches on a PC using MATLAB language.
Monte Carlo | Taylor second approx. | Chebyshev, n=15 |
2.048 sec. | 0.0076 sec. | 0.012 sec. |
Due to the steepness of the function C(y), the Chebyshev approximation is sensible to the truncation interval [a,b]. In our numerical computations, we have used a=-4 and b=0.25. Within the range of parameter considered, most values of YT(2) lie on the selected interval [a,b]; hence, truncation does not affect the mixed exponential-power function by much. Otherwise, the approximation might be improved by taking C(a) if y<a and by C(b) if y>b.
The selection of a and b is a delicate point. A priori, there is no theoretical limit in the length of the interval. Moreover, a large interval will result in a better approximation. On the flip side, the effort in the estimation of coefficients ck in the polynomial expansion will increase for larger intervals. An intuitive criterion to select the truncation values might consist in looking at points where the function C(y) flattens out.
On the other hand, if the second asset has a high volatility, the length of the truncating interval should increase to avoid having too many points outside the interval.
The method is stable for the number of points considered in the trapezoidal rule. Also, the approximation gets close to the actual price after a fairly moderate number of polynomials. For n=10, the method shows a good approximation within an error in the order of a penny. For n=15 and n=20 the approximation improves even more. For approximations of larger orders, the gain in precision does not compensate the increase in computational time.
Figure 3 shows prices of a spread contract based on a Chebyshev approximation of order n=15. Maturity times range from one month to one year, while strike prices go from zero (exchange option) to 3 dollars. Results are consistent with an increase in the contract prices with higher maturity and their decrease with the increase of the strike price. Figure 3(c) at the bottom, conveniently scaled, shows the difference between both prices. The average relative error is 0.0075%.
Figure 3: Prices of a spread contract with maturity between one month and one year and strike price between 0 and 3 dollars under the benchmark parametric set. (c) Difference between prices obtained for both methods.
(a) Price computed by a Monte Carlo approach
[figure omitted; refer to PDF]
(b) Price computed by a Chebyshev expansion
[figure omitted; refer to PDF]
(c) Difference of prices between Monte Carlo and Chebyshev
[figure omitted; refer to PDF]
Figure 4 shows the price of the contract as a function of volatilities of both underlying assets. For symmetry, we only have considered σ1 ≥σ2 with volatilities ranging within 10%-50%. Again, we can observe a remarkable agreement between prices obtained by Monte Carlo and Chebyshev approximation. As expected, the price of the spread increases with the increment in the volatility. The average relative error is 0.0023%.
Figure 4: Prices of a spread contract with log-return volatilities between 10% and 50% under the benchmark parametric set. (c) Difference between prices obtained for both methods. We consider, without loss of generality, the case σ1 ≥σ2 .
(a) Price computed by a Monte Carlo approach, n=107
[figure omitted; refer to PDF]
(b) Price computed by a Chebyshev expansion, order 15
[figure omitted; refer to PDF]
(c) Difference of prices between Monte Carlo and Chebyshev
[figure omitted; refer to PDF]
6. Conclusions
We compare three methods to price spreads options under a bivariate Black-Scholes model with correlated Brownian motions versus a standard Monte Carlo approach. Our results show that Taylor approximation does not offer a uniform convergence, and hence a poor result when values are far from the point around the expansion is taken. For some values in the parameter set, it may affect the corresponding expected value.
The approximation based on Chebyshev polynomials seems to be appropriate in terms of the balance offered between accuracy and computational cost. Moreover, the method is suitable to be implemented in more general models provided the conditional distribution is available.
Acknowledgments
Both authors would like to recognize NSERC for its support.
[1] N. D. Pearson, "An Efficient Approach for Pricing Spread Options," The Journal of Derivatives , vol. 3, no. 1, pp. 76-91, 1995.
[2] M. Li, S.-J. Deng, J. Zhou, "Closed-form approximations for spread option prices and Greeks," Journal of Derivatives , vol. 15, no. 3, pp. 58-80, 2008.
[3] R. Pachon, "Numerical pricing of european options with arbitrary payoffs," SSRN Electronic Journal , 2016.
[4] M. Gass, K. Glau, M. Mahlstedt, M. Mair, "Chebyshev interpolation for parametric option pricing," https://arxiv.org/abs/1505.04648
[5] E. Kirk Correlation in the Energy Markets, in Managing Energy Price Risk , Risk Publications and Enron, London, UK, 1995.
[6] R. Carmona, V. Durrleman, "Pricing and hedging spread options," SIAM Review , vol. 45, no. 4, pp. 627-685, 2003.
[7] S. Borokova, F. J. Permana, H. Weide, "A closed form approach to the valuation and hedging of basket and spread options," The Journal of Derivatives , vol. 14, no. 4, pp. 8-24, 2009.
[8] A. Alvarez, M. Escobar, P. Olivares, "Spread Options under Stochastic Covariance and Jumps. Pricing two dimensional derivatives under stochastic correlation," International Journal of Financial Markets and Derivatives , vol. 2, no. 4, pp. 265-287, 2011.
[9] A. Venkatramanan, C. Alexander, "Closed form approximations for spread options," Applied Mathematical Finance , vol. 18, no. 5, pp. 447-472, 2011.
[10] M. Li, J. Zhou, S.-J. Deng, "Multi-asset spread option pricing and hedging," Quantitative Finance , vol. 10, no. 3, pp. 305-324, 2010.
[11] J. Hull, A. White, "The pricing of options on assets with stochastic volatilities," The Journal of Finance , vol. 42, no. 2, pp. 281-300, 1987.
[12] N. Ju, "Pricing Asian and basket options via Taylor expansion," The Journal of Computational Finance , vol. 5, no. 3, pp. 79-103, 2002.
[13] Q.-J. Meng, D. Ding, "An efficient pricing method for rainbow options based on two-dimensional modified sine-sine series expansions," International Journal of Computer Mathematics , vol. 90, no. 5, pp. 1096-1113, 2013.
[14] F. Fang, C. W. Oosterlee, "A novel pricing method for European options based on Fourier-cosine series expansions," SIAM Journal on Scientific Computing , vol. 31, no. 2, pp. 826-848, 2008.
[15] W. Margrabe, "The value of an option to exchange one asset for another," The Journal of Finance , vol. 33, no. 1, pp. 177-186, 1978.
[16] Y. L. Tong The Multivariate Normal Distribution , Springer, Berlin, Germany, 1989.
[17] J. C. Mason, D. C. Handscomb Chebyshev Polynomials , CRC Press Company, Boca Raton, Fla, USA, 2003.
[18] A. Gil, J. Segura, N. Temme Numerical Methods for Special Functions , SIAM, Philadelphia, Pa, USA, 2007.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
Copyright © 2016 Pablo Olivares and Alexander Alvarez. 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.
Abstract
We propose a closed-form approximation for the price of basket options under a multivariate Black-Scholes model. The method is based on Taylor and Chebyshev expansions and involves mixed exponential-power moments of a Gaussian distribution. Our numerical results show that both approaches are comparable in accuracy to a standard Monte Carlo method, with a lesser computational effort.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer