Existence of Limit Cycles in the Solow Model with Delayed-Logistic Population Growth

Abstract

This paper is devoted to the existence and stability analysis of limit cycles in a delayed mathematical model for the economy growth. Specifically the Solow model is further improved by inserting the time delay into the logistic population growth rate. Moreover, by choosing the time delay as a bifurcation parameter, we prove that the system loses its stability and a Hopf bifurcation occurs when time delay passes through critical values. Finally, numerical simulations are carried out for supporting the analytical results.

(ProQuest: ... denotes non-US-ASCII text omitted.)

Carlo Bianca 1 and Luca Guerrini 2

Academic Editor:A. Agouzal and Academic Editor:J. Banas

1, Dipartimento di Scienze Matematiche, Politecnico, Corso Duca degli Abruzzi 24, 10129 Torino, Italy
2, Department of Management, Polytechnic University of Marche, 60121 Ancona, Italy

Received 26 August 2013; Accepted 1 December 2013; 28 January 2014

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

1. Introduction

Investigations on the existence of limit cycles and the related stability analysis in nonlinear autonomous differential equations go back to Poincaré [1], who has proved that if differential equations admit a limit cycle then the difference between the number of nodes, centers, and foci enclosed by a limit cycle and the number of enclosed saddle points are equal to one. In 1924 Ivar Otto Bendixson established a sufficient condition for nonexistence of limit cycles [2]. Thereafter, the existence of limit cycles in systems of nonlinear differential equations has been of great interest in the pure and applied mathematics; see the survey [3].

This paper is concerned with the existence and stability analysis of a differential equations system, with time delay, that can be proposed for the modeling of the economic growth.

The origin of the economic growth theory goes back to Solow [4] and Swan's [5] neoclassical growth theory. A standard assumption in the economic growth theory is that population always grows at a constant and positive rate. Recently, the idea that such a specification is unrealistic has been arisen. Indeed this assumption implies, as in the Malthus model [6], that the population size goes to infinity as time goes to infinity.

Assuming a nonconstant population growth rate but variable and bounded over time, Guerrini proposed in [7] a generalization of the Solow model. Specifically a standard neoclassical technology and a Verhulst [8] logistic population growth law are assumed. The mathematical analysis performed in [7] shows the existence of a unique globally stable steady state to which the economy adjusts, so sharing the technical and qualitative properties with the Solow model.

It is well known that, in the logistic law models, the time evolution of the population density at time t depends on the population density at the same time. Obviously, this is a roughly approximation of the reality, since the real process of reproduction is not instantaneous in time and is strictly dependent on the previous instants of time. Therefore, the logistic law appears to be inadequate for the time description of the population density growth; the interested reader is referred to Hutchinson [9] for a complete description in the case of the ecology.

Bearing all the above in mind, the logistic law has been modified and a time delay has been inserted. The resulting mathematical model is known as the delayed logistic equation or Hutchinson's equation. It is worth precising that the introduction of time delay has been carried out also in biological systems, see paper [10], and in economics systems (see [11-15]).

The present paper aims at investigating how the delayed logistic law may affect the dynamics of the Solow model and consequently provokes the existence and the stability of limit cycles and Hopf bifurcations. The corresponding mathematical model is a system of two first-order nonlinear delayed differential equations, whose linear stability is discussed by analyzing the associated characteristic transcendental equation. The analytical results show that, as the magnitude of the time delay increases, the system loses its stability and a Hopf bifurcation occurs when the time delay passes through a sequence of critical values. Numerical simulations refer to a sensitive analysis on the time delay and support the analytical results.

It is worth stressing that the problem whether time delay may induce cyclic behavior in the Solow economic growth model has been already discussed in the pertinent literature; see, among others, paper [16] and the references cited therein.

The contents of this paper are outlined as follows. After this introduction, Section 2 is meant to be the description of the economy where the delayed Solow model is derived. Analytical investigations on the existence of steady states, stability and Hopf bifurcation analysis are dealt with in Section 3. Section 4 is concerned with further analysis on the Hopf bifurcation, namely, the direction of the bifurcation (when it is subcritical or supercritical) and the bifurcating period. Numerical simulations obtained by sensitive analysis on the magnitude of the time delay are performed in Section 5. A critical analysis of the results and discussions of future research perspectives are postponed to Section 6.

2. The Delayed Mathematical Model

