`Journal of Applied MathematicsVolume 2012 (2012), Article ID 417942, 14 pageshttp://dx.doi.org/10.1155/2012/417942`
Research Article

## Numerical Solutions of a Variable-Order Fractional Financial System

1School of Business, Central South University, Hunan, Changsha 410083, China
2Department of Applied Mathematics, Central South University, Hunan, Changsha 410083, China

Received 10 May 2012; Revised 21 July 2012; Accepted 6 August 2012

Copyright © 2012 Shichang Ma et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

The numerical solution of a variable-order fractional financial system is calculated by using the Adams-Bashforth-Moulton method. The derivative is defined in the Caputo variable-order fractional sense. Numerical examples show that the Adams-Bashforth-Moulton method can be applied to solve such variable-order fractional differential equations simply and effectively. The convergent order of the method is also estimated numerically. Moreover, the stable equilibrium point, quasiperiodic trajectory, and chaotic attractor are found in the variable-order fractional financial system with proper order functions.

#### 1. Introduction

Variable-order fractional calculus (i.e., the fractional differentiation and integration of variable order) is the generalization of classical calculus and fractional calculus, which were invented by Newton and Leibnitz hundreds of years ago. Now the study on it becomes a hotpot in recent ten years [16]. It has turned out that many problems in physics, biology, engineering, and finance can be described excellently by models using mathematical tools from variable-order fractional calculus, such as mechanical applications [4], diffusion process [5], multifractional Gaussian noise [7], and FIR filters [8]. For more details, see [6, 9] and references therein.

As the generalized form of fractional differential equations (see [1013] for a comprehensive review of fractional calculus and fractional differential equations), the variable-order fractional differential equations mean the differential equations with variable-order fractional derivatives [2, 46, 9]. For lack of direct viewing physical interpretation, the variable-order fractional calculus developed at a very low pace during the foregone several decades. In 1995, the fractional integration and differentiation of variable-order fractional are studied by Samko in [2], which contains mathematical analysis of variable-order fractional differential and integral operators but not variable-order fractional models. The later research results have shown that many complex physical phenomena can be better described by using variable-order fractional differential equations. For instance, an experimental investigation of temperature-dependent variable-order fractional integrator and differentiator is presented in [9]. Based on the investigation, it is easier to understand the physical meaning of variable-order fractional calculus and better to know how the variable-order fractional operators are applied in physical processes. As a concrete example, the variable-order fractional calculus through the description of a simple problem in mechanics is presented in [4]. A mathematical definition for the variable-order fractional differential operator is given and suitable for mechanical modelling therein. But to the best of our knowledge, there is no financial system described in variable-order fractional sense. We eagerly attempt to study the numerical methods for variable-order fractional differential equations and find their high-order numerical solutions.

In this paper, we consider a variable-order fractional financial system: where is the Caputo variable-order fractional differential operator whose definition is given in the next section and , , are the order functions which are continuous and bounded. Variable denotes the interest rate, is the investment demand, and represents the price index. The parameters , , and denote the savings amount, cost per investment, and the elasticity of demand of commercial markets, respectively.

Unfortunately, most variable-order fractional differential equations do not have an exact analytical solution, even with simply constant coefficients, since the order of equations is a function but not an integer. Thus the effective and applicable numerical techniques for solving such equations are always needed. The Adams-Bashforth-Moulton method, studied and discussed penetratingly in [14, 15], is a type of predictor-corrector methods and a relatively new approach to provide numerical approximation to various nonlinear problems, including fractional differential equations. In recent years, the application of the method is extended to more concrete physical and mathematical models (for more details, see [1626]). But the work on numerical methods for variable-order fractional differential equations is still less common. Motivated by the work mentioned previously, we want to employ the Adams-Bashforth-Moulton method for solving a variable-order fractional financial system and obtain its high-order numerical results.

The remainder of the paper is organized as follows. In Section 2, we give some definitions and mathematical background of variable-order fractional calculus which are essential for understanding the financial system. In Section 3, we introduce the known Adams-Bashforth-Moulton method. In Section 4, numerical experiments are given. The effectiveness of the Adams-Bashforth-Moulton method for solving variable-order fractional differential equation is illustrated, and the convergent order is also estimated. Finally, we propose our conclusions in Section 5.

#### 2. Preliminaries

In this part, we give some definitions of fractional derivative and variable-order fractional derivative. There exist different approaches for defining the fractional derivatives [11]; the common used cases are called Grünwald-Letnikov definition, Riemann-Liouville definition, and Caputo definition.

