arxiv: v3 [q-fin.cp] 13 Nov 2017

Transcription

1 April, 8 Quantitative Finane SingCFSEuroV3reformatted To appear in Quantitative Finane, Vol., No., Month XX, 35 Singular Fourier-Padé Series Expansion of European Option Pries arxiv:76.679v3 q-fin.cp] 3 Nov 7 Tat Lung Ron Chan University of East London, Water Lane, Stratford, UK, E5 4LZ v. released May 6 We apply a new numerial method, the singular Fourier-Padé SFP method invented by Drisoll and Fornberg,, to prie European-type options in Lévy and affine proesses. The motivation behind this appliation is to redue the ineffiieny of urrent Fourier tehniques when they are used to approximate pieewise ontinuous non-smooth probability density funtions. When tehniques suh as fast Fourier transforms and Fourier series are applied to prie and hedge options with non-smooth probability density funtions, they ause the Gibbs phenomenon; aordingly, the tehniques onverge slowly for density funtions with jumps in value or derivatives. This seriously adversely affets the effiieny and auray of these tehniques. In this paper, we derive priing formulae and their option Greeks using the SFP method to resolve the Gibbs phenomenon and restore the global spetral onvergene rate. Moreover, we show that our method requires a small number of terms to yield fast error onvergene, and it is able to aurately prie any European-type option deep in/out of the money and with very long/short maturities. Furthermore, we ondut an error-bound analysis of the SFP method in option priing. This new method performs favourably in numerial experiments ompared with existing tehniques. Keywords: singular Fourier-Padé series, European-type options, Lévy proesses, affine proesses. Introdution The appliation of Fourier tehniques to priing and hedging finanial derivatives has flourished in omputational and mathematial finane over the last deade beause most of the underlying proesses driving derivatives have a harateristi funtion i.e., the Fourier transform of the probability density funtion PDF. In general, two Fourier tehniques are applied to option priing: fast Fourier transforms FFTs and Fourier series, suh as Fourier osine COS series. To explain the first tehnique, FFTs in option priing, we first onsider a European option priing formula that is a risk-neutral expetation of the disounted payoff of the option driven by an underlying proess. This formation naturally implies the integration of the disounted payoff and the PDF of the underlying proess. As the underlying proess has a harateristi funtion, one an apply an inverse ontinuous Fourier transform to reover the PDF, disretise the ontinuous Fourier integral and use an FFT algorithm that omputes disrete Fourier transforms in ON log N operations to obtain the option pries f. Carr and Madan 999, Lewis, Lipton, Chourdakis 4, Itkin 5, Lord et al. 8. The main disadvantage of using these methods is that they require thousands of grid points and onsiderable omputational time to reah an aeptable level of auray. Another approah is to use Fourier series to expand a PDF in terms of a partial sum of orthogonal basis funtions. The basis funtion an be omplex an exponential, e.g., Chan 6, or COS, e.g. Fang and Oosterlee 9a, funtion. Eah oeffiient in the series is approximated Corresponding author. t.l.han@uel.a.uk

2 April, 8 Quantitative Finane SingCFSEuroV3reformatted via the harateristi funtion of the PDF. By integrating eah basis funtion with the disounted payoff funtion of the option, we an obtain an option prie. This tehnique originated in Fang and Oosterlee 9a to prie European-style options with COS basis funtions. The suess of Fourier COS series in European option priing has lead to the rapid expansion of priing and hedging options with early-exerise features and exoti options, like Asian, multi-asset or barrier options e.g., Leentvaar and Oosterlee 8, Fang and Oosterlee 9b,, Zhang and Oosterlee 3. Compare with FFTs, the advantage of using Fourier series for option priing is that they an ahieve a global spetral exponential onvergene rate and require fewer summation terms as long as the governing PDF is suffiiently smooth. However, using any type of Fourier series to represent a C ν pieewise ontinuous non-smooth funtion, e.g., a non-smooth PDF, is notoriously fraught. Disontinuities ause the Gibbs phenomenon, whih has two important onsequenes for the Fourier partial sum of length N: i failure to onverge at the jump, and ii pointwise onvergene elsewhere at the rate ON. More generally, if the funtion f and its derivatives up to order ν are ontinuous but f ν is disontinuous i.e., f has a jump of order ν, then the global onvergene rate is ON ν. The impat of the Gibbs phenomenon an lead to inaurate priing and hedging when the approximate option pries are generated via FFT or Fourier series methods at or around the jumps. Ruijter et al. 3 attempt to resolve the Gibbs phenomenon and improve the algebrai-index onvergene rate of the plain COS method by adding filters. They all this the filter COS method. In their paper, they intensively test six different filters, e.g., the Fejér filter and the exponential filter, and disover that the exponential filter gives a better algebrai-index onvergene rate for priing European options. The method is implemented to represent the non-smooth PDF of the Variane Gamma VG model. Although they provide some improvements over the urrent COS method, filtering only works well away from the peak of the VG PDF; at the peak, the approximation performs worse. The method also requires almost a thousand COS partial summation terms to reah aeptable auray when the PDF is non-smooth. Aording Ortiz-Graia and Oosterlee 3, 6, another disadvantage of using the COS method is that the COS exhibits periodiity in the viinity of the integration boundaries, and longmaturity option round-off errors may aumulate near the domain boundaries. A series of papers of using wavelets, e.g., B-spline and Shannon wavelets, have addressed this problem. Wavelet methods are similar to Fourier series methods; they use wavelets to represent a PDF and then integrate the disounted payoff funtion with wavelet basis funtions to reah the losed-form representation of an option priing formula. Compared with the COS method, the wavelet methods are flexible and aurate for priing options with long maturities. However, the hoie of the wavelet basis funtions determines whether the method an ahieve a spetral onvergene rate. In some ases, the B-spline and Haar wavelets annot ahieve spetral onvergene f. Table in Ortiz-Graia and Oosterlee 3 and Ortiz-Graia and Oosterlee 6. Although Ortiz-Graia and Oosterlee use Shannon wavelets Shannon Wavelets Inverse Fourier Tehnique SWIFT to orret this problem and ahieve spetral onvergene in some ases where the PDF is smooth f. Table and Table 4 in Ortiz-Graia and Oosterlee 6, the auray of SWIFT heavily relies on a sale parameter. The sale parameter is determined by the speed of deay of a harateristi funtion, whih in turn, determines the auray of the approximation of the orresponding PDF. Although, in general, Ortiz-Graia and Oosterlee 6 onsider or the best hoie of the sale, there is no theoretial reason for this hoie, and different sale parameters are hosen for different PDFs throughout their paper. Finally, in our numerial tests, the auray of the COS method is very sensitive to small and A vetor spae in whih funtions are ν ontinuously differentiable. A funtion is alled pieewise ontinuous on an interval if the funtion is made up of a finite number of ν times differentiable ontinuous piees.