This section deals with the derivation of the Solow mathematical model with a time delay. Specifically a closed economy is considered, which consists of a single good Y = Y ( t ) that is used either for consumption or investment. The good is produced by labor L = L ( t ) and physical capital K = K ( t ) accordingly to a neoclassical production function Y = F ( K , L ) ; see [17]. Each worker has a unit of time available each period that is supplied inelastically in the labor market. There is a full employment in the economy, so that population size and workforce can be used interchangeably. Since the economy is assumed to be closed, change in the capital stock equals gross investment less depreciation, namely, [figure omitted; refer to PDF] where s denotes the constant fraction of the gross income Y saved and δ is the constant fraction of capital stock that disappears as a result of depreciation. Population growth rate evolves according to the following delayed logistic equation [9]: [figure omitted; refer to PDF] where _{L d} : = L ( t - τ ) , τ > 0 is the time delay, and a > b > 0 are positive constants such that a - b L ( 0 ) > 0 , with L ( 0 ) being the density of the population at the initial time. Setting k = K / L and writing the production function in intensive form f ( k ) = F ( K / L , 1 ) , we obtain the following mathematical model described by two nonlinear delayed differential equations: [figure omitted; refer to PDF] for some initial function L ( t ) = [straight phi] ( t ) , for t ∈ [ - τ , 0 ] . In contrast to classical dynamical systems with zero delay, the definition of the initial function [straight phi] , usually called history function, is required. The domain of [straight phi] is defined over the range of time delimited by the time delay.

It is worth stressing that with respect to the paper [18], the mathematical model considered in this paper (to which we refer as the Solow model with delayed-logistic population rate) includes the time delay in the logistic term. Moreover, the mathematical analysis will be performed with respect to an arbitrary neoclassical production function f .

3. Asymptotic Analysis: Steady States and Hopf Bifurcation

This section is devoted to the study of the asymptotic behavior of the mathematical model (3). The main goal of the present paper is to investigate the existence of equilibrium points (steady states), their stability analysis, and the possible existence of Hopf bifurcations.

The steady states of the mathematical model (3) coincide with the steady states of the corresponding model with zero delay. Hence, there exists a unique nontrivial steady state ( _{k *} , _{L *} ) , such that s f ( _{k *} ) = δ _{k *} and _{L *} = a / b . Let x = k - _{k *} and y = L - _{L *} , so that the equilibrium in (3) is shifted to the origin. By applying in system (3) the Taylor expansion with center the origin, the model (3) recasts into the following system: [figure omitted; refer to PDF] where the nonlinear parts _{g 1} ( x , _{y d} ) , _{g 2} ( y , _{y d} ) read: [figure omitted; refer to PDF] where _{y d} = y ( t - τ ) . Writing the linear part of system (4) in matrix form [figure omitted; refer to PDF] it is immediately seen that the resulting characteristic equation is [figure omitted; refer to PDF] Equation (7) is a quasi-polynomial equation, which has, in general, an infinite number of (complex) roots. When there is no time delay, that is, τ = 0 , the characteristic equation reduces to a simple quadratic equation with two negative real characteristic roots, λ = - a and λ = s ^{f [variant prime]} ( _{k *} ) - δ = - s [ f ( _{k *} ) - _{k *} ^{f [variant prime]} ( _{k *} ) ] / _{k *} < 0 . Thus, the steady state of system (4) is locally asymptotically stable.

Let τ > 0 be fixed. Taking the time delay as a bifurcation parameter, we investigate the effects of the time delay on the dynamics of the model (3). Accordingly, we look for values of τ so that the steady state changes its asymptotic stability from local stability to instability.

It is well known that the steady state of system (4) is locally asymptotically stable if each of the characteristic roots of (7) has negative real part. Since λ = 0 is not a zero of (7), we need to examine only when this equation has pure imaginary roots λ = ± i ω , where ω is a positive real number.

Lemma 1.

The characteristic equation (7) has a pair of simple conjugate pure imaginary roots λ = ± i _{ω 0} ( _{ω 0} > 0 ) at τ = _{τ j} , where [figure omitted; refer to PDF] Furthermore, the following transversality condition: [figure omitted; refer to PDF] holds true.

Proof.

