Abstract and Applied Analysis
Volume 2015 (2015), Article ID 760136, 12 pages
http://dx.doi.org/10.1155/2015/760136

Research Article

The Existence and Uniqueness of Global Solutions to the Initial Value Problem for the System of Nonlinear Integropartial Differential Equations in Spatial Economics: The Dynamic Continuous Dixit-Stiglitz-Krugman Model in an Urban-Rural Setting

Minoru Tabata¹ and Nobuoki Eshima²

¹Department of Mathematical Sciences, Osaka Prefecture University, Sakai, Osaka 599-8531, Japan
²Department of Statistics, Oita University, Oita 879-5593, Japan

Received 18 May 2014; Accepted 29 July 2014

Academic Editor: Juan J. Nieto

Copyright © 2015 Minoru Tabata and Nobuoki Eshima. 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.

Abstract.

Assume that economic activities are conducted in a bounded continuous domain where workers move toward regions that offer higher real wages and away from regions that offer below-average real wages. The density of real wages is calculated by solving the nominal wage equation of the continuous Dixit-Stiglitz-Krugman model in an urban-rural setting. The evolution of the density of workers is described by an unknown function of the replicator equation whose growth rate is equal to the difference between the density of real wages and the average real wage. Hence, the evolution of the densities of workers and real wages is described by the system of the nominal wage equation and the replicator equation. This system of equations is an essentially new kind of system of nonlinear integropartial differential equations in the theory of functional equations. The purpose of this paper is to obtain a sufficient condition for the initial value problem for this system to have a unique global solution.

1. Introduction

The new economic geography (NEG) is a new branch of spatial economics that was initiated by Krugman in the early 1990s. This new branch has attracted many social scientists and becomes one of the most important major branches of spatial economics at present. In 2008 Krugman received the Nobel Memorial Prize in Economic Sciences (officially Sveriges Riksbank Prize in Economic Sciences in Memory of Alfred Nobel) for his great contribution to the NEG (see [1–6]). A large number of mathematical models have been built in the NEG. In particular there are many models described by nonlinear integropartial differential equations that are new and important in the theory of functional equations (see [1–3]). Hence the NEG is regarded as a new frontier of the theory of nonlinear integropartial differential equations.

The Krugman core-periphery model (the CP model) is the origin of the NEG (see [1, Chapter 5]) since various models are constructed as its extension. The CP model is a discrete model. In this model economic activities are conducted at two points. These two points represent a core region and a periphery region, respectively. An extension to the case of a finite set of points has been studied in [1]. This model is called the Dixit-Stiglitz-Krugman model (DSK model). Its mathematical foundation is studied in [7–13]. Moreover, in [14], we consider an extension of the CP model to the case of a bounded continuous domain where economic activities are conducted continuously in space. This model is called the continuous DSK model (cDSK model).

These models are static models with no population dynamics. It is very important in spatial economics to build population dynamics into them. Hence Krugman constructs the dynamic DSK model (dDSK model) by combining the DSK model with the replicator dynamics that is one of the most fundamental dynamics in evolutionary game theory (see [15] and [16, Chapter 3]). His dynamic model is very important in spatial economics since it describes economies of agglomeration in the case where workers move from one point to another to seek higher real wages within a finite set of points at which economic activities are conducted (see [1, p. 62, p. 77]). Hence, by following this line, we consider the dynamic cDSK model (dcDSK model) in this paper; that is, we combine the cDSK model with the replicator dynamics. This dynamic model is regarded as a continuous version of the dDSK model and explains agglomeration of capital and concentration of workers when workers move in a bounded continuous domain where economic activities are conducted continuously in space.

Let us discuss the dcDSK model from the viewpoint of the theory of functional equations. The dcDSK model is described by the system of the nominal wage equation and the replicator equation. We refer to this system of equations as the dcDSK system. The nominal wage equation is a nonlinear integral equation that contains the density of nominal wages as an unknown function and the density of workers as a known function (see [14, ]). Hence, if we solve this equation under the condition that the density of workers is a given function, then we can obtain the density of nominal wages. However, the integral kernel of the nominal wage equation contains not only these densities but also the price index. We must note that the price index itself is a nonlinear integral operator acting on the density of nominal wages and the density of workers (see [14, ]). Therefore, we can say that the nominal wage equation is a double nonlinear integral equation (see [14, Remark 2.3] for mathematical difficulties caused by the double nonlinearity).