3 April, 8 Quantitative Finane SingCFSEuroV3reformatted large option pries. The COS method tends to be less aurate when we measure small option pries f. Table 8 in Setion 8. To irumvent the aforementioned disadvantages of FFT, Fourier series and wavelet methods, we propose a novel method, the singular Fourier-Padé SFP method, whih exhibits the following harateristis: i global spetral onvergene rate for pieewise ontinuous PDFs, ii fast error onvergene with fewer partial summation terms required, iii aurately pries any European-type option with the features of deep in/out of the money and very long/short maturities, iv onsistently aurate for approximating large or small option pries throughout, and v does not require a sale parameter to adjust its auray. Why do we hoose the SFP method? Compared with the SFP method, the Fourier-Padé tehnique is not a better hoie. The Fourier-Padé tehnique f. Chisholm and Common 98, Small and Charron 988, Geer 995 is a famous tehnique for approximating a non-smooth funtion suh that spetral onvergene an be ahieved away from the jumps and onvergene is not globally degraded. Nevertheless, Drisoll and Fornberg, note that the fundamental limitation of the Fourier-Padé method is that the use of poles to approximate branh uts is ineffiient, and logarithmi singularities translated by the jumps on the unit irle are very diffiult for the Padé approximants to simulate. Aordingly, the method fails to onverge at the jump. To address this problem, Drisoll and Fornberg, add proper logarithmi branh singularity terms into the Fourier-Padé approximation proess. They all this method the singular Fourier-Padé method. In all numerial tests f. Drisoll and Fornberg,, the SFP an aelerate error onvergene to the true solutions for analyti periodi funtions or non-smooth funtions, and the global effets of the Gibbs phenomenon due to a jump are largely eliminated. If the funtion is very diffiult, when jump loations are known in advane, the method still appears to onverge globally and exponentially and offer 4 6 digits of auray at the jumps. Most importantly, in addition to omparing the Fourier-Padé tehnique, Drisoll and Fornberg ompare the SFP method with other two famous methods of overoming the Gibbs phenomenon the Gegenbauer method Gottlieb and Shu 997, a method to projet the Fourier partial sum onto a spae spanned by a Gegenbauer polynomial, and singularity removal Ekhoff 997, a method of singularity removal for funtions with jump disontinuities. In the numerial tests omparing these three methods f. Drisoll and Fornberg,, the SFP method yields better global spetral onvergene and faster error onvergene than do the other methods. Based on all these benefits of the SFP method, we hose it over the other available numerial methods. The remainder of this paper is strutured as follows. Setion provides an introdution. Setion desribes the SFP method. Setion 3 introdues the finanial stohasti models we examine in this paper. Setion 4 desribes and proves the formulation of the SFP option priing/greeks formulae for different styles of European options. In Setion 5, we desribe the SFP algorithm and the Fourier-Padé method to find the jump loations on a non-smooth funtion. Setion 6 desribes the hoie of trunated integration intervals. Setion 7 provides the error analysis of the SFP method in option priing and hedging. Setion 8 disusses, analyses and ompares the numerial results of the SFP method with those of other numerial methods. Finally, we onlude and disuss possible future developments in Setion 9.

4 April, 8 Quantitative Finane SingCFSEuroV3reformatted. Singular Fourier-Padé Interpretation and Corretion of the Gibbs Phenomenon If we onsider a funtion f with a formal power series representation k= b kx k, and a rational funtion defined by R N,M = P N /Q M, where P N and Q M are the polynomials of P N x = N p n x n and Q M x = n= M q m x m, respetively, then we say that R N,M = P N /Q M is the linear Padé approximant of order N, M of the formal series satisfying the ondition N M p n x n n= m= M+N q m x m k= m= b k x k = Ox N+M+. Here, f is approximated by M+N k= b k x k, To obtain the approximant RN, M, we simply alulate the oeffiients of polynomials P N and Q M by solving a system of linear equations. To obtain {q m } M m=, we first normalise q = to ensure that the system is well-determined and has a unique solution in. Then, we onsider the oeffiients for x N+,..., x M+N, and we an yield a Toeplitz* linear system: b N+ b N b N b N+ M q. b N+ b N+ b.. N bn+ M q = b N+M b N+ b N+ b N q M One {q m } M m= is known, {p n} N n= is found through the terms of order N and less in. This yields p = Bq, where b ij = b i j. For example, if N = M, one obtains p p p N b q q q M. = b b b N b b Computing Padé approximants through linear algebra, as we have done here, is simple but not neessarily the most effiient or stable numerial method f. Drisoll and Fornberg,. Now, suppose that f be a pieewise analyti funtion defined on the interval π, π, with s jump loations in f at t = ζ s π, π, s =,..., S. We define a jump as an atual disontinuous point on f or a disontinuous point appearing after its derivatives. Then the omplex Fourier series CFS representation is given by ft = k= b k e ikt, with b k = π fte ikt dt. 5 π π The transformation z = e it, whih maps the interval π, π onto the unit irle in the omplex A Toeplitz matrix or diagonal-onstant matrix is an invertible matrix in whih eah desending diagonal from left to right is onstant.

5 April, 8 Quantitative Finane SingCFSEuroV3reformatted plane, transforms the Fourier series into the following Laurent series in z, whih an be split into fz = k= b k e ikt = bk z k + k= k= b k z k = f + z + f z, 6 where the prime sums indiate that the zeroth term should be halved. The Fourier-Padé approximation of f ± is omprised of polynomials The resulting approximant is then defined as P ± N z = Q± M zf ± z + Oz N+M+, z. 7 P + N z Q + M z + P N z Q M z. 8 However, Drisoll and Fornberg, note that this approximant does not reprodue very well at/around the jump loations of the funtion and make the approximation inaurate. Therefore, they suggest that every jump in value of f at t = ζ an be attributed to a logarithm of the form log z e iζ. 9 This logarithmi singularity in f ±, whih is diffiult for the Padé approximant to simulate, an be exploited to enhane the approximation proess. This is the rationale behind the SFP method introdued in Drisoll and Fornberg,. We modify the Fourier-Padé approximant 7 to obtain the following ondition: where P ± N z + L = Q± M zf ± z + Oz U+, L = S s= L ± s z log z e iζs for some polynomials L s, s =,..., S, and U is determined by S and the degrees of P N, Q M and L s. If we extend the SFP method to support any pieewise analyti real funtion f in a finite interval a, b] with a set of jump loations {ζ s } S s= a, b] appearing in f, the CFS representation of the funtion is defined as fx = Re k= b k e i π b a kx ], with b k = b a b a π i fxe kx b a dx. Here, Re represents the real part of the funtion. As we fous on approximating a real funtion, we an further obtain ] fx = Re b k e i π kx b a + b. 3