Let λ = i ω ( ω > 0 ) be a root of the characteristic equation (7). Then λ must be a root of [figure omitted; refer to PDF] Substituting λ = i ω into (10), we have [figure omitted; refer to PDF] Separating the real and imaginary parts yields [figure omitted; refer to PDF] Since ω and a are positive constants, ω τ = π / 2 + 2 j π , j = 0,1 , 2 , ... , and ω = a = _{ω 0} . To ensure that the roots λ = ± i _{ω 0} are simple roots, we differentiate (10) with respect to λ and we obtain d D ( λ , τ ) / d λ = - 1 + a τ ^{e - λ τ} . It is easy to see that the conditions [figure omitted; refer to PDF] hold true only if - 1 ± i _{ω 0} τ = 0 . This concludes the first part of the statement.

Next, let λ ( τ ) = μ ( τ ) + i ω ( τ ) be the root of (7) such that μ ( _{τ j} ) = 0 and ω ( _{τ j} ) = _{ω 0} . Differentiating (10) with respect to τ and using (10), one has [figure omitted; refer to PDF] On the other hand, [figure omitted; refer to PDF] Then the proof of the theorem has been gained.

Since the sign of d [ Re λ ( _{τ j} ) ] / d τ is positive, then each crossing of the real part of the characteristic roots at _{τ j} must be from left to right. With being the crossing direction, always toward instability, the following result on the distribution of roots of (7) is obtained.

Proposition 2.

If τ ∈ [ 0 , _{τ 0} ) , all roots of (7) have negative real parts. If τ = _{τ 0} , all roots of (7), except for ± i _{ω 0} , have negative real parts. If τ ∈ ( _{τ j} , _{τ j + 1} ) , for j = 0,1 , 2 , ... , (7) has 2 ( j + 1 ) roots with positive real parts.

Bearing all the above in mind, we have that the hypotheses for Hopf bifurcation are satisfied at _{τ j} . These allow us to prove the following theorem.

Theorem 3.

The positive steady state ( _{k *} , _{L *} ) of system (3) is locally asymptotically stable for τ ∈ [ 0 , _{τ 0} ) and unstable for τ > _{τ 0} . Moreover, system (3) undergoes a Hopf bifurcation at ( _{k *} , _{L *} ) , when τ = _{τ j} , for j = 0,1 , 2 , ... .

Proof.

Spectral properties stated in Proposition 2 lead immediately to stability properties of the zero equilibrium of system (4) and of the positive equilibrium ( _{k *} , _{L *} ) of system (3). Hence, we have the proof.

4. Analysis of the Hopf Bifurcations

This section is meant to the study of the direction of the Hopf bifurcation and the stability of bifurcating periodic solutions. The results are obtained by applying the normal form theory and the center manifold theorem stated in [19].

As already shown, the system (3) undergoes a Hopf bifurcation at the positive equilibrium point ( _{k *} , _{L *} ) when τ = _{τ j} and λ = i _{ω 0} is the corresponding purely imaginary root of the characteristic equation at this point. Setting τ = _{τ j} + μ , μ ∈ ... , then μ = 0 is the Hopf bifurcation value of system (3) for the new bifurcation parameter μ .

