Dynamic Keynesian Model of Economic Growth with Memory and Lag

A mathematical model of economic growth with fading memory and continuous distribution of delay time is suggested. This model can be considered as a generalization of the standard Keynesian macroeconomic model. To take into account the memory and gamma-distributed lag we use the Abel-type integral and integro-differential operators with the confluent hypergeometric Kummer function in the kernel. These operators allow us to propose an economic accelerator, in which the memory and lag are taken into account. The fractional differential equation, which describes the dynamics of national income in this generalized model, is suggested. The solution of this fractional differential equation is obtained in the form of series of the confluent hypergeometric Kummer functions. The asymptotic behavior of national income, which is described by this solution, is considered.


Introduction
Advanced mathematical methods of fractional calculus [1][2][3][4][5] are a powerful tool for describing the fading memory and spatial non-locality.Fractional derivatives and integrals of non-integer order have different applications in natural and social sciences [6,7].
In this article, we suggest a generalization of one of the most famous models of economic growth, which is associated with the founder of modern macroeconomic theory, John M. Keynes [8][9][10].In the suggested generalization, we take into account two types of phenomena: (I) long memory with power-law fading and (II) continuously distributed lag with gamma distribution of delay time.
The continuously distributed lag has been considered in economics starting with the works of Michal A. Kalecki [11] and Alban W.H. Phillips [12,13].The macrodynamic models of business cycles, where the continuous uniform distribution of delay time is used, were considered by Michal A. Kalecki in 1935 [11], (see also Section 8.4 of [14], (pp.251-254)).The economic growth models with continuously distributed lag were proposed by Alban W.H. Phillips [12,13] in 1954.In his works, the distribution of delay time has been described by exponential distribution.Operators with continuously distributed lag were considered by Roy G.D. Allen [14] (pp.[23][24][25][26][27][28][29], in 1956. The time delay (lag) is caused by finite speeds of processes, i.e., the change of one variable does not lead to instant changes of another variable.Therefore, the distributed lag (time delay) cannot be considered as a long memory in processes.For example, in physics, the retarded potential of the electromagnetic field is well known.The change in the value of the electromagnetic field at the point of observation r 1 is delayed with respect to the change in the sources of the field located at the point r 2 at the time t = |r 1 r 2 | /c, where c is the speed of propagation of disturbances.It is known that the processes of propagation of the electromagnetic field in a vacuum do not mean the presence of memory in this process.
An application of advanced mathematical methods of fractional calculus in Keynesian economic models with continuous time was proposed by authors [30,31] in 2016 (see also [32][33][34]).The fractional differential equations of the dynamic Keynesian model with power-law memory and their solutions have been considered in [30][31][32][33][34]. Continuously distributed lag was not discussed in these works.
For macroeconomics, it is important to simultaneously take into account lagging and memory phenomena.In this article, we consider memory with power-law fading and lag with gamma distribution of delay time.The memory is described by the Riemann-Liouville fractional integrals and the Caputo fractional derivatives.The distributed lag is described by the translation T τ , in which the delay time τ > 0 is considered as a random variable that is distributed by probability law (distribution) on positive semiaxis.The composition of these operators is represented as the Abel-type integral and integro-differential operators with the confluent hypergeometric Kummer function in the kernel.Using these operators, we propose the fractional differential equation for the generalized dynamic Keynesian model that describes the fractional dynamics of national income.We obtain a solution for this equation that describes the macroeconomic growth with power-law fading memory and gamma distribution of delay time.

Standard Dynamic Keynesian Model
In macroeconomic growth models, two types of variables are used [14,35,36].First, exogenous variables are considered as independent quantities that are external to the considered economic model.Secondly, endogenous variables are internal variables that are formed within the model.The endogenous variables are described as functions of exogenous variables.In models with continuous time, all these variables are considered as functions of time t.
Let us consider the standard dynamical Keynesian model with continuous time.In the Keynesian model, the following variables are used to describe the dynamics of the revenue and expenditure parts of the economy: Y(t) is a national income; G(t) is the government expenditure; C(t) describes the consumption expenditure; I(t) describes the investment expenditure, E(t) is a total expenditure, i.e., E(t) is defined as the sum of all expenditures: In dynamic equilibrium, we have In this case, the balance equation establishes the equality of the national income to the sum of all expenditures In the Keynesian model, it is assumed that the consumption expenditure in period t depends on the income level in the same period.The consumption expenditure C(t) is regarded as an endogenous variable equal to the amount of domestic consumption of some part of the national income and final consumption independent of income.As a result, the consumption expenditure C(t) is described by the linear equation of the economic multiplier where m(t) is the multiplier factor that describes the marginal propensity to consume (0 < m(t) < 1), and the function b(t) > 0 describes the autonomous consumption that does not depend on income.The expression m(t)Y(t) describes the part of consumption that depends on income.
In the static model, the investment expenditure and government expenditure are considered as exogenous variables.In the dynamic Keynesian model, the investment expenditure I(t) is treated as endogenous and it is assumed to depend on the level of income [35], (pp.95-97).The investment expenditure I(t) is determined by the rate of change of the national income.This assumption is described by the equation of the economic accelerator where v(t) is the rate of acceleration, which characterizes the level of technology and state infrastructure, and Y (1) (t) = dY(t)/dt is the first-order derivative of the income function Y(t) with respect to the time variable.
In the Keynesian model, government expenditure G(t), the propensity to consume m(t), the rate of acceleration v(t), and the autonomous consumption b(t) are exogenous variables that are specified as external to the model and characterize the functioning and development of the economy.These variables, as functions of time, are assumed to be given.
The purpose of the dynamic Keynesian model is to describe the behavior of the national income.For this, it is necessary to find the national income Y(t), as a function of time t.Substituting the multiplier Equation (4) and the accelerator Equation (5) into the balance Equation (3), we obtain This equation can be written in the form Equation (7) of the dynamic Keynesian model is a non-homogeneous linear differential equation with a first-order derivative.We see that the functions G(t) and b(t) are included in Equation ( 7) as a sum.This could be expected since G(t) is an independent expenditure on investment, that is independent of income, and b(t) is an independent expenditure on consumption, also not dependent on national income.From the point of view of the main purpose of the Keynesian model, which is to describe the dynamics of the national income, these two types of expenditure simply complement each other [35], (pp.95-98).Therefore, it is convenient to use the sum which describes the independent expenditure.In this case, the consumption function is that part of consumption that depends on income.All this allows us to write down the equation of the standard dynamic Keynesian model in the form where 0 < 1 − m(t) < 1, and v(t) > 0. Equation ( 9) is a non-homogeneous first-order differential equation that describes the standard dynamic Keynesian model, which does not take into account the effects of memory and delay.

Dynamic Keynesian Model with Memory
In the standard Keynesian model, Equation ( 9) implies an instantaneous change in the investment expenditure, when the rate of growth of national income changes.This means that the equation of this model does not take into account the effects of memory and delay.Mathematically, this is due to the fact that the standard model equation is a first-order differential equation.The derivative of the first order, which is used in accelerator Equations (5), implies an instantaneous change of the investment expenditure I(t), when changing the rate of the national income Y(t).Because of this, accelerator Equation ( 5) does not take into account memory and lag.Multiplier Equation (4) also assumes that the consumption expenditure C(t) changes instantly when the national income changes.As a result, model Equation ( 9) can describe an economy, in which agents have no memory.This fact greatly limits the applicability of the standard model to describe the real processes in the economy.To expand the scope of the model, we should take into account that economic agents can remember the history of changes of the national income and the investment expenditure, because it affects the behavior of these agents.
Generalization of the standard Keynesian model, in which memory [37,38] is taken into account, was proposed by the authors [30,31].Let us briefly describe this generalized model with memory.The equation of investment accelerator with memory [37,39] can be written as where M(t − τ) is the memory function.Note that Equation (10) can also be used to describe the distributed lag.In this case, M(t − τ) is called the weighting function, which is interpreted as the probability density function.Equation ( 5) can be obtained from (10) in the case M(t − τ) = δ(t − τ) and n = 1.Substitution of the investment I(t) in the form of expression (10), and the consumption expenditure (4) into balance Equation (3), we get where G b (t) is defined by (8).Equation ( 9) of the standard Keynesian model without memory and lag can be obtained from ( 11) by using M(t − τ) = δ(t − τ) and n = 1.Equation ( 11) describes the fractional dynamics of the national income within the framework of the Keynesian model of growth with memory.If the parameters m(t) and v(t) are given, the growth of national income Y(t) is conditioned by the behavior of the independent expenditure (8).
The memory with one-parameter power-law fading is described [37,39] by the function where Γ(α) is the gamma function and n − 1 < α ≤ n.Using ( 12), the accelerator with memory ( 10) is represented [37,39] as where is the Caputo fractional derivative [4,5].In general, the rate of acceleration v(t) depends on the parameter of memory fading, i.e., v(t) = v(t, α).The parameter α > 0 is interpreted as a fading parameter of power-law memory [37].The concept of an accelerator with memory [39] allows us to get the equation of the Keynesian model with power-law memory in the form of the fractional differential equation where the Caputo fractional derivative can be represented by the Laplace convolution where n = [α] + 1 for α / ∈ N and n = α for α ∈ N, and the function Y(τ) has integer-order derivatives Y (j) (τ), j = 1, . . ., (n − 1), that are absolutely continuous.

Memory and Lag by Abel-Type Integral and Derivative with Kummer Functions
The economic accelerator and multiplier with continuously distributed lag were proposed by Alban W.H. Phillips [12,13] in 1954 (see also Sections 3.4, 3.5 and 8.7 in [14]).The distribution of delay time has been described by the exponential distribution.In 1956, the operators with continuously distributed lag were considered by Roy G.D. Allen in the book [14], (pp.[23][24][25][26][27][28][29].In the general case, the distribution of delay time can be described by other probability distributions [40], not just exponential distributions.Note that the time delay is caused by finite speeds of processes.Therefore, the distributed lag (time delay) cannot be interpreted as a memory.

Fractional Integral with Memory and Lag
The translation operator T τ is defined [1], (pp.95-96), by the expression (T τ Y)(t) = Y(t − τ ), where τ > 0 is the delay time.In the general case, the delay time τ > 0 can be considered as a random variable, which is distributed by probability law (distribution) on positive semiaxis [40].The translation operator with the continuously distributed delay time can be defined [40] by the equation where we assume that Y(t) and M T (t) are piecewise continuous functions on R and the integral In Equation (17), the kernel M T (τ) is the weighting function that satisfies the condition To take into account the distributed time delay and power-law fading memory, we can use a composition of the translation operator (17) and integration of non-integer order.The Riemann-Liouville fractional integral with a distributed time delay can be defined [40] by the equation where (I α RL,t 0 Y)(t) is the Riemann-Liouville fractional integral of the order α ∈ R + , and M T (τ) is the weighting function that satisfies the conditions (18).Here we can assume that (I α RL,t 0 Y)(t) and M T (t) are piecewise continuous functions on R such that RL,t 0 is defined [4], (p.92), by the expression where α ≥ 0 is the order of the fractional integral, Γ(α) is the gamma function, and τ ∈ [t 0 , t].However, in order for definition (19) to be correct, it needs to take into account that the fractional integral ( 20) is defined for the case t > t 0 [4], (p.69).By virtue of this, we can consider the two following cases.
For t 0 = −∞, Equation (20) defines the Liouville fractional integral (I α L,+ Y)(t), where τ ∈ (−∞, t].Then the Liouville fractional integration with continuously distributed lag is defined in the form Considering that integral (20) is defined for the case t > t 0 = 0 [4], (p.69), we can assume that (I α RL,0+ Y)(t) = 0 for t < t 0 = 0 in expression (20) with t 0 = 0.Here the function Y(τ) is defined on the finite interval [t 0 , t 1 ] with 0 = t 0 < t 1 < ∞, [4], (p.69).This allows us to use the upper limit t > 0 insteat of infinity in Equation (19), such that As a result, we can define the fractional integration with the gamma-distributed lag in the form where M T (τ) = M λ,a T (τ) is the probability density function of the gamma distribution with the shape parameter a > 0 and the rate parameter λ > 0. If a = 1, the function (24) describes the exponential distribution.
In economics, the gamma distribution ( 24) is applied to take into account waiting times, when there is a sharp increase in the average delay time.For example, this distribution is used to describe delays orders in queues, delays in payments, and to take into account the likelihood of risk events.The distribution also describes the time of receipt of the order for the enterprise, the service life of device components, and time between store visits.
Using Equation ( 29), the memory-and-lag function ( 27) can be written as that defines the kernel of operator (28).In general, the function (31) can be interpreted as a new memory function.Note that equality F 1,1 (a; c; z) = Γ(c)E a 1,c (z) (see Equation 5.1.18 of [43]) allows us to represent the memory kernel (31) through the three parameter Mittag-Leffler functions [43].
As a result, the Riemann-Liouville fractional integral with gamma distribution of delay time can represented [40] by the equation where α > 0 is the order of integration and the parameters a > 0, λ > 0 describe the shape and rate of the gamma distribution, respectively.It is known (see Section 37 of [2] and [44]), the integral operators of the form are called the Abel-type integral operators.For example, we can consider the confluent hypergeometric Kummer function F 1,1 (β; α; λ(t − τ)) as the kernel K(t, τ) of the operator (33).It is known that the Abel-type (AT) fractional integral operator with Kummer function in the kernel (see equation 37.1 in [1], (p.731), and [45]) is defined by the equation In paper [45], the integral operator ( 34) is denoted as K 0 (β, α, λ).As a result, the Riemann-Liouville fractional integral with gamma-distributed lag (32) can be expressed through the AT fractional integral (34) by the equation The AT fractional integral (34) can be represented as an infinite series of the Riemann-Liouville fractional integrals where (β Since the Riemann-Liouville fractional integrals (20) are bounded in L p (t 0 , t 1 ), where p ≥ 1, t 1 < ∞ (see the proof of Theorem 2.6 in [1], (pp.48-51)), then the series (37) may be summed for |λ| < I 1 RL,a+ −1 L p (t 0 ,t 1 ) .After evaluating the sums one may remove this restriction on λ and on the sum (37), since the suggested fractional integral operators (32) are analytic functions with respect to λ (see [1], (p.732)).
Using Equation (37) and the semigroup property of the Riemann-Liouville fractional integrals, we can get the semigroup property for the Riemann-Liouville fractional integral with gamma-distributed lag (32) where is the beta function.Equality (38) directly follows from Equation 37.14 of [1], (p.733).Using Theorem 37.1 of [1], (p.733), and Equation (35), we can state that the suggested fractional integral with lag (32) has the same range in L p (t 0 , t 1 ) as the Riemann-Liouville fractional integrals and it is bounden from L p (t 0 , t 1 ) onto I α RL,0+ [L p (t 0 , t 1 )] ⊂ L p (t 0 , t 1 ).
Using the condition of the invertibility of the AT operators (34), which is described by Theorem 37.2 of [1], (p.736), and Equation 37.32 of [1], (p.735), we get that the solution of the fractional integral Equation can be represented in the form where D α RL,0+ and D a RL,0+ are the Riemann-Liouville fractional derivatives of orders α > 0 and a > 0 respectively.These derivatives can be defined by the Laplace convolution as , where D n t = d n /dt n , n ∈ N. Note that Equation ( 40) can be interpreted as an equation of economic multiplier with power-law memory and distributed lag [39], which is a generalization of the multiplier Equation ( 4), where m(t) is the multiplier factor.In this case, Equation ( 41) can be considered as the equation of an economic accelerator with memory [39].The equation of multiplier with memory is a reversible, such that the dual (inverse) equation describes an accelerator with memory (see Section 4 of [39]).

Fractional Derivative with Memory and Lag
Using the integral operator (32), we can define the fractional derivatives with continuously distributed lag [40].For example, the Caputo fractional derivative with gamma-distributed lag is defined by the Laplace convolution in the form The convolution is an associative operation that allows us to write where M λ,a;n−α TRL (t) is defined by Equation (31).This allows us to represent operators (42) in the form where n − 1 < α ≤ n.Using Equation ( 31), we have the representation of the kernel M λ,a;n−α As a result, the Caputo fractional derivative with gamma-distributed lag is represented [40] by the Equation where n − 1 < α ≤ n.Fractional differential operator (46) can be expressed through the Riemann-Liouville fractional integral (32) with gamma-distributed lag in the form where n − 1 < α ≤ n.
Using the Laplace transform of the Caputo fractional derivative and the gamma distribution function, we get [40] the Laplace transform of the Caputo fractional derivative with gamma-distributed lag in the form where n − 1 < α ≤ n.
Let us note that the operator ( 42) with a = 1 describes the Caputo fractional derivative with exponentially distributed lag [40].Using the fact that the Caputo fractional derivatives with integer values α = n ∈ N are the integer-order derivatives (D α C,0+ Y)(t) = d n Y(t)/dt n , the operator ( 42) with a = 1 and α = n describes the integer-order derivatives with the exponential distribution [40].Note that these operators with exponentially distributed lag were defined in the works of Caputo and Fabrizio [46,47], where they have been misinterpreted as fractional derivatives of non-integer orders.We can state [40] that the derivative of integer order with exponentially distributed lag coincides with the Caputo-Fabrizio operator of the order β = n − 1/(λ + 1), where λ is the rate parameter of the distribution and n = [β] + 1.As a result, the Caputo-Fabrizio operator can be interpreted as an integer-order derivative with the exponentially distributed lag.
The proposed operator (42) in the form ( 46) can be interpreted as a new generalized operator with the memory function given by the confluent hypergeometric function.The generalized fractional derivative (46) can be used to simultaneously account of long memory with power-law fading and distributed lag with the gamma distribution of delay time.In the next section, we describe a macroeconomic model with memory and distributed lag by using this proposed operator.

Fractional Differential Equation of a Keynesian Model with Memory and Lag
Let us take into account that the relationship between the investment expenditure I(t) and the national income Y(t) depends on memory and lag effects.For the case of the power-law memory and gamma distribution of the delay time, we can use the generalized accelerator equation where D λ,a;α T;C;0+ is the Caputo fractional derivative with distributed lag given by (46).Equation ( 49) describes the economic accelerator that takes into account the power-law fading memory and the gamma-distributed lag.
Substituting expressions (49) and ( 4) into balance Equation (3), we obtain the fractional differential equation of the Keynesian model with power-law memory and gamma distribution of delay time in the form where α is the parameter of memory fading, a > 0 is the shape parameter and λ > 0 is the rate parameter of the gamma distribution of the delay time.
Let us consider the case where v(t) and m(t) are constant quantities.Then the Keynesian model with one-parameter power-law memory and gamma-distributed lag is described by the fractional differential equation where The general solution of Equation ( 51) can be written as where Y 0 (t) is the solution of Equation ( 51) with F(t) = 0, i.e., the homogeneous equation and Y F (t) is particular solution of (51) that can be represented in the form where G α [t − τ] is the generalized Green function [4], (p.281,295).Equation ( 54) yields the solution Y F (t) for Equation (51) with initial conditions, Y (j) (0) = 0 for all j = 0, . . ., (n − 1).
Theorem 1.The fractional differential Equation where D λ,a;α T;C;0+ is the fractional derivative of order α > 0 with gamma-distributed lag, in which a > 0 and λ > 0 are the shape and rate parameters of the gamma distribution respectively, has the solution where n = [α] + 1, and S γ α,δ [µ, λ|t] is the special function that is defined by the expression where F 1,1 (a; b; z) is the confluent hypergeometric Kummer function (30).
Proof.The first step is to find a solution for the homogeneous Equation (53).Using the Laplace transform of Equation ( 53), we get Then we can write where µ = ω λ −a .Using Equation 5.4.9 of [48] in the form where Re(c − a) > 0, we get [40] the Laplace transform of the function (57) (Theorem 1) as Using Equation ( 61) the solution of the homogenous fractional differential Equation ( 53) has the form where [ω λ −a , λ t] is defined by Equation (57).The second step is to find a particular solution (54) of Equation (55).The Laplace transform of Equation ( 55) with conditions Y (j) (0) = 0 for all j = 0, . . .(n − 1) gives the expression that can be rewritten in the from The equality where µ = ω λ −a , gives Using ( 61) with δ = a, γ = α, we have As a result, we obtain that describes the particular solution of Equation (55).Substitution of ( 62) and ( 68) into (52) gives (56).This ends the proof.
As a result, the solution of Equation ( 55) of the Keynesian model with one-parameter power-law memory and gamma-distributed lag for constant v(t) = v and m(t) = m is described by the expression where G b (τ) = G(τ) + g(τ) and ω = (1 − m)/v.

Asymptotic Behavior of National Income Growth with Memory and Lag
In economic theory, the important concept of the Harrod's warranted rate of growth [14] is used, which is also called the technological growth rate.The warranted growth rate describes the growth when the following two conditions are satisfied.The first condition is the constancy of the structure of the economy.This condition means that the parameters of the model do not change over time.In the Keynesian model, we should consider the parameters (t), b(t), and m(t) as constant quantities.The second condition is the absence of external influences.This condition means the absence of exogenous variables.In the Keynesian model, we should consider the case of the absence of the independent expenditure, i.e., G b (t) = 0. Mathematically the warranted growth rate can be obtained from the asymptotic expression of the solution of the homogeneous differential equation of the macroeconomic model.In the standard Keynesian model, the solution of Equation ( 9) with G b (t) = 0 has the form Y(t) = Y(0) exp(ω t).Therefore, the warranted growth rate of this model is described by the value ω = (1 − m)/v, when G b (t) = 0.
The Keynesian model with memory was suggested by authors [30,31] in 2016 (see also [32][33][34]).For this growth model, the fractional differential equation, its solution, and properties have been described.We proved [32,33] that the warranted growth rate with memory is equal to the value ω eff (α) = ω 1/α , where α > 0 is a parameter of power-law memory fading and ω is the rate of growth without memory (α = 1).The warranted growth rates of models with memory do not coincide with the growth rates of the standard Keynesian model.The memory effects can significantly change the growth rates of the economy [32][33][34] and lead to new types of behavior for the same parameters of the economic model.The principles of changing growth rates by memory have been proposed in [32,33].The memory effects can both increase and decrease the warranted growth rates in comparison with the standard Keynesian model.For the memory fading with α < 1, we get a slowdown in the growth and decline of the economy.We can state that memory with α < 1 leads to inhibition of economic growth or decline, i.e., we have stagnation of the economy if α < 1.For the memory fading with α > 1 we have an improvement in the economy.In this case, the memory effect leads either to the slowdown in the decline rate or to the replacement of the decline with growth, or to the increase in the growth rate.
To consider the warranted growth rate of national income for the Keynesian model with memory and distributed lag, we should obtain an asymptotic behavior of the solution (62) of homogenous model Equation ( 53).This solution is expressed by the function S γ α,δ [µ, λ|t] that is represented as an infinite sum (57) of the confluent hypergeometric Kummer function F 1,1 (a; c; z).We can use the asymptotic expression of the function F 1,1 (a; c; z) at infinity z → −∞ that is given in [4] (p. 29), in the form where z = −λt < 0. Therefore the asymptotic expression for t → ∞ is Equation ( 71) leads to the asymptotic expression at infinity ( t → ∞ ) of the function (57) in the form where αk + γ = 0, 1, 2, . . .for all integer k.Solution ( 62) is expressed by the function S α−j−1 α,a [µ, λ|t], where µ = ω λ −a .Then using Equation (72), we get the asymptotic expression For the gamma distribution with the integer shape parameter (a = m ∈ N), which is called the Erlang distribution, expression (73) gives where e −iπm = (−1) m is used.Using Equation 5.1.18 of [43], (p.99), in the form F 1,1 (a; c; z) = Γ(c)E a 1,c (z), we can get asymptotic expressions ( 72) and (74) by using asymptotic expressions of the three parameter Mittag-Leffler functions E δ(k+1) 1,δ(k+1)−αk−γ (−λt) (for example, see [49]).Using Equation (74) we can see that the series (74) can be represented through the two-parameter Mittag-Leffler function with the negative first parameter [50] as This allows us to state that the function S α−j−1 α,m [ω λ −a , λ|t] in the long time limit leads to a series, which can be interpreted as a two-parameter Mittag-Leffler function, which is analogous to the property of the three-parameter Mittag-Leffler function described in [51].
Equation (74) allows us to give the asymptotic expression As a result, the asymptotic behavior of solution (69) with G b (t) = 0 can be described by the Equation that characterizes warranted growth of national income, which is represented by solution (69) with G b (τ) = 0 and a = m ∈ N. Equation (77) allows us to state that the warranted growth has the power-law form with the power −α + j, where j ∈ {0, . . . ,n − 1} is the smallest value at which Y (j) (0) = 0.As a result, the warranted growth of the national income with memory and distributed lag has the power-law type instead of the exponential type of growth with memory without time delay [32,33], where warranted growth rate is ω eff (α) = ω 1/α .Therefore we can state that the distributed lag (time delay) suppresses the effects of fading memory.

Conclusions
The standard Keynesian model [8][9][10]14] describes the dynamics of national income in the absence of long memory and distributed time delay.The Keynesian model with power-law memory has been suggested by authors [30,31].The effects of continuously distributed lag are not considered in [30,31].In this paper, we generalize the Keynesian model with memory by taking into account gamma distribution of delay time.To take into account the distributed lag, we use the operators that are compositions of the translation operator with distributed delay time and the fractional derivatives or integrals.These operators allow us to take into account the memory and lag in the economic accelerator.These operators are the Abel-type integro-differential operators with the confluent hypergeometric Kummer function in the kernel.The solution of the suggested fractional differential Equation, which describes the fractional dynamics of national income, has been suggested.The asymptotic behavior of economic processes with memory and distributed lag demonstrates power-law growth.In the absence of delays, the processes with fading memory demonstrate exponential growth.The warranted growth rate with memory [30][31][32][33][34] is equal to the value ω eff (α) = ω 1/α , where α > 0 is a memory fading parameter.Therefore effects of long memory can significantly accelerate the growth rate of the economy by several orders of magnitude [30][31][32][33].Fading memory can lead to an increase in the growth rate in processes without lag [30][31][32][33][34].The appearance of distributed lag does not accelerate growth due to the memory effect (see also [52,53]).Moreover, we can state that the lag can suppress the effect of fading memory.The distributed lag leads to slower growth.We assume that the suggested approach and model can be used for economic growth modeling by analogy with the computer simulation of the economy in [54][55][56][57][58].It should also be noted that fractional differential Equations have been applied to describe power-law memory in continuous-time finance [59][60][61][62][63][64][65][66][67][68][69][70][71][72].This fact allows us to assume that the proposed approach can be used to take into account the continuously distributed time delay and memory in financial processes with the waiting-time distribution [60,62].