Perturbation Expansion for Option Pricing with Stochastic Volatility

We fit the volatility fluctuations of the S&P 500 index well by a Chi distribution, and the distribution of log-returns by a corresponding superposition of Gaussian distributions. The Fourier transform of this is, remarkably, of the Tsallis type. An option pricing formula is derived from the same superposition of Black-Scholes expressions. An explicit analytic formula is deduced from a perturbation expansion around a Black-Scholes formula with the mean volatility. The expansion has two parts. Th...