Let C = C ( [ - _{τ j} , 0 ] , ^{... 2} ) be the space of continuous real-valued functions. For any [straight phi] = ( _{[straight phi] 1} , _{[straight phi] 2} ) ∈ C , let [figure omitted; refer to PDF] where [figure omitted; refer to PDF] By Riesz's representation theorem, there exists a matrix whose components are bounded variation functions η ( θ , μ ) in θ ∈ [ - _{τ j} , 0 ] such that [figure omitted; refer to PDF] In fact, we can choose [figure omitted; refer to PDF] where Γ ( θ ) is the Dirac delta function. For [straight phi] ∈ C , we define [figure omitted; refer to PDF] Then system (4) can be rewritten in the following form: [figure omitted; refer to PDF] where _{u t} = u ( t + θ ) , for θ ∈ [ - _{τ j} , 0 ] . For ψ ∈ C ~ = C ( [ 0 , _{τ j} ] , ^{... 2} ) , we define [figure omitted; refer to PDF] For [straight phi] ∈ C and ψ ∈ C ~ , we define the following bilinear form: [figure omitted; refer to PDF] Therefore, A ( 0 ) and ^{A *} ( 0 ) are adjoint operators. Let q ( θ ) and ^{q *} ( r ) be the eigenvectors of A ( 0 ) and ^{A *} ( 0 ) corresponding to the eigenvalues i _{ω 0} and - i _{ω 0} , respectively, and normalized so that Y9; ^{q *} ( r ) , q ( θ ) YA; = 1 . It is well known that q ( θ ) can be computed from A ( 0 ) q ( θ ) = i _{ω 0} q ( θ ) , where q ( θ ) = q ( 0 ) ^{e i _{ω 0} θ} , similarly for ^{q *} ( r ) . For a solution _{u t} of (17) at μ = 0 , we define [figure omitted; refer to PDF] where [figure omitted; refer to PDF] Then (17) valued in the points of the center manifold ...9E; is described by [figure omitted; refer to PDF] where _{g 0} = g ( 0 , _{u t} ) , with g defined by (17). Now (26) can be written as [figure omitted; refer to PDF] with [figure omitted; refer to PDF] Next, substitute _{u t} ( θ ) = W ( t , θ ) + 2 Re [ z ( t ) q ( θ ) ] into _{g 0} and denote [figure omitted; refer to PDF] Hence, from g ( z , z - ) = ^{q - *} ( 0 ) _{g 0} , we obtain [figure omitted; refer to PDF] The next step is to compute _{W 20} ( θ ) and _{W 11} ( θ ) since the term _{g 21} is dependent on them. By (22) and (28) we have [figure omitted; refer to PDF] where [figure omitted; refer to PDF] Substituting (26) and (28) into W . = _{W z} z . + _{W z -} z - . and comparing the coefficients of the resulting equation with those of (32), we get [figure omitted; refer to PDF] It follows that [figure omitted; refer to PDF] Comparing the coefficients of (33) with (32), for θ ∈ [ - _{τ j} , 0 ) , we obtain [figure omitted; refer to PDF] and for θ = 0 we have [figure omitted; refer to PDF] Solving (35) for _{W 20} ( θ ) and _{W 11} ( θ ) , one has [figure omitted; refer to PDF] where _{E 1} , _{E 2} ∈ ^{... 2} can be determined by setting θ = 0 in H ( z , z - , θ ) . Substituting (38) into (34), using (37), and noticing that [figure omitted; refer to PDF] we arrive at [figure omitted; refer to PDF] As a result, _{E 1} and _{E 2} are calculated. Based on the foregoing analysis, we can see that each _{g i j} is determined. Therefore, we can compute the following quantities: [figure omitted; refer to PDF] which give the properties of bifurcating periodic solutions.

Theorem 4.

Let ( _{k *} , _{L *} ) be the steady state of the model (3). Then one has the following.

(1) The direction of the Hopf bifurcation of the system (3), at the equilibrium ( _{k *} , _{L *} ) for τ = _{τ j} , is subcritical (resp., supercritical) if _{μ 2} < 0 (resp., _{μ 2} > 0 ).

(2) The bifurcating periodic solution on the center manifold is unstable (resp., locally asymptotically stable) if _{β 2} > 0 (resp., _{β 2} < 0 ).

(3) The period of the bifurcating periodic solution decreases (resp., increases) if _{T 2} < 0 (resp., _{T 2} > 0 ).

5. Numerical Simulations

This section is concerned with some numerical simulations of the mathematical model (3). Specifically we depict the behavior of the solutions when the time delay is varied, thus performing a sensitive analysis on the parameter τ .

The mathematical model (3) is characterized by six nonnegative parameters which have an economic meaning and a function f . The set of simulations of this section is obtained by choosing f as the Cobb-Douglas function, namely, [figure omitted; refer to PDF] fixing the following magnitude of the following five parameters: [figure omitted; refer to PDF] and performing a sensitive analysis on the time delay τ .

The first set of simulations refer to τ = 9 . As Figure 1 shows, oscillations of the solutions occur (Figure 1(a)) for the all length of the simulation. A cycle limit is reached; see Figure 1(b). Further numerical simulations show that, when we increase the time delay, namely, for τ ...5; 9 , we have a cycle limit whose period increases; see Figure 2(b). Moreover, looking at Figure 2(a), we can see that also the magnitude of the oscillations is increased.

The time evolution of the functions k ( t ) (dotted line) and L ( t ) for α = 0.8 , s = 0.5 , a = 0.2 , b = 0.5 , δ = 0.9 , and τ = 9 (a). The phase space of L ( t ) versus k ( t ) for t ∈ [ 0,500 ] (b).

(a) [figure omitted; refer to PDF]

(b) [figure omitted; refer to PDF]

