# A Generic Decomposition Formula for Pricing Vanilla Options under Stochastic Volatility Models

Authors: Raúl Merino; Josep Vives
Journal: International Journal of Stochastic Analysis
 (2015)

---

## Abstract

We obtain a decomposition of the call option price for a very general stochastic volatility diffusion model, extending a previous decomposition formula for the Heston model. We realize that a new term arises when the stock price does not follow an exponential model. The techniques used for this purpose are nonanticipative. In particular, we also see that equivalent results can be obtained by using Functional Itô Calculus. Using the same generalizing ideas, we also extend to nonexponential models the alternative call option price decomposition formula written in terms of the Malliavin derivative of the volatility process. Finally, we give a general expression for the derivative of the implied volatility under both the anticipative and the nonanticipative cases.


---

## Body

## 1. Introduction

Stochastic volatility models are a natural extension of the Black-Scholes model in order to manage the skew and the smile observed in real data. It is well known that in these models the average of future volatilities is a relevant quantity. See, for example, [1, Chapter 2]. Unfortunately, adding a stochastic volatility structure makes pricing and calibration more complicated, as closed formulas do not always exist. See, for example, [2], for a first reference on this topic. Moreover, even when these formulas exist, like for the Heston model (see [3]), in general, they do not allow a fast calibration of the parameters.During last years, different developments for finding approximations to the closed-form option pricing formulas have been published. Malliavin techniques are naturally used to solve this problem in [4, 5] as the average future volatility is an anticipative quantity. Otherwise, a nonanticipative method to obtain an approximation of the pricing formula is developed for the Heston model in [6]. The method is based on the use of the adapted projection of the average future volatility. As a result, the model allows obtaining a decomposition of the call option price in terms of such future volatility.In the present paper, we generalize [6] to general stochastic volatility diffusion models. Similarly, following the same kind of ideas, we extend the expansion based on Malliavin calculus obtained in [4, 5]. This is important because Heston model is not the unique stochastic volatility model currently used in practice, and some of them, like SABR model, are not of exponential type. For a general discussion about stochastic volatility models in practice, see [7].The main ideas developed in this paper are the following:(i)
A generic call option price decomposition is found without having to specify the volatility structure.
(ii)
A new term emerges when the stock option prices do not follow an exponential model, as, for example, in the SABR case.
(iii)
The Feynman-Kac formula is a key element in the decomposition. It allows expressing the new terms that emerge under the new framework (i.e., stochastic volatility) as corrections of the Black-Scholes formula.
(iv)
The decomposition found using Functional Itô calculus appears to be the same as the decomposition obtained through our techniques.
(v)
A general expression of the derivative of the implied volatility, both for nonanticipative and anticipative cases, is given.


## 2. Notations

LetS
=
{
S
(
t
)
,
t
∈
[
0
,
T
]
} be a strictly positive price process under a market chosen risk neutral probability that follows the model(1)
d
S
t
=
μ
t
,
S
t
d
t
+
θ
t
,
S
t
,
σ
t
ρ
d
W
t
+
1
-
ρ
2
d
B
t
,where W and B are independent Brownian motions, ρ
∈
(
-
1,1
), μ
:
[
0
,
T
]
×
R
+
→
R, θ
:
[
0
,
T
]
×
R
+
2
→
R
+, and σ
(
t
) is a positive square-integrable process adapted to the filtration of W. We assume on μ and σ sufficient conditions to ensure the existence and uniqueness of the solution of (1). Notice that we do not assume any concrete volatility structure. Thus, our decompositions can be adapted to many different models. In particular, we cover the following models:(i)
Black-Scholes model:μ
(
t
,
S
(
t
)
)
≔
r
S
(
t
), θ
(
t
,
S
(
t
)
,
σ
(
t
)
)
≔
σ
S
(
t
), ρ
=
0, r
>
0, and σ
>
0.
(ii)
CEV model:μ
(
t
,
S
(
t
)
)
≔
r
S
(
t
), θ
(
t
,
S
(
t
)
,
σ
(
t
)
)
≔
σ
S
(
t
)
β with β
∈
(
0,1
], ρ
=
0, r
>
0, and σ
>
0.
(iii)
Heston model:μ
(
t
,
S
(
t
)
)
≔
r
S
(
t
), θ
(
t
,
S
(
t
)
,
σ
(
t
)
)
≔
σ
(
t
)
S
(
t
), r
>
0, σ
>
0, and(2)
d
σ
2
t
=
k
θ
-
σ
2
t
d
t
+
ν
σ
2
t
d
W
t
,where k, θ, and ν are positive constants satisfying the Feller condition 2
k
θ
>
ν
2.
(iv)
SABR model:μ
(
t
,
S
(
t
)
)
≔
r
S
(
t
), θ
(
t
,
S
(
t
)
,
σ
(
t
)
)
≔
σ
(
t
)
S
(
t
)
β with β
∈
(
0,1
], r
>
0, σ
>
0, and (3)
d
σ
t
=
α
σ
t
d
W
twith α
>
0
.For existence and unicity of the solution in the Heston case, see, for example, [8, Section 2.2]. For the CEV and SABR models, see [9] and the references therein.The following notation will be used in all the paper:(i)
We will denote byB
S
(
t
,
S
,
σ
) the price of a plain vanilla European call option under the classical Black-Scholes model with constant volatility σ, current stock price S, time to maturity τ
=
T
-
t, strike price K, and interest rate r. In this case, (4)
B
S
t
,
S
,
σ
=
S
Φ
d
+
-
K
e
-
r
τ
Φ
d
-
,where Φ
(
·
) denotes the cumulative probability function of the standard normal law and (5)
d
±
=
l
n
S
/
K
+
r
±
σ
2
/
2
τ
σ
τ
.
(ii)
We use in all the paper the notationE
t
[
·
]
≔
E
[
·
∣
F
t
], where F
t
,
t
≥
0 is the natural filtration of S.
(iii)
In our setting, the call option price is given by(6)
V
t
=
e
-
r
τ
E
t
S
T
-
K
+
.
(iv)
Recalling that from the Feynman-Kac formula, the operator(7)
L
θ
≔
∂
t
+
1
2
θ
t
,
S
t
,
σ
t
2
∂
S
2
+
μ
t
,
S
t
∂
S
-
rsatisfies L
θ
B
S
(
t
,
S
(
t
)
,
θ
(
t
,
S
(
t
)
,
σ
(
t
)
)
)
=
0.
(v)
We will also use the following definitions fory
≥
0:(8)
G
t
,
S
t
,
y
≔
S
2
t
∂
S
2
B
S
t
,
S
t
,
y
,
H
t
,
S
t
,
y
≔
S
t
∂
S
G
t
,
S
t
,
y
,
K
t
,
S
t
,
y
≔
S
2
t
∂
S
2
G
t
,
S
t
,
y
,
L
t
,
S
t
,
y
≔
θ
t
,
S
t
,
y
S
t
.


