Non-Linear Macroeconomic Models of Growth with Memory

In this article, two well-known standard models with continuous time, which are proposed by two Nobel laureates in economics, Robert M. Solow and Robert E. Lucas, are generalized. The continuous time standard models of economic growth do not account for memory effects. Mathematically, this is due to the fact that these models describe equations with derivatives of integer orders. These derivatives are determined by the properties of the function in an infinitely small neighborhood of the considered time. In this article, we proposed two non-linear models of economic growth with memory, for which equations are derived and solutions of these equations are obtained. In the differential equations of these models, instead of the derivative of integer order, fractional derivatives of non-integer order are used, which allow describing long memory with power-law fading. Exact solutions for these non-linear fractional differential equations are obtained. The purpose of this article is to study the influence of memory effects on the rate of economic growth using the proposed simple models with memory as examples. As the methods of this study, exact solutions of fractional differential equations of the proposed models are used. We prove that the effects of memory can significantly (several times) change the growth rate, when other parameters of the model are unchanged.


Introduction
The Solow model is a dynamic single-sector model of economic growth (see Solow articles [1,2] and books [3][4][5]). In this model, the economy is considered without structural subdivisions, i.e., as a one-sector model. It is assumed that the economy produces only universal products, which can be consumed both in the non-production and production consumptions, i.e., can be consumed, C(t), or invested, I(t), (see [4] p. 25). Exports and imports are not taken into account. This model describes the capital accumulation, labor or population growth, and increases in productivity, which is commonly called the technological progress. The Solow model can be used to estimate the separate effects on economic growth of capital, labor, and technological change.
The Solow model is a generalization of the Harrod-Domar model, which includes a productivity growth as new effect. This relatively simple growth model was independently proposed by Robert M. Solow in 1956 [1,2]. In 1987 Solow was awarded the Nobel Memorial Prize in Economic A more detailed overview of economic models with memory is presented in work [39] (see also book [40] pp. 5-32 and references therein).
However, non-linear economic models with fading memory and exact solutions of fractional differential equations of such models have not been investigated. This article proposes two non-linear models of economic growth with memory. The suggested models are generalizations of two well-known models, which are proposed by two Nobel laureates in economics, Robert M. Solow and Robert E. Lucas in works [1,41], respectively. The first proposed model is a generalization of the model of long-run growth, which is considered Robert M. Solow in paper "A contribution to the theory of economic growth" [1]. In the standard model, the capital depreciation is neglected. In the proposed generalization, we take into account the power-law memory. Moreover, we assume the power-law form of the dynamics of the labor and knowledge. The exact analytical solution of the non-linear fractional differential equation, which describes the proposed model, is suggested. The second proposed model is the standard growth model for closed economy without capital depreciation, which is considered by Robert E. Lucas in paper "Making a Miracle" [41], (see also [42][43][44][45]). In the generalized model, we take into account the memory with power-law fading [46]. The suggested model with memory is described by the non-linear fractional differential equation. The explicit expression of analytical solution of this non-linear equation is obtained.
The aim of the work is not to search for memory effects in economics, but to provide mathematical proof of the possibility of the existence of (at least) acceleration in growth rates, when the parameter of memory fading changes, while other parameters of the process remain unchanged.
The main hypothesis used in this article is the assumption of the power-law character of memory fading, which is described by a single fading parameter. The assumption allows us to use fractional calculus and the fractional differential equations [27][28][29][30][31][32], where the non-integer order of the derivatives and integrals is described by this fading parameter.
We should note that the power-law form of memory fading is distinguished by the fact that this form can be considered as a first approximation for wide class of memory functions. In paper [47], using the generalized Taylor series in the Trujillo-Rivero-Bonilla form, we proved that the memory function for a wide class of functions can be represented through the sum of one or several power-law memory functions. As a result, the integr-differential equations of economic models with general memory functions approximately can be written [47] through the Riemann-Liouville fractional integrals and the Caputo fractional derivatives of non-integer orders [29][30][31][32]48].
In this article, simple non-linear economic models with memory are investigated by using the obtained exact solutions of the fractional differential equations of the proposed models. The main aim of this article is to study the possible influence of memory effects on the rate of economic growth using the proposed simple models. As the methods of this study, we use exact solutions of equations, which describe economic growth with power-law memory. Expressions of the warranted rates of growth with memory are suggested and compared with the warranted rates of standard models. We prove that the effects of memory can significantly (several times) change the growth rate, when other parameters of the model are unchanged.
Note that non-linear growth models with continuous time, which take into account memory effects, have not previously been proposed in the modern literature. The exact solutions of the equations of these models and expressions for growth rates which take into account the influence of memory effects have not been proposed before this article.
The article is organized as follows. In Section 2, an economic model with memory is proposed, which is a generalization of the well-known Solow model, considered by Robert M. Solow in paper "A contribution to the theory of economic growth" [1]. In Section 3, we proposed an economic model with memory, which is a generalization of the well-known Solow model, considered by Robert E. Lucas in paper "Making a Miracle" [41]. In each section, the standard models of economic growth without memory are first considered, then the equations of the generalized models with memory are derived, solutions of these equations are obtained, expressions of the warranted rate of growth with memory are derived, and then compared with the growth rate without memory. In the last paragraphs of Sections 2 and 3, three principles are formulated that make it possible to qualitatively describe the changes arising in the behavior of the economic process, when taking into account the effects of power memory. Finally, brief conclusions of the article are drawn.