The time evolution of the functions k ( t ) (dotted line) and L ( t ) for α = 0.8 , s = 0.5 , a = 0.2 , b = 0.5 , δ = 0.9 , and τ = 11 (a).

(a) [figure omitted; refer to PDF]

(b) [figure omitted; refer to PDF]

When the magnitude of the time delay decreases, also the oscillations and the cycle limit period decrease; see Figure 3. Finally, when the magnitude of the time delay is less than 8 , then the oscillations decrease in time and the dimension of the cycle limit is reduced; see Figure 4.

The time evolution of the functions k ( t ) (dotted line) and L ( t ) for α = 0.8 , s = 0.5 , a = 0.2 , b = 0.5 , δ = 0.9 , and τ = 8 (a). The phase space of L ( t ) versus k ( t ) for t ∈ [ 0,1000 ] (b).

(a) [figure omitted; refer to PDF]

(b) [figure omitted; refer to PDF]

The time evolution of the functions k ( t ) (b) for α = 0.8 , s = 0.5 , a = 0.2 , b = 0.5 , δ = 0.9 , and τ = 7 (a). The phase space of L ( t ) versus k ( t ) for t ∈ [ 0,500 ] (b).

(a) [figure omitted; refer to PDF]

(b) [figure omitted; refer to PDF]

It is worth stressing that if the time delay is less than 5 then the cycle limit reduces to the point ( 0.05,0.4 ) and oscillations disappear for τ ...4; 2 .

The above set of simulations support also the direction of the Hopf bifurcation. Indeed, a supercritical Hopf bifurcation leads from a decaying oscillation to growth and saturation of a sustained oscillation.

6. Conclusions and Research Perspectives

The mathematical model proposed in this paper is a generalization of the classical Solow model. Specifically the population growth rate is not assumed to be constant and a time delay is inserted in order to better approximate the real-world economy. Analytical results have shown the existence of a steady state and the possibility of a Hopf bifurcation. Moreover, the existence of limit cycles has been proved and numerical simulations, obtained by performing a sensitive analysis on the time delay parameter, have plotted the limit cycles in some cases.

The existence of limit cycles is an important topic in the economic systems. Indeed, limit cycle is a trajectory for which the economy of the system would be constant over a cycle; namely, on an average there is no loss or gain of economy. Limit cycle is an outcome of delicate energy balance due to the presence of nonlinear term into the mathematical model. Therefore, the cycles can be interpreted as economic fluctuations.

The mathematical model (3) can be further generalized by including the time delay also into the k function. Moreover, the time evolution of the population density can be modeled by a kinetic approach; see the review paper [20] and the references cited therein. The kinetic approach proposed in papers [21-24] allows taking into account the interactions that occur among the individuals of a population. Moreover, these models, even if they refer to nonequilibrium systems (namely, systems subjected to external force fields), can attain a nonequilibrium stationary state [25] and a bifurcation analysis can be performed. Finally, by asymptotic methods we can derive macroscopic evolution equations (see, among others, papers [26, 27]) that are of great interest in the pure [28] and applied mathematics.

It is worth mentioning that an important research perspective is the comparison of the generalized model introduced in the present paper with the experimentally measurable quantities. This is a work in progress and results will be presented to due course.

Acknowledgment

The first author acknowledges the support by the FIRB project RBID08PP3J-Metodi matematici e relativi strumenti per la modellizzazione e la simulazione della formazione di tumori, competizione con il sistema immunitario, e conseguenti suggerimenti terapeutici.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

References

[1] H. Poincaré, "Mémoire sur les courbes définies par une équation différentialle. I, II," Journal de Mathématiques Pures et Appliquées , vol. 7, pp. 375-422, 1881.

[2] I. Bendixson, "Sur les courbes définies par des équations différentielles," Acta Mathematica , vol. 24, no. 1, pp. 1-88, 1901.

[3] Yu. S. Ilyashenko, "Finiteness theorems for limit cycles," Uspekhi Matematicheskikh Nauk , vol. 45, no. 2(272), pp. 143-200, 1990.

[4] R. M. Solow, "A contribution to the theory of economic growth," Quarterly Journal of Economics , vol. 70, pp. 65-94, 1956.

[5] T. W. Swan, "Economic growth and capital accumulation," Economic Record , vol. 32, pp. 334-361, 1956.