Before we state the definitions of fractional derivative and variable-order fractional derivative, let us recall the difference approximation of integer-order derivative, which reads and so on. Similarly, the th order integer derivative could be defined as follows.

Definition 2.1 (see [11, 14]). Let be a positive integer; then the th order derivative of a given differentiable function is

Similarly, replacing integer with noninteger , one can have the following Grünwald-Letnikov fractional differential derivative.

Definition 2.2 (see [11, 14]). Let be a positive real number, , and . Then with being called the Grünwald-Letnikov fractional derivative of order of the function .

For describing the definitions of fractional differential operator and variable-order fractional differential operator of other types, we need the following Euler's Gamma function.

Definition 2.3 (see [11, 14]). The function , defined by is called Euler's Gamma function (or Euler's integral of the second kind). Particularly, we have for is a positive integer.

Definition 2.4 (see [11]). The Riemann-Liouville fractional order derivative of is defined as where , provided the right side is pointwise defined on .

Definition 2.5 (see [11, 14]). Let be a real number, and . Then the Caputo fractional derivative of is defined as where is a positive integer number and is the Euler's Gamma function.

According to the previous definitions of fractional derivative, we can propose the variable-order fractional differential operator now. Replacing the constant order with a given function, the fractional derivative is indeed extended to the variable-order fractional sense.

Definition 2.6. Let be a positive real number, , , and . Then with being called the Grünwald-Letnikov variable-order fractional derivative of order of the function .

Definition 2.7. The Riemann-Liouville variable-order fractional derivative of is defined as where , provided the right side is pointwise defined on .

Definition 2.8 (see [1]). Let be a continuous and bounded function, , and . Then are called the Caputo variable-order fractional derivative of , where , is the integer part of , and is the Euler's Gamma function.

Remark 2.9. In most modelling applications, we note that the Caputo version is usually preferred in the descriptions of financial and physical models, because the initial conditions of the system have a straightforward financial and physical meaning, and therefore it is in general possible to provide these data, for example, by suitable measurements. Here, we consider the objective function in the interval and suppose that satisfies some smoothness conditions in every subinterval of with . For the sake of simplification, our fractional operator can be simplified as .

Since the variable-order fractional differential operators are related to integration or summation, it may lead to the various definitions of such variable-order fractional differential operators: to let the argument of be the current time instant ; to let the argument of be the same as that of , which is the integration variable . There is no need to list all the possible definitions of Grünwald-Letnikov variable-order fractional derivative, Riemann-Liouville variable-order fractional derivative, and Caputo variable-order fractional derivative here. For their concrete definitions, we refer to [6, Definition 6–14].

There is one more fundamental difference between differential operator of integer order and variable-order fractional (including fractional order) derivatives [14]. The former are local operators, and the latter are not. The meaning of “local” here is as follows. In order to calculate for , it is sufficient to know in an arbitrarily small neighborhood of . This follows from the classical representation of as a limit of a difference quotient, and it also tells us a truth that the classical differential operator does not have memory property. However, to calculate for (i.e., is a function related to , or is not an integer), the definitions tell us that we need to know throughout the entire interval . All the historical data are related to the current state value, and therefore we say that the fractional order differential operators and variable-order fractional differential operators have memory property. Specifically, since the variable-order fractional differential operator employs function as order, it does have time-varying long memory property when the order function is related to time.

Most of fractional differential equations, even common ordinary differential equations, such as Riccati equations, do not have exact analytic solutions, not to mention the differential equations with variable-order fractional derivatives, which are much more complex. Therefore, the numerical methods are necessary and deserve comprehensive study. Several numerical methods have been proposed to solve the fractional differential equations; for more details, see [15, 27, 28]. Hence, we can employ them for solving variable-order fractional differential equations after some necessary modifications. In this paper, we use the Adams-Bashforth-Moulton method to solve variable-order fractional differential equations in simulations. To illustrate the numerical results of nonlinear variable-order fractional differential equations by means of this algorithm, consider the following differential equation:

Similar to the fractional order sense, the differential equation (3.1) is equivalent to the following Volterra integral equation [29]: where denotes the integer part of . For better describing the discretization of (3.2), we select the following equispaced grid: where . Then (3.2) becomes where the predicted value is determined by in which we can deduce the coefficients as

The estimation error of this approximation is described as follows [14]: where , . Now we apply the mentioned method to determine the numerical solution of variable-order fractional system. Consider the following variable-order fractional system: with   () and initial condition . According to the previous method, the system (3.8) can be discretized as follows: where

