Abstract
Multivariate probability density functions of returns are constructed in order to model the empirical behavior of returns in a financial time series. They describe the well-established deviations from the Gaussian random walk, such as an approximate scaling and heavy tails of the return distributions, long-ranged volatility-volatility correlations (volatility clustering) and return-volatility correlations (leverage effect). Free parameters of the model are fixed over the long term by fitting 100+ years of daily prices of the Dow Jones 30 Industrial Average. The multivariate probability density functions which we have constructed can be used for pricing derivative securities and risk management.