Solow Model of Long-Run Growth with Memory
In this section, we consider a generalization of the model of long-run growth, which is considered in Solow's paper "A contribution to the theory of economic growth" [1]. In this model, the capital depreciation is neglected. In the proposed generalization, we take into account the power-law memory. Moreover, we assume the power-law form of the dynamics of the labor and knowledge. The exact analytical solution of the non-linear fractional differential equation, which describes the proposed model, is suggested.

Long-Run Growth without Memory
Let us describe the standard economic model of growth, which is considered in Solow's paper [1] pp. 66-67. The Solow model uses four variables: Output Y(t), capital K(t), labor L(t), and knowledge A(t). At any time, the economy has some amounts of capital, labor, and knowledge, which are combined to produce output. The production function is considered in the following form: The Solow model uses assumptions about how the stocks of labor L(t), knowledge A(t), and capital K(t) change over time. The initial levels of capital, labor, and knowledge are taken as given, and are assumed to be strictly positive. In the standard model, the labor and knowledge growth are described [3] p. 13, by the following equations: where ρ and g are exogenous parameters and f (1) (t) denotes a first derivative with respect to time. The solutions of Equations (2) and (3) can be represented in the following form: Let us consider a power-law generalization of Equations (2) and (3) in the following form: where q and p are exogenous parameters. For q = 1 and p = 1, Equations (5) and (6) give the standard Equations (2) and (3).
Let us obtain solutions of Equations (5) and (6) for the case q 1 and p 1. Equation (5) can be written as: Using the standard chain rule, we get: The solution of Equation (8) has the following form: where c = L 1−q (0). For simplicity, we will assume that L(0) = 0. Then, c = 0 and Equation (9) is given as: where t > 0 and ρ (1 − q) > 0. For t = t 0 > 0, Equation (10) has the following form: Using expression (11), we can represent solution (10) by the following expression: Analogously, we get: where we assume that A(0) = 0. Then, we represent solution (13) in the following form: Output Y(t) is divided between the consumption C(t) and investment I(t). In the standard model, the fraction s ∈ [0, 1] of the output devoted to investment is exogenous and constant such that: Equation (15) describes the economic multiplier without memory and lag. It is known that depreciation (and disposal of fixed assets) can be described as a truncated exponential distribution of the time lag (up to a numerical factor), which can be represented as a differential equation of the first order ( [49] p. 20, [50] p. 86, [51] pp. 252-253 (see also [52] pp. [25][26][27]).
In a general case, the depreciation can be described by the following equation: where D(t − τ) is the kernel of the integral operator that describes the depreciation. The kernel D(t − τ) characterizes the share of fixed assets put into the operation at the time and continuing to operate at time t > τ.
In economics, depreciation is usually described by the kernel in the following exponential form: Therefore, the standard form of depreciation and the exponential distributed lag up to a numerical factor can be considered as a truncated mathematical expectation value with an exponential distributed delay time in the following form: Equation (18) can be represented in the form of the ordinary differential equation of the first order. Let us consider the derivative of the first order for Equation (18). Then, we get: As a result, Equation (18) gives: Substituting expression (20) into Equation (21), that is, taking into account that the capital depreciates at a rate of δ ∈ [0, 1], we obtain for the standard model the following differential equation: Let us consider the model without capital depreciation, i.e., δ = 0. Then, Equation (21) takes the following form: The substitution of (1) into Equation (22) gives: where L(t) and A(t) are described by Equations (12) and (14). The substitution of (12) and (14) into Equation (23) gives: where q 1 and p 1. The solution of Equation (24), which describes a model without memory, will be obtained as a special case of the generalized model, when the memory fading parameter is equal to one which means the absence of memory [7,33].

Long-Run Growth with Power-Law Memory
We should note that Equation (22) cannot take into account the memory, since the derivative K (1) (t) is determined by the behavior of Y(τ) only at the same time instant τ = t, and does not take into account the history of changes of Y(τ) in the past, when τ ∈ (0, t).
To take into account the changes of Y(τ) in the past, instead of Equation (15), we can consider the equation of the multiplier with memory in the following form: where the function s(t − τ) allows us to take into account the changes of output Y(τ) in the past τ ∈ (0, t).
The function s(t − τ) is interpreted as a memory function. In order to have standard dimensions of economic quantities, the time variable t is considered as a dimensionless parameter, i.e., we change the variable t old → t new = t old /t c , where t c is a characteristic time of processes.
The sequential action of the memory effects, which are described by Equation (25) and the depreciation, which is described by Equation (17) can be represented in the following form: where we use the associativity of the Laplace convolution * and the kernel: This representation allows us to interpret Equation (26) as an equation describing some generalized form of memory and/or depreciation.
If the integral operator in (26) is a fractional integral operator, then for such an operator a left inverse operator, which is a fractional derivative must exist [48]. We should emphasize that the main property of any generalized (fractional) derivative is to be a left-inverse operator to the corresponding generalized (fractional) integral operator. This requirement is important for a self-consistent mathematical theory of the fractional operators to have a general fractional calculus of these operators.
For example, if the integral operator in (26) is the Riemann-Liouville fractional integral I α RL,0+ of the order α > 0, then Equation (26) can be rewritten in the following form: Considering the action of the Caputo fractional derivative D α C,0+ and taking into account that this derivative is the left inverse to the integral operator I α RL,0+ , we obtain the following equation: The parameter α > 0 is the order of the fractional derivative or integral. This parameter is interpreted as a memory fading parameter. For integer values of α, these fractional operators take the form of standard integrals and derivatives of integer orders (for example, D 1 C,0+ K (t) = K (1) (t)). The path described briefly opens up broad prospects for accounting for sequential and simultaneous actions of two and three effects, namely, fading memory, depreciation, and distributed lag.
To simplify the description, this article will consider a more simple way of describing the effect of power memory based on exact analytical expressions for solutions. The substitution of expression (25) into Equation (21) with δ = 0 gives: Let us consider the power-law memory. In this case, the memory function can be described by the following equation: where s is the savings parameter, and Γ(α) is the gamma function. For function (31), Equation (30) can be written in the following form: where I µ RL;0+ is the Riemann-Liouville fractional integral [29] of the order µ > 0. Note that Equation (32) with the Riemann-Liouville fractional integral can be considered [33,47] as an approximation of Equation (30) with generalized memory functions s(t, τ).
Using the property that the Caputo fractional derivative is the left inverse operator for the Riemann-Liouville fractional integral (see Lemma 2.21 of [29] p. 95), the action of the Caputo derivative on Equation (32) where D µ C;0+ is the Caputo fractional derivative [29] of the order µ > 0. As a result, to take into account the power-law memory, we can use the following equation: where D α C;0+ is the Caputo fractional derivative of the order α = µ + 1 > 1 with µ > 0, i.e., α > 1. In the general case, we can consider Equation (34) for the order α > 0 including α ∈ (0, 1).
Note that Equation (34) coincides with Equation (29), which is based on the approach described at the beginning of this subsection.
Let us consider the specific example of a production function. We will use the Cobb-Douglas function. In this case, Equation (1) takes the following form: where we can consider b = 1 − a. This production function is easy to analyze, and it appears to be a good first approximation to actual production functions. Therefore, the Cobb-Douglas function is useful for modeling. The substitution of expression (35) into Equation (34) gives: Then, substituting expressions (12) and (14) into Equation (36), we obtain the following equation: where Equation (37) is the non-linear fractional differential equation that describes the generalization of the Solow model without capital depreciation, where we take into account the power-law fading memory.

Remark 1. Using Equation
with t > 0 should be satisfied for parameters of Equation (40) with the Caputo derivative instead of α − γ(α) ≥ 0 that is used in [29].
Remark 2. Note that in equations 3.5.48 and 3.5.49 of [29] p. 209, there are typos: The signs of the parameters γ(α) and α in the gamma functions must be opposite. In Equation (40), these signs are correct.
Taking into account these Remarks and Propositions 3.8 and 3.9 of [29] pp. 209-210, we can formulate the following conditions for the existence of solution (40) for Equation (37).
Conditions (46) for the parameters, under which solution (40) exists, are important for applications in economic models.

Rate of Growth with Power-Law Memory
For solution (40), the rate of growth with power-law memory at time t > 0 is approximately equal to: where a 1. For case α = 1, which corresponds to the absence of memory, Equation (47) can be written in the following form: Let us consider the case when we have growth with memory, i.e., the increasing function K(t) for non-integer values of α. In this case, the growth rate with memory, R(α) should be positive for non-integer values of α ∈ (n − 1, α), n ∈ N, i.e., we have the following inequality: Solutions (40) and (41) exist only if the following conditions are satisfied: Obviously, using Equations (47) and (50), we get that condition (49) has the following form: which will be satisfied for any solutions for t > 0. Mathematically, this statement follows from the fact that if α − γ(α) > n − 1 with n ∈ N, then α − γ(α) > 0.
As a result, we can formulate the following principle.

Principle 1. Principle of Inevitability of Growth with Memory:
For economic models with memory, which are described by Equations (5), (6), and (36) with a 1, q 1, and p 1, the rate of growth with power-law memory (47) is positive.
for non-integer values of α. As a result, the inclusion of power-law memory effects leads to the inevitability of capital growth K(t).
In addition, we have the lower boundary of the growth rates for processes with memory, which is described by the following statement.

Statement 2.
The rate of growth with power-law memory for processes, which are described by Equation (37), has the lower boundary that is defined by the inequality: where n = [α] + 1 for non-integer values of α.
This statement is based on the fact that the solution exists only if the following condition holds: Using definition (47) of the rate of growth with memory, we get (53).

Remark 3.
Solutions, in which the decline (recession) is realized, do not exist within the framework of the suggested model with memory. In this model, the processes with memory cannot have negative values of the rate of growth with a power-law memory. Note that for processes without memory, which are described by Equation (47) with α = 1, the negative growth rates (R(1) < 0) can be realized.
Let us consider the condition under which the growth rate with memory, R(α), is greater than the growth rate without memory, R(1). This condition is represented by the inequality: Using expressions (47) and (48), inequality (55) gives: where we assume that a ∈ (0, 1). Taking into account conditions (46) for the existence of solution (40) and that n = 1 for α ∈ (0, 1), we can formulate the following statement. Statement 3. The rate R(α) of growth with memory for the capital K(t) is greater than the rate R(1) of growth without memory R(α) > R(1), if the following conditions are satisfied: Let us consider the question of how much growth with memory can be greater than growth without memory. To describe this, we give the ratio of growth rates R(α) and R(1). This ratio has the following form: where we assume that R(1) > 0, conditions (57) are satisfied, and β is given by (39). Note that α ∈ (0, 1), in inequality R(α) > R(1) is also satisfied since a > 1 and α + β < 1 + β < 0.

Dynamics of Capital Per Unit of Effective Labor
The amount of capital per unit of effective labor [3] p. 10, is defined by the following equation: The solution of capital in form (43) and expressions (12) and (14) allow us to describe the amount of capital per unit of effective labor (59) by the following equation: Remark 4. Note that the memory model described by the equation with a fractional derivative cannot be reduced to changing the parameterization of the standard model. This is due to the qualitatively different properties of the derivative of non-integer orders. For example, such standard properties of the first-order derivative, as the product rule and the chain rule, do not hold for fractional derivatives. For example, the expression k(t), which is given by solutions of equations for k(t), cannot coincide (for processes with memory) with the expression of k(t) = K(t)/A(t)L(t), which is given by solutions of equations for K(t), and L(t), A(t). For standard models, both expressions are the same. We can state that the variable per capita (per unit of effective labor) constructed from the solutions of original variables is not a solution of the equation for the variable per capita (per unit of effective labor). Let us give a brief explanation of this statement (for a detailed description, see article [46]. Let us assume that we have the solution L sol (t) and A sol (t) of differential equations for L(t) and A(t), respectively. We can consider two variables K(t) and k(t) and two equations K (1) (t) = F(t, K(t)) and k (1) (t) = f (t, k(t). In the standard model, solutions K sol (t) and k sol (t) of these equations can be derived by the formula K sol (t) = k sol (t)A sol (t)L sol (t). The fractional differential equations for the variables K(t) and k(t) cannot be related by such formulas, if the memory effect is taken into account.
For the capital per unit of effective labor k(t), the rate of growth with power-law memory at t > t 0 is approximately equal to: where q 1, p 1, a 1, and β is defined by Equation (38). Expression (61) can be represented in the following form: For α = 1, which describes the case of absence memory, expression (62) takes the following form: Expression (63) describes the rate of growth without memory for a 1. Let us consider the condition under which the growth rate with memory, r(α), is greater than the growth rate without memory, r(1). This condition is represented by the following inequality: Using Equations (62) and (63), inequality (64) takes the following form: Inequality (65) coincides with inequality (56). Therefore, Statement 3 and the system of inequalities (57) also should be satisfied for dynamics of the capital per unit of effective labor.
For the capital per unit of effective labor k(t), the rate r(α) of growth with memory is greater than the rate r(1) of growth without memory r(α) > r(1) if conditions (57) are satisfied.
Let us consider the question of how much growth with memory can be greater than growth without memory. Using Equations (62) and (63), the ratio of r(α) and r(1) is expressed by the following equation: where we assume that r(1) > 0, and the parameter β is defined by Equation (38). For b = 1 − a with a ∈ (0, 1), we have the linearly homogeneous production functions Cobb-Douglas. In this case, the rate of growth with power-law memory (62) takes the following form: Using that α ∈ (n − 1, n), we get: for non-integer values of α. In this case, the ratio of r(α) and r(1) is expressed by the following equation: Using that α ∈ (n − 1, n), we obtain the following inequality: for non-integer values of α.

Example 1.
Using expression (69), we see that the rate of growth with memory can be several times (in α > 1 time) higher than the standard rate of growth without memory. For example, if α = 3.2, then the rate of growth with memory is 220% more (is 3.2 times more) than the standard rate of growth without memory, when other initial parameters of the models (with memory and without memory) are the same.
As a result, we find that the growth rate with memory is greater than the growth rate without memory by α times. This allows us to formulate the following principle, where we assume that conditions (57) are satisfied.

Principle 2. Principle of Changing Growth Rates by Memory:
In the case b = 1 − a with a ∈ (0, 1) and conditions (57) are satisfied, the rate of growth with a power-law memory for capital per unit effective labor can be greater than the rate of growth without memory (r(α) > r(1)) in α times, if the parameter of memory fading α is more than one (α > 1). The growth rate with memory is less than the growth rate without memory, i.e., r(α) < r(1), if the parameter α is less than one (0 < α < 1).
As a result, we can conclude that the fading memory can significantly (several times) change the growth rate, when other parameters of the model are unchanged. Therefore, we should not neglect the memory in economic models.
As a prospect for further studies of non-linear economic growth models with memory, it is possible to emphasize the need to take into account both simultaneous and sequential actions of the effects of memory and depreciation.

Solow-Lucas Model of Closed Economy with Memory
In this section, we consider a generalization of the standard growth model for closed economy without capital depreciation, which is considered by Robert E. Lucas in paper "Making a Miracle" [41]. In the generalized model, we take into account the memory with power-law fading. The suggested model with memory is described by the non-linear fractional differential equation. The explicit expression of analytical solution of this non-linear equation is obtained.

Solow-Lucas Model for Closed Economy without Memory
Let us consider the Solow model for closed economy that is described by Robert E. Lucas in [41] pp. 253-254 (see also Section 3.2 in [42] pp. 73-75). We will call this economic model the Solow-Lucas model. Note that Robert E. Lucas Jr received the Nobel Memorial Prize in Economic Sciences in 1995 "for having developed and applied the hypothesis of rational expectations, and thereby having transformed macroeconomic analysis and deepened our understanding of economic policy" [43,44].
Let us consider an economy that uses physical capital k(t) and the human capital h(t) to produce a single good y(t), the dynamics of which is described by the following equation: Here, the human capital input is multiplied by u, where u is the fraction of time that people spend on producing goods.
The growth of human capital h(t) depends on the amount of time spent on production, adjusted for quality (see Equation (13) in [45]) in the following form: The solution of Equation (72) can be written as: where h(0) is the human capital at time t = 0.
The growth of physical capital depends on the savings rate s, such that we have the following equation: which is interpreted as an economic accelerator without memory and lag. The substitution of expressions (71) and (73) into Equation (74) gives: In this model, it is assumed that variables s and u are considered as given constants. The model, which is described by Equations (71)-(75), will turn out to be a Solow model rewritten for a closed economy [41,42].
In this model, the rate of technological changes (the average Solow balance) is described by the following expression: and the initial technology level is equal to Ah 1−a (0). The long run growth rate of both capital and production per worker is δ (1 − u). The growth rate of human capital, and the ratio of physical and human capital converge to constant values. In the long run, the level of income is proportional to the economy's initial stock of human capital.

Remark 5.
Note that in the Solow-Lucas model, the average values of k(t), h(t), and y(t) are taken as original variables to build a model, in contrast to the Solow model that use K(t), L(t), Y(t) as original variables. In standard models, the use of these two types of variables is equivalent. Unlike standard models, in models with memory the use of different types of variables leads to different models that are not equivalent in the general case [46]. This non-equivalence is associated with the violation of the standard form of the product rule for the fractional derivative of non-integer orders. As a consequence, we have a violation of the standard relationship of the rates of change of absolute (volume) and relative (specific) indicators.

Solow-Lucas Model for Closed Economy with Memory
In the proposed generalization of the standard Solow-Lucas model, we will use the generalizations of Equations (71), (72), and (74).
(1) We assume that the economy, which uses physical capital k(t) and the human capital h(t) to produce a single good y(t), is described by the following equation: where u is the fraction of time people spend producing goods. If b = 1 − a, then Equation (77) gives Equation (71).
(2) The growth of human capital is assumed to be the power law in the following form: If θ = 1, then Equation (78) gives Equation (72) of the standard model.
(3) Equation (74) does not take into account the memory effects, since the derivative dk(t)/dt is determined by the behavior of y(τ) only at the same time of instant τ = t, and does not take into account the history of changes of y(τ) in the past, when τ ∈ (0, t). To take into account the changes of y(τ) in the past, we can use the following equation: where function s(t, τ) allows us to take into account that the growth of physical capital k(t) depends on the changes of y(τ) in the past at τ ∈ (0, t). Function s(t, τ) can be interpreted as a memory function. In order to have standard dimensions of economic quantities, the time variable t is considered as a dimensionless parameter, i.e., we changed the variable t old → t new = t old /t c , where t c is a characteristic time of economic process. An important property of memory is described by the principle of memory fading, which states that the increasing of the time interval leads to a decrease in the corresponding contribution [7,33]. Let us consider the power-law form of the memory fading. In this case, the memory function can be described by the following equation: where s is the savings rate and Γ(α) is the gamma function. For function (80), Equation (79) can be written in the following form: where I µ RL;0+ is the Riemann-Liouville fractional integral of the order µ > 0, [29] pp. 69-70. We should note that (81) with the Riemann-Liouville fractional integral can be considered as an approximation of the equations with generalized memory functions (79). In paper [47], using the generalized Taylor series in the Trujillo-Rivero-Bonilla form for the memory function, we proved that Equation (79) for a wide class of memory functions can be represented through the Riemann-Liouville fractional integrals (and the Caputo fractional derivatives) of non-integer orders.
The action of the Caputo derivative on Equation (81) gives: where D µ C;0+ is the Caputo fractional derivative of the order µ > 0, [29] pp. 90-99. In Equation (82), we can use the property that the Caputo fractional derivative is the left inverse operator for the Riemann-Liouville fractional integral (see the equation of Lemma 2.21 in [29] p. 95), in the following form: for y(t) ∈ L ∞ (0, t) or y(t) ∈ C[0, t], and the property of the Caputo fractional derivative: Using expressions (83) and (84), Equation (82) can be written as: with α = µ + 1. If α = 1, then Equation (85) takes the form of (74). The substitution of (77) into (85) gives the equation which describes the dynamics with memory for physical capital k(t) in the following form: where the behavior of h(t) is given by Equation (78). Equation (86) describes the proposed generalization of the Solow-Lucas model for a closed economy with power-law memory. For α = 1, Equation (86) describes the Solow-Lucas model without memory, where b 1 − a and θ 1.

Growth Rates of Closed Economy with Memory
Let us obtain the solution of the non-linear fractional differential Equation (78) which describes the Solow-Lucas model of closed economy with power-law memory.
We first solve Equation (78) with θ 1 for the human capital. Equation (78) can be written as: Therefore, the solution of Equation (87) with θ 1 has the following form: where we can use c = h 1−θ (0). For simplicity, we will assume that h(0) = 0 and h(t 0 ) = 1. Then, we get c = 0 and solution (88) can be written in the following form: Equation (89) with t = t 0 > 0 has the following form: where δ (1 − u)(1 − θ) > 0. Assuming that h(t 0 ) = 1, we get: and solution (89) can be represented by the following equation: The substitution of (92) into (86) gives: where α = µ + 1 with µ > 0. Equation (93) can be written in the following form: where θ 1 and the constant λ is defined as: Using equations 3.5.47 and 3.5.48 of Section 3.5.3 in [29] p. 209, the following equation: with λ 0, a 1, α = µ + 1 > 1, and has the following solution: where λ is defined by Equation (95), and The condition for the existence of a solution can be written as: As a result, the rate of growth with power-law memory (α 1) at time t > 0 is approximately equal to: for t > t 0 > 0 and θ 1, a 1, where we use (97). For α = 1, expression (101) takes the following form: Expression (102) describes the rate of growth without memory for the case θ 1, a 1.
Using that solution (98) exists if conditions (100) are satisfied, i.e., α − γ(α) > n − 1, we get: for non-integer values of α ∈ (n − 1, n) and t > 0. As a result, we get that R(α) > 0 for non-integer values of α. Therefore, the rate of growth with power-law memory is positive for t > 0 and non-integer values of α.
We should note that the generalized model, which is described by Equation (79) assumes that α = µ + 1, where µ > 0. Therefore, we may not consider the case with 0 < α < 1.
Let us consider the condition, under which the rate of growth with memory, R(α), is greater than the rate of growth without memory, R(1), that is: Using expressions (101) and (102), inequality (104) gives: where we assume that a ∈ (0, 1). Using Equation (105), we can state that for a ∈ (0, 1), we have R(α) > R(1) for all t > 0. As a result, we can formulate the following principle.

Principle 3. Principle of Changing Growth Rates by Memory:
If a solution of Equation (96) for the physical capital k(t) exists and a < 1, then the rate of growth with power-law memory, R(α), is greater than the rate of growth without memory, R(1), i.e., R(α) > R(1).

Statement 4.
The rate R(α) of growth with memory for the physical capital k(t) is greater than the rate R(1) of growth without memory (R(α) > R(1)), if the following conditions are satisfied: The ratio of R(α) and R(1), which are given by expressions (101) and (102), is expressed by the following equation: For the special case 0 < 1 − θ = b < 1, Equation (107) gives: where we use α = µ + 1.

Example 2.
Using expression (108), we see that the rate of growth with memory can be several times (in µ > 0 time) higher than the standard rate of growth without memory. For example, if α = 3.2 and 1 − θ = b ∈ (0, 1), then the rate of growth with memory is 110% more (is 2.1 times more) than the standard rate of growth without memory, when other initial parameters of the models (with memory and without memory) are the same.
In this case, we see that the effects of memory can increase the growth rate by more than two times in comparison with the standard Solow-Lucas model without memory.
As a result, we can conclude that the memory effects can significantly (several times) change the growth rate, when other parameters of the model are unchanged, and we should not neglect the memory in economic models.

Conclusions
This article proposes generalizations of two well-known standard models with continuous time, which are described by Nobel laureates in economics, Robert M. Solow and Robert E. Lucas. In the suggested generalizations of the standard models, the memory with power-law fading is taken into account. As a mathematical tool for accounting for fading memory, the article uses integral derivatives of non-integer order. Note that the violation of the standard product (Leibniz) rule for fractional derivatives of non-integer orders leads to the non-equivalence of models that are equivalent in the absence of memory [46]. For the suggested model with memory, we derive non-linear fractional differential equations and its explicit analytical solutions. We prove that the rate of growth with power-law memory can be greater than the rate of growth without memory. For non-linear growth models with power-law memory, the growth rates can be changed by memory effects.
In this article, we proposed the following three principles, which can be used to describe the changes arising in the economic dynamics, when taking into account the power-law memory: (I) The Principle of Inevitability of Growth with Memory states that in the first proposed model, which are described by Equations (5), (6) and (36), with a 1, q 1, and p 1, the rate of growth with power-law memory (47) is positive (52) for non-integer values of α. As a result, the inclusion of power-law memory effects leads to the inevitability of capital growth K(t). (II) The Principle of Changing Growth Rates by Memory states that in the first proposed model, the rate of growth with power-law memory for capital per unit of effective labor can be greater than the rate of growth without memory in α times, if the parameter of memory fading α is more than one (α > 1). (III) The Principle of Changing Growth Rates by Memory states that in the second proposed model, the rate of growth with power-law memory is greater than the rate of growth without memory, if a solution of Equation (96) for the physical capital k(t) exists and a < 1.
Therefore, we should emphasize that the memory effects can significantly (several times) change the growth rate, when other parameters of the model are unchanged, and we should not neglect the memory in economic models. For example, if α = 3.2, then the rate of growth with memory is 220% more (3.2 times more) and 110% more (2.1 times more) for first and second proposed models, respectively than the standard rate of growth without memory.
Let us note here a possible direction for further mathematical research on non-linear models of economic growth with memory. An interesting and important part of future research is to account for both simultaneous and sequential actions of memory effects and depreciation. Standard economic models often use the exponential kernel of the integral operator to describe depreciation. In Equation (17), kernel D(t − τ) = e −δ (t−τ) describes the standard (exponentially distributed) depreciation (up to a numerical factor) in form (18). The following types of integral operator kernels can be used to account for memory and depreciation effects in economic growth models: (a) In Equation (17), kernel D(t − τ) = (t − τ) α−1 /Γ(α) describes the interconnection of variables with power-law memory. In this case, integral operator (17) is called the Riemann-Liouville fractional integral. Economic models of this type were considered, for example, in works [7,[33][34][35][36]38]. (b) We can use kernel D(t − τ) = e −δ (t−τ) (t − τ) α−1 /Γ(α), which describes the gamma distributed lag by Equation (17). The economic model with this type of distributed lag was considered in work [53]. Note that kernels of exponential amortization and power-law memory are special cases of this kernel for the case δ = 0 and α = 1, respectively. (c) To describe the sequential action of two such effects as fading memory and continuously distributed delay (or continuously distributed amortization), it is necessary to use integral fractional operators with distributed lag, proposed in works [54] and then applying them in linear economic models [55][56][57][58]. (d) Kernel D(t − τ) = (t − τ) α−1 E γ ρ,α [ω(t − τ) ρ ] describes the connection of variables with power-law memory and depreciation in special cases. Here, E γ ρ,α [ω(t − τ) ρ ] is the three-parameter Mittag-Leffler function [59,60], which is also called the Prabhakar function [61][62][63]. Integral operator (17) with this kernel is called the Prabhakar fractional integral [61][62][63]. The economic model with this type was considered in works [64,65], where a new economic model through learning-by-doing with memory was proposed. For the first time, the fractional derivative of Riemann-Liouville type, which has the Prabhakar function in the kernel, was proposed by Anatoly A. Kilbas, Megumi Saigo, and Ram K. Saxena in work [61] in 2004. This operator can be known as the Kilbas-Saigo-Saxena (KSS) fractional derivative. Note that the KSS operators are left-inverse for the Prabhakar fractional integrals [61]. Then, the fractional derivatives and integrals began to be studied in various papers (for example, see [62][63][64] and references therein). An application of these operators in economics with memory is proposed in papers [64,65]. We should emphasize that the main property of any generalized (fractional) derivative is to be a left-inverse operator to the corresponding generalized (fractional) integral operator. This requirement is important for a self-consistent mathematical theory of the fractional operators to have a general fractional calculus of these operators.
Funding: This research received no external funding.

Conflicts of Interest:
The author declares no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.