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 [1–6]. 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 [10–13] 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, 4–6, 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: π·π‘ž1(𝑑)𝐷π‘₯(𝑑)=𝑧(𝑑)+(𝑦(𝑑)βˆ’π‘Ž)π‘₯(𝑑),π‘ž2(𝑑)𝑦(𝑑)=1βˆ’π‘π‘¦(𝑑)βˆ’π‘₯(𝑑)2,π·π‘ž3(𝑑)𝑧(𝑑)=βˆ’π‘₯(𝑑)βˆ’π‘π‘§(𝑑),(1.1) where π·π‘žπ‘–(𝑑) is the Caputo variable-order fractional differential operator whose definition is given in the next section and π‘ž1(𝑑), π‘ž2(𝑑), π‘ž3(𝑑) 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 [16–26]). 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 π‘“ξ…ž(𝑑)=limβ„Žβ†’0𝑓(𝑑+β„Ž)βˆ’π‘“(𝑑)β„Žπ‘“(firstorderdifferenceforderivative),ξ…žξ…ž(𝑑)=limβ„Žβ†’0𝑓(𝑑+β„Ž)βˆ’2𝑓(𝑑)+𝑓(π‘‘βˆ’β„Ž)β„Ž2(secondorderdifferenceforderivative),(2.1) 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 𝑓(𝑛)(𝑑)=limβ„Žβ†’01β„Žπ‘›π‘›ξ“π‘˜=0(βˆ’1)π‘˜βŽ›βŽœβŽœβŽπ‘›π‘˜βŽžβŽŸβŽŸβŽ π‘“(π‘‘βˆ’π‘˜β„Ž).(2.2)

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 π‘š=[π‘ž]+1. Then π·π‘žπ‘“(𝑑)=lim𝑁→+∞1β„Žπ‘žπ‘π‘ξ“π‘˜=0(βˆ’1)π‘˜βŽ›βŽœβŽœβŽπ‘žπ‘˜βŽžβŽŸβŽŸβŽ π‘“ξ€·π‘‘βˆ’π‘˜β„Žπ‘ξ€Έ,(2.3) 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 Ξ“βˆΆ(0,+∞)→ℝ, defined by Ξ“ξ€œ(π‘₯)∢=0+βˆžπ‘‘π‘₯βˆ’1π‘’βˆ’π‘‘π‘‘π‘‘,(2.4) is called Euler's Gamma function (or Euler's integral of the second kind). Particularly, we have Ξ“(𝑛)=(π‘›βˆ’1)! for 𝑛 is a positive integer.

Definition 2.4 (see [11]). The Riemann-Liouville fractional order derivative of 𝑓(𝑑) is defined as π‘Žπ·π‘žπ‘‘1𝑓(𝑑)=𝑑Γ(π‘šβˆ’π‘ž)π‘šπ‘‘π‘‘π‘šξ‚Άξ€œπ‘‘π‘Ž(π‘‘βˆ’πœ)π‘šβˆ’1βˆ’π‘žπ‘“(𝜏)π‘‘πœ,(2.5) where π‘šβˆ’1β‰€π‘ž<π‘š, π‘šβˆˆβ„• provided the right side is pointwise defined on 𝑑>π‘Ž.

Definition 2.5 (see [11, 14]). Let π‘ž>0 be a real number, π‘“βˆˆπ’žπ‘š[0,𝑑] and 0β‰€πœβ‰€π‘‘. Then the Caputo fractional derivative of 𝑓(𝑑) is defined as π·π‘žβŽ§βŽͺ⎨βŽͺ⎩1𝑓(𝑑)=βˆ«Ξ“(π‘šβˆ’π‘ž)𝑑0𝑓(π‘š)(𝜏)(π‘‘βˆ’πœ)π‘ž+1βˆ’π‘šπ‘‘π‘‘πœ,π‘šβˆ’1β‰€π‘ž<π‘š,π‘šπ‘“(𝑑)π‘‘π‘‘π‘š,π‘ž=π‘š,(2.6) 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, π‘“βˆˆπ’žπ‘š[0,𝑇], 𝑑≀𝑇, and π‘š=[max0≀𝑑≀𝑇{π‘ž(𝑑)}]+1. Then π·π‘ž(𝑑)𝑓(𝑑)=limβ„Žπ‘β†’01β„Žπ‘π‘π‘ž(𝑑)ξ“π‘˜=0(βˆ’1)π‘˜βŽ›βŽœβŽœβŽπ‘˜βŽžβŽŸβŽŸβŽ π‘“ξ€·π‘ž(𝑑)π‘‘βˆ’π‘˜β„Žπ‘ξ€Έ,(2.7) with β„Žπ‘=(π‘‘βˆ’0)/𝑁 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 π‘Žπ·π‘‘π‘ž(𝑑)1𝑓(𝑑)=𝑑Γ(π‘šβˆ’π‘ž(𝑑))π‘šπ‘‘π‘‘π‘šξ‚Άξ€œπ‘‘π‘Ž(π‘‘βˆ’πœ)π‘šβˆ’1βˆ’π‘ž(𝜏)𝑓(𝜏)π‘‘πœ,(2.8) where π‘š=[max0≀𝑑≀𝑇{π‘ž(𝑑)}]+1, π‘šβˆˆβ„• provided the right side is pointwise defined on 𝑑>π‘Ž.