The replicator equation is a nonlinear integropartial differential equation whose unknown function denotes the density of workers. Its coefficient denotes the growth rate of worker population (see [1, , ] and [16, p. 73]). The coefficient is equal to the difference between the density of real wages and the average real wage, where the density of real wages is defined as the density of nominal wages deflated by a fractional power of the price index, and the average real wage is defined as the integral of the product of the density of workers and the density of real wages (see [1, (5,1), ]). Hence, the coefficient is regarded as a double nonlinear integral operator acting on the density of workers and the density of nominal wages.

Moreover, the coefficient of the replicator equation of the dcDSK system contains an unknown function implicitly in the sense that the coefficient is determined by solving the nominal wage equation under the condition that the unknown function is given, in contrast to the replicator equation whose coefficient explicitly contains an unknown function in evolutionary game theory (see [16, ]). If we can define an operator that maps the density of workers to the density of real wages by solving the nominal wage equation under the condition that the density of workers is a given function, then the replicator equation is regarded as a nonlinear integropartial differential equation whose coefficient contains the operator that acts on an unknown function.

For these reasons we deduce that the dcDSK system is an essentially new kind of system of nonlinear integropartial differential equations. Therefore, it is important to study this system not only in spatial economics but also in the theory of functional equations. In this paper we prove a sufficient condition for the initial value problem for the dcDSK system to have a unique global solution and obtain estimates of the solution. The main result is Theorem 4.

2. The System of Equations

Let us introduce the notations. Let be a domain of a Euclidean space. By we denote the Banach space of all Lebesgue summable functions of . By we denote the Banach space of all essentially bounded functions of . By we denote the Banach space of all uniformly bounded continuous functions of .

We assume that economic activities are conducted continuously in a bounded domain of an -dimensional Euclidean space, where is a positive integer. If , then the model is unrealistic from the viewpoint of economics, but we accept such a case for mathematical generality in this paper. We denote the norms of the Banach spaces and by and , respectively; that is, we define

Let . By we denote the Banach space of all functions such that

By , , , and we denote the set of all positive-valued functions of , , , and , respectively. By , , , and we denote the set of all nonnegative-valued functions of , , , and , respectively.

Let us introduce the dcDSK system. The nominal wage equation is a nonlinear integral equation of the following form (see [1, ] and [14, ]):

where is an unknown function that denotes the density of nominal wages at time and at point . By we denote the price index, which is a nonlinear integral operator of the following form (see [1, ] and [14, ]):

where represents the density of workers at time and at point . By we denote the elasticity of substitution among varieties of manufactured goods. We assume that

When increases, varieties of manufactured goods are close to perfect substitutes; as decreases toward 1, the desire to consume a greater variety of manufactured goods increases (see [1, p. 46] and [2, p. 308]). We denote the income at time and at point by . The dcDSK model consists of a monopolistically competitive sector (manufacturing) and a perfectly competitive sector (agriculture) (see [1, p. 61]). Hence the income consists of agricultural income and manufacturing income; that is, it has the following form (see [1, (5.3)]):

where and denote the share of manufacturing expenditure and the share of agricultural expenditure, respectively, and we denote the density of farmers at point by . We assume that is a given function such that (see [14, , ])

Note that this function is independent of the time variable .

The function represents the iceberg form of transport costs (see [17]). We refer to this function as the transport cost function. We reasonably accept the following assumption (see [14, Assumption 2.1]).

Assumption 1. The transport cost function is a nonnegative-valued continuous function of such that

Considering (6), and noting that the right-hand side of (4) is a nonlinear integral operator acting on the density of workers and the density of nominal wages , we see that the right-hand side of (3) is a double nonlinear integral operator acting on these densities.

We define the density of real wages at time and at point by deflating the density of nominal wages by a fractional power of the price index as follows (see [1, ], (4), and (7)):

The density of real wages can be regarded as a nonlinear integral operator acting on and .

The replicator equation is a nonlinear integropartial differential equation of the following form (see [1, ] and [13, ]):

where denotes a positive constant. We define

where is the average real wage defined as follows in the same way as [1, ] and [13, ]:

It follows from (11), (12), and (13) that workers move toward regions that offer higher real wages and away from regions that offer below-average real wages (see [1, p. 62]). Considering (4) and (10), we see that (12) is a double nonlinear integral operator acting on and . In this paper for simplicity we assume that

Hence (12) is the growth rate. The dcDSK system consists of (3), (10), and (11). In Section 4 we define an operator that maps to by solving (3) under the condition that is a given function. Substituting this operator in (11), we can transform the dcDSK system into the replicator equation whose coefficient contains the operator that maps to .