## 3. A Decomposition Formula Using Itô Calculus

In this section, following the ideas in [6], we extend the decomposition formula to a generic stochastic volatility diffusion process. We note that the new formula can be extended without having to specify the underlying volatility process, obtaining a more flexible decomposition formula. When the stock price does not follow an exponential process, a new term emerges. The formula proved in [6] is a particular case.It is well known that if the stochastic volatility process is independent from the price process, then the pricing formula of a plain vanilla European call is given by(9)
V
t
=
E
t
B
S
t
,
S
t
,
σ
-
t
,where σ
-
2
(
t
) is the so-called average future variance and it is defined by(10)
σ
-
2
t
≔
1
T
-
t
∫
t
T
σ
2
s
d
s
.Naturally, σ
-
(
t
) is called the average future volatility. See [1, page 51].The idea used in [6] consists in considering the adapted projection of the average future variance(11)
v
2
t
≔
E
t
σ
-
2
t
=
1
T
-
t
∫
t
T
E
t
σ
2
s
d
sto obtain a decomposition of V
(
t
) in terms of v
(
t
). This idea switches an anticipative problem related with the anticipative process σ
-
(
t
) into a nonanticipative one with the adapted process v
(
t
). We apply this technique to our generic stochastic differential equation (1).Theorem 1 (decomposition formula).
For allt
∈
[
0
,
T
), we have (12)
V
t
=
B
S
t
,
S
t
,
v
t
+
1
2
·
E
t
∫
t
T
e
-
r
u
-
t
G
u
,
S
u
,
v
u
·
L
2
u
,
S
u
,
σ
u
-
σ
2
u
d
u
+
1
8
·
E
t
∫
t
T
e
-
r
u
-
t
K
u
,
S
u
,
v
u
d
M
,
M
u
+
ρ
2
E
t
∫
t
T
e
-
r
u
-
t
L
u
,
S
u
,
σ
u
·
H
u
,
S
u
,
v
u
d
W
,
M
u
,where M
(
t
)
≔
∫
0
T
‍
E
t
σ
2
(
s
)
d
s
=
∫
0
t
‍
σ
2
(
s
)
d
s
+
(
T
-
t
)
v
(
t
)
2.Proof.
Notice thate
-
r
T
B
S
(
T
,
S
(
T
)
,
v
(
T
)
)
=
e
-
r
T
V
(
T
). As e
-
r
t
V
(
t
) is a martingale, we can write (13)
e
-
r
t
V
t
=
E
t
e
-
r
T
V
T
=
E
t
e
-
r
T
B
S
T
,
S
T
,
v
T
.Our idea is to apply the Itô formula to the process e
-
r
t
B
S
(
t
,
S
(
t
)
,
v
(
t
)
)
.
As the derivatives ofB
S are not bounded, we have to use an approximation to the identity changing B
S
(
t
,
S
,
σ
) by (14)
B
S
n
t
,
S
,
σ
≔
B
S
t
,
S
,
σ
ψ
n
S
,where ψ
n
(
S
)
=
ϕ
(
1
/
n
S
) for some ϕ
∈
C
b
2 such that ϕ
(
S
)
=
1 for all |
S
|
<
1 and ϕ
(
S
)
=
0 for all |
S
|
>
2, and v
(
t
) by v
ε
(
t
)
=
1
/
(
T
-
t
)
δ
+
∫
t
T
‍
E
σ
2
(
s
)
d
s, where ε
>
0, and finally apply the dominated convergence theorem. For simplicity, we skip this mollifying argument across the paper.
So, applying Itô formula, using the fact that(15)
∂
σ
B
S
t
,
S
,
σ
=
S
2
σ
τ
∂
S
2
B
S
t
,
S
,
σand the Feynman-Kac operator (7), we deduce (16)
e
-
r
T
B
S
T
,
S
T
,
v
T
-
e
-
r
t
B
S
t
,
S
t
,
v
t
=
∫
t
T
e
-
r
u
L
v
S
B
S
u
,
S
u
,
v
u
d
u
+
∫
t
T
e
-
r
u
∂
S
B
S
u
,
S
u
,
v
u
θ
u
,
S
u
,
σ
u
·
ρ
d
W
u
+
1
-
ρ
2
d
B
u
+
1
2
∫
t
T
e
-
r
u
S
2
u
·
∂
S
2
B
S
u
,
S
u
,
v
u
d
M
u
+
1
2
∫
t
T
e
-
r
u
S
2
u
·
∂
S
2
B
S
u
,
S
u
,
v
u
L
2
u
,
S
u
,
σ
u
d
u
-
σ
2
u
d
u
+
1
8
∫
t
T
e
-
r
u
S
2
u
·
∂
S
2
S
2
u
∂
S
2
B
S
u
,
S
u
,
v
u
d
M
,
M
u
+
ρ
2
∫
t
T
e
-
r
u
θ
u
,
S
u
,
σ
u
·
∂
S
S
2
u
∂
S
2
B
S
u
,
S
u
,
v
u
d
W
,
M
u
.
Taking conditional expectation and multiplying bye
r
t, we have(17)
E
t
e
-
r
T
-
t
B
S
T
,
S
T
,
v
T
=
B
S
t
,
S
t
,
v
t
+
1
2
E
t
∫
t
T
e
-
r
u
-
t
S
2
u
∂
S
2
B
S
u
,
S
u
,
v
u
·
L
2
u
,
S
u
,
σ
u
d
u
-
σ
2
u
d
u
+
1
8
·
E
t
∫
t
T
e
-
r
u
-
t
S
2
u
·
∂
S
2
S
2
u
∂
S
2
B
S
u
,
S
u
,
v
u
d
M
,
M
u
+
ρ
2
E
t
∫
t
T
e
-
r
u
-
t
θ
u
,
S
u
,
σ
u
∂
S
S
2
u
·
∂
S
2
B
S
u
,
S
u
,
v
u
d
W
,
M
u
.Remark 2.
In [6], the following operators are defined for X
(
t
)
=
log
⁡
S
(
t
):(i)
G
~
(
t
,
X
(
t
)
,
σ
(
t
)
)
:
=
∂
x
2
-
∂
x
B
S
(
t
,
X
(
t
)
,
σ
(
t
)
)
.
(ii)
H
~
(
t
,
X
(
t
)
,
σ
(
t
)
)
:
=
∂
x
3
-
∂
x
2
B
S
(
t
,
X
(
t
)
,
σ
(
t
)
)
.
(iii)
K
~
(
t
,
X
(
t
)
,
σ
(
t
)
)
:
=
∂
x
4
-
2
∂
x
3
+
∂
x
2
B
S
(
t
,
X
(
t
)
,
σ
(
t
)
)
.
We observe the following:(i)
G
~
(
t
,
X
(
t
)
,
σ
(
t
)
)
=
G
(
t
,
S
(
t
)
,
σ
(
t
)
)
.
(ii)
K
~
(
t
,
X
(
t
)
,
σ
(
t
)
)
=
K
(
t
,
S
(
t
)
,
σ
(
t
)
)
.
(iii)
H
~
(
t
,
X
(
t
)
,
σ
(
t
)
)
=
H
(
t
,
S
(
t
)
,
σ
(
t
)
).Remark 3.
We have extended the decomposition formula in [6] to the generic SDE (1). When we apply Itô Calculus, we realize that Feynman-Kac formula absorbs some of the terms that emerge. It is important to note that this technique works for any payoff or any diffusion model satisfying Feynman-Kac formula.Remark 4.
Note that whenθ
(
t
,
S
(
T
)
,
σ
(
t
)
)
=
σ
(
t
)
S
(
t
) (i.e., the stock price follows an exponential process), then (18)
V
t
=
B
S
t
,
S
t
,
v
t
+
1
8
·
E
t
∫
t
T
e
-
r
u
-
t
K
u
,
S
u
,
v
u
d
M
,
M
u
+
ρ
2
E
t
∫
t
T
e
-
r
u
-
t
σ
u
H
u
,
S
u
,
v
u
d
W
,
M
·
u
,and the term(19)
1
2
E
t
∫
t
T
e
-
r
u
-
t
S
2
u
∂
S
2
B
S
u
,
S
u
,
v
u
·
L
2
u
,
S
u
,
σ
u
-
σ
2
u
d
uvanishes.
Indeed, we show that, due to the use of Feynman-Kac formula,(20)
Movement  of  the  asset
+
Movement  of  the  volatility
=
1
2
θ
u
,
S
u
2
·
∂
S
2
B
S
u
,
S
u
,
v
u
d
u
+
∂
σ
B
S
u
,
S
u
,
v
u
d
v
u
=
1
2
σ
2
u
S
2
u
·
∂
S
2
B
S
u
,
S
u
,
v
u
d
u
+
1
2
S
2
u
·
∂
S
2
B
S
u
,
S
u
,
v
u
d
M
+
v
2
d
u
-
σ
2
d
u
=
1
2
·
S
2
u
∂
S
2
B
S
u
,
S
u
,
v
u
d
M
+
v
2
d
u
,where 1
/
2
S
2
(
u
)
∂
S
2
B
S
(
u
,
S
(
u
)
,
v
(
u
)
)
v
2 is used into the Feynman-Kac formula and (21)
E
t
∫
t
T
1
2
S
2
∂
S
2
B
S
u
,
S
u
,
v
u
d
M
=
0
.