Definition 2.8 (see [1]). Let π‘ž(𝑑)>0 be a continuous and bounded function, 𝑓(𝜏)βˆˆπ’žπ‘š[0,𝑑], and 0β‰€πœβ‰€π‘‘. Then π·π‘ž(𝑑)⎧βŽͺ⎨βŽͺβŽ©βˆ«π‘“(𝑑)=𝑑01⋅𝑓Γ(π‘šβˆ’π‘ž(𝜏))(π‘š)(𝜏)(π‘‘βˆ’πœ)π‘ž(𝜏)+1βˆ’π‘šξ‚Όπ‘‘π‘‘πœ,π‘šβˆ’1β‰€π‘ž(𝑑)<π‘š,π‘šπ‘“(𝑑)π‘‘π‘‘π‘š,π‘ž(𝑑)=π‘š(2.9) are called the Caputo variable-order fractional derivative of 𝑓(𝑑), where π‘š=[max0≀𝑑≀𝑇{π‘ž(𝑑)}]+1, [π‘₯] 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 [0,𝑇] and suppose that 𝑓(β‹…) satisfies some smoothness conditions in every subinterval of (0,𝑑) 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: (1) to let the argument of π‘ž be the current time instant 𝑑; (2) 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 [0,π‘₯]. 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.

3. The Adams-Bashforth-Moulton Method

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: π·π‘ž(𝑑)𝑦𝑦(𝑑)=𝑓(𝑑,𝑦(𝑑)),0≀𝑑≀𝑇,(π‘˜)(0)=𝑦0(π‘˜),π‘˜=0,1,…,π‘šβˆ’1.(3.1)

Similar to the fractional order sense, the differential equation (3.1) is equivalent to the following Volterra integral equation [29]: 𝑦(𝑑)=[π‘ž]ξ“π‘˜=0𝑦0(π‘˜)π‘‘π‘˜+ξ€œπ‘˜!𝑑0(π‘‘βˆ’π‘ )π‘ž(𝑠)βˆ’1Ξ“(π‘ž(𝑠))𝑓(𝑠,𝑦(𝑠))𝑑𝑠,(3.2) where [π‘ž] denotes the integer part of π‘ž. For better describing the discretization of (3.2), we select the following equispaced grid: π‘‡β„Ž=𝑁,𝑑𝑛=π‘›β„Ž,(3.3) where 𝑛=0,1,2,…,𝑁. Then (3.2) becomes π‘¦β„Žξ€·π‘‘π‘›+1ξ€Έ=[π‘ž(𝑑𝑛+1)]ξ“π‘˜=0𝑦0(π‘˜)π‘‘π‘˜π‘›+1+β„Žπ‘˜!π‘ž(𝑑𝑛+1)Ξ“ξ€·π‘žξ€·π‘‘π‘›+1𝑓𝑑+2𝑛+1,π‘¦π‘ƒβ„Žξ€·π‘‘π‘›+1+𝑛𝑗=0β„Žπ‘ž(𝑑𝑛+1)π‘Žπ‘—,𝑛+1Ξ“ξ€·π‘žξ€·π‘‘π‘›+1𝑓𝑑+2𝑗,π‘¦β„Žξ€·π‘‘π‘—,ξ€Έξ€Έ(3.4) where the predicted value π‘¦π‘ƒβ„Ž(𝑑𝑗) is determined by π‘¦π‘ƒβ„Žξ€·π‘‘π‘›+1ξ€Έ=[π‘ž(𝑑𝑛+1)]ξ“π‘˜=0𝑦0(π‘˜)π‘‘π‘˜+1π‘˜!Ξ“ξ€·π‘žξ€·π‘‘π‘›+1𝑛𝑗=0𝑏𝑗,𝑛+1𝑓𝑑𝑗,π‘¦β„Žξ€·π‘‘π‘—,ξ€Έξ€Έ(3.5) in which we can deduce the coefficients as π‘Žπ‘—,𝑛+1=⎧βŽͺ⎨βŽͺβŽ©π‘›π‘ž(𝑑𝑛+1)+1βˆ’ξ€·ξ€·π‘‘π‘›βˆ’π‘žπ‘›+1ξ€Έξ€Έ(𝑛+1)π‘ž(𝑑𝑛+1),𝑗=0,(π‘›βˆ’π‘—βˆ’2)π‘ž(𝑑𝑛+1)+1+(π‘›βˆ’π‘—)π‘ž(𝑑𝑛+1)+1βˆ’2(π‘›βˆ’π‘—+1)π‘ž(𝑑𝑛+1)+1𝑏,1≀𝑗≀𝑛,1,𝑗=𝑛+1,𝑗,𝑛+1=β„Žπ‘ž(𝑑𝑛+1)π‘žξ€·π‘‘π‘›+1ξ€Έξ€·(π‘›βˆ’π‘—+1)π‘ž(𝑑𝑛+1)βˆ’(π‘›βˆ’π‘—)π‘ž(𝑑𝑛+1)ξ€Έ.(3.6)

