Probability distribution of returns in the Heston model with stochastic volatility*

Abstract

We study the Heston model, where the stock price dynamics is governed by a geometrical (multiplicative) Brownian motion with stochastic variance. We solve the corresponding Fokker-Planck equation exactly and, after integrating out the variance, find an analytic formula for the timedependent probability distribution of stock price changes (returns). The formula is in excellent agreement with the Dow-Jones index for the time lags from 1 to 250 trading days. For large returns, the distribution is exponential in log-returns with a time-dependent exponent, whereas for small returns it is Gaussian. For time lags longer than the relaxation time of variance, the probability distribution can be expressed in a scaling form using a Bessel function. The Dow-Jones data for 1982-2001 follow the scaling function for seven orders of magnitude.