## 4. Basic Elements of Functional Itô Calculus

In this section, we give the insights of the Functional Itô Calculus developed in [10–13].LetX
:
[
0
,
T
]
×
Ω
↦
R be an Itô process, that is, a continuous semimartingale defined on a filtered probability space (
Ω
,
F
,
(
F
t
)
t
∈
[
0
,
T
]
,
P
) which admits the stochastic integral representation (22)
X
t
=
x
0
+
∫
0
t
μ
u
d
u
+
∫
0
t
σ
u
d
W
u
,where W is a Brownian motion and μ
(
t
) and σ
(
t
) are continuous processes, respectively, in L
1
(
Ω
×
[
0
,
T
]
) and L
2
(
Ω
×
[
0
,
T
]
).We defineD
(
[
0
,
T
]
,
R
) the space of cadlag functions. Given a path x
∈
D
(
[
0
,
T
]
,
R
), we will denote by x
t its restriction to [
0
,
t
]. For h
≥
0, the horizontal extension x
t
,
h is defined as (23)
x
t
,
h
u
=
x
t
u
=
x
u
,
u
∈
0
,
t
;
x
t
,
h
u
=
x
t
,
u
∈
t
,
t
+
hand the vertical extension is defined as (24)
x
t
h
u
=
x
t
u
=
x
u
,
u
∈
0
,
t
;
x
t
h
t
=
x
t
+
h
,
that  is
x
t
h
u
=
x
u
+
h
1
t
=
u
.A process Y
:
[
0
,
T
]
×
R
→
R, progressively measurable with respect to the natural filtration of X, may be represented as (25)
Y
t
=
F
t
,
X
s
,
0
≤
s
≤
t
=
F
t
,
X
tfor a certain measurable functional F
:
[
0
,
t
]
×
D
(
[
0
,
t
]
,
R
)
→
R. Let F
∞ be the space of local Lipschitz functionals with respect to the norm of the supremum on D
(
[
0
,
t
+
h
]
,
R
); that is, there exists a constant C
>
0 such that for any compact K and for any x
∈
D
(
[
0
,
t
]
,
K
) and y
∈
D
(
[
0
,
t
+
h
]
,
K
) we have(26)
F
t
,
x
t
-
F
t
+
h
,
y
t
+
h
≤
C
x
t
,
h
-
y
t
+
h
∞
.Under this framework, we have the next definitions of derivative.Definition 5 (horizontal derivative).
The horizontal derivative of a functionalF
∈
F
∞ at t is defined as (27)
D
t
F
t
,
x
t
=
l
i
m
h
→
0
+
F
t
+
h
,
x
t
,
h
-
F
t
,
x
t
h
.Definition 6 (vertical derivative).
The vertical derivative of a functionalF
∈
F
∞ at t is defined as (28)
∇
x
F
t
,
x
t
=
l
i
m
h
→
0
+
F
t
,
x
t
h
-
F
t
,
x
t
h
.Of course we can consider iterated derivatives as ∇
x
x.We also have the following Itô formula that works for nonanticipative functionals:Theorem 7 (Functional Itô Formula).
For any nonanticipative functionalF
∈
F
∞ and any t
∈
[
0
,
T
], we have(29)
F
t
,
X
t
-
F
0
,
X
0
=
∫
0
t
D
u
F
u
,
X
u
d
u
+
∫
0
t
∇
x
F
u
,
X
u
d
X
u
+
1
2
∫
0
t
∇
x
x
F
u
,
X
u
d
X
,
X
u
,provided that D
t
F, ∇
x
F, and ∇
x
x
F belong to F
∞.Proof.
See [11, 12].


## 5. A General Decomposition Using Functional Itô Calculus