3. Result and Discussion

We consider the initial value problem by imposing the following initial condition on the dcDSK system (3), (10), and (11):

where is a given function of . This function denotes the initial density of workers. The following assumption is imposed on this function in [1, pp. 61–63] and [14, , ].

Assumption 2. (i) .
(ii) .

Hence, we accept this assumption in this paper also. Let be a constant. If a function

belongs to

and satisfies the dcDSK system for a.e. and the initial condition (15), then we say that the function (16) is a solution to the initial value problem in . If a function (16) is defined for a.e. and is a solution to the initial value problem in for each , then we say that the solution is global. No boundary condition needs to be imposed on the density of workers, since the evolution of the density of workers can be determined uniquely in (17) by the initial condition (15) as done in Section 5.

We define a function , in order to state a sufficient condition for the initial value problem to have a unique global solution. Consider the quadratic equation

where denotes an unknown quantity. It follows from (5) that this equation has a positive solution and a negative solution. We denote the positive solution by . We see easily that

By making use of this positive solution, we define the following quadratic equation:

where denotes an unknown quantity. This quadratic equation has a positive solution and a negative solution, since . We denote the positive solution by . We see easily that

The following lemma is proved in [12, Lemma 3.2] (see [12, ]).

Lemma 3. (i)    for each .
(ii)    is a strictly monotone-decreasing smooth function of for each .
(iii)    is a strictly monotone-decreasing smooth function of for each .
(iv)   , , for each .
(v)   , , for each .

The following theorem is the main result of this paper, which is proved in Sections 4 and 5 (see Assumption 1).

Theorem 4. If , , and satisfy (5), (7), and the following inequality:

then the initial value problem has a unique global solution , and this solution satisfies the following:

where , , are positive-valued functions of and

Let us discuss this theorem. The inequalities (5), (7), and (22) are a sufficient condition for the initial value problem to have a unique global solution. Note that is independent of . Making use of Lemma 3, (i), (ii), (iv), we deduce that

Hence we can say that (22) holds for a much larger set of . By Lemma 3, (i), we see that (5), (7), and (22) imply the following inequality:

Applying this inequality, Assumption 1, (5), and (7) to (33) and (34), we see that

Making use of (5), (7), (8), and the above theorem, we see easily that the right-hand sides of (3), (4), and (12) belong to for each and that the right-hand side of (11) belongs to for each . The conservation law of workers follows from (17) and (28). It follows from (27) that

Hence, no worker moves toward a point where no worker lives. Combining (25) and (26), we see that

The functions , , are defined in Lemma 17.

Remark 5. We impose (7) on the dcDSK model; that is, we consider the dcDSK model in an urban-rural setting. In this paper we cannot treat the case

that is, we cannot consider the dcDSK model in an urban setting (see [1, p. 331]) because it follows from Lemma 3, (i), that (41) cannot be substituted in (22). The DSK model with (41) is studied in [13]. The cDSK model with (41) is studied in [14, Theorem 3.2].

4. Solutions of the Nominal Wage Equation

Let us solve the nominal wage equation (3) under the condition that the density of workers is a given function. In this section we do not deal with the replicator equation. Hence, we have no need to consider the time evolution of the density of workers, the density of nominal wages, and the density of real wages. Therefore, for simplicity of symbols, we omit the time variable from these densities, and we denote them by , , and in (3), (4), (6), and (10) in this section. We refer to these equations with the same numbers. No confusion should arise. We assume that is a given function that satisfies the same condition as Assumption 2 as follows:

Proposition 6. If , , and satisfy (5), (7), and (36), then the following statements (i) and (ii) hold.

(i)The nominal wage equation (3) has a solution in .(ii)If (3) has a solution in , then the solution satisfies the following:

Proof. In [14, Theorem 3.1, (i), (ii)], from (5), (7), and (36), we prove (i), and we prove that if (3) has a solution in , then the solution satisfies (44)–(47). Applying (42), (43), and (47) to (12) and (13), we obtain (48) and (49).

Lemma 7. Let and be positive constants. Let be a constant. If , , satisfy the following inequality:

then

where

Proof. By the mean value theorem we see easily that if , then there exists a constant such that

Substituting or in this equality and making use of (50), we see that

Applying (52) to the right-hand side of this inequality, we obtain the present lemma.

