Received 30 October 2015; accepted 22 December 2015; published 25 December 2015

ABSTRACT

The study analyses some problems arising in stochastic volatility models by using Ito’s lemma and its applications to boundary Cauchy problem by giving the solution of vanilla option pricing models satisfying the partial differential equation obtained by assuming stochastic volatility in replication problems and risk neutral probability.

Keywords:

Contingent Claim, Stochastic Volatility, Ito’s Lemma, Cauchy problem, Bivariate

1. Introduction

In finance Wiener process and geometric Brown process are largely used. The name came from George Brown in the 1827 noted that the volatility of a small particle suspended in a liquid increased with the time. Wiener gave a mathematical formal assumption on the phenomena from the term of Wiener process. The main property of the Wiener process is that it is a forward process such that we may integrate it although it is a function of infinite variation; the main idea is that the process converges to the discrete process because the limit tends toward the discrete process when it is shared in sub intervals. From this we may approximate the Wiener process in the instant as. This permits to give the proof of Ito’s lemma by using Taylor series. Indeed, the problem is more complicated because the diffusion process is a standardized normal distribution. This permits to give the solution to the Cauchy problem subject to boundary although the problem is parabolic. The case of stochastic volatility may be viewed as Cauchy problem where the diffusion process of Ito’s lemma is a bivariate standardized normal distribution. Thus we may solve easily the problem of option pricing with stochastic volatility in risk neutral world by using the integrant factor.

2. The Model and Its Assumptions

The geometric Brown process is used in finance to indicate a formal assumption for the dynamic of the prices that does not permit to assume negative value, formally we have:

(1)

where denotes the drift of the distribution and it is the average in the dt, σ denotes the volatility of the distribution and denotes a Wiener process such that it may be decomposed by the following:

(2)

We may assume the following for the Wiener process:

(3)

This means that a Wiener process is a forward process, the uncertainty is to the end of the process in T + dt. From this we may obtain explanation of Ito’s lemma by using Taylor series, if we take a function of S as F(S) we may write Ito’s lemma in the following way:

(4)

We may note that:

(5)

(6)

From this we obtain as dt tends to zero:

(7)

By substituting dS we obtain Ito’s lemma:

(8)

We may see now as to obtain the expected value of a normal distribution as such we have the following:

(9)

As such we have the following:

(10)

This may be rewritten by:

(11)

From this we may obtain explanation for Ito’s lemma, if we take a function of S as we may write Ito’s lemma in the following way:

(12)

where:

(13)

As result:

(14)

Because:

(15)

Now we may analyze the following parabolic problem:

(16)

Subject to the following constraint:

(17)

The solution it is easy to solve, because if we take Ito’s lemma and we take the expectation we obtain that the solution to the parabolic problem is given by:

(18)

It is interesting to introduce the concept of stochastic volatility, as such we may write Ito’s lemma in the following form:

(19)

where denotes a standardized bivariate normal distribution with the following form:

(20)

where:

(21)

(22)

(23)

(24)

The PDE that an option must satisfy by assuming stochastic volatility is given by the following:

(25)

where:

(26)

(27)

(28)

The solution it is easy to solve, because if we take Ito’s lemma and we take the expectation we obtain that the solution to the parabolic problem is given by the following by using the integrant factor:

(29)

The final pay off of a Call and Put option is given by the following:

(30)

(31)

The prices of the options are given by the expectation of the final pay off discounted for the Call options, instead, for the Put options we assume that to replicate the value of options we have to have available the amount of money K that will produce a risk free earnings:

(32)

(33)

So from 33 we may note that for the Put options will be discounted only the process. Now to compute the value of the option is a problem because we have stochastic interest rate so the solution is to take the default free zero coupon bond as forward measure. We assume the following process for the default free zero coupon bond:

(34)

(35)

As result we obtain the following pricing formula for the options by using the respective numeraire:

(36)

where:

(37)

(38)

(39)

(40)

(41)

(42)

where:

(43)

(44)

3. Numerical Results

We may compare the model with [1] and [2] , as results we have the following figures for Call options (Tables 1-3).

As results we have the following figures for Put options (Tables 4-6).

We may note from the numerical results that for rational value of parameters [1] converges to [2] , instead, the bivariate approach permits to capture the skew for options deep in the money, so it is a good candidate in period of crisis. We may analyses further the formulations obtained by using Monte Carlo simulations (Appendix), but the approach has a drawback, for the Put option the simulation gives always the same results so we have to use the fair Put Call parity, as such we have the following prospects (Tables 7-9).

Indeed, with the simulations the normal distribution is more skewed, so we used a calibrated parameters drift, to compare the two approaches [3] -[7] .

4. Conclusion

The study permits to obtain closed form solution for option pricing with stochastic volatility by assuming normal distribution obtained by the properties of the bivariate standardized normal distribution that is the solution of the Cauchy problem with stochastic volatility due to the properties of Ito’s lemma.

Cite this paper

RossanoGiandomenico, (2015) Option Pricing with Stochastic Volatility. Journal of Applied Mathematics and Physics,03,1645-1653. doi: 10.4236/jamp.2015.312189

References

Appendix

To run the simulation we used the following VBA code.