The estimation error of this approximation is described as follows [14]: max𝑗=0,1,…,𝑁||π‘¦ξ€·π‘‘π‘—ξ€Έβˆ’π‘¦β„Žξ€·π‘‘π‘—ξ€Έ||ξ€·β„Ž=𝑂𝛼,(3.7) where 𝛼=min(2,1+π‘žπ‘š), π‘žπ‘šβˆΆ=min0β‰€πœβ‰€π‘‡{π‘ž(𝜏)}. Now we apply the mentioned method to determine the numerical solution of variable-order fractional system. Consider the following variable-order fractional system: π·π‘ž1(𝑑)π‘₯=𝑓1𝐷(π‘₯,𝑦,𝑧),π‘ž2(𝑑)𝑦=𝑓2𝐷(π‘₯,𝑦,𝑧),π‘ž3(𝑑)𝑧=𝑓3(π‘₯,𝑦,𝑧),(3.8) with 0<π‘žπ‘–(𝑑)≀1  (𝑖=1,2,3) and initial condition (π‘₯0,𝑦0,𝑧0). According to the previous method, the system (3.8) can be discretized as follows: π‘₯𝑛+1=π‘₯0+β„Žπ‘ž1(𝑑𝑛+1)Ξ“ξ€·π‘ž1𝑑𝑛+1𝑓+21ξ€·π‘₯𝑃𝑛+1,𝑦𝑃𝑛+1,𝑧𝑃𝑛+1ξ€Έ+𝑛𝑗=0β„Žπ‘ž1(𝑑𝑛+1)𝛾1,𝑗,𝑛+1Ξ“ξ€·π‘ž1𝑑𝑛+1𝑓+21ξ€·π‘₯𝑗,𝑦𝑗,𝑧𝑗,𝑦𝑛+1=𝑦0+β„Žπ‘ž2(𝑑𝑛+1)Ξ“ξ€·π‘ž2𝑑𝑛+1𝑓+22ξ€·π‘₯𝑃𝑛+1,𝑦𝑃𝑛+1,𝑧𝑃𝑛+1ξ€Έ+𝑛𝑗=0β„Žπ‘ž2(𝑑𝑛+1)𝛾2,𝑗,𝑛+1Ξ“ξ€·π‘ž2𝑑𝑛+1𝑓+22ξ€·π‘₯𝑗,𝑦𝑗,𝑧𝑗,𝑧𝑛+1=𝑧0+β„Žπ‘ž3(𝑑𝑛+1)Ξ“ξ€·π‘ž3𝑑𝑛+1𝑓+23ξ€·π‘₯𝑃𝑛+1,𝑦𝑃𝑛+1,𝑧𝑃𝑛+1ξ€Έ+𝑛𝑗=0β„Žπ‘ž3(𝑑𝑛+1)𝛾3,𝑗,𝑛+1Ξ“ξ€·π‘ž3𝑑𝑛+1𝑓+23ξ€·π‘₯𝑗,𝑦𝑗,𝑧𝑗,(3.9) where π‘₯𝑃𝑛+1=π‘₯0+𝑛𝑗=0𝛽1,𝑗,𝑛+1Ξ“ξ€·π‘ž1𝑑𝑛+1𝑓1ξ€·π‘₯𝑗,𝑦𝑗,𝑧𝑗,𝑦𝑃𝑛+1=𝑦0+𝑛𝑗=0𝛽2,𝑗,𝑛+1Ξ“ξ€·π‘ž2𝑑𝑛+1𝑓2ξ€·π‘₯𝑗,𝑦𝑗,𝑧𝑗,𝑧𝑃𝑛+1=𝑧0+𝑛𝑗=0𝛽3,𝑗,𝑛+1Ξ“ξ€·π‘ž3𝑑𝑛+1𝑓3ξ€·π‘₯𝑗,𝑦𝑗,𝑧𝑗,𝛾𝑖,𝑗,𝑛+1=⎧βŽͺ⎨βŽͺβŽ©π‘›π‘žπ‘–(𝑑𝑛+1)+1βˆ’ξ€·π‘›βˆ’π‘žπ‘–ξ€·π‘‘π‘›+1(𝑛+1)π‘žπ‘–(𝑑𝑛+1),𝑗=0,(π‘›βˆ’π‘—βˆ’2)π‘žπ‘–(𝑑𝑛+1)+1+(π‘›βˆ’π‘—)π‘žπ‘–(𝑑𝑛+1)+1βˆ’2(π‘›βˆ’π‘—+1)π‘žπ‘–(𝑑𝑛+1)+1𝛽,1≀𝑗≀𝑛,1,𝑗=𝑛+1,𝑖,𝑗,𝑛+1=β„Žπ‘žπ‘–(𝑑𝑛+1)π‘žπ‘–ξ€·π‘‘π‘›+1ξ€Έξ€·(π‘›βˆ’π‘—+1)π‘žπ‘–(𝑑𝑛+1)βˆ’(π‘›βˆ’π‘—)π‘žπ‘–(𝑑𝑛+1)ξ€Έ.(3.10)

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]: π·π‘ž(𝑑)[],π‘₯(𝑑)=𝑓(𝑑),π‘‘βˆˆ0,1π‘₯(0)=0.005,(4.1) where π‘ž(𝑑)=(2/3)+(1/5)𝑑 and 𝑓(𝑑)≑1.

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.

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.