In this section, we apply Functional Itô Calculus to the problem of finding a decomposition for the call option price. The decomposition problem is an anticipative path-dependent problem. Using a smart choice of the volatility process into the Black-Scholes formula, we can convert it into a nonanticipative one. It is natural to wonder whether the Functional Itô Calculus brings some new insides into the problem.We consider the functional(30)
F
t
,
S
t
,
σ
t
2
=
e
-
r
t
B
S
t
,
S
t
,
f
t
,
σ
t
2
,where σ
2 is the path-dependent process and f
∈
F
∞ is a nonanticipative functional.Under this framework, we calculate the derivatives using Functional Itô Calculus with respect to the variance. Then we write them in terms of the classical Black-Scholes derivatives. We must realize that, for simplicity, the new derivatives are calculated with respect to the variance instead of the volatility of the process.Remark 8.
If∂ denotes the classical derivative, we have(i)
Alternative Vega:(31)
∇
σ
2
F
=
e
-
r
t
∂
f
B
S
t
,
S
t
,
f
t
,
σ
t
2
∇
σ
2
f
t
,
σ
t
2
.
(ii)
Alternative Vanna:(32)
∇
σ
2
,
S
F
=
∂
f
,
S
B
S
t
,
S
t
,
f
t
,
σ
t
2
∇
σ
2
f
t
,
σ
t
2
.
(iii)
Alternative Vomma:(33)
∇
σ
2
,
σ
2
F
=
e
-
r
t
∂
f
,
f
B
S
t
,
S
t
,
f
t
,
σ
t
2
∇
σ
2
f
t
,
σ
t
2
2
-
e
-
r
t
∂
f
B
S
t
,
S
t
,
f
t
,
σ
t
2
∇
σ
2
2
f
t
,
σ
t
2
.
(iv)
Alternative Theta:(34)
D
t
F
=
-
r
e
-
r
t
B
S
t
,
S
t
,
f
t
,
σ
t
2
+
e
-
r
t
∂
t
B
S
t
,
S
t
,
f
t
,
σ
t
2
+
e
-
r
t
∂
f
B
S
t
,
S
t
,
f
t
,
σ
t
2
D
t
f
t
,
σ
t
2
.Theorem 9 (decomposition formula).
For allt
∈
[
0
,
T
), S
(
t
), and f
(
t
,
σ
t
2
)
>
0, we have (35)
V
t
=
B
S
t
,
S
t
,
f
u
,
σ
t
2
+
E
t
∫
t
T
e
-
r
u
-
t
f
u
,
σ
u
2
τ
G
u
,
S
u
,
f
u
,
σ
u
2
·
D
u
f
u
,
σ
u
2
d
u
+
1
2
·
E
t
∫
t
T
e
-
r
u
-
t
G
u
,
S
u
,
f
u
,
σ
u
2
·
L
2
u
,
S
u
,
σ
u
-
f
2
u
,
σ
u
2
d
u
+
1
2
·
E
t
∫
t
T
e
-
r
u
-
t
f
u
,
σ
u
2
2
·
τ
2
K
u
,
S
u
,
f
u
,
σ
u
2
d
f
u
,
σ
u
2
,
f
u
,
σ
u
2
+
ρ
E
t
∫
t
T
e
-
r
u
-
t
L
2
u
,
S
u
,
σ
u
·
H
u
,
S
u
,
f
u
,
σ
u
2
f
u
,
σ
u
2
·
τ
d
W
u
,
f
t
,
σ
u
2
.Proof.
Notice thatF
(
T
,
X
(
T
)
,
σ
T
2
)
=
e
-
r
T
B
S
(
T
,
S
(
T
)
,
f
(
T
,
σ
T
2
)
)
=
e
-
r
T
V
T. As e
-
r
t
V
(
t
) is a martingale, we can write (36)
e
-
r
t
V
t
=
E
t
e
-
r
T
V
T
=
E
t
e
-
r
T
B
S
T
,
S
T
,
f
T
,
σ
T
2
=
E
t
F
T
,
S
T
,
σ
T
2
.Our idea is to apply an approximation to the identity argument as in Theorem 1 and then use the Functional Itô formula to (37)
F
t
,
S
t
,
σ
t
2
=
e
-
r
t
B
S
t
,
S
,
f
t
,
σ
t
2
.
We deduce that(38)
F
T
,
S
T
,
σ
T
2
-
F
t
,
S
t
,
σ
t
2
=
∫
t
T
D
u
F
u
,
S
u
,
σ
u
2
d
u
+
∫
t
T
∇
S
F
u
,
S
u
,
σ
u
2
d
S
u
+
∫
t
T
∇
σ
2
F
u
,
S
u
,
σ
u
2
d
u
+
1
2
∫
t
T
∇
S
2
F
t
,
S
t
,
σ
t
2
d
S
,
S
u
+
1
2
∫
t
T
∇
σ
2
2
F
u
,
S
u
,
σ
u
2
d
σ
u
2
,
σ
u
2
+
1
2
∫
t
T
D
u
2
F
u
,
S
u
,
σ
u
2
d
u
+
∫
t
T
∇
S
,
σ
2
2
F
u
,
S
u
,
σ
u
2
d
S
u
,
σ
u
2
+
∫
t
T
∂
f
∇
S
F
u
,
S
u
,
σ
u
2
d
S
u
,
f
u
,
σ
u
2
+
∫
t
T
∂
f
∇
σ
2
F
u
,
S
u
,
σ
u
2
d
σ
2
,
f
u
,
σ
u
2
.
Note that(i)
asS
(
t
) is not path-dependent, we have that ∇
S
(
·
)
=
∂
S
(
·
);
(ii)
asu
>
t and f is a nonanticipative functional, then ∇
σ
2
(
u
)
f
(
t
,
σ
t
2
)
=
0.
So, we have(39)
F
T
,
S
T
,
σ
T
2
-
F
t
,
S
t
,
σ
t
2
=
∫
t
T
D
u
F
u
,
S
u
,
σ
u
2
d
u
+
∫
t
T
∂
S
F
u
,
S
u
,
σ
u
2
d
S
u
+
1
2
∫
t
T
∂
S
2
F
u
,
S
u
,
σ
u
2
d
S
,
S
u
+
1
2
∫
t
T
D
u
2
F
u
,
S
u
,
σ
u
2
d
u
+
∫
t
T
∂
f
,
S
2
F
u
,
S
u
,
σ
u
2
d
S
u
,
f
u
,
σ
u
2
.We deduce that (40)
F
T
,
S
T
,
σ
T
2
-
F
t
,
S
t
,
σ
t
2
=
∫
t
T
L
f
u
,
σ
u
2
B
S
d
u
+
∫
t
T
e
-
r
u
∂
f
B
S
u
,
S
u
,
f
u
,
σ
u
2
D
u
f
u
,
σ
u
2
d
u
+
1
2
∫
t
T
e
-
r
u
∂
S
2
B
S
u
,
S
u
,
f
u
,
σ
u
2
·
θ
2
u
,
S
u
,
σ
u
-
S
2
f
2
u
,
σ
u
2
d
u
+
∫
t
T
∂
S
B
S
u
,
S
u
,
f
u
,
σ
u
2
θ
u
,
S
u
,
σ
u
·
ρ
d
W
u
+
1
-
ρ
2
d
B
u
+
1
2
·
∫
t
T
e
-
r
u
∂
f
2
B
S
u
,
S
u
,
f
u
,
σ
u
2
d
f
u
,
σ
u
2
,
f
u
,
σ
u
2
+
ρ
∫
t
T
e
-
r
u
∂
f
,
S
2
B
S
u
,
S
u
,
f
u
,
σ
u
2
·
θ
u
,
S
u
,
σ
u
d
W
u
,
f
u
,
σ
u
2
.Taking now conditional expectations, using (15), and multiplying by e
r
t, we obtain that (41)
e
-
r
T
-
t
E
t
F
T
,
S
T
,
σ
T
2
=
B
S
t
,
S
t
,
f
u
,
σ
t
2
+
E
t
∫
t
T
e
-
r
u
f
u
,
σ
u
2
T
-
t
·
S
2
∂
S
2
B
S
u
,
S
u
,
f
u
,
σ
u
2
D
u
f
u
,
σ
u
2
d
u
+
1
2
·
E
t
∫
t
T
e
-
r
u
G
u
,
S
u
,
f
u
,
σ
u
2
·
L
2
u
,
S
u
,
σ
u
-
f
2
u
,
σ
u
2
d
u
+
1
2
·
E
t
∫
t
T
e
-
r
u
f
u
,
σ
u
2
2
·
τ
2
K
u
,
S
u
,
f
u
,
σ
u
2
d
f
u
,
σ
u
2
,
f
u
,
σ
u
2
+
ρ
E
t
∫
t
T
e
-
r
u
H
u
,
S
u
,
f
u
,
σ
u
2
f
u
,
σ
u
2
·
τ
θ
u
,
S
u
,
σ
u
d
W
u
,
f
u
,
σ
u
2
.Remark 10.
Note that Functional Itô formula proved in [11] holds for semimartingales but in [12] is also proved for Dirichlet process. In both cases, the hypothesis hold by definition of f and differentiability of the derivatives of Black-Scholes function when τ
,
S
,
σ
>
0. Therefore, this technique can be applied to these models.Remark 11.
Note that Theorem9 coincides with Theorem 1 when we choose the volatility function as f
(
t
,
σ
t
2
)
=
v
(
t
). We found an equivalence of the ideas developed by [6, 10–13] in the decomposition problem. Both formulas come from different points of view; the ideas under [10–13] are based on an extension to functionals of the work [14], while the main idea of [6] is to change a process by its expectation. Realize that standard Itô Calculus also can be applied to Dirichlet processes (for more information see [14]).Remark 12.
Realize that Theorem9 holds for any nonanticipative f
(
t
,
σ
t
2
). It is not trivial to find a different nonanticipative process f
(
t
,
σ
t
2
) different from the one chosen in [6].