[6] T. R. Malthus An Essay on the Principle of Population , Oxford University Press, Oxford, UK, 1999.

[7] L. Guerrini, "The Solow-Swan model with a bounded population growth rate," Journal of Mathematical Economics , vol. 42, no. 1, pp. 14-21, 2006.

[8] P. F. Verhulst, "Notice sur la loi que la population suit dans son accroissement," Correspondance Mathématique et Physique , vol. 10, pp. 113-121, 1838.

[9] G. E. Hutchinson, "Circular causal systems in ecology," Annals of the New York Academy of Sciences , vol. 50, no. 4, pp. 221-246, 1948.

[10] C. T. H. Baker, G. A. Bocharov, C. A. H. Paul, "Mathematical modelling of the interleukin-2 T-cell system: a comparative study of approaches based on ordinary and delay differential equations," Journal of Theoretical Medicine , vol. 2, pp. 117-128, 1997.

[11] C. Bianca, M. Ferrara, L. Guerrini, "Hopf bifurcations in a delayed-energy-based model of capital accumulation," Applied Mathematics & Information Sciences , vol. 7, no. 1, pp. 139-143, 2013.

[12] C. Bianca, M. Ferrara, L. Guerrini, "The Cai model with time delay: existence of periodic solutions and asymptotic analysis," Applied Mathematics & Information Sciences , vol. 7, no. 1, pp. 21-27, 2013.

[13] L. V. Ballestra, L. Guerrini, G. Pacelli, "Stability switches and Hopf bifurcation in a Kaleckian model of business cycle," Abstract and Applied Analysis , vol. 2013, 2013.

[14] L. Guerrini, M. Sodini, "Nonlinear dynamics in the Solow model with bounded population growth and time-to-build technology," Abstract and Applied Analysis , vol. 2013, 2013.

[15] C. Bianca, M. Ferrara, L. Guerrini, "Qualitative analysis of a retarded mathematical framework with applications to living systems," Abstract and Applied Analysis , vol. 2013, 2013.

[16] P. J. Zak, "Kaleckian lags in general equilibrium," Review of Political Economy , vol. 11, pp. 321-330, 1999.

[17] L. Gallaway, V. ShuklaThe, "The neoclassical production function," The American Economic Review , vol. 64, pp. 348-358, 1974.

[18] C. Bianca, L. Guerrini, "On the Dalgaard-Strulik model with logistic population growth rate and delayed-carrying capacity," Acta Applicandae Mathematicae , vol. 128, pp. 39-48, 2013.

[19] B. Hassard, D. Kazarino, Y. Wan Theory and Application of Hopf Bifurcation , Cambridge University Press, Cambridge, UK, 1981.

[20] C. Bianca, "Thermostatted kinetic equations as models for complex systems in physics and life sciences," Physics of Life Reviews , vol. 9, pp. 359-399, 2012.

[21] C. Bianca, "Modeling complex systems by functional subsystems representation and thermostatted-KTAP methods," Applied Mathematics & Information Sciences , vol. 6, pp. 495-499, 2012.

[22] C. Bianca, M. Ferrara, G. Guerrini, "High-order moments conservation in thermostatted kinetic models," Journal of Global Optimization , 2013.

[23] C. Bianca, "An existence and uniqueness theorem to the Cauchy problem for thermostatted-KTAP models," International Journal of Mathematical Analysis , vol. 6, no. 17-20, pp. 813-824, 2012.

[24] C. Bianca, "Modeling complex systems with particles refuge by thermostatted kinetic theory methods," Abstract and Applied Analysis , vol. 2013, 2013.

[25] C. Bianca, "Existence of stationary solutions in kinetic models with Gaussian themostats," Mathematical Methods in the Applied Sciences , vol. 36, no. 13, pp. 1768-1775, 2013.

[26] A. Bellouquid, C. Bianca, "Modelling aggregation-fragmentation phenomena from kinetic to macroscopic scales," Mathematical and Computer Modelling , vol. 52, no. 5-6, pp. 802-813, 2010.

[27] P. Degond, B. Wennberg, "Mass and energy balance laws derived from high-field limits of thermostated Boltzmann equations," Communications in Mathematical Sciences , vol. 5, no. 2, pp. 355-382, 2007.

[28] M. A. Ragusa, "Commutators of fractional integral operators on Vanishing-Morrey spaces," Journal of Global Optimization , vol. 40, no. 1-3, pp. 361-368, 2008.