An Analysis on the Fractional Asset Flow Differential Equations

Transcription

1 Article An Analysis on the Fractional Asset Flow Differential Equations Din Prathumwan 1, Wannika Sawangtong 1,2 and Panumart Sawangtong 3, * 1 Department of Mathematics, Faculty of Science, Mahidol University, Bangkok 14, Thailand; din.pra@student.mahidol.ac.th D.P.; wannika.saw@mahidol.ac.th W.S. 2 Centre of Excellence in Mathematics, Commission on Higher Education, Ministry of Education, Si Ayutthaya Road, Bangkok 14, Thailand 3 Department of Mathematics, Faculty of Applied Science, King Mongkut s University of Technology North Bangkok, Bangkok 18, Thailand * Correspondence: panumart.s@sci.kmutnb.ac.th; Tel.: Academic Editor: Hari M. Srivastava Received: 3 March 217; Accepted: 12 June 217; Published: 16 June 217 Abstract: The asset flow differential equation AFDE is the mathematical model that plays an essential role for planning to predict the financial behavior in the market. In this paper, we introduce the fractional asset flow differential equations FAFDEs based on the Liouville-Caputo derivative. We prove the existence and uniqueness of a solution for the FAFDEs. Furthermore, the stability analysis of the model is investigated and the numerical simulation is accordingly performed to support the proposed model. Keywords: asset flow differential equations AFDEs; Liouville-Caputo derivative; fixed point theorems; locally asymptotically stable 1. Introduction In recent years, the important thing for financial behavior has been understanding and classifying the dynamics of the financial system such as price crashes, bubbles, momentum and liquidity. Over the last three decades, the mathematical model has played an essential role in preparing planning to predict the financial behavior in the market in the future. Thus, a number of mathematical models have been used to predict financial behavior [1 12. Moreover, in order to obtain better simulation of future phenomena, nonlinear continuous models have also developed. Researchers proposed several nonlinear continuous models for describing the complicated dynamics of the financial system. For example, in 28, Fanti et al. [1 proposed the behavior of capital stock by using a dynamic Investment Saving-Liquidity Preference Money Supply IS-LM model with delay, and Chatterjee et al. [11 proposed a dynamic Hecker-Ohlin model and considered the equilibrium of the model. Moreover, Yang et al. [12 proposed a herd mechanism of an open financial market for considering volatility and parameters of bubbles and collapses in the market. In this paper, we study the analysis of a nonlinear system that is the asset flow differential equations AFDEs. AFDEs were proposed and developed by Gaginalp et al. [13. They investigated the solution of autonomous equations with dimension n 3. This model may describe many behaviors in real financial markets [ The model variables are explained by the demand, the supply, the market price, the transition rate, the trend-based and value-based of the investor preferences at sudden times [13,17,18. Fractional calculus was introduced for over three centuries. At the beginning, fractional calculus was only in a theoretical sense. The traditional system of integer-order of differential equations may fail to explain complicated incidents in real phenomena. However, nowadays, fractional order of differential equations has been applied to model the complicated real situations involving many Mathematics 217, 5, 33; doi:1.339/math5233

2 Mathematics 217, 5, 33 2 of 17 areas for example, physics, engineering, epidemic, finance and sciences [ The concept of a fractional model helps with expressing the data structure for real life situations more than the integer model. The main dominant point of fractional order equations compared to integer order equations is its memory, which is the history of the system [ Since variables in the financial system have long memories, the fractional order system is suitable for studying the behavior in financial systems. There are many papers proposing fractional-order differential equations for the financial system [25 28 consisting of the parts of products, money, bond and labor force. Here, we propose a system of time-fractional differential equations modified from the Gaginalp s model [13 to describe the behavior of asset flow in the financial market. A number of researchers studied about the existence and uniqueness for the solution of fractional differential equations [ However, there has been a lack of research focusing on the global existence and uniqueness for the solution of the fractional asset flow differential equations FAFDEs. In this paper, we propose an alternative way to prove the local existence solution of the system of fractional-order differential equations by the Banach fixed point theorem, which helps to explain the solution of the fractional model. Moreover, the global existence for the solution of the fractional model can be express by continuation theorem. The remainder of this work is organized as follows. In Section 2, the definitions and lemma involving the fractional-order differential equations are presented. In addition, we present the fractionalorder model in this section. In Section 3, the local existence and uniqueness theorem of the model s solution is presented. In Section 4, we propose the continuation theorems for the model to extend the local solution into the global solution. The stability analysis of the fractional model is presented and is followed by the numerical simulation in Section 6. Finally, section 7 concludes and discusses this work. 2. Model Formulation 2.1. Fractional Calculus Definitions of fractional integral and derivatives [19,35 are as follows : Definition 1. The Riemann Liouville integral of order < α < 1 for a function f is defined as: where Γ is the gamma function. RL I α f t = 1 t s α 1 f sds, 1 Definition 2. The Riemann Liouville derivative of order < α < 1 for a function f :, R is defined as: RL Dt α 1 d t f s f t = ds. 2 Γ1 α dt t s α The regularization definition of the Riemann Liouville derivative was introduced by Liouville Caputo in 1967 [36. Definition 3. The Liouville Caputo type fractional derivative of order < α < 1 for a function f :, R is defined by: Dt α 1 t f f t = s ds. 3 Γ1 α t s α The most regular definition is the Liouville Caputo-type definition because the initial condition for fractional order differential equations with Liouville Caputo-type derivatives is similar in form to the integer-order differential equation [19 22.

3 Mathematics 217, 5, 33 3 of 17 Lemma 1. Banach fixed point theorem Let A be a closed convex bounded subset of a Banach space X and V be a contraction mapping from A into A. Then, there exists a unique z A such that Vz = z Fractional Model The behavior of price dynamic in the market based on: x 1 t the market price of an asset at time t, x 2 t the fraction of total asset invested in the stock at time t, x 3 t the trend-based component of the investor preference at time t, x 4 t the value-based component of the investor preference at time t, P a D S the fundamental value, the demand, the supply. Gaginalp et al. [13 proposed the dynamical system of AFDEs to describe those four state variables in the market price : dx 1 dt dx 2 dt dx 3 dt dx 4 dt = x 1 t F D, S = kt1 x 2 t 1 ktx 2 t + x 2 t1 x 2 tf D = a 1 a 2 F x 3 t, S P a x = a 3 a 1 t 4 x P 4 t, a D, S where x 1 = x 1,, x 2 = x 2,, x 3 = x 3,, x 4 = x 4, and x i, are positive constants for i = 1, 2, 3, 4. The value of the market price of an asset at time t can be only a positive number and cannot be approached to infinity, the value of the fraction of total assets is defined between and 1. The demand D and the supply S can be determined as D = kt1 x 2 t and S = 1 ktx 2 t, respectively. kt is the transition rate that is dependent on x 3 t and x 4 t defined by kt = tanhx 3 t + x 4 t, Ft is an increasing continuous function with F1 =, and P a, a 1, a 2, a 3 and a 4 are positive constants. In this article, the FAFDEs based on the Liouville Caputo fractional derivative are introduced as follows: 4 Dt α x 1 = x 1 tf Gt, x 1 t, x 2 t, x 3 t, x 4 t, Dt α x 2 = Kt, x 1 t, x 2 t, x 3 t, x 4 t x 2 t + x 2 t1 x 2 tf Gt, x 1 t, x 2 t, x 3 t, x 4 t, Dt α x 3 = b 1 FGt, x 1 t, x 2 t, x 3 t, x 4 t b 2 x 3 t, Dt α x 4 = b 3 b 4 x 1 t b 5 x 4 t, 5 where D α t is the Liouville Caputo fractional derivative of order < α < 1, b j are positive constants for j = 1, 2, 3, 4, 5 and x i = x i, are positive constants for i = 1, 2, 3, 4. We need the following assumptions; let X be a Banach space with its norm X. A1 F : [, R is increasing and continuously differentiable where F1 =, and satisfies that for any p, q [,, Fp Fq L 1 p q for some positive constant L 1.

4 Mathematics 217, 5, 33 4 of 17 A2 G :, X 4 [, is continuous and satisfies there being two constants L 2 and L 3 such that: 4 Gt, x 1, x 2, x 3, x 4 Gt, y 1, y 2, y 3, y 4 L 2 x i t y i t i=1 A3 and sup t, Gt,,,, L 3 for any t >. K :, X 4, is continuously differentiable and has the following : there exists two positive constants L 4 and L 5 such that: 4 Kt, x 1, x 2, x 3, x 4 Kt, y 1, y 2, y 3, y 4 L 4 x i t y i t i=1 and sup t, Kt,,,, L 5 for any t >. 3. Local Existence and Uniqueness of a Solution In this section, we show that system 5 has a unique continuous solution by the Banach fixed point theorem. To do so, we introduce a Banach space: C[, T = {u : [, T R such that u is continuous} and its norm is defined by: with T being a positive constant. u C[,T = sup ut, t [,T Theorem 1. There exists a finite time h such that system 5 has a unique continuous solution on [, h. Proof. Let M be a constant with M 4 max{ x i,, i = 1, 2, 3, 4} + 1. Choose a positive constant h where: h < min { 1 1 Γα + 1 α Γα + 1 α,, 4M FL 3 + ML 1 L 2 4L L 4 M + M1 + M[ FL 3 + ML 1 L Γα + 1 α Γα + 1 α,, 4b 1 [ FL 3 + ML 1 L 2 + b 2 M b 3 + b 4 + b 5 M 1 } Γα + 1 α 1 + b 2 + b 4 + b 5 + L M FL 3 + 3M M + b 1 L 1 L 2 6 Let X 4 = {x = x 1, x 2, x 3, x 4 x i X for i = 1, 2, 3, 4} with X = C[, h. Then, X 4 is a Banach space equipped with the norm x X 4 = x 1 X + x 2 X + x 3 X + x 4 X for any x = x 1, x 2, x 3, x 4. We construct the set U = {x X 4 : x X 4 M}. Clearly, U X 4 is closed, bounded and convex.

5 Mathematics 217, 5, 33 5 of 17 We also define the operator B : U X 4 by Bx = B 1 x, B 2 x, B 3 x, B 4 x where B i : U X with i = 1, 2, 3, 4 and B i are given by: B 1 x = x 1, + 1 t s α 1 [x 1 sfgs, x 1 s, x 2 s, x 3 s, x 4 s ds, B 2 x = x 2, + 1 t s α 1 [Ks, x 1 s, x 2 s, x 3 s, x 4 s x 2 s + x 2 s1 x 2 sfgs, x 1 s, x 2 s, x 3 s, x 4 sds, B 3 x = x 3, + 1 B 4 x = x 4, + 1 t s α 1 [b 1 FGs, x 1 s, x 2 s, x 3 s, x 4 s b 2 x 3 s ds, t s α 1 [b 3 b 4 x 1 s b 5 x 4 s ds. Let x, y U. We want to show that B : U U. Consider B 1 xt x 1, + 1 t s α 1 x 1 s FGs, x 1 s, x 2 s, x 3 s, x 4 s ds x 1, + M x 1, + M x 1, + t s α 1 [ 4 FGs,,,, + L 1 L 2 x i s i= t s α 1 [ FL 3 + L 1 L 2 M ds M [ FL 3 + L 1 L 2 M hα α. It follows from the definition of h 6 that: Consider B 1 x X x 1, + 1/4. 7 B 2 xt x 2, + 1 t s α 1 [L 5 + L 4 M + M + M1 + M FGs, x 1 s, x 2 s, x 3 s, x 4 s ds x 2, + 1 [L 5 + L 4 M + M + M1 + M [ FL 3 + L 1 L 2 M hα α. ds By the inequality 6, B 2 x X x 2, + 1/4. 8 Consider B 3 xt x 3, + 1 t s α 1 [b 1 [ FL 3 + L 1 L 2 M + b 2 Mds x 3, + 1 [b 1 [ FL 3 + L 1 L 2 M + b 2 M hα α. Then, B 3 x X x 3, + 1/4. 9

6 Mathematics 217, 5, 33 6 of 17 Consider B 4 xt x 4, + 1 t s α 1 [b 3 + b 4 M + b 5 Mds x 4, + 1 [b 3 + b 4 + b 5 M hα α. Then, B 4 x X x 4, + 1/4. 1 Therefore, from the inequalities 7 1 and the definition of h as in the inequality 6: Bx X 4 = B 1 x X + B 2 x X + B 3 x X + B 4 x X 4max{ x i,, i = 1, 2, 3, 4} + 1 M. Then, this means that B maps U into itself. Next, we want to show that B is a contraction mapping. Consider B 1 xt B 1 yt 1 t s α 1 x 1 sfgs, x 1 s, x 2 s, x 3 s, x 4 s y 1 sfgs, y 1 s, y 2 s, y 3 s, y 4 s ds 1 t s α 1 [ x 1 s FGs, x 1 s, x 2 s, x 3 s, x 4 s FGs, y 1 s, y 2 s, y 3 s, y 4 s + x 1 s y 1 s FGs, y 1 s, y 2 s, y 3 s, y 4 s ds 1 t s α 1 [ML 1 L 2 x y X 4 + x y X 4 [ FL 3 + L 1 L 2 Mds 1 [2ML 1L 2 + FL 3 hα α x y X 4 = L 6 x y X 4, where L 6 = 2ML 1L 2 + FL 3 h α. Γα + 1 Thus, B 1 xt B 1 yt L 6 x y X 4. 11

7 Mathematics 217, 5, 33 7 of 17 Consider B 2 xt B 2 yt 1 t s α 1 Ks, x 1 s, x 2 s, x 3 s, x 4 s Ks, y 1 s, y 2 s, y 3 s, y 4 s + y 2 s x 2 s + x 2 sfgs, x 1 s, x 2 s, x 3 s, x 4 s x2 2sFGs, x 1s, x 2 s, x 3 s, x 4 s y 2 sfgs, y 1 s, y 2 s, y 3 s, y 4 s + y 2 2 sfgs, y 1s, y 2 s, y 3 s, y 4 s ds 1 t s [ Ks, α 1 x 1 s, x 2 s, x 3 s, x 4 s Ks, y 1 s, y 2 s, y 3 s, y 4 s y 2 s x 2 s + x 2 s y 2 s FGs, x 1 s, x 2 s, x 3 s, x 4 s + x 2 s FGs, x 1 s, x 2 s, x 3 s, x 4 s FGs, y 1 s, y 2 s, y 3 s, y 4 s + x2 2s FGs, x 1s, x 2 s, x 3 s, x 4 s FGs, y 1 s, y 2 s, y 3 s, y 4 s + x2 2s y2 2 s FGs, y 1s, y 2 s, y 3 s, y 4 s ds 1 t s [L α 1 4 x y X 4 + y 2 s x 2 s + x 2 s y 2 s [ FL 3 + L 1 L 2 M + x 2 s L 1 L 2 x y X 4 + x 2 2 s [L 1L 2 x y X 4 + x 2 s + y 2 s x 2 s y 2 s [ FL 3 + L 1 L 2 M 1 t s [L α 1 4 x y X 4 + x y X 4 + x y X 4 [ FL 3 + L 1 L 2 M + ML 1 L 2 x y X 4 + M 2 L 1 L 2 x y X 4 + 2M x y X 4 [ FL 3 + L 1 L 2 M ds 1 1 t s α 1 x y X 4 ds [ L ML 1 L 2 + 3M 2 L 1 L M FL 3 ds [ L ML 1 L 2 + 3M 2 L 1 L M FL 3 = L 7 x y X 4, h α α x y X 4 where L 7 = 1 + L 4 + 2ML 1 L 2 + 3M 2 L 1 L M FL 3 h α. Γα + 1 Thus, B 2 xt B 2 yt L 7 x y X Consider B 3 xt B 3 yt 1 t s [b α 1 1 FGs, x 1 s, x 2 s, x 3 s, x 4 s FGs, y 1 s, y 2 s, y 3 s, y 4 s + b 2 x 3 s y 3 s ds where L 8 = b 1L 1 L 2 + b 2 h α. Γα t s [b α 1 1 L 1 L 2 x y X 4 + b 2 x y X 4 ds b 1L 1 L 2 + b 2 = L 8 x y X 4, h α α x y X 4

8 Mathematics 217, 5, 33 8 of 17 Thus, B 3 xt B 3 yt L 8 x y X Consider B 4 xt B 4 yt 1 t s b α 1 4 x 1 s y 1 s + b 5 x 4 s y 4 s ds b 4 + b 5 h α α x y X 4 = L 9 x y X 4, where L 9 = b 4 + b 5 h α Γα + 1. Thus, B 4 xt B 4 yt L 9 x y X We have that from the inequalities and the definition of h as in the inequality 6: Bx By X 4 L 6 + L 7 + L 8 + L 9 x y X 4, 15 with < L 6 + L 7 + L 8 + L 9 < 1. This shows that B is a contraction mapping. Therefore, by Banach fixed point theorem, system 5 has a unique continuous solution on [, h. 4. Continuation Theorems In this section, we extend the continuation theorem for the system of FAFDEs. Firstly, we give the following lemma. Lemma 2. Gronwall-Bellman inequality [37 Suppose γ >, at is a nonnegative function that is locally integrable on t < T, and gt is a nonnegative, nondecreasing continuous bounded function defined on t < T. If ut is nonnegative and locally integrable on t < T with: on this interval, then: ut at + gt ut at + gt t s γ 1 usds [ n=1 Furthermore, if at is nondecreasing on t < T, then: where E γ z is the Mittag Leffler function defined by: gtγγ n t s nγ 1 as ds. Γnγ ut ate γ gtγγt γ, E γ z = k= z k Γkγ + 1 for z >. Theorem 2. If x 1 and x 2 are bounded on [,, then system 5 has a solution xt for t [,. Proof. Suppose that solution xt of system 5 has the maximal existence interval defined by [, β with β <. Let us consider: for t [, β,

9 Mathematics 217, 5, 33 9 of 17 x 4 t x 4, + 1 t s α 1 b 3 b 4 x 1 s b 5 x 4 s x 4, + 1 ds t s [b α b 4 x 1 s + b 5 x 4 s ds. By the assumption of theorem, x 1 t is bounded on [,, we have that x 1 t C 1 for t [, β and then for t [, β, x 4 t x 4, + 1 t s [b α b 4 C 1 + b 5 x 4 s ds x 4, + b 3 + b 4 C 1 Γα + 1 βα + b 5 It follows from the Lemma 2 that, for t [,, x 4 t t s α 1 x 4 s ds. [ x 4, + b 3 + b 4 C 1 b5 Γα + 1 βα E α tα. It is well-known that the E α is an analytic function [38 and then we can conclude the x 4 is bounded on [, β. Next, we consider x 3 on [, β. Then, for t [, β, x 3 t x 3, + 1 t s α 1 b 1 FGs, x 1 s, x 2 s, x 3 s, x 4 s b 2 x 3 s ds [ x 3, + 1 t s α 1 4 b 1 FL 3 + L 1 L 2 x i s i= + b 2 x 3 s x 3, + 1 t s [b α 1 1 FL 3 + L 1 L 2 [ x 1 s + x 2 s + x 4 s + b 1 FL 3 + L 1 L 2 x 3 s + b 2 x 3 s ds. ds Since x 1, x 2 and x 4 are bounded on [, β, C 3 > such that sup t [,β x 1 t + x 2 t + x 3 t C 3. Then; x 3 t x 3, t s [b α FL 3 + L 1 L 2 C 3 + b 1 L 1 L 2 + b 2 x 3 s ds t s α 1 b 1 L 1 L 2 + b 2 x 3 s ds x 3, + b 1 2 FL 3 + L 1 L 2 C 3 β α + 1 Γα + 1 By Lemma 2, we obtain: [ x 3 t x 3, + b 1 2 FL 3 + L 1 L 2 C 3 β α Γα + 1 t s α 1 b 1 L 1 L 2 + b 2 x 3 s ds. b1 L E 1 L 2 + b 2 α t α. Thus, x 3 is bounded on [, β. Hence, xt is bounded on [, β. Since x are continuous and bounded on [, β,lim t β xt is finite. Let lim t β xt = xβ with xβ = x 1 β, x 2 β, x 3 β, x 4 β. Then, xt C 4 [, β. Consider t [β, β + 1. Let M > 4max{x i β, i = 1, 2, 3, 4} + 1.

10 Mathematics 217, 5, 33 1 of 17 Similar to Theorem 1, we can construct the positive constant h depending on M and h < 1 such that our problem has a solution on [β, β + h. This part is shown in the below. Let X = C[β, β + h be a Banach space with the norm y X = sup yt and t [β,β+ h X 4 = {y = y 1, y 2, y 3, y 4 y i X for i = 1, 2, 3, 4}. Then, X 4 is a Banach space equipped with the norm y X 4 = y 1 X + y 2 X + y 3 X + y 4 X for any y = y 1, y 2, y 3, y 4. { } We construct the set Ũ = y X 4 : y X 4 M. We define the operator B : Ũ X 4 by By = B 1 y, B 2 y, B 3 y, B 4 y for y Ũ and for any y Ũ, B i is constructed by: B 1 y = x 1 β + 1 t s α 1 y 1 sfgs, y 1 s, y 2 s, y 3 s, y 4 sds, β where x 1 β = x 1, + 1 β β s α 1 y 1 sfgs, y 1 s, y 2 s, y 3 s, y 4 sds, B 2 y = x 2 β + 1 t s α 1 [Ky 1 s, y 2 s, y 3 s, y 4 s y 2 s β + y 2 s1 y 2 sfgs, y 1 s, y 2 s, y 3 s, y 4 sds, β and x 2 β = x 2, + 1 β s α 1 [Ky 1 s, y 2 s, y 3 s, y 4 s y 2 s +y 2 s1 y 2 sfgs, y 1 s, y 2 s, y 3 s, y 4 sds, B 3 y = x 3 β + 1 t s α 1 [b 1 FGs, y 1 s, y 2 s, y 3 s, y 4 s b 2 y 3 sds, β and x 3 β = x 3, + 1 β β s α 1 [b 1 FGs, y 1 s, y 2 s, y 3 s, y 4 s b 2 y 3 sds, B 4 y = x 4 β + 1 t s α 1 [b 3 b 4 y 1 s b 5 y 4 sds, β β and x 4 β = x 4, + 1 β s α 1 [b 3 b 4 y 1 s b 5 y 4 sds. To make it clear about the construction of h, without loss of generality, we only show that B 1 maps Ũ into Ũ. Let us consider that for y Ũ: B 1 y x 1 β + 1 t s α 1 y 1 s FGs, y 1 s, y 2 s, y 3 s, y 4 s ds x 1 β + 1 M x 1 β + β β t s α 1 [ FGs,,,, + L Mds [ β + h α βα. α α M [ FL 3 + L M In this case, the constant h is determined by: [ 1 α α h < 4 M[ FL 3 + L M + βα β. We can construct h in this way and follow method as in Theorem 1.

11 Mathematics 217, 5, of 17 By considering as the Theorem 1, B has a fixed point x, that is xt = B xt for t [β, β + h. We next define the function x on [, β + h by: x t = { xt, xt, t [, β, t [β, β + h. We can see that x is a continuous solution of our problem on [, β + h with β + h > β. We thus get a contradiction. Therefore, system 5 has a solution on [,. Remark 1. Because x 1 and x 2 are the market price and the fraction of total asset invested, respectively. By assumptions of the model [13, it is possible to assume that x 1 and x 2 are bounded. 5. Stability Analysis In order to determine the stability of system 5, we firstly consider the equilibrium point. To find the equilibrium point C = x1, x 2, x 3, x 4, we set all equations of system 5 to zero : 16 D α t x 1 =, D α t x 2 =, D α t x 3 =, D α t x 4 =. 17 where: Then: C 1 = x 1,1, x 1,2, x 1,3, x 1,4, x1,2 x = K 1, x 2,, b 3 b 4 x1, b 5 x 1,3 =, x 1,4 = b 3 b 4 x 1 b 5. C 2 = x 2,1, x 2,2, x 2,3, x 2,4, where: x2,1 =, 1 F + 4F K x2,2 = x 2,3 = b 1 b 2 F, x 2,4 = b 3 b 5., x2, b 1 b F, b 3 2 b F 2 2F, C 3 = x 3,1, x 3,2, x 3,3, x 3,4,

12 Mathematics 217, 5, of 17 where: x3,1 =, 1 F 4F K x3,2 = x 3,3 = b 1 b 2 F, x 3,4 = b 3 b 5,, x2, b 1 b F, b 3 2 b F 2 2F, with 4F K, x2, b 1 b F, b 3 2 b F 2 and F = FGt, x j,1, x j,2, x j,3, x j,4 for j = 1, 2, 3. It is easy to see that the Jacobian matrix of system 5 is: where: F x j,1 F x 2 x j,1 F x 3 x j,1 F x 4 JCj = K x1 x j,2 2 F x1 µ j K x3 x j,2 2 F x3 K x4 x j,2 2 F x4 b 1 F x1 b 1 F x2 b 1 F x3 b 2 b 1 F x4, 18 b 4 b 5 µ j = K x2 1 + x j,2 F x 2 + F x j,2 2 F x2 2x j,2 F, F xi = K xi = [ FGt, x1, x x i t 2, x 3, x 4 Cj, [ Kt, x1, x x i t 2, x 3, x 4 Cj, for i = 1, 2, 3, 4 and j = 1, 2, 3. The characteristic equation of JCj for j = 1, 2, 3 is: where: λ 4 + ρ 1 λ 3 + ρ 2 λ 2 + ρ 3 λ + ρ 4 =, 19 ρ 1 = x j,2 2 F x2 + 2x j,2 1F x j,2 F x 2 b 1 F x3 K x2 + b 2 + b 5 + 1, ρ 2 = x j,2 2 F x2 x 1 F x1 F + b 2 + b 5 + x j,2 b 1 F x2 F x3 2b 1 F x3 F 2[F 2 + F x2 F + 2b 2 F + b 5 F b 2 + b 5 F x2 + x j,1 b 4 F x4 F x2 K x1 b 1 F x1 F x3 + K x2 2b 2 + b 5 b 1 F x3 + F 1 F + b 1 F x3 K x2 b 1 b 5 F x3 b 1 F x2 K x3 b 1 F x3 b 2 K x2 b 5 K x2 + b 2 b 5 + b 2 + b 5,

13 Mathematics 217, 5, of 17 ρ 3 = x j,2 2 b 2 + b 5 x j,1 F x 1 F x2 b 2 + b 5 F x2 F + b 2 b 5 F x2 + x j,1 x j,2 2b 4 F x4 [2F F x2 + b 1 F x1 F x3 [F x2 2F + x j,1 b 1 F x1 F x3 F + b 1 F x1 F x3 [K x2 b 5 b 4 F x4 F b 1 F x1 F x2 K x3 b 2 + b 5 F x2 K x1 b 1 F x1 F x3 + b 4 F x2 K x4 + b 4 F x4 [b 2 K x x j,2 [2F F x2 b 1 F x3 F b 1 b 5 F x3 2F + b 1 b 5 F x3 F x2 2b 2 + b 5 F 2 + b 2 + b 5 F x2 F + 2b 2 b 5 [F F x2 + b 1 b 5 [F x3 K x2 F x2 K x3 b 1 b 5 F x3 + b 2 b 5 [1 + K x2 + b 1 [F x2 K x3 + b 5 K x2 F x3 F x3 F F + b 2 + b 5 [F 2 + [b 1 F x3 + b 2 + b 5 K x2 2b 2 b 5 b 2 b 5 F, ρ 4 = x j,1 x j,2 2[b 2 b 4 F x4 b 1 b 5 F x1 F x3 F + b 5 F x1 F x2 [b 1 F x3 + b 2 x j,2 b 2b 4 F x2 F x4 + x j,1 b 1 b 5 [F x1 F x3 K x2 + F F x1 F x2 K x3 F x3 + b 2 b 4 [F x2 K x4 F x4 K x2 F x4 F + b 2 [b 4 F x4 b 5 F x2 K x1 + x j,2 b 5 F x1 F x2 [b 2 b 2 x 2 b 1 F x3 + 2b 5 [b 1 F x3 b 2 F 2 + b 1 b 5 [F x2 K x3 + F x3 1 K x2 + F x3 F F + b 2 b 5 [K x2 1 + F F. In order to obtain main results, we need the following proposition, referred to in [39,4. Proposition 1. The equilibrium points are locally asymptotically stable if all the eigenvalues λ of the Jacobian matrix evaluated at the equilibrium points satisfy: argλ > απ 2. 2 Theorem 3. If ρ 4 >, ρ 1 >, ρ 1 ρ 2 > ρ 3 and ρ 1 ρ 2 ρ 3 = ρ 2 1 ρ 4 + ρ 2 3, then the equilibrium point is locally asymptotically stable. Proof. From the equations 18 and 19, we can find eigenvalues of the Jacobian matrix corresponding to the equilibrium points from the fourth-order characteristic equation. Under assumptions of this theorem, Ahmed et al. [41 show that argλ > απ/2 for all eigenvalues λ. From Proposition 1, we can conclude that the equilibrium points of system 5 is locally asymptotically stable. 6. Numerical Simulation The exact solution of some fractional order differential equations can not be found. Thus, the numerical solutions for fractional order differential equation model 5 can be approached by using an algorithm based on an Adams Bashforth type predictor-corrector method [42,43. The numerical solutions of system 5 are carried out by setting functions Ft = γ log t and Kt, x 1 t, x 2 t, x 3 t, x 4 t = tanh x 3 t + x 4 t. The model parameters are chosen to be γ = 1.1, b 1 =.87, b 2 =.1, b 3 =.27, b 4 =.4 and b 5 =.45, with the initial conditions x 1 = 6.5, x 2 =.5, x 3 = and x 4 =.15. The approximate solutions of market price x 1, fraction of total asset x 2, trend-based component x 3 and value-based component x 4 of the FAFDEs with the different values of order, α = 1,.9 and.8 are shown in Figures 1 4, respectively.

14 Mathematics 217, 5, of = 1. =.9 = x 1 t t Figure 1. Market price x 1 t for order α = 1.,.9 and = 1. =.9 = x 2 t t Figure 2. Fraction of total asset x 2 t for order α = 1.,.9 and = 1. =.9 =.8.25 x 3 t t Figure 3. Trend-based component x 3 t for order α = 1.,.9 and.8.

15 Mathematics 217, 5, of = 1. =.9 =.8.8 x 4 t t Figure 4. Value-based component x 4 t for order α = 1.,.9 and.8. The results in Figures 1 3 indicate that the behavior for the market price of an asset, the fraction of total asset invested in the stock, and the trend-based component of the investor preference have the same pattern at each time for an integer order system with α = 1. and a fractional order system with α =.9 and α =.8. To investigate the effect of fractional derivative order on market behavior, the numerical results demonstrate that a decrease in the derivative order is associated with a decrease in the maximum value of x 1, x 2 and x 3. Figure 4 shows that decreasing of derivative order leads to a decrease in the value-based component at each time. In addition, the market price as shown in Figure 1 represents the pattern of head and shoulders or rounding top, in which a major high appears by a local maximum. In order to execute planning in the future, investors can buy or sell when the pattern completes. By Theorem 3, the equilibrium point C = x1, x 2, x 3, x 4 is locally stable, where x1 = , x 2 =.56, x 3 = and x 4 =.121. In other words, the values of x 1, x 2, x 3 and x 4 converge to equilibrium point C at some long-term behaviors. 7. Discussion and Conclusions In this paper, we have proposed a fractional order asset flow differential equations model, which is a generalization of an integer order model proposed by Gaginalp [13. Then, we propose and prove the local existence and uniqueness theorems for the solution of the model by applying the Banach fixed point theorem. After that, we extend the local solution, t, h, to a global solution, t,, by proposing the continuation theorem for the model. The stability of the model is also investigated and we obtain the sufficient condition of the parameters. The numerical solutions of the model are illustrated by using an Adams Bashforth type predictor-corrector method. Acknowledgments: This research was funded by King Mongkut s University of Technology North Bangkok, contract no. KMUTMB-GEN Appreciation is extended to Faculty of Science, Mahidol University, Bangkok, Thailand. Author Contributions: Din Prathumwan and Panumart Sawangtong developed the mathematical model. Wannika Sawangtong and Panumart Sawangtong proposed the research idea for proving theorems. Wannika Sawangtong and Din Prathumwan were responsible for simulating the results by the MATLAB programming. Din Prathumwan wrote this paper. All authors contributed in editing and revising the manuscript. Conflicts of Interest: The authors declare no conflict of interest. References 1. Ma, J.H.; Chen, Y.S. Study for the bifurcation topological structure and the global complicated character of a kind of nonlinear finance system I. Appl. Math. Mech. 21, 22,

16 Mathematics 217, 5, of Ma, J.H.; Chen, Y.S. Study for the bifurcation topological structure and the global complicated character of a kind of nonlinear finance system II. Appl. Math. Mech. 21, 22, Cesare, L.D.; Sportelli, M. A dynamic IS-LM model with delayed taxation revenues. Chaos Solitons Fractals 25, 25, Fanti, L.; Manfredi, P. Chaotic business cycles and fiscal policy: An IS-LM model with distributed tax collection lags. Chaos Solitons Fractals 27, 32, Frisch, R.; Holme, H. The Characteristic Solutions of a Mixed Difference and Differential Equation Occuring in Economic Dynamics. Econometrica 1935, 3, He, X.Z.; Li, K.; Wei, J.; Zheng, M. Market stability switches in a continuous-time financial market with heterogeneous beliefs. Econ. Model. 29, 26, Chiarella, C.; Dieci, R.; Gardini, L. Asset price dynamics in a financial market with fundamentalists and chartists. Discret. Dyn. Nat. Soc. 21, 2, He, X.Z.; Zheng, M. Dynamics of moving average rules in a continuous-time financial market model. J. Econ. Behav. Organ. 21, 76, Chiarella, C.; He, X.Z.; Zheng, M. An analysis of the effect of noise in a heterogeneous agent financial market model. J. Econ. Dyn. Control 211, 35, Zhou, L.; Li, Y. A generalized dynamic IS-LM model with delayed time in investment processes. Appl. Math. Comput. 28, 196, Chatterjee, P.; Shukayev, M. A stochastic dynamic model of trade and growth: Convergence and diversification. J. Econ. Dyn. Control 212, 36, Yang, H.; Li, L.; Wang, D. Research on the stability of open financial system. Entropy 215, 17, Caginalp, G.; Balenovich, D. Asset flow and momentum: Deterministic and stochastic equations. Philos. Trans. R. Soc. 1999, 357, Caginalp, G.; Ermentrout, G.B. A kinetic thermodynamics approach to the psychology of dluctuations in financial market. Appl. Math. Lett. 199, 3, Merdan, H.; Alisen, M. A mathematical model for asset pricing. Appl. Math. Comput. 211, 218, Sullivan, A.; Sheffrin, S.M. Economics: Principles in Action; Prentice Hall: Upper Saddle River, NJ, USA, Caginalp, G.; Balenovich, D. Market oscillations induced by the competition between value-based and trend-based investment strategies. Appl. Math. Financ. 1994, 1, Caginalp, G.; Ermentrout, B. Numerical studies of differential equations related to theoretical financial markets. Appl. Math. Lett. 1991, 4, Podlubny, I. Fractional Differential Equations. Mathematics in Science and Engineering; Academic Press: San Diego, CA, USA, 1999; Volume Area, I.; Batarfi, H.; Losada, J.; Nieto, J.J.; Shammakh, W.; Torres, A. On a fractional order Ebola epidemic model. Adv. Differ. Equ. 215, 215, Sofuoglu, Y.; Ozalp, N. A fractional rrder model on bilingualism. Commun. Fac. Sci. Univ. Ank. Ser. A 214, 214, Sofuoglu, Y.; Ozalp, N. Fractional order bilingualism model without conversion from dominant unilingual group to bilingual group. Differ. Equ. Dyn. Syst. 215, 215, Hilfer, R. Applications of Fractional Calculus in Physics; World Scientific: Hackensack, NJ, USA, Mainardi, F. Fractional relaxation-oscillation and fractional diffusion-wave phenomena. Chaos Solitons Fractals 1996, 7, Dadras, S.; Momeni, H.R. Control of a fractional-order economical system via sliding mode. Phys. A Stat. Mech. Appl. 21, 389, Wang, Z.; Huang, X.; Shi, G. Analysis of nonlinear dynamics and chaos in a fractional order financial system with time delay. Comput. Math. Appl. 211, 62, Chen, W.C. Nonlinear dynamics and chaos in a fractional-order financial system. Chaos Solitons Fractals 28, 36, Li, Q.; Zhou, Y.; Zhao, X.; Ge, X. Fractional order stochastic differential equation with application in european option pricing. Discret. Dyn. Nat. Soc. 214, 214, Li, C.; Sarwar, S. Existence and contiuation of solution for Caputo type fractional differential equations. Electron. J. Differ. Equ. 216, 216, 1 14.

17 Mathematics 217, 5, of Babakhani, D.; Gejji, V.D. Existence of positive solutions of nonlinear fractional differential equations. J. Math. Anal. Appl. 23, 278, Diethlm, K.; Ford, N.J. Analysis of fractional differential equations. J. Math. Anal. Appl. 22, 265, Yu, C.; Gao, G. Existence of fractional differential equations. J. Math. Anal. Appl. 25, 31, Kou, C.; Zhou, H.; Li, C.P. Existence and continuation theorems of fractional Riemann Liouville type fractional differential equations. Int. J. Bifurc. Chaos 212, 22, Deng, W.H.; Li, C.P.; Guo, Q. Analysis of fractional differential equations with multi orders. Fractals 27, 15, Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations; Elsevier Science: Amsterdam, The Netherlands, Caputo, M. Linear model of dissipation whose Q is almost frequency independent-ii. Geophys. J. Int. 1967, 13, Gronwall, T.H. Note on the derivatives with respect to a parameter of the solutions of a system of differential equations. Ann. Math. 1919, 2, Gorenflo, R.; Kilbas, A.A.; Mainardi, F.; Rogosin, S.V. Mittag-Leffer Function, Related Topics and Applications; Springer: New York, NY, USA, Tavazoei, M.S.; Haeri, M. Unreliability of frequency-domain approximation in recognising chaos in fractional-order systems. IET Signal Proc. 27, 1, Petráš, I. Stability of fractional-order systems with rational orders: A survey. Fract. Calc. Appl. Anal. 29, 12, Ahmed, E.; El-Sayed, A.M.A.; El-Saka, H.A.A. On some Routh Hurwitz conditions for fractional order differential equations and their applications in Lorenz, Rössler, Chua and Chen systems. Phys. Lett. A 26, 358, Diethelm, K.; Ford, N.J.; Freed, A.D. Detailed error analysis for a fractional Adams method. Numer. Algorithms 24, 36, Diethelm, K. An algorithm for the numerical solution of differential equations of fractional order. Electron. Trans. Numer. Anal. 1997, 5, 1 6. c 217 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution CC BY license