Lemma 8. Assume that , , and satisfy (5), (7), and (22). If , , satisfy (42) and (43), and is a solution of (3) with , , respectively, then

where is a positive-valued function of .

Proof. Let us transform the nominal wage equation (3). Define the following new unknown function:

Defining the following nonlinear integral operator:

we rewrite (4) as follows:

Making use of this equality and defining the following nonlinear integral operator (see (6)):

we rewrite (3) equivalently as follows:

Hence, we see that

where is defined by in the same way as (56), . Subtract both sides of (61) with from those of (61) with . The right-hand side of the equality thus obtained is transformed as follows (see (6)):

where

We denote the th term of the right-hand side of this equality by , .
Substituting (58) in (46), we see that

where

It follows from (5) and Assumption 1 that

Applying this inequality, (7), (45), (56), and (64) to , we see that

where

Applying (7), (64), and (66) to , we see that

where

Making use of (5), (45), (56), and Lemma 7 with and , we obtain the following inequality:

where

Applying this inequality and (43) to (69), we see that

where

Applying (5), (43), (45), (56), (66), and Lemma 7 with and to (57) and performing the same calculations as done in proving (73), we see that

where

Combining this inequality and (64), we see that

where

Applying (7), (42), (45), and (56) to (6), we see that

Integrating both sides of these inequalities with respect to and making use of (7), (8), (42), and (43), we see that

Applying (7) and (37) to this inequality, we see that

where

Applying (66), (77), and (81) to , we obtain the following inequality:

where

Performing calculations similar to, but easier than, those done in proving (75), we see that

where

We obtain the following inequality by combining this inequality and (64) in the same way as (77):

where

Applying this inequality, (66), and (81) to , we obtain the following inequality:

where

Applying this inequality, (67), (73), and (83) to (62), we obtain the following inequality:

where

Substituting (34) in the definitions of , , we rewrite , , as follows in the same way as [12, , , , ]:

where

It is proved in [12, Lemma 7.2, ] that (5), (7), and (22) imply the following inequality:

Making use of (5), (7), and (36), we see easily that

Hence, it follows from (91) and (95) that

where

Making use of this inequality, (56), and (71) and defining

we obtain Lemma 8.

Proposition 9. If , , and satisfy (5), (7), and (22), then the nominal wage equation (3) has a unique solution .

Proof. Recalling that (36) follows from (5), (7), and (22) and combining Proposition 6 and Lemma 8 with , we obtain this proposition.

Remark 10. (i) If is not a continuous domain but a finite set of points, then the nominal wage equation is not a nonlinear integral equation but a nonlinear equation in an Euclidean space whose dimension is equal to the number of points of (see [12, (2.5)]). This subject is treated in [12]. In [12] we prove a result similar to Propositions 6 and 9 by analyzing this nonlinear equation in the Euclidean space. However, in this paper, Propositions 6 and 9 are proved in the Banach spaces in contrast to the finite dimensional proof done in [12]. Propositions 6 and 9 are similar to, but essentially different from, the results obtained in [12].
(ii) The inequalities (5), (7), and (36) are a sufficient condition for the nominal wage equation (3) to have a solution in (see Proposition 6). The inequalities (5), (7), and (22) are a sufficient condition for this solution to be unique. We make use of (22) in order to obtain (95) in the proof of Lemma 8.
(iii) It is proved in [14, Theorem 3.1, (iii)] that if and are sufficiently small, then (3) has a unique solution. However, the condition (22) is not accepted in [14]. Recalling (35), we see that the result of [14] can be regarded as a corollary of Proposition 9.

By Proposition 9 we can define an operator that maps to , where satisfies (42) and (43). We denote this operator by

that is, satisfies (3). Applying this operator to (10), we can define an operator that maps to . We denote this operator by ; that is, we define

Proposition 11. If , , satisfy (42) and (43), then

where is a positive-valued function of (see (99) for ).

Proof. For the result follows form Lemma 8. Next, we prove the case . Let us define

We define , , by (103) in the same way as (56). Making use of (75), (85), and (97), we see that

where

Substitute (58) and , , in (10). Subtracting both sides of the equalities thus obtained from each other, we obtain

Apply (5), (7), (42), (43), (45), (64), (105), and Lemmas 8 and 7 with and , , to the right-hand side of this equality. Defining

we obtain the present lemma when .

Remark 12. (i) We can prove Propositions 6, 9, and 11 with no boundary condition on , , and . We make use of (42) and (43) only. See [14, Remark 2.3, (i)].
(ii) For simplicity we omit the time variable in this section. By replacing , , and by , , and respectively, we make use of Propositions 6, 9, and 11 in the next section.