## 6. Basic Elements of Malliavin Calculus

In the next section, we present a brief introduction to the basic facts of Malliavin calculus. For more information, see [15].Let us consider a Brownian motionW
=
W
(
t
)
,
t
∈
0
,
T defined on a complete probability space Ω
,
F
,
P. Set H
=
L
2
(
0
,
T
), and denote by W
(
h
) the Wiener integral of a function h
∈
H. Let S be the set of random variables of the form F
=
f
(
W
(
h
1
)
,
…
,
W
(
h
n
)
), where n
≥
1, f
∈
C
b
∞, and h
1
,
…
,
h
n
∈
H. Given a random variable F of this form, we define its derivative as the stochastic process D
t
W
F
,
t
∈
[
0
,
T
] given by(42)
D
t
W
F
=
∑
i
=
1
n
∂
x
i
f
W
h
1
,
…
,
W
h
n
h
i
x
,
t
∈
0
,
T
.The operator D
W and the iterated operators D
W
,
n are closable and unbounded from L
2
(
Ω
) into L
2
(
[
0
,
T
]
n
×
Ω
), for all n
≥
1. We denote the closure of S with respect to the norm (43)
F
n
,
2
2
≔
F
L
2
Ω
2
+
∑
k
=
1
n
D
W
,
k
F
L
2
0
,
T
k
×
Ω
2
.We denote by δ
W the adjoint of the derivative operator D
W. Note that δ
W is an extension of the Itô integral in the sense that the set L
a
2
(
[
0
,
T
]
×
Ω
) of square integrable and adapted processes is included in D
o
m
δ and the operator δ restricted to L
a
2
(
[
0
,
T
]
×
Ω
) coincides with the Itô stochastic integral. We use the notation δ
(
u
)
=
∫
0
T
‍
u
(
t
)
d
W
(
t
). We recall that L
W
n
,
2
≔
L
2
[
0
,
T
]
;
D
W
n
,
2 is contained in the domain of δ for all n
≥
1.We will use the next Itô formula for anticipative processes.Proposition 13.
Let us consider the processesX
(
t
)
=
x
(
0
)
+
∫
0
t
‍
u
(
s
)
d
W
(
s
)
+
∫
0
t
‍
v
(
s
)
d
s, where u
,
v
∈
L
a
2
(
[
0
,
T
]
×
Ω
). Furthermore, consider also a process Y
(
t
)
=
∫
t
T
‍
θ
(
s
)
d
s, for some θ
∈
L
1,2. Let F
:
R
3
→
R be a twice continuously differentiable function such that there exists a positive constant C such that, for all t
∈
[
0
,
T
], F and its derivatives evaluated in (
t
,
X
(
t
)
,
Y
(
t
)
) are bounded by C. Then it follows that (44)
F
t
,
X
t
,
Y
t
=
F
0
,
X
0
,
Y
0
+
∫
0
t
∂
s
F
s
,
X
s
,
Y
s
d
s
+
∫
0
t
∂
x
F
s
,
X
s
,
Y
s
d
X
s
+
∫
0
t
∂
y
F
s
,
X
s
,
Y
s
d
Y
s
+
∫
0
t
∂
x
,
y
2
F
s
,
X
s
,
Y
s
D
-
Y
s
u
s
d
s
+
1
2
·
∫
0
t
∂
x
2
F
s
,
X
s
,
Y
s
u
2
s
d
s
,where (
D
-
Y
)
(
s
)
≔
∫
s
T
‍
D
s
W
Y
(
r
)
d
r.Proof.
See [4].The next proposition is useful when we want to calculate the Malliavin derivative.Proposition 14.
Letσ and b be continuously differential functions on R with bounded derivatives. Consider the solution X
=
X
t
,
t
∈
[
0
,
T
] of the stochastic differential equation: (45)
X
t
=
x
0
+
∫
0
t
σ
X
s
d
W
s
+
∫
0
t
b
X
s
d
s
.Then, we have (46)
D
s
X
t
=
σ
X
s
·
exp
⁡
∫
s
t
σ
′
X
s
d
W
s
+
∫
s
t
λ
s
d
s
1
0
,
t
s
,where λ
(
s
)
=
[
b
′
-
1
/
2
σ
′
2
]
(
X
(
s
)
).Proof.
See [15, Section 2.2].


## 7. Decomposition Formula Using Malliavin Calculus