6 April, 8 Quantitative Finane SingCFSEuroV3reformatted Based on this representation, we denote z as exp i π b a x power series of f equal to Re and then transform f into a trunated ] U b k z k + b. 4 The transformation z = exp i π b a x also suggests that the jump loation ζ translates into log z, where ε = e i π ζ b a 5 ε in f. Finally, applying 4 and 5, we an reah the final form of the SFP approximant given by where S P + N z + U L + N s z log z/ε s = b k z k + b Q + M z + OzU+, 6 s= P + N z = N n= p nz n, Q + M z = M m= q mz m, L + N s z = N s n l s= n s z ns, s =,..., S, U = N + M + S + S s= N s. 7 One we have determined the unknown oeffiients of {p n } N n=, {q m} M m= and {l n s } Ns 5, the SFP representation of fx an be formulated as fx = Re P + N z + S s= L+ N s z log z/ε s Q + M z n s= f. Setion, with z = e i π b a x Stohasti Models of the Asset Dynamis We assume fritionless finanial equity markets and no arbitrage, and take as given an equivalent martingale measure Q hosen by the market. All stohasti proesses defined in the following are assumed to live on the omplete filtered probability spae Ω, F, {F t } t, Q. Standard referenes for exponential Lévy proesses an be found in Shoutens 3 and Cont and Tankov 4; for affine proesses, see Duffie et al Exponential Lévy Proesses The stok prie proess S t t under Q driven by an exponential Lévy proess an be defined: S T = S t e r qt t+xt Xt+ωT t 9 = S t e r q+ωt t+xt t, where X T, X t and X T X t are all Lévy proesses. As Lévy proesses have independent stationary inrements, we an say that X T X t = X T t. Throughout the paper, r and q denote the onstant risk-free interest rate and the onstant dividend yield, respetively; S t represents the

7 April, 8 Quantitative Finane SingCFSEuroV3reformatted known stok prie at time t; and S T represents the random stok prie at time T. The ondition that S T e r qt t t is a martingale will be guaranteed by an appropriate hoie of the meanorreting ompensator ω as follows: ω = T t E e XT t, where E e XT t is assumed to be finite for all t T. Given a Lévy proess XT t t, define the orresponding harateristi funtion as follows: ϕu = Ee iuxt t ] = e T tφu, u R. Here, φu is a ontinuous funtion with the Lévy-Khinthine representation given by φu = A u iγu + e iuχ + iuχ χ νdχ, χ X T t, with harateristi triplet A, γ, ν. There is substantial onsideration of exponential Lévy proesses in stok proess modelling. Due to spae limitations, we annot show all the proesses in this paper; rather, we seletively investigate the proesses due to their popularity or their numerial implementation hallenge. However, this does not mean that our method only works on the proesses we seleted. As the paper evolves, we an see that as long as a proess has a harateristi funtion, the option prie and its risk driven by the proess an be properly approximated via our method. 3.. Exponential Affine Proesses An exponential affine harateristi funtion is given by ϕu = Ee iuxt ] = expcu, t + X Du, t. 3 Here, the funtion Cu, t is fully haraterised by dcu, t/dt = Du, t and the funtion Du, t satisfies a Riati equation Duffie et al. 3. More preisely, we an find funtions Cu, t and Du, t with Cu, t = and Du, t = u, suh that M t = expcu, T t + X t Du, T t 4 is a martingale proess. In the family of affine proesses, we fous on investigating the Heston model in this paper. The stohasti differential equation SDE of the Heston model is written as ds t = r qs t dt + y t S t dw,t, 5 dy t = λȳ y t dt + η y t dw,t, 6 where L t and y t denote the stohasti log-asset prie variable and the variane of the asset prie proess, respetively. In this proess, the speed of mean reversion λ, the mean level of variane ȳ and the volatility of volatility η are onstant values greater than or equal to zero. Additionally, the Brownian motions W,t and W,t are orrelated with the orrelation oeffiient ρ s, and ω = Ee i ixt ] is the mean-orreting ompensator. The model harateristi funtion fits in the

8 April, 8 Quantitative Finane SingCFSEuroV3reformatted general affine harateristi funtion framework and is given by with ϕu = exp iu r qt + ω + y λȳ tλ η iρ s ηu E log F e Et F η e Et F e Et λ iρ sηu E +, u R, 7 E = λ iρ s ηu + u + izη, F =λ iρ s ηu E/λ iρ s ηu + E. This harateristi funtion is uniquely speified beause we take x + yi suh that its real part is nonnegative, and we restrit the omplex logarithm to its prinipal branh. In this ase, as Lord and Kahl prove, the resulting harateristi funtion is the orret one for all omplex numbers z in the analyti strip of the harateristi funtion. In the SDE, we have two possible onditions with respet to λ, ȳ and η: λȳ η, 8 λȳ < η. 9 The model satisfies the Feller property if 8 holds; otherwise, 9 holds. If a proess fulfils the property, the proess never hits zero; onversely, if it does not, the proess an reah. Condition 9 is a very important property for the Heston SDEs beause they an only have a unique solution when we speify a boundary ondition at. In mathematial finane, the hosen boundary ondition is that the proess remains at. We define this as an absorbing boundary ondition. When the proess reahes and is allowed to leave, we all it a refleting boundary. These two boundary onditions are ruial for priing early-exerise options, inluding Amerian options and barrier options. 4. Singular Fourier-Padé Representation of European Option Pries and their Option Greeks 4.. Option Priing Formulae In this setion, we derive losed-form formulae for European-style options using the SFP method. Considering the PDF f of a stohasti proess, the urrent log-prie x := log S, the strike prie of K and maturity T t, we an express the option prie V x, K, t starting at time t with its ontingent laim paying out GS T, K as follows: V x, K, t = e rt t EGS T, K S t = e x 3 = e rt t EGS t e XT Xt, K 3 = e rt t Ge x+χ, Kfzdz, χ X T X t. 3