5. The Iteration Scheme

The purpose of this section is to prove Theorem 4. Let us construct an iteration scheme to obtain a solution to the initial value problem for the dcDSK system (3), (10), and (11). Let be a constant such that (see (33))

Decompose the time interval into an infinite number of intervals , , where

Let us define by the following iteration scheme (see (15)):

where we define (see (12)–(14))

Lemma 13. For each , the following (i–viii) statements hold.

(i).(ii)  for each .(iii),    for each .(iv)  for each .(v)  for each .(vi)  for each .(vii)  for each .(viii)  for each .

Proof. By (111) and Assumption 2, we can consider that (i) and (ii) hold when . Hence, making use of (115) and Propositions 6 and 9, we obtain (iii)–(viii) when . Applying (i), (ii), and (viii) with , (109), and (110) to (112) and (113) when , we obtain (i) with . Integrating both sides of (112) with respect to , recalling (12) and (13), and making use of (i) with , we obtain (ii) with . Making use of (i) and (ii) with , (115), and Propositions 6 and 9, we obtain (iii)–(viii) when . Assume that (i)–(viii) hold when , . We can prove (i)–(viii) with in the same way as in obtaining (i)–(viii) with from (i)–(viii) with . Therefore we prove this lemma.

Note that Propositions 6 and 9 are proved on the basis of (42) and (43). Hence, we need Lemma 13, (i), (ii), in order to make use of Propositions 6 and 9 in the proof of Lemma 13, (iii)–(viii).

Differentiating both sides of (112) with respect to , we obtain the following equality:

where

Lemma 14. for each , , where

Proof. Subtract both sides of (114) with from those of (114). Applying (12) and (13) to the equality thus obtained, we obtain

Applying Lemma 13, (ii), (vi), and Proposition 11 with and to the right-hand side of this equality, we see that

Applying Lemma 13, (ii), (viii), and (110) to (112), we see that

Applying this inequality, Lemma 13, (ii), (viii), and (120) to (117), we obtain the present lemma.

Lemma 15. The function

converges in (17) as for each .

Proof. Replace by in (116), where is a positive constant that satisfies the same inequality as (109). Subtracting both sides of the equality thus obtained from those of (116), we obtain

where

Performing the same calculations as done in proving (120) and making use of Lemma 13, (ii), we obtain

where

Applying this inequality and Lemma 14 to (125), we deduce that

Let us solve (123) with respect to (124) by considering (125) as a perturbation term. We perform the same calculations as done in proving Gronwall’s lemma. Replace by in (123), multiply both sides by

and integrate both sides with respect to . Applying Lemma 13, (viii), and (129) to the equality thus obtained and noting that (see (111)), we can easily obtain an integral inequality for . Consider this integral inequality as a differential inequality whose unknown function is

Solve this differential inequality with respect to . Applying the inequality thus solved to the integral inequality for , we obtain

where

Applying this inequality and Proposition 11 with to (115), we deduce that

From (132)–(134), we obtain the present lemma.

Lemma 16. The limit

is a global solution and satisfies (24), (25), and (27)–(32).

Proof. Making use of Lemmas 13 and 15, we deduce that (135) satisfies (25) and (28)–(32). Recalling (115), we can substitute (122) in (3) and (10). Let in the equalities thus obtained. By Lemma 15, we see easily that (135) satisfies (3) and (10). Replace by in (116). Multiply both sides by (130) and integrate both sides with respect to . We obtain an integral equation whose unknown function is . Apply (111) and Lemmas 14 and 15 to this integral equation, recall (12) and (13), and let . We deduce that (135) satisfies the following equality:

Hence, (135) satisfies (11) and (15) (see (14)). Applying (32) to (136), we see that (135) satisfies (24) and (27).

The following lemma gives necessary conditions for the initial value problem to have a solution that belongs to (17).

Lemma 17. If the initial value problem has a solution

in (17) for some , then this solution satisfies (23), (25), (26), and (28)–(32) in , where , (see Proposition 11).

Proof. Recalling (12), (13), and the definition of (17), we see easily that

where

Hence, integrating both sides of (11) with respect to and recalling (12) and (13), we see easily that

where

Making use of (138) and (139), we can transform (11) into (136). Applying (138) and Assumption 2, (ii), to (136), we obtain