Remark 3.1. Since the Adams-Bashforth-Moulton method has been successfully applied in differential equations with fixed fractional order, and the stability of the Adams-Bashforth-Moulton method has been verified in [14], then it is unnecessary to repeat them here. For variable-order fractional case, the fractional order is bounded, so that every iteration can be treated as one single step iteration for different fixed fractional orders. In what follows, we provide a test example to illustrate the effectiveness of such method in solving differential equations of variable-order fractional derivatives.

#### 4. Numerical Experiments

In this section, two numerical examples are considered. In the first simple example, we compute its numerical solution with different step sizes and estimate the convergent order. While in the second example, it is found that with some proper parameters and order functions, the variable-order fractional financial system exhibits the stable equilibrium point, quasiperiodic trajectory, and chaotic motion. Three different cases, the integer order case, the fractional order case, and the variable-order case, are studied separately. Therefore we can see how different behaviors appear.

Example 4.1. Consider the following differential equation [29]: where and .

Applying the Adams-Bashforth-Moulton method, we evaluate the numerical solutions with different step sizes. Figure 1 shows the approximate solution for Example 4.1. From the figure, it is clear to see that the Adams-Bashforth-Moulton method can solve variable-order fractional differential equation simply and effectively. Moreover, we estimate the posteriori error and convergent order of the Adams-Bashforth-Moulton method in calculation which are shown in Table 1.

Table 1: Maximum error and convergent order of Example 4.1 with different step sizes (computed once via the Adams-Bashforth-Moulton method).
Figure 1: The numerical solutions of Example 4.1 with different step sizes.

Furthermore, we can apply the corrector step in Adams-Bashforth-Moulton method for more times when solving fractional differential equations, since it has been shown that the number of corrector steps may be usefully varied to improve the convergence order, based on the fractional differential equations [14]. In what follows, applying the computing procedures once again, we treat the corrected value as predicted value and obtain more accurate numerical solutions. The study results are shown in Table 2.

Table 2: Maximum error and convergent order of Example 4.1 with different step sizes (computed twice via the Adams-Bashforth-Moulton method).

Example 4.2. Consider the following variable-order fractional financial system:

We will solve the financial system (4.2) above by using the Adams-Bashforth-Moulton method. It is found that the variable-order fractional financial system (4.2) can exhibit the stable equilibrium point, quasiperiodic trajectory, and chaotic motion with different order functions. The phase portrait and largest Lyapunov exponent are used to identify the dynamics of the financial system (4.2). The positive largest Lyapunov exponent implies that the financial system (4.2) will generate chaotic motion. The dynamics of the financial system (4.2) is somewhat different from classical financial system whose order functions remain one, which is discussed in [30, 31]. In what follows, we start to study this example in three cases, and then see how different behaviors appear.

Case 1 (the integer order case). Let . Now the financial system (4.2) becomes a classical ordinary differential system. By the Adams-Bashforth-Moulton method, we can easily compute its numerical results and observe that the financial system (4.2) exhibits chaotic motion. The largest Lyapunov exponent of the system is . The following Figure 2 shows our numerical simulations.

Figure 2: The numerical solutions of Example 4.2 with integer order derivatives and step size ((a) 3D portrait; (b) - projection).

Case 2 (the fractional order case). First we choose the commensurate fractional order case: . Now the financial system (4.2) becomes a fractional order differential system. By the Adams-Bashforth-Moulton method, we can easily compute its numerical results and observe that the financial system (4.2) exhibits chaotic motion. The largest Lyapunov exponent of the system is . Figure 3 shows our numerical simulations.
Furthermore, we consider the incommensurate fractional order case: , , . Now the financial system (4.2) also becomes a fractional order differential system. By the Adams-Bashforth-Moulton method, we can easily compute its numerical results and observe that the financial system (4.2) exhibits chaotic motion. The largest Lyapunov exponent of the system is . Figure 4 shows our numerical simulations.

Figure 3: The numerical solutions of Example 4.2 with commensurate fractional order derivatives and step size ((a) 3D portrait; (b) - projection).
Figure 4: The numerical solutions of Example 4.2 with incommensurate fractional order derivatives and step size ((a) 3D portrait; (b) - projection).