9 April, 8 Quantitative Finane SingCFSEuroV3reformatted Furthermore, if we hoose an interval, d] satisfying the ondition fχe iuχ dχ Ee iuxt Xt ] = ϕu, 33 where ϕu is a harateristi funtion of X T X t, we an approximate the priing formula: V x, K, t e rt t Ge x+χ, Kfχdχ. 34 Theorem When a dividend-paying risky asset prie proess S t t with a traeable, analytial harateristi funtion ϕ has a urrent asset prie e x = S, risk-free interest rate r and ompounded ontinuous dividend q, the SFP priing formula of a European vanilla all option driven by this proess with maturity T and strike prie K is with P + V x, K, t = e rt t N Re z + ] S s= L+ N s z log z/ε s Q + M z, 35 z = e i π x+log K d ε s = e i π d ζs P + N z = N n= p nz n, Q + M z = M L + N s z = N s n l s= n s z ns, s =,..., S. m= q mz m, 36 Here, ζ s is the jump in V x, K, t. Proof: First, the payoff funtion of a European all option is Ge x+χ, K = max e x+χ K, in 34. To onform to the SFP approximation framework, we first transform the payoff: max e x+χ K, = K max e x+χ log K,. 37 Aordingly, by replaing x + χ log K with y, we have a new form of V x, K, t denoted as e rt t K maxe y, f y x + log K dy. 38 As we intend to make the SFP method more effiient, we define a trunated omputational interval, d] f. Setion 6, whih satisfies ondition 33, to replae, ]. Then, V x, K, t is reformulated as e rt t K maxe y, f y x + log K dy =e rt t Ke y f y x + log K dy. 39 Using the Fourier transform shift theorem and the CFS expansion shown in 3, we express f y x + log K as Re b k e i π b a ky + b ], 4

10 April, 8 Quantitative Finane SingCFSEuroV3reformatted where b k = d π i fye ky d dy e i π k x+log K d and b = d We substitute 4 into 39 and apply Fubini s theorem; V x, K, t an be omputed as ] d e rt t Ke y Re b k e i π ky d + b dy = e rt t Re d fydy d Ke y dy π i fye ky d dy e i π In the equation above, basi alulus implies that Ke y e i π d ky dy = d k x+log K fydy. 4 Ke y e i π d ky dy+. 4 Kd π i e k+d Kd π i d e kd d, 43 iπk + d iπk Ke y dy = Ke d Kd. 44 In the meantime, beause of ondition 33, we an also see that π i fye ky d dy ϕ π d k and fydy ϕ =. 45 For the sake of simpliity in 4, we set B k = d ϕ Ĝ k = π Kd iπk+d d k, B = d ϕ = d e, π i k+d π i d e kd d, Ĝ = Ke d Kd Kd iπk 46 to obtain a simplified form: V x, K, t = e rt t Re B k Ĝ k e i π d k x+log K + B Ĝ ]. 47 To express our final priing formula with the SFP representation, we approximate B k Ĝ k e i π d k x+log K + B Ĝ 48 in 47 with 6. To do so, we first replae x + log K with a variable y in 48, and then, this allows us to set exp i π d y equal to z. The transformation z = exp i π d y maps the interval, d] onto the unit irle in z. This hange also transforms the jumps ζ along V x, K, t into z with

11 April, 8 Quantitative Finane SingCFSEuroV3reformatted the form of ε = exp i π d ζ. Finally, expressing 48 with a new variable of z, we have B k Ĝ k z k + B Ĝ. 49 Substituting the equation above with f z in 6, we obtain the approximant given by S P + N z + U L + N s z log z/ε s = B k Ĝ k z k + B Ĝ Q + M z + OzU+ 5 s= P + N z = N n= p nz n, Q + M z = M m= q mz m, L + N s z = N s n l s= n s z ns, s =,..., S, ε s = e i π d ζs, U = N + M + S s= N s. One we an determine the unknown oeffiients of {p n } N n=, {q m} M m= and {l n s } Ns the algorithm shown in Setion 5 and replae n s= 5 in 5 via B k Ĝ k e i π d k x+log K + B Ĝ with P + N z + S s= L+ N s z log z/ε s Q + M z, z = exp i π x + log K d in 47, we reah our first SFP priing formula for European vanilla all options, as in 35. Q.E.D. We an apply the same tehnique to seek the SFP priing formula of vanilla put options. To ahieve this, we first transform the put payoff funtion Ge x+χ, K = max K e x+χ, into K max e x+χ log K,. 5 Then, we follow the proof of Theorem to reah the CFS representation of V x, K, t given by with Ĝ k = Ĝ = e rt t Re Ke y e i π d ky dy = B k Ĝ k e i π d k x+log K + B Ĝ ], 53 Kd π i e k+ Kd π i d e k d, iπk + d iπk 54 Ke y dy = Ke K. 55 In the same manner shown in the proof of Theorem, we approximate the CFS expansion in 53 with the SFP approximant 6. We an therefore yield the same expression of 35 with different oeffiient values of {p n } N n=, {q m} M m= and {l n s } Ns n s=.

