Modern financial engineering is a part of applied mathematics that studies market models. Each model is characterized by several parameters. Some of them are familiar to a wide audience, for example, the price of a risky security, or the risk free interest rate. Other parameters are less known, for example, the volatility of the security. This parameter determines the rate of change of security prices and is determined by several factors. For example, during the periods of stable economic growth the prices are changing slowly, and the volatility is small. During the crisis periods, the volatility significantly increases. Classical market models, in particular, the celebrated Nobel Prize awarded Black–Scholes–Merton model (1973), suppose that the volatility remains constant during the lifetime of a financial instrument. Nowadays, in most cases, this assumption cannot adequately describe reality. We consider a model where both the security price and the volatility are described by random functions of time, or stochastic processes. Moreover, the volatility process is modelled as a sum of two independent stochastic processes. Both of them are mean reverting in the sense that they randomly oscillate around their average values and never escape neither to very small nor to very big values. One is changing slowly and describes low frequency, for example, seasonal effects, another is changing fast and describes various high frequency effects. We formulate the model in the form of a system of a special kind of equations called stochastic differential equations. Our system includes three stochastic processes, four independent factors, and depends on two small parameters. We calculate the price of a particular financial instrument called European call option. This financial contract gives its holder the right (but not the obligation) to buy a predefined number of units of the risky security on a predefined date and pay a predefined price. To solve this problem, we use the classical result of Feynman (1948) and Kac (1949). The price of the instrument is the solution to another kind of problem called boundary value problem for a partial differential equation. The resulting equation cannot be solved analytically. Instead we represent the solution in the form of an expansion in the integer and half-integer powers of the two small parameters mentioned above. We calculate the coefficients of the expansion up to the second order, find their financial sense, perform numerical studies, and validate our results by comparing them to known verified models from the literature. The results of our investigation can be used by both financial institutions and individual investors for optimization of their incomes.