Case 3 (the variable-order case). Firstly, we choose , , and with , which are rigid monotony decrease functions. Now the financial system (4.2) becomes a variable-order fractional financial system. Figure 5 shows the approximate solutions for Example 4.2. From the figure, it is clear to see that the Adams-Bashforth-Moulton method can solve variable-order fractional differential equation simply and effectively, and financial system (4.2) tends to the stable equilibrium point.
Secondly, we choose , , and with , which are rigid monotony decrease functions. Now the financial system (4.2) becomes a variable-order fractional financial system. Figure 6 shows the approximate solutions for Example 4.2. From the figure, it is clear to see that the Adams-Bashforth-Moulton method can solve variable-order fractional differential equation simply and effectively, and financial system (4.2) tends to the quasiperiodic trajectory.
Furthermore, we study variable-order fractional financial system (4.2) with nonlinear order functions. In simulation, by choosing , , and with , we see that the chaotic attractor exists in such financial system. The largest Lyapunov exponent of the system is . The following Figures 7 and 8 illustrate the order functions and numerical solution we obtain.
In practice, the stable equilibrium point and quasiperiodic trajectory mean that the financial system is stable during evolution, while the chaotic motion implies that the financial system is unpredictable.

Figure 5: The numerical solutions of Example 4.2 with initial value and step size : stable equilibrium point ((a) 3D portrait; (b) - projection).
Figure 6: The numerical solutions of Example 4.2 with initial value and step size : quasiperiodic trajectory ((a) 3D portrait; (b) - projection).
Figure 7: The nonlinear order functions in simulation.
Figure 8: The numerical solution of Example 4.2 with step size and nonlinear order functions ((a) 3D portrait; (b) - projection).

#### 5. Conclusions

Variable-order fractional calculus has been highly neglected since it was proposed. Nevertheless, the scientific community has found a large variety of applications which can be modeled and more clearly understood by using this branch of mathematics. Proposed as a generalization of a fractional financial system, a variable-order fractional financial system, which consists of interest rate, investment demand, and price index, has been studied in this paper. By using the Adams-Bashforth-Moulton method, we obtain the numerical solution of the variable-order fractional financial system. We estimate the convergent order of the Adams-Bashforth-Moulton method, which satisfies our theoretical analysis. Moreover, it is found that the variable-order fractional financial system exhibits chaotic attractor when proper order functions are chosen. According to the numerical examples, we can easily conclude that the Adams-Bashforth-Moulton method is a powerful technique to calculate approximate solutions of variable-order fractional equations.

Last but not the least, we hope our work about variable-order fractional calculus would generate interest from related scholars in the future and also hope that their work may result in significant contributions to this field.

#### Acknowledgments

The authors would like to thank the main editor and anonymous referees for their valuable comments and suggestions leading to improvement of this paper. The first author is partly supported by the National Natural Science Foundation of China (Grant no. 70921001), Research Fund for the Humanities and Social Sciences of Higher Education of China (Grant no. 06JD82002), and the Philosophy and Social Science Foundation Commissioned Project of Hunan Province, China (Grant no. 2010JD02). The second author is supported by the National Foundation for Studying Abroad (no. 201206370079).

#### References