12 April, 8 Quantitative Finane SingCFSEuroV3reformatted From the vanilla all and put priing formulae, we an see that they share the same format as the SFP approximant of 35. It is no different for other European-type options disussed in this paper. The only differene is that the values of {p n } N n=, {q m} M m= and {l n s } Ns n s= in the SFP formula of eah option are ompletely different from one another. This is attributed to the jumps ζ in eah option and their payoff funtions, whih lead to different Fourier transformations of Ĝk and Ĝ in 47. For the reader s information, we enumerate all the options we investigate in this paper and their orresponding payoff funtions Ge x+χ, K and Fourier transformations G k and G in Tables and. 4.. Option Greeks Now, we turn our attention to deriving the option Greeks. Although aurately valuing finanial laims plays a key role in finanial modelling, the risk management hedging of these derivative instruments is equally important. Finanial institutions manage option risk when they sell options to their lients through the analysis of the Greeks, whih are defined as the sensitivities of the option prie to single-unit hanges in the values of state variables or parameters. Suh sensitivities represent the different dimensions of the risk assoiated with an option. In this paper, we will fous on deriving three option Greeks Delta, Gamma, and Vega. Delta,, is defined as the rate of hange in the option value with respet to hanges in the underlying asset prie; Gamma, Γ is the rate of hange of with respet to hanges in the underlying prie; and finally, Vega is the measurement of an option s sensitivity to hanges in the volatility of the underlying asset prie. In general, volatility measures the amount and speed at whih the prie moves up and down, and it is often based on hanges in the reent, historial pries of a trading instrument. Other Greeks, suh as Theta, an be derived in a similar fashion; however, depending on the harateristi funtion, the derivation expression might be rather lengthy. We omit them here, as many terms are repeated. Delta is the first derivative of the value of V of the option with respet to the underlying instrument prie S. Hene, differentiating the CFS expansion of V 47 with respet to S, we have t = V x, K, t S = V x, K, t x x S = e rt t x Re i π d k B k Ĝ k e i π k x+log K d ]. 56 In a similar fashion, we an obtain Γ t by differentiating t with respet to S suh that Γ t = V x, K, t S = t S = t x x S, 57 and eventually, Γ t = e rt t x Re i π d k i π d k + B k Ĝ k e i π k x+log K d ]. 58 It is also easy to obtain the formula for Vega, V y t, where y t is the initial value of the volatility at time t. For example, for the Heston model, as y is the initial value of the volatility in 7, we derive Vega as follows: V x, K, y, t y rt t = e Re B k y Ĝ k e i π d k x+log K ], 59

13 April, 8 Quantitative Finane SingCFSEuroV3reformatted with B k y = ϕ π d k, y y, 6 where ϕ ontains the parameter y. To obtain our first SFP representation of, we first let y = x + log K, z = exp then transform all the jumps ζ in t into ε = exp the CFS representation in 56 into the form f z = U i π d y and i π d ζ in 56. Aordingly, this transforms i π d k B k Ĝ k z k. 6 Based on the equation above, using 6, we an eventually obtain the SFP approximant given by S P + N z + L + N s z log z/ε s = f zq + M z + OzU+. 6 s= Applying the approximation algorithm in Setion 5 to determine the oeffiients of P + N, Q+ M, and L + N s, we an obtain the SPF formula for t with the form P + e rt t x N Re z + S s= L+ N s z log z/ε s Q + M z. 63 To determine the SFP approximant of Γ t and Vega, we follow the same idea of approximating t but replae f z for U i π d k i π d k + B k Ĝ k z k and B k y Ĝ k z k Singular Fourier-Padé Algorithm and Loating Singularities The approah to omputing the polynomial oeffiients needed in the SFP method is fairly straightforward. To demonstrate the algorithm, we fous on a simple ase where the option priing and Greeks formulae are infinitely smooth apart from jumps loated at the endpoints and d. As we onsider z = exp i π d y in either the option priing formula or the Greeks formula, the jump of and d in the z-plane is -. For sake of simpliity, we denote f z as the CFS representation of any European-style priing formula or its option Greeks one. With some supersripts dropped for larity and knowing that s =, in 6, we have P N z + L N z log zε = f zq M z + Oz U+, 65

14 April, 8 Quantitative Finane SingCFSEuroV3reformatted Table.: Payoff funtions and their transforms for a variety of finanial ontingeny laims Finanial Contingeny Claim Payoff Funtion Transformed Payoff Funtion GS T, K Ge x+χ, K Call maxs T K, K maxe x+χ log K, Put maxk S T, K maxe x+χ log K, Covered Call mins T, K K mine x+χ log K, + K Cash-or-Nothing Call ST K e x+χ log K Cash-or-Nothing Put ST K e x+χ log K Asset-or-Nothing Call S T ST K e x+χ e x+χ log K Asset-or-Nothing Put S T ST K e x+χ e x+χ log K Asymmetri Call S n T Kn ST K K n e nx+χ log K e x+χ log K Asymmetri Put K n S n T S T K K n e nx+χ log K e x+χ log K Symmetri Call Symmetri Put S T K n ST K K S T n ST K n j= n j= n j n j e jx+χ+n j log K e x+χ log K n j n j e n jx+χ+j log K e x+χ log K where N + M + N = U. Both L N and f z have Taylor series and CFS expansions, respetively, to determine U; therefore, their expansions are log zεs = f z = U zk ε k + 66 U B k Ĝ k z k + B Ĝ. 67 Our goal is to derive a linear system for the unknown polynomial oeffiients. Note that Q M z and L N z are determined only by terms of order greater than N. Aordingly, we seek a linear solution to BĜ L ] ] q =. 68 l Here, BĜ is the M + N + M + Toeplitz matrix B U +Ĝ B U + U Ĝ U B U +Ĝ U + B U. B U Ĝ U B Ĝ +Ĝ... U + B Ĝ..... BU Ĝ U B U Ĝ U, 69

15 April, 8 Quantitative Finane SingCFSEuroV3reformatted Table.: Complex Fourier transforms for a variety of finanial ontingeny laims Finanial Contingeny Claim Fourier Transform Fourier Transform Call Put Kd iπk+d Kd iπk+d Covered Call Kd iπk+d Cash-or-Nothing Call Cash-or-Nothing Put Asset-or-Nothing Call π i e d k+d π i e d k+ π i e d k+ Asset-or-Nothing Put d iπk+d Asymmetri Call Asymmetri Put Symmetri Call Symmetri Put K n d iπk+nd K n d iπk+nd n j= n j= Ĝ k Ĝ Kd iπk Kd iπk + Kd iπk π i e d kd π i e d k π i e d k d iπk π i e d kd d iπk π i e d k d π i iπk+d e b a k+d π i e d k+ π i e d k+nd Kn d iπk π i e d kd π i e d k+ π i e d k n j n j K n d iπk+jd Kn d iπk π i e d k+jd n j n j+ K n d π i iπk+n jd e d k+n j n K n j Ke d Kd Ke K Ke + K d e d e K n n end K n d K n n en K n j= n n j n j ejd + n K n d j= n j n j+ K n n j en j K n and L is the M +N +N + matrix defined similarly using the Taylor oeffiients of log+z. The vetors q = {q m } M m= and l = {l n } N n = hold the unknown polynomial oeffiients in order of inreasing degree. As the olumn dimension of the matrix in 68 is one greater than its row dimension, we an onlude that there is one nonzero solution to 68. In many ases, this an be made into a square system by hoosing, say, q =. However, if one does not want to assume that any partiular oeffiient is nonzero, one an solve 68 by a singular value deomposition f. Gonnet et al. 3. Finally, the unknown oeffiients of p = {p n } N n= an be obtained by multipliation through the following matrix system: B Ĝ l B Ĝ B Ĝ p = B q l l l. 7 Ĝ l B U Ĝ U If there is more than one jump loation in the option priing/greeks urve 65, this suggests the following modifiation of the equation: P N z + L N z log zε L Ns z log z = f zq M z + Oz U+. 7 ε S Aordingly, we have to modify 68 to produe a new L matrix and a vetor of oeffiients for eah loation to reflet the hanges. Aording to Drisoll and Fornberg,, there is no l U