Solving (141) with this initial condition, we see that satisfies (28) for each . Making use of this result and (139), we can apply Propositions 6, 9, and 11 to (137). Applying Propositions 6 and 9 to (137), we see that (137) satisfies (25) and (29)–(32) for each . Integrate both sides of (11) with respect to , calculate the absolute values of both sides, and integrate them with respect to . Applying (28) and (32) to the right-hand side of the equality thus obtained, we see that (137) satisfies (23) for each . Combining this result and Proposition 11, we see that (137) satisfies (26) for each .

Proof of Theorem 4. Let

be a solution that belongs to (17). Substitute this solution and (135) in (11). Subtracting both sides of the equalities thus obtained from each other, we obtain the following equation in the same way as (123):

where

Note that both (135) and (144) satisfy (25) and (28)–(32) (see Lemmas 16 and 17). Making use of this result, we can perform the same calculations as done in proving (127). Hence we obtain

Making use of this inequality in place of (129) and performing the same calculations as done in proving (132)–(134), we prove that (144) is the same as (135). By Lemmas 16 and 17 we prove Theorem 4.

Conflict of Interests

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

Acknowledgments

This work was supported by Grant-in-aid for Scientific Research (22530242), Ministry of Education, Culture, Sports, Science and Technology of Japan. The authors would like to express their deepest gratitude to the anonymous reviewer for his/her valuable suggestions.

References

M. Fujita, P. Krugman, and A. J. Venables, The Spatial Economy, The MIT Press, Boston, Mass, USA, 1999.

M. Fujita and J.-F. Thisse, Economics of Agglomeration, Cities, Industrial Location, and Regional Growth, Cambridge University Press, 2002.

M. Fujita, “The evolution of spatial economics: from Thünen to the new economic geography,” The Japanese Economic Review, vol. 61, no. 1, pp. 1–32, 2010.

P. Krugman, “Increasing returns and economic geography,” The Journal of Political Economy, vol. 99, no. 3, pp. 483–499, 1991.

P. Krugman, The Self-Organizing Economy, Blackwell, London, UK, 1996.

P. Krugman, “The official web site of the Nobel Prize in Economic Sciences,” http://www.nobelprize.org/nobel_prizes/economics/laureates/2008/

S. B. S. D. Castro, J. Correia-da-Silva, and P. Mossay, “The core-periphery model with three regions and more,” Papers in Regional Science, vol. 91, no. 2, pp. 401–418, 2012.

P. Mossay, “Increasing returns and heterogeneity in a spatial economy,” Regional Science and Urban Economics, vol. 33, no. 4, pp. 419–444, 2003.

P. Mossay, “The core-periphery model: a note on the existence and uniqueness of short-run equilibrium,” Journal of Urban Economics, vol. 59, no. 3, pp. 389–393, 2006.

P. Mossay and T. Tabuchi, “Preferential trade agreements harm third countries,” CIRJE Discussion Paper No. F-811, University of Tokyo, 2011.

N. G. Pavlidis, M. N. Vrahatis, and P. Mossay, “Existence and computation of short-run equilibria in economic geography,” Applied Mathematics and Computation, vol. 184, no. 1, pp. 93–103, 2007.

M. Tabata, N. Eshima, Y. Kiyonari, and I. Takagi, “The existence and uniqueness of short-run equilibrium of the Dixit–Stiglitz–Krugman model in an urban–rural setting,” IMA Journal of Applied Mathematics, 2015.

M. Tabata, N. Eshima, and Y. Sakai, “Existence, uniqueness, and computation of short-run and long-run equilibria of the Dixit-Stiglitz-Krugman model in an urban setting,” Applied Mathematics and Computation, vol. 234, pp. 339–355, 2014.

M. Tabata, N. Eshima, Y. Sakai, and I. Takagi, “An extension of Krugman's core-periphery model to the case of a continuous domain: existence and uniqueness of solutions of a system of nonlinear integral equations in spatial economics,” Nonlinear Analysis B: Real World Applications, vol. 14, no. 6, pp. 2116–2132, 2013.

J. Hofbauer and K. Sigmund, The Theory of Evolution and Dynamical Systems: Mathematical Aspects of Selection, Cambridge University Press, 1988.

J. W. Weibull, Evolutionary Game Theory, The MIT Press, Cambridge, Mass, USA, 1995.

P. A. Samuelson, “The transfer problem and transport costs: the terms of trade when impediments are absent,” The Economic Journal, vol. 62, pp. 278–304, 1952.