1. S. Umarov and S. Steinberg, “Variable order differential equations and diffusion processes with changing modes,” submitted, http://www.arxiv.org/pdf/0903.2524.pdf.
2. S. G. Samko, “Fractional integration and differentiation of variable order,” Analysis Mathematica, vol. 21, no. 3, pp. 213–236, 1995.
3. C. F. Lorenzo and T. T. Hartley, “Variable order and distributed order fractional operators,” Nonlinear Dynamics, vol. 29, no. 1–4, pp. 57–98, 2002.
4. C. F. M. Coimbra, “Mechanics with variable-order differential operators,” Annalen der Physik, vol. 12, no. 11-12, pp. 692–703, 2003.
5. H. G. Sun, W. Chen, and Y. Q. Chen, “Variable-order fractional differential operators in anomalous diffusion modeling,” Physica A, vol. 388, no. 21, pp. 4586–4592, 2009.
6. D. Valério and J. S. Costa, “Variable-order fractional derivatives and their numerical approximations,” Signal Processing, vol. 91, no. 3, pp. 470–483, 2011.
7. H. Sheng, H. G. Sun, Y. Q. Chen, and T. S. Qiu, “Synthesis of multifractional Gaussian noises based on variable-order fractional operators,” Signal Processing, vol. 91, no. 7, pp. 1645–1650, 2011.
8. C. C. Tseng, “Design of variable and adaptive fractional order FIR differentiators,” Signal Processing, vol. 86, no. 10, pp. 2554–2566, 2006.
9. H. Sheng, H. G. Sun, C. Coopmans, Y. Q. Chen, and G. W. Bohannan, “A physical experimental study of variable-order fractional integrator and differentiator,” The European Physical Journal, vol. 193, no. 1, pp. 93–104, 2011.
10. K. B. Oldham and J. Spanier, The Fractional Calculus: Theory and Applications of Differentiation and Integration to Arbitrary Order, Academic Press, 1974.
11. I. Podlubny, Fractional Differential Equations, Academic Press, San Diego, Calif, USA, 1999.
12. P. L. Butzer and U. Westphal, An Introduction to Fractional Calculus, World Scientific, Singapore, 2000.
13. K. S. Miller and B. Ross, An Introduction to the Fractional Calculus and Fractional Differential Equations, Wiley-Interscience, New York, NY, USA, 1993.
14. K. Diethelm, The Analysis of Fractional Differential Equations, Springer, Berlin, Germany, 2010.
15. K. Diethelm, N. J. Ford, and A. D. Freed, “A predictor-corrector approach for the numerical solution of fractional differential equations,” Nonlinear Dynamics, vol. 29, no. 1–4, pp. 3–22, 2002.
16. K. Yang and J. M. Wang, “A temperature prediction-correction method for estimating surface soil heat flux from soil temperature and moisture data,” Science in China D, vol. 51, no. 5, pp. 721–729, 2008.
17. A. Bnouhachem and M. A. Noor, “Numerical comparison between prediction-correction methods for general variational inequalities,” Applied Mathematics and Computation, vol. 186, no. 1, pp. 496–505, 2007.
18. X. L. Fu, “A two-stage prediction-correction method for solving monotone variational inequalities,” Journal of Computational and Applied Mathematics, vol. 214, no. 2, pp. 345–355, 2008.
19. P. Midya, B. Roeckner, P. Rakers, and P. Wagh, “Prediction correction algorithm for natural pulse width modulation,” in Proceedings of the 109th AES Convention, September 2000.
20. X. M. Yuan, “The prediction-correction approach to nonlinear complementarity problems,” European Journal of Operational Research, vol. 176, no. 3, pp. 1357–1370, 2007.
21. L. J. Wicker, “A two-step Adams-Bashforth-Moulton split-explicit integrator for compressible atmospheric models,” Monthly Weather Review, vol. 137, no. 10, pp. 3588–3595, 2009.
22. K. H. Sun, X. Wang, and J. C. Sprott, “Bifurcations and chaos in fractional-order simplified Lorenz system,” International Journal of Bifurcation and Chaos, vol. 20, no. 4, pp. 1209–1219, 2010.
23. K. Diethelm, N. J. Ford, and A. D. Freed, “Detailed error analysis for a fractional Adams method,” Numerical Algorithms, vol. 36, no. 1, pp. 31–52, 2004.
24. E. Misirli and Y. Gurefe, “Multiplicative Adams Bashforth-Moulton methods,” Numerical Algorithms, vol. 57, no. 4, pp. 425–439, 2011.
25. S. K. Agrawal, M. Srivastava, and S. Das, “Synchronization of fractional order chaotic systems using active control method,” Chaos, Solitons & Fractals, vol. 45, no. 6, pp. 737–752, 2012.
26. G. Psihoyios and T. E. Simos, “Trigonometrically fitted predictor-corrector methods for IVPs with oscillating solutions,” Journal of Computational and Applied Mathematics, vol. 158, no. 1, pp. 135–144, 2003.
27. P. Zhuang, F. Liu, V. Anh, and I. Turner, “Numerical methods for the variable-order fractional advection-diffusion equation with a nonlinear source term,” SIAM Journal on Numerical Analysis, vol. 47, no. 3, pp. 1760–1781, 2009.
28. Y. F. Xu and Z. M. He, “The short memory principle for solving Abel differential equation of fractional order,” Computer & Mathematics with Applications, vol. 62, no. 12, pp. 4796–4805, 2011.
29. A. Razminia, A. F. Dizaji, and V. J. Majd, “Solution existence for non-autonomous variable-order fractional differential equations,” Mathematical and Computer Modelling, vol. 55, no. 3-4, pp. 1106–1117, 2012.
30. W. C. Chen, “Nonlinear dynamics and chaos in a fractional-order financial system,” Chaos, Solitons & Fractals, vol. 36, no. 5, pp. 1305–1314, 2008.
31. X. S. Zhao, Z. B. Li, and S. Li, “Synchronization of a chaotic finance system,” Applied Mathematics and Computation, vol. 217, no. 13, pp. 6031–6039, 2011.