16 April, 8 Quantitative Finane SingCFSEuroV3reformatted 4 T= 5 T=. st differentiation of the denisty Density of the VG model y - - T= y Density of the VG model st differentiation of the denisty y # 4 T= y Figure.: Density funtions top of the VG model and their first derivative bottom. A pieewise ontinuous non-smooth PDF with a jump after/before differentiation is shown in the graphs on the top and bottom right. The parameters are taken from VG Para. rigorous optimal formula for hoosing the degrees M, N, and N,...,N s. Beause the denominator polynomial Q M is shared, we allow M to be the largest, with the others being equal as far as possible. For the ase of just one jump loation, taking N at roughly 4% of the total available degrees of freedom seems to work well. Experiments suggest that these hoies an affet the observed auray, oasionally by as muh as an order of magnitude, but on average, there is little variation within a broad range of hoies. The methods desribed above all assume that the loations of all jumps are known in the option priing/greeks urve. For jumps whose loations are unknown, we follow the suggestion of Drisoll and Fornberg, to use the Fourier-Padé algorithm to estimate their loations. As we approximate a PDF rather than a payoff funtion with the CFS method see the proof of Theorem, we only onsider jumps existing in the PDF. The existene of a jump is attributed to a ombination parameters and/or a very short maturity that auses a sharp-peaked PDF urve. Using the Fourier-Padé algorithm f. the Fourier-Padé approximant in Setion to estimate the jumps in a PDF is fairly simple. We first express the PDF as the CFS representation: Re Then, we an differentiate 7 with respet to y to obtain Re ϕ π ] d k e i π ky b a + ϕ. 7 i π b a k ϕ π ] d k e i π ky b a. 73 Finally, we let z = exp i π d y in the two equations above, and they are ready for the Fourier-Padé approximation. In general, when the PDF has a jump, the sharp-peaked jump point will have an enormously large value after differentiation. Then, based on this result, we find that point in the original PDF. In other words, Figure is a graphial illustration of the outlooks of two PDFs top under the VG model f. Setion 8.. and the first derivative bottom after the Fourier-Padé approximation. In the figure, we an see that the PDFs in the top and bottom panels are smooth without any jumps, as they do not show any point with an extremely large value. However, in the graph on the bottom right, we an see that there is a jump that produes a value of 5 4.

17 April, 8 Quantitative Finane SingCFSEuroV3reformatted 6. Choie of Trunated Intervals As we will show in Setion 7, the hoie of the interval, d] is affeted by the auray of the SFP method. A minimum and substantial interval, d] an apture most of the mass of a PDF suh that the SFP method an in turns produe sensible global spetral onvergene rate. In this short setion, we show how to onstrut an interval related to the losed-form formulas of stohasti proess umulants. The idea of using the umulants is first proposed by Fang and Oosterlee 9a to onstrut the definite interval, d] in 33. Based on their ideas, we have the following expression for, d]: d = + L log S 74 K = d, 75 where,, and 4 are the first, seond and fourth umulants, respetively, of the stohasti proess and L, ]. For simple, less-ompliated finanial models, we also obtain losed-form formulas for,, and 4, whih are shown in Table A of Appendix A. However, in the Heston model, we use the absolute value of and ignore the value of 4 due to the negative value of and the lengthy representation of 4 f. Fang and Oosterlee 9a. We therefore have + L rather than + L + 4. The trunated intervals only work for smooth PDFs without any jumps. In general, by trial and error, if there is a jump in a PDF, we add.5 to 74 to allow the better onvergene of our method. Aordingly, the formulae an be transformed as: d = + L log S K = d. 77 Remark Without a doubt, in most ases, the trunated interval of 76 works for any normal non-smooth PDF. However, when the onditions beome very extreme, for example, T = e 6 in BSM Para4 in Setion 8.., the PDF is extremely narrow, thin and spiky at one point. The interval of 76 an still perform but is not optimal. To improve it, we replae.5 in the equation with.. This is attributed to the desire for higher auray and onvergene with fewer terms required over a short interval range. 7. Error Analysis In this setion, we demonstrate that the total error from priing European-style options an be made very small by hoosing a suitably large interval, d]. Furthermore, we show that global spetral onvergene an be ahieved at any point along the option priing urve even though the input PDF is a C ν pieewise ontinuous non-smooth funtion. In this paper, there are three types of approximation errors in any all/put option. i Integration trunation error: ɛ := Ge x+χ, Kfχdχ Ge x+χ, Kfχdχ 78

18 April, 8 Quantitative Finane SingCFSEuroV3reformatted ii Error related to approximating 4 with 47: ɛ := Re d Re = Re d π i fye ky d dy Ĝk e i π k x+log K d + d B k Ĝ k e i π d k x+log K + B Ĝ ] d π i fye ky d dy B k Ĝ k e i π k x+log K d + fydy Ĝ fydy B Ĝ ] 79 ] iii Trunated SFP series error: ] P ɛ 3 := Re B k Ĝ k z k + B + N Ĝ Re z + ] S s= L+ N s z log z/ε s Q + M z, z = e i π d x+log K and ε s = e i π d ζs 8 If we introdue the onept of the umulative probability density funtion CDF F χ suh that fχdχ = df χ, we an simplify the integration trunation error as follows: ɛ = Ge x+χ, Kfχdχ Ge x+χ, Kfχdχ = Ge x+χ, Kfχdχ + Ge fχdχ x+χ d Ge x+χ, K F χdχ + Ge x+χ, K χ F χdχ χ 8 d : if χ =, d,,. 8 We an see that ɛ is bounded and approahes zero as long as, d] is hosen reasonably suh that F d when d < or F when >. We are also able to adapt the same idea to investigate the bound of ɛ. Aordingly, taking into aount expi πk d y, we first investigate the error ɛ := d π i fye ky d dy B k