In this section, we use the Malliavin calculus to extend the call option price decomposition in an anticipative framework. This time, the decomposition formula has one term less than in the Itô formula’s setup.We recall the definition of the future average volatility as(47)
σ
-
t
≔
1
T
-
t
∫
t
T
σ
2
s
d
s
.Theorem 15 (decomposition formula).
For allt
∈
[
0
,
T
), we have (48)
V
t
=
E
t
B
S
t
,
S
t
,
σ
-
t
+
1
2
·
E
t
∫
t
T
e
-
r
u
G
u
,
S
u
,
σ
-
u
·
L
2
u
,
S
u
,
σ
u
-
σ
2
u
d
u
+
ρ
2
·
E
t
∫
t
T
e
-
r
u
-
t
L
u
,
S
u
,
σ
u
H
u
,
S
u
,
σ
-
u
·
∫
u
T
D
u
W
σ
2
r
d
r
d
u
,where(49)
G
t
,
S
t
,
σ
t
≔
S
2
t
∂
S
2
B
S
t
,
S
t
,
σ
t
,
H
t
,
S
t
,
σ
t
≔
S
t
∂
S
G
t
,
S
t
,
σ
t
,
L
t
,
S
t
,
σ
t
≔
θ
t
,
S
t
,
σ
t
S
t
.Proof.
Notice thate
-
r
T
B
S
(
T
,
S
(
T
)
,
σ
-
(
T
)
)
=
e
-
r
T
V
T. As e
-
r
t
V
(
t
) is a martingale, we can write (50)
e
-
r
t
V
t
=
E
t
e
-
r
T
V
T
=
E
t
e
-
r
T
B
S
T
,
S
T
,
σ
-
T
.So, using the approximation to the identity argument and applying the Itô formula presented in Proposition 13 to(51)
e
-
r
t
B
S
t
,
S
t
,
σ
-
t
,we deduce by using (15) and (7) that (52)
e
-
r
T
B
S
T
,
S
T
,
σ
-
T
-
e
-
r
t
B
S
t
,
S
t
,
σ
-
t
=
∫
t
T
e
-
r
u
L
σ
-
S
B
S
u
,
S
u
,
σ
-
u
d
u
+
1
2
·
∫
t
T
e
-
r
u
S
2
u
∂
S
2
B
S
u
,
S
u
,
σ
-
u
·
θ
u
,
S
u
,
σ
u
S
u
2
d
u
-
1
2
∫
t
T
e
-
r
u
S
2
u
·
∂
S
2
B
S
u
,
S
u
,
σ
-
u
σ
2
u
d
u
+
∫
t
T
e
-
r
u
∂
S
B
S
u
,
S
u
,
σ
-
u
θ
u
,
S
u
,
σ
u
·
ρ
d
W
u
+
1
-
ρ
2
d
B
u
+
ρ
2
·
∫
t
T
e
-
r
u
θ
u
,
S
u
,
σ
u
·
∂
S
S
2
u
∂
S
2
B
S
u
,
S
u
,
σ
-
u
·
∫
u
T
D
u
W
σ
2
r
d
r
d
u
.
Taking conditional expectation and multiplying bye
r
t, we have (53)
E
t
e
-
r
T
-
t
B
S
T
,
S
T
,
σ
-
T
=
E
t
B
S
t
,
S
t
,
σ
-
t
+
1
2
·
E
t
∫
t
T
e
-
r
u
-
t
G
u
,
S
u
,
σ
-
u
·
L
2
u
,
S
u
,
σ
u
-
σ
2
u
d
u
+
ρ
2
·
E
t
∫
t
T
e
-
r
u
-
t
L
u
,
S
u
,
σ
u
·
H
u
,
S
u
,
σ
-
u
∫
u
T
D
u
W
σ
2
r
d
r
d
u
.Remark 16.
As it is expected, a new term emerges when it is considered (1) like it happens in Theorem 1.Remark 17.
In particular, whenθ
(
t
,
S
(
T
)
,
σ
(
t
)
)
=
σ
(
t
)
S
(
t
), (54)
V
t
=
B
S
t
,
S
t
,
σ
-
t
+
ρ
2
E
t
∫
t
T
e
-
r
u
-
t
σ
u
·
H
u
,
S
u
,
σ
-
u
∫
u
T
D
u
W
σ
2
r
d
r
d
u
.Additionally, the gamma effect is cancelled as we have seen in the Itô formula section.Remark 18.
Note that whenv
(
t
) is a deterministic function, we have that all decomposition formulas are equal.Remark 19.
Whenρ
=
0, we have (55)
E
t
B
S
t
,
S
t
,
σ
-
t
-
B
S
t
,
S
t
,
v
t
=
1
2
·
E
t
∫
t
T
e
-
r
u
-
t
G
u
,
S
u
,
σ
-
u
-
G
u
,
S
u
,
v
u
L
2
u
,
S
u
,
σ
u
d
u
-
1
2
·
E
t
∫
t
T
e
-
r
u
-
t
G
u
,
S
u
,
σ
-
u
-
G
u
,
S
u
,
v
u
σ
2
u
d
u
-
1
8
·
E
t
∫
t
T
e
-
r
u
-
t
K
u
,
S
u
,
v
u
d
M
,
M
u
.In particular, when θ
(
t
,
S
(
t
)
,
σ
(
t
)
)
=
σ
(
t
)
S
(
t
), (56)
E
t
B
S
t
,
S
t
,
σ
-
t
-
B
S
t
,
S
t
,
v
t
=
-
1
8
·
E
t
∫
t
T
e
-
r
u
-
t
K
u
,
S
u
,
v
u
d
M
,
M
u
.The difference between the two approaches is given by the vol-vol of the option.


## 8. An Expression for the Derivative of the Implied Volatility