Example 4.2. Consider the following variable-order fractional financial system: π·π‘ž1(𝑑)𝐷π‘₯(𝑑)=𝑧(𝑑)+(𝑦(𝑑)βˆ’3)π‘₯(𝑑),π‘ž2(𝑑)𝑦(𝑑)=1βˆ’0.1𝑦(𝑑)βˆ’π‘₯(𝑑)2,π·π‘ž3(𝑑)𝑧(𝑑)=βˆ’π‘₯(𝑑)βˆ’π‘§(𝑑).(4.2)

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 π‘ž1=π‘ž2=π‘ž3=1. 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 0.2292. The following Figure 2 shows our numerical simulations.

Case 2 (the fractional order case). First we choose the commensurate fractional order case: π‘ž1=π‘ž2=π‘ž3=0.95. 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 0.1734. Figure 3 shows our numerical simulations.
Furthermore, we consider the incommensurate fractional order case: π‘ž1=0.79, π‘ž2=0.95, π‘ž3=0.99. 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 0.1551. Figure 4 shows our numerical simulations.

Case 3 (the variable-order case). Firstly, we choose π‘ž1(𝑑)=0.76βˆ’(0.01/100)𝑑, π‘ž2(𝑑)=0.88βˆ’(0.06/100)𝑑, and π‘ž3(𝑑)=0.82βˆ’(0.36/100)𝑑 with π‘‘βˆˆ[0,100], 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 π‘ž1(𝑑)=0.99βˆ’(0.01/100)𝑑, π‘ž2(𝑑)=0.85βˆ’(0.01/100)𝑑, and π‘ž3(𝑑)=0.89βˆ’(0.01/100)𝑑 with π‘‘βˆˆ[0,200], 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 π‘ž1(𝑑)=0.99βˆ’(0.01/100)𝑑, π‘ž2(𝑑)=0.99βˆ’(0.02/100)𝑑, and π‘ž3(𝑑)=(1/12)𝑒sin((1/25)𝑑)+0.76 with π‘‘βˆˆ[0,300], we see that the chaotic attractor exists in such financial system. The largest Lyapunov exponent of the system is 0.1340. 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.

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).