19 April, 8 Quantitative Finane SingCFSEuroV3reformatted in ɛ. If we expand the equation above, we obtain ɛ := d π i fye ky d dy d d ϕ i π d k 83 = π i fye ky d dy π i fye ky d dy d d 84 fydy + fydy 85 d d = F F d + F F d 86 : if y =, d,,. 87 Based on the result above, ɛ := Re d Re. d π i fye ky d dy B k Ĝ k e i π k x+log K d + ] fydy B Ĝ ɛ Ĝ k e i π d k x+log K + d fydy d ] fydy Ĝ 88 To onlude that ɛ is approahing zero, we first note that there is no approximate error of G k and G beause of their losed-form expressions. Then, one ɛ tends to zero, the first term of the equation will also diminish to zero. The last term of the equation tends to zero beause B equals d d fydy f. 45 and 46. Finally, the SPF series trunation error is also bounded f. Drisoll and Fornberg, and an be formulated as follows: ɛ 3 := Re = Oz N+M+ Ns s= Ns+. B k Ĝ k e i π d k x+log K + B Ĝ ] Re P + N z + ] S s= L+ N s z log z/ε s Q + M z 89 The error term of ɛ 3 tends to zero with a global spetral rate of Oz N+M+ Ns s= Ns+ even though the input PDF is C ν pieewise ontinuous. Before we illustrate the total error bound when approximating any true European-type option prie V x, K, t defined as P + ɛ := V x, K, t e rt t N Re z + ] S s= L+ N s z log z/ε s Q + M z, 9 we first summarise the whole approximation proedure of European-type option pries and note where ɛ, ɛ, and ɛ 3 lie. We start off by seeking a definite interval, d] that allows us to approximate

20 April, 8 Quantitative Finane SingCFSEuroV3reformatted V x, K, t defined on, ] in 3 with the form V x, K, t e rt t Ge x+χ, Kfχdχ. The interval, d] we proposed satisfies ondition 33. As a result, we obtain our first approximation error ɛ. As V x, K, t is now approximated in, d], this implies that we an onstrut a CFS expansion of V x, K, t, like the one in 4. Then, beause inluding a harateristi funtion ϕ in the CFS expansion allows for a more aurate approximation, we have another CFS expansion of V x, K, t, given as in 47. Aordingly, we have ɛ, an approximation error of 4, being approximated by 47. Finally, ɛ 3 is the error of 5, whih is the formula for approximating 47 with the SFP approximant 6. By ombining the results of ɛ, ɛ and ɛ 3, we an determine the total error bound ɛ; hene, we have an inequality of P + ɛ = V x, K, t e rt t N Re z + ] S s= L+ N s z log z/ε s Q + M z t P + = e rt Ge x+χ N, Kfχdχ Re z + ] S s= L+ N s z log z/ε s Q + M z e rt t Ge x+χ, Kfχdχ Ge x+χ, Kfχdχ + Re Re d π i fye ky d dy Ĝk e i π k x+log K d + d B k Ĝ k e i π d k x+log K + B Ĝ ] + fydy Ĝ ] P Re B k Ĝ k z k + B + N Ĝ Re z + ] S s= L+ N s z log z/ε s Q + M z e rt t ɛ + ɛ + ɛ 3 < e rt t ɛ + ɛ + Oz N+M+ S s= Ns+. ] 9 9 Remark 3 Aording to Drisoll and Fornberg,, the rate of Oz N+M+ Ns Ns+ s= is not spatially uniform beause onvergene at a jump is somewhat limited by the well-known numerial ill-onditioning of the straightforward Padé problem. Moreover, when a urve has a disontinuous point zeroth-order jump, it also has a onvergene rate that is the average of the one-sided limit values. Through the numerial experiments in Drisoll and Fornberg, we an see that the SFP method an yield spetral onvergene rate at jumps. This is not neessarily the ase if the jumps are very diffiult to interpolate. Nevertheless, as Drisoll and Fornberg suggest, we an still have 4 6 digits of auray at the jumps when they are known in advane. When the jumps are not known in advane, we an use the Fourier-Padé algorithm to find them see Setion 5. By inorporating this tehnique, the SFP method still an perform and yield spetral onvergene or 4th-order onvergene when the jumps are very diffiult to interpolate.

21 April, 8 Quantitative Finane SingCFSEuroV3reformatted 8. Numerial Results In this setion, we demonstrate the performane of the SFP method through various numerial tests. The purpose of this setion is first to test whether the error onvergene analysis presented in Setion 7 is in line with the numerial findings in this setion. Seond, we test the ability of the SFP method to prie any European-style option that are deep in/out of the money and have long/short maturities. Third, we analyse whether the SFP method an provide onsistent auray when approximating small or large values of option pries. Finally, in a major development of the method, we test whether the SFP method an retain global spetral onvergene even when the PDF is pieewise ontinuous. A number of popular numerial methods are implemented to test the SFP method in terms of the error onvergene, onvergene rate and omputational time. These methods inlude the COS method a Fourier COS series method, Fang and Oosterlee 9a, the filter-cos method a COS method with an exponential filter to resolve the Gibbs phenomenon; see Ruijter et al. 3, the CONV method an FFT method, Lord et al. 8, the Lewis-FRFT method a frational FFT method, Lewis, Chourdakis 4, and the B-spline, Haar wavelet and the SWIFT methods a wavelet-based method; see Ortiz-Graia and Oosterlee 3, 6. When we implement the CONV and Lewis-FRFT methods, we use Simpson s rule for the Fourier integrals to ahieve fourth-order auray. In the filter-cos method, we use an exponential filter and set the auray parameter to as Ruijter et al. 3 report that this filter provides better algebrai onvergene than the other options. We also set the damping fators of the CONV and Lewis-FRFT to and any value greater than zero, respetively. As the SFP method requests approximating jumps in logarithmi series, we onsider and apply the endpoints and d as our two known jumps for all non-smooth/smooth PDFs. Only the jump of the non-smooth PDF of the Blak-Sholes-Merton BSM model is known as its mean value. For the rest of the non-smooth PDFs, we use the Fourier-Padé algorithm f. Setion 5 to find their loations. In all numerial experiments, we use the parameter U to denote the number of terms of the SFP method and N to denote the number of terms/grid points of the others. When we measure the approximation errors of the numerial methods, we use absolute errors, the infinity norm errors R and the L norm errors R as the measurement units. Moreover, to improve the auray of our method, we use the all-put parity V all x, K, t = V put x, K, t + S exp qt K exp rt to approximate all pries one we have put pries ready. Finally, all the CPU times presented in seonds are determined after averaging the omputational time over experiments. A MaBook Pro with a.8 GHz Intel Core i7 CPU and two 8 GB DDR SDRAM ahe memory is used for all experiments. The ode is written in MATLAB Rb. Finally, the MATLAB ode of implementing the COS method and the FFT method, suh as the CONV method and the like, is retrieved from von Sydow et al Exponential Lévy Proesses 8... The Blak-Sholes-Merton Model. The first numerial experiments are performed using the BSM model f. Blak and Sholes 973, Merton 973. The stok dynamis driven by the BSM model a geometri Brownian proess are given by S T = S e r q σ T +σw T, 93 where W T is a risk-neutral Brownian motion, and σ is the volatility. The harateristi funtion of the model is also defined as ϕu = exp T iur q σ σ u, u R. 94