In this section, we give a general expression for the derivative of the implied volatility under the framework of Itô Calculus and Malliavin calculus. A previous calculation of this derivative in the case of exponential models by using Malliavin calculus is given in [5].LetI
(
S
(
t
)
) denote the implied volatility process, which satisfies by definition V
(
t
)
=
B
S
(
t
,
S
(
t
)
,
I
(
S
(
t
)
)
). We calculate the derivative of the implied volatility in the standard Itô case.Proposition 20.
Under (1), for every fixed t
∈
[
0
,
T
) and assuming that (
v
(
t
)
)
-
1
<
∞ a.s., we have that (57)
∂
S
I
S
∗
t
=
E
t
∫
t
T
∂
S
F
2
u
,
S
∗
u
,
v
u
d
u
∂
σ
B
S
t
,
S
∗
t
,
I
S
∗
t
-
E
t
∫
t
T
F
1
u
,
S
∗
u
,
v
u
+
∂
S
F
3
u
,
S
∗
u
,
v
u
d
u
2
S
∂
σ
B
S
t
,
S
∗
t
,
I
S
∗
t
,where (58)
E
t
∫
t
T
F
1
u
,
S
u
,
v
u
d
u
=
1
2
·
E
t
∫
t
T
e
-
r
u
-
t
G
u
,
S
u
,
v
u
·
L
2
u
,
S
u
,
σ
u
-
σ
2
u
d
u
+
1
8
·
E
t
∫
t
T
e
-
r
u
-
t
K
u
,
S
u
,
v
u
d
M
,
M
u
+
ρ
2
E
t
∫
t
T
e
-
r
u
-
t
·
θ
u
,
S
u
,
σ
u
S
u
H
u
,
S
u
,
v
u
d
W
,
M
·
u
,
E
t
∫
t
T
F
2
u
,
S
u
,
v
u
d
u
=
ρ
2
E
t
∫
t
T
e
-
r
u
-
t
·
θ
u
,
S
u
,
σ
u
S
u
H
u
,
S
u
,
v
u
d
W
,
M
·
u
,
E
t
∫
t
T
F
3
u
,
S
u
,
v
u
d
u
=
1
2
·
E
t
∫
t
T
e
-
r
u
-
t
G
u
,
S
u
,
v
u
·
L
2
u
,
S
u
,
σ
u
-
σ
2
u
d
u
+
1
8
·
E
t
∫
t
T
e
-
r
u
-
t
K
u
,
S
u
,
v
u
d
M
,
M
u
.Proof.
Taking partial derivatives with respect toS
(
t
) on the expression V
(
t
)
=
B
S
(
t
,
S
(
t
)
,
I
(
S
(
T
)
)
), we obtain (59)
∂
S
V
t
=
∂
S
B
S
t
,
S
t
,
I
S
T
+
∂
σ
B
S
t
,
S
t
,
I
S
T
∂
S
I
S
t
.On the other hand, from Theorem 1, we deduce that (60)
V
t
=
B
S
t
,
S
t
,
v
t
+
E
t
∫
t
T
F
1
u
,
S
u
,
v
u
d
u
,which implies that (61)
∂
S
V
t
=
∂
S
B
S
t
,
S
t
,
v
t
+
E
t
∫
t
T
∂
S
F
1
u
,
S
u
,
v
u
d
u
.Using the fact that (
v
(
t
)
)
-
1
<
∞, we can check that ∂
S
V
(
t
) is well defined and finite a.s. Thus, using the fact that S
∗
t
=
K
exp
⁡
(
r
(
T
-
t
)
), (59), and (61), we obtain (62)
∂
S
I
S
∗
t
=
∂
S
B
S
t
,
S
∗
t
,
v
t
-
∂
S
B
S
t
,
S
∗
t
,
I
S
t
∂
σ
B
S
t
,
S
∗
t
,
I
S
t
+
E
t
∫
t
T
∂
S
F
1
u
,
S
∗
u
,
v
u
d
u
∂
σ
B
S
t
,
S
∗
t
,
I
S
t
.From [16], we know that ∂
S
I
0
(
t
)
=
0, where I
0
(
t
) is the implied volatility in the case ρ
=
0, so (63)
∂
S
B
S
t
,
S
∗
t
,
v
t
=
∂
S
B
S
t
,
S
∗
t
,
I
0
S
t
-
E
t
∫
t
T
∂
S
F
3
u
,
S
∗
u
,
v
u
d
u
.So, we have that (64)
∂
S
I
S
∗
t
=
∂
S
B
S
t
,
S
∗
t
,
I
0
t
-
∂
S
B
S
t
,
S
∗
t
,
I
S
∗
t
∂
σ
B
S
t
,
S
∗
t
,
I
S
∗
t
+
E
t
∫
t
T
∂
S
F
2
u
,
S
∗
u
,
v
u
d
u
∂
σ
B
S
t
,
S
∗
t
,
I
S
∗
t
.On the other hand, we have that (65)
∂
S
B
S
t
,
S
∗
t
,
v
t
=
ϕ
d
,
B
S
t
,
S
∗
t
,
v
t
=
S
ϕ
d
-
ϕ
-
d
,where ϕ is the standard Gaussian density. Then (66)
∂
S
B
S
t
,
S
∗
t
,
v
t
=
B
S
t
,
S
∗
t
,
v
t
+
S
2
S
,
∂
S
B
S
t
,
S
∗
t
,
I
0
t
-
∂
S
B
S
t
,
S
∗
t
,
I
S
∗
t
=
1
2
S
B
S
t
,
S
∗
t
,
I
0
t
-
B
S
t
,
S
∗
t
,
I
S
∗
t
=
-
1
2
S
E
t
∫
t
T
F
1
u
,
S
∗
u
,
v
u
+
∂
S
F
3
u
,
S
∗
u
,
v
u
d
u
.Now, we derive the implied volatility using Malliavin calculus. This is done in [5] for the case θ
(
t
,
S
(
T
)
,
σ
(
t
)
)
=
σ
(
t
)
S
(
t
).Proposition 21.
Under (1), for every fixed t
∈
[
0
,
T
) and assuming that (
σ
~
(
t
)
)
-
1
<
∞ a.s., we have that (67)
∂
S
I
S
∗
t
=
E
t
∫
t
T
∂
S
F
2
u
,
S
∗
u
,
σ
-
u
d
u
∂
σ
B
S
t
,
S
∗
t
,
I
S
∗
t
-
E
t
∫
t
T
F
1
u
,
S
∗
u
,
σ
-
u
+
∂
S
F
3
u
,
S
∗
u
,
σ
-
u
d
u
2
S
∂
σ
B
S
t
,
S
∗
t
,
I
S
∗
t
,where (68)
E
t
∫
t
T
F
1
u
,
S
u
,
σ
-
u
d
u
=
1
2
·
E
t
∫
t
T
e
-
r
u
-
t
G
u
,
S
u
,
σ
-
u
·
L
2
u
,
S
u
,
σ
u
-
σ
2
u
d
u
+
ρ
2
·
E
t
∫
t
T
e
-
r
u
-
t
L
u
,
S
u
,
σ
u
·
H
u
,
S
u
,
σ
-
u
∫
u
T
D
u
W
σ
2
r
d
r
d
W
,
M
·
u
,
E
t
∫
t
T
F
2
u
,
S
u
,
σ
-
u
d
u
=
ρ
2
·
E
t
∫
t
T
e
-
r
u
-
t
L
2
u
,
S
u
,
σ
u
·
∫
u
T
D
u
W
σ
2
r
d
r
d
W
,
M
u
,
E
t
∫
t
T
F
3
u
,
S
u
,
σ
-
u
d
u
=
1
2
·
E
t
∫
t
T
e
-
r
u
-
t
G
u
,
S
u
,
σ
-
u
·
L
2
u
,
S
u
,
σ
u
-
σ
2
u
d
u
.Proof.
See [5] or the previous proof.Remark 22.
Note that this is generalization of the formula proved in [5]. In that case, F
1
=
F
2 and F
3
=
0.


## 9. Examples

In this section, we provide some applications of the decomposition formula to well-known models in Finance.


### 9.1. Heston Model

We consider that the stock price follows the Heston Model (1). Using Theorem 1 or Theorem 9, we have(69)
V
t
=
B
S
t
,
X
t
,
v
t
+
ρ
2
·
E
t
∫
t
T
e
-
r
u
-
t
H
u
,
X
u
,
v
u
·
∫
u
T
e
-
k
r
-
s
d
r
σ
2
u
ν
d
u
+
1
8
·
E
t
∫
t
T
e
-
r
u
-
t
K
u
,
X
u
,
v
u
·
∫
u
T
e
-
k
r
-
s
d
r
2
ν
2
σ
2
u
d
u
.Using Theorem15, we have that (70)
V
t
=
B
S
t
,
S
t
,
σ
-
t
+
ρ
2
·
E
t
∫
t
T
e
-
r
u
-
t
H
u
,
S
u
,
σ
-
u
·
∫
u
T
D
u
W
σ
2
r
d
r
σ
u
d
u
,where D
u
W
σ
2
(
r
)
=
ν
σ
(
u
)
e
x
p
ν
/
2
∫
u
r
‍
(
1
/
σ
(
s
)
)
d
W
(
s
)
+
∫
u
r
‍
-
k
-
ν
2
/
8
σ
2
(
s
)
d
s
.


### 9.2. SABR Model

We consider that the stock price follows the SABR model (3). Using Theorem 1 or Theorem 9, we have(71)
V
t
=
B
S
t
,
S
t
,
v
t
+
1
2
·
E
t
∫
t
T
e
-
r
u
-
t
G
u
,
S
u
,
v
u
σ
2
u
·
S
2
β
-
1
u
-
1
d
u
+
1
8
·
E
t
∫
t
T
e
-
r
u
-
t
K
u
,
S
u
,
v
u
d
M
,
M
u
+
ρ
2
E
t
∫
t
T
e
-
r
u
-
t
H
u
,
S
u
,
v
u
σ
u
d
W
,
M
·
u
,where(72)
d
M
,
M
=
4
α
2
σ
4
t
∫
t
T
e
α
2
s
-
t
d
s
2
d
t
,
d
M
,
W
=
2
α
σ
2
t
∫
t
T
e
α
2
s
-
t
d
s
d
t
.Using Theorem15, we have that (73)
V
t
=
E
t
B
S
t
,
S
t
,
σ
-
t
+
1
2
·
E
t
∫
t
T
e
-
r
u
-
t
G
u
,
S
u
,
σ
-
u
σ
2
u
·
S
2
β
-
1
u
-
1
d
u
+
ρ
2
·
E
t
∫
t
T
e
-
r
u
-
t
H
u
,
S
u
,
σ
-
u
·
∫
u
T
D
u
W
σ
2
r
d
r
σ
u
d
u
,where D
u
W
σ
2
(
r
)
=
2
α
σ
2
(
u
)
1
[
0
,
r
]
(
u
).


## 9.1. Heston Model

We consider that the stock price follows the Heston Model (1). Using Theorem 1 or Theorem 9, we have(69)
V
t
=
B
S
t
,
X
t
,
v
t
+
ρ
2
·
E
t
∫
t
T
e
-
r
u
-
t
H
u
,
X
u
,
v
u
·
∫
u
T
e
-
k
r
-
s
d
r
σ
2
u
ν
d
u
+
1
8
·
E
t
∫
t
T
e
-
r
u
-
t
K
u
,
X
u
,
v
u
·
∫
u
T
e
-
k
r
-
s
d
r
2
ν
2
σ
2
u
d
u
.Using Theorem15, we have that (70)
V
t
=
B
S
t
,
S
t
,
σ
-
t
+
ρ
2
·
E
t
∫
t
T
e
-
r
u
-
t
H
u
,
S
u
,
σ
-
u
·
∫
u
T
D
u
W
σ
2
r
d
r
σ
u
d
u
,where D
u
W
σ
2
(
r
)
=
ν
σ
(
u
)
e
x
p
ν
/
2
∫
u
r
‍
(
1
/
σ
(
s
)
)
d
W
(
s
)
+
∫
u
r
‍
-
k
-
ν
2
/
8
σ
2
(
s
)
d
s
.


## 9.2. SABR Model

We consider that the stock price follows the SABR model (3). Using Theorem 1 or Theorem 9, we have(71)
V
t
=
B
S
t
,
S
t
,
v
t
+
1
2
·
E
t
∫
t
T
e
-
r
u
-
t
G
u
,
S
u
,
v
u
σ
2
u
·
S
2
β
-
1
u
-
1
d
u
+
1
8
·
E
t
∫
t
T
e
-
r
u
-
t
K
u
,
S
u
,
v
u
d
M
,
M
u
+
ρ
2
E
t
∫
t
T
e
-
r
u
-
t
H
u
,
S
u
,
v
u
σ
u
d
W
,
M
·
u
,where(72)
d
M
,
M
=
4
α
2
σ
4
t
∫
t
T
e
α
2
s
-
t
d
s
2
d
t
,
d
M
,
W
=
2
α
σ
2
t
∫
t
T
e
α
2
s
-
t
d
s
d
t
.Using Theorem15, we have that (73)
V
t
=
E
t
B
S
t
,
S
t
,
σ
-
t
+
1
2
·
E
t
∫
t
T
e
-
r
u
-
t
G
u
,
S
u
,
σ
-
u
σ
2
u
·
S
2
β
-
1
u
-
1
d
u
+
ρ
2
·
E
t
∫
t
T
e
-
r
u
-
t
H
u
,
S
u
,
σ
-
u
·
∫
u
T
D
u
W
σ
2
r
d
r
σ
u
d
u
,where D
u
W
σ
2
(
r
)
=
2
α
σ
2
(
u
)
1
[
0
,
r
]
(
u
).


## 10. Conclusion

In this paper, we notice that the idea used in [6] can be used for a generic Stochastic Differential Equation (SDE). There is no need to specify the volatility process. Only existence and uniqueness of the solution of the SDE are needed, allowing much more flexibility in the decomposition formula. We observe what is the effect of assuming that the stock price follows an exponential process and how a new term arises in a general framework. Additionally, we have computed the decomposition using three different methods: Itô formula, Functional Itô Calculus, and Malliavin calculus. In the case of call options, the idea used in [6] is equivalent to the use of the Functional Itô formula developed in [10–13] but without the need of the theory behind the Functional Itô Calculus. Both formulas can be applied to Dirichlet process, in particular, to the fractional Brownian motion with Hurst parameter equal or bigger than 1
/
2. Furthermore, we realize that the Feynman-Kac formula has a key role in the decomposition process.



---

*Source: 103647-2015-06-17.xml*