22 April, 8 Quantitative Finane SingCFSEuroV3reformatted The parameters for the experiments are seleted from the following: BSM Para : S =, σ =.5, r =.3, T =., q =., 95 BSM Para : S =, σ =.5, r =., T = 5 or, q =, K =, 96 BSM Para3 : K =, σ =., r =.6, T = e 6, q =. 97 In the first numerial test BSM Para, we first hek for onvergene behaviour against a range of strikes K from to for deep in/out-of-the money and at-the-money vanilla put options, respetively. The parameters are retrieved from von Sydow et al. 5. Using the same parameters, we hek for onvergene behaviour against a range of strikes K from 8 to for ash-or-nothing put options rather than vanilla options in our seond numerial test. We ompare our method with the COS method, the CONV method and the Lewis-FRFT method in the first test but only with the COS method in the seond test. The third numerial experiment BSM Para is devoted to omparing the performane of the COS method, the SFP method, the SWIFT method for long maturity all options. As we sometimes enounter these options in the insurane and pension industry, it is worth testing our method against them. The parameters are retrieved from Ortiz-Graia and Oosterlee 6 for the test. Finally, for the last numerial test BSM Para3, we use the values of σ, r and q from Andriopoulos et al. 3 and hek for onvergene behaviour against a range of stok pries S from 8 to for very short maturity vanilla all options and their option Delta and Gamma. Our method is ompared with the COS, filter-cos, Lewis-FRFT and CONV methods. Referene values for all the tests are based on the BSM analytial formula. In eah numerial test, exept for the third one, we delare 5 different option pries within the range of either K or S to test the effiieny of our method and the others. We represent all reovered PDFs via the SFP method in Figure for all sets of parameters. We an see that the non-smooth reovered PDF top right an be obtained via BSM Para3. We set L = in the trunated interval proposed in Setion 6. One should note that for the trunated interval in the last numerial test BSM Para3, we replae.5 with. in 76 see Remark, as we want our method to have higher auray and a better onvergene rate. It is quite rare for the maturity time to be infinitesimal as in BSM Para3, but to demonstrate the effiieny of our method, it is worth performing suh a test. All the tables in this setion and the others suggest that the differene in the omputational time aross methods is not large. It takes less than. seonds to approximate 5 option pries for any method when N and U equal 64. Without onsidering the auray of the methods, this is a quite reasonable time frame in whih to produe a substantial number of option pries at one. However, when we onsider the error onvergene and onvergene rate of the methods, there are sizeable differenes among them. Figure 3 and its numerial presentation in Table 3 are dediated to the first numerial test BSM Para. In this test, the SFP method has the smallest R and R errors, with R and R less than 5.8e 3 as U = 64 when the strike inreases from to. It also yields global spetral onvergene against whih none of the other methods an ompete. This suggests that the trunated range proposed in this paper an work over longer range than that proposed in Fang and Oosterlee 9a. In Table 4, we use the same parameters BSM Para to examine the ability of the COS and SFP methods to prie ash-or-nothing put options. Again, the two methods an ahieve very low onvergene error and spetral onvergene. Table 5 refers to the third test BSM Para and repliates Table 3 in Ortiz-Graia and Oosterlee 6. In this test, the SFP method impressively provides high auray over the SWIFT and COS methods with fewer terms required. For example, when T = 5, we only need 3 terms to obtain the true solution with 7-digit auray in the SFP method, whereas the other two require more terms to obtain -digit auray. Finally, in the last test BSM Para4, as the umulants are too small to reate a meaningful trunated range for the COS and filter-cos methods, we use 76 to allow the method to generate relevant option pries. In Figure 4, the SFP method dominates the other four methods to yield the global spetral onvergene rate away from the jump, while the input

23 April, 8 Quantitative Finane SingCFSEuroV3reformatted Density of the GBM model y Density of the GBM model y Density of the GBM model T=5 T= y Figure.: Probability density funtions are generated by three sets of parameters: BSM Para top left, BSM Para3 top right and BSM Para bottom. CONV Lewis-FRFT log jerrorj - -4 log jerrorj log jerrorj Strike COS log jerrorj Strike SFP Strike Strike Figure 3.: Error envelopes of the CONV, Lewis-FRFT, COS and SFP methods for priing European deep at/in/out-of-the-money put options under the BSM model with strikes inreasing from to. The number of terms/grid points inreases in a sequene of 8 red line, 6 blue line, 3 green line and 64 blak line. Spetral onvergene is observed in the SFP method. The parameters are taken from BSM Para. non-smooth PDF see top right, Figure behaves like a Dira Delta funtion and has an almost infinitely thin spike near the origin. In Table 6, when we test the auray of the option pries around/at the jump S = the SFP method dominate the other methods. It has zero absolute error in the ase of S = 95 and 6.69e 5 error from the true solution in the ase of S = when U = 64. This is in line with the findings of Drisoll and Fornberg as the SFP offers 4-6 digit auray at the jumps when the funtion is very diffiult. Obviously, the non-smooth PDF in the last test is very diffiult to approximate, as it behaves like a Dira Delta. Finally, using the same parameters of BSM Para3, we reover the all Delta and Gamma via the SFP approximant in Figure 5. In the graphs, the SFP provides a solution of extremely high auray and global spetral onvergene apart from the jump The Variane Gamma Model. A VG proess Madan and Seneta 99, Madan and Milne 99, Madan et al. 998 is an infinite ativity Lévy proess and is a subordinate version of Brownian motion f. Cont and Tankov 4. The most important feature of this model is that the Brownian motion is evaluated in random time t determined by an independent inreasing Lévy The Dira Delta funtion is a generalised funtion or distribution on the real number line that is zero everywhere exept at zero.