On Complex Dynamics of Differentiated Products: Cournot Duopoly Model under Average Profit Maximization

In this paper, some important dynamic characteristics such as multistability and synchronization phenomena are investigated for a game of an economic Cournot duopoly whose time evolution is received by the iteration of a noninvertible map in the plane. In the asymmetric case, the equilibrium points of game’s map are calculated, and their stability conditions are obtained.(e obtained results show that the Nash equilibrium point loses its stability through flip bifurcation. Under some restrictions, the map’s coordinate axes form an invariant manifold, and hence their dynamics are studied based on a one-dimensional discrete dynamic map. In the symmetric case where both firms are identical, the map has the property of symmetry, and this implies that the diagonal q1 � q2 forms an invariant manifold and therefore synchronization phenomena occur. Global analysis of the behavior of the noninvertible map is carried out through studying critical manifolds of the map that categorize it as Z4 − Z2 − Z0 type. Furthermore, global bifurcation of the basins of attraction is confirmed through contact between the critical curves and the boundaries of escaping domain.


Introduction
Dynamic economic games whose evolutions has been depended on discrete time maps have attracted many researchers due to the interesting dynamic behaviors and bifurcations types occurring in such maps. Analyzing the dynamic characteristics of such maps has been dated at least to [1] by Puu who passed away in 2020 after enriching this direction with many useful books and papers. In the current paper, we focus on a special type of these maps that are known as duopoly maps. In duopoly game, the market possesses only two competing firms (or players) whose target is seeking the optimal quantities in the case of Cournot by maximizing their profits, or seeking the optimal prices in the case of Bertrand. Cournot, the French economist who started the core of such kind of games in 1838, introduced a duopoly model adopting quantities as decision variables. Later in 1880s and based on Cournot model, Bertrand suggested his model known as Bertrand model and considered prices as competitive means. In 1950, the Nash equilibrium point detected by the famous American mathematician John Nash has been raised in the equilibrium theory of noncooperative games. In general, this theory provided research with significant tools in studying oligopoly models. As a result, Nash point became important aspect in studying such models and was found in several Cournot and duopoly games. Indeed, one has to highlight game theory as an important tool of applied mathematics that was used and adopted to study the interaction among players in those economic games. It has been used in studying the strategic interactions among rational competitors and decisionmakers. Game theory was extensively developed in the 1950s by many scholars and was explicitly applied to all fields of social science, as well as logic, biology, economics systems, and computer science. Recently, the theory of game has extensively been applied in economics for studying competitions among competing firms. Even though static equilibrium analysis of such games is limited and has not provided valuable information on games' evolutions, it has been investigated by many researchers (see, for example, [2]). Studying the game evolution through dynamics gives more information on the stability/ instability of Nash equilibrium, the types of bifurcations by which it loses its stability, and the topological properties of basin of attraction in the phase plane. Rand [3] was the first to introduce and analyze some of these dynamic characteristics.
Firm in duopoly competition simultaneously chooses its decision or strategy taking into account two important factors: its own actions and actions by its competitor. In literature, several works have been reported for studying the dynamic characteristics of duopoly games. ese dynamic characteristics included bifurcation and its types, either flip or Neimark-Sacker; peculiar shapes of attractive basins for attracting sets and chaotic attractors; and multistability phenomena. Furthermore, studying such games requires adopting specific kinds of utility consumer preferences. e most popular one that is well-known as Cobb-Douglas utility was adopted in many studies in literature. e reason is that Cobb-Douglas utility is a particular functional form of the production function. It has been widely used to represent the technological relationship between the inputs and outputs that can be produced by those inputs. Moreover, literature reported many other utility functions such as constant elasticity of substitution (CES) and Singh and Vives utility. One can get more information and properties of those utilities from literature ( [4][5][6][7][8][9][10]). e dynamic characteristics of duopoly models are complex and reveal interesting information on their evolutions. erefore, many interesting results on these complex dynamics have reported many interesting results on the stability of Nash equilibrium point of these games and opened the routes to attract many researchers. Some related works and results for these models are reported here. For example, a game of two competing firms whose decision variables are different (Cournot duopoly game) has been analyzed in [11]. Other studies on duopoly and triopoly in which firms seek the optimality of production have been introduced and investigated in ( [12][13][14]). In [15], a Cournot duopoly game whose players are rational, adopting bounded rationality rule and seeking maximization of their relative profits, has been introduced and analyzed. Using a theoretical framework, a Cournot-Bertrand game has been formalized in [8]. In [16], Cournot-Bertrand model has been studied by Naimzada et al. using a two-dimensional discrete linear map. ey focused on investigating the dynamic characteristics of Nash equilibrium points and its instability through two types of bifurcations. Additionally, analyzing the game's evolution requires representing it by a discrete time map so that firms can update of their productions.
is requires some kinds of adjustment rules. Literature has presented and reported many of these rules. e first rule among those rules which has been extensively adopted is known as bounded rationality. It has been deeply adopted in many studies for process of modeling discrete time maps representing the evolutions of these games. It is described as a gradient rule because it depends on firms' marginal profits and requires firms to carry out an estimation of their marginal profits to determine whether they increase or decrease, so they can update their outputs next period of time. ere are also other mechanisms which have been used in the modeling process such as tit-for-tat rule and local monopolistic approximation mechanism ( [17][18][19][20][21][22][23][24][25][26][27]).
Based on a linear inverse demand functions and adopting the bounded rationality mechanism, a Cournot duopoly game with players seeking the maximization of their average profits is introduced and analyzed in this manuscript. Most firms do not rank profits as the major goal. e work of modem firms is so complex that they do not think merely about profit maximization. e basis of the difference between the objectives of the neoclassical firm and the modem corporation arises from the fact that the profit maximization objective relates to the entrepreneurial behavior while modem corporations are motivated by different objectives because of the separate roles of shareholders and managers. Due to the above discussion, we have introduced a general case with different objectives. e game is described by a nonlinear discrete dynamic map that possesses four equilibrium points, three of which are boundary points, and the fourth is an interior point and coincides with Nash equilibrium point. We use in this manuscript local and global analysis to investigate the dynamics of these equilibrium points. e analysis shows that Nash point loses its local stability through flip bifurcation. e structure of the game's map shows that it is being attracted to the origin, and hence its coordinate axes form an invariant manifold. Consequently, the map's dynamics are studied by a onedimensional map conjugating the famous logistic map. It is also shown that under some restrictions, the synchronization phenomena are studied and show that the diagonal in map's phase plane forms an invariant manifold. Furthermore, our findings show that the game's map is noninvertible and its phase plane is divided into three zones that make it belong to Z 4 − Z 2 − Z 0 type. e sections in the current paper are summarized as follows. In Section 2, the discrete form of the game's map is introduced. e equilibrium point and its stability are investigated both analytically and numerically through local analysis in Section 3. In Section 4, the invariant manifold is studied. e basin of attraction and critical curves are analyzed in Section 5. In Section 6, the global analysis and synchronization phenomena are investigated based on some numerical experiments.

The Model
Suppose an economic market consisting of two competing firms (or players) whose decision variables are quantities, q 1 and q 2 . Both firms adopt linear inverse demand functions given by e horizontal differentiated products mainly refer to difference in extrinsic features between the products, like color and shape. Such differentiation mostly adapts to differences in consumer preferences. In (1), the parameter b refers to the degree of horizontal product differentiation that is taken in the interval [0, 1]. For b � 0, each firm behaves as a monopoly player due to the independent varieties between firms' products. b � 1 means that both products are perfect substitutes or 2 Discrete Dynamics in Nature and Society are homogenous, and hence p 1 � p 2 . In this paper, we focus on the case of b ∈ (0, 1). e parameter a denotes maximum price in market. Suppose that the competing firms adopt different nonnegative marginal costs and hence their total costs are taken as C i (q i ) � c i q i , i � 1, 2 (c 1 ≠ c 2 ). Now, the total profit for each firm takes the following form.
Unlike many studies in literature, the current game assumes that both firms focus on maximization of their average profits rather than profits only. e average profits then become (3) e parameter r i , i � 1, 2, where 0 < r i < 1, is directly used to measure altruism. at is to say, it explains the firms' care of average performance. When r i � 0, i � 1, 2, it means firms give care to their own profits only. Substituting (4) into (5), one gets and then the marginal average profits (zφ i /zq i , i � 1, 2) become Maximizing the average profits requires firms to know some characteristics of market and the nature of competition. Knowing the characteristics of market needs complete information on the market and the strategies of firms' rivals, and these become impossible. erefore, firms believe that updating output productions requires watching the changes taking place in the marginal average profits. at is to say, if zφ i /zq i > 0, this means firms will increase their production next period of time; otherwise, they decrease productions or become naive. Accordingly, firms estimate production based on bounded rationality rule that has the following form.
Adopting bounded rationality rule makes us assume that which makes relative production directly proportional to the firms marginal average profits. Equation (6) possesses some important characteristics that can be found in literature ( [25][26][27][28]). Substituting (5) into (6), one simply has a two-dimensional discrete map that will be used throughout the manuscript to discuss the game's repetition (or evolution).

Local Stability
In order to study the stability of the fixed points (8), we recall the Jacobian matrix at these points, and hence the following propositions are raised (their proofs are given in Appendix).

Proposition 1.
e point e o is an unstable repelling node.

Proposition 2.
e point e 1 is a stable point provided that Discrete Dynamics in Nature and Society

Proposition 3.
e point e 2 is a stable point provided that

Proposition 4.
e point e * that is called Nash point is locally stable provided that Proposition 5. e point e * loses its stability due to flip bifurcation only.
e obtained results show that the stability condition of Nash point depends on the parameters a, b, c 1 , c 2 , r 1 , r 2 and the speed of adjustment parameters v 1 and v 2 . Some numerical experiments are carried out here to validate the theoretical results given above. As in many studies in literature, we keep our focus on studying the influences of speed parameters on the stability of Nash while fixing the other parameters. Let us suppose the set of parameters a � 1, b � 0.35, c 1 � 0.4, c 2 � 0.3, r 1 � 0.7, r 2 � 0.5, and v 1 � 1 (or v 2 � 3.852) on varying v 2 (or on varying v 1 ). Figures 1(a) and 1(b) show the one-dimensional period-doubling bifurcation diagram on varying v 1 with respect to the two decisional variables q 1 and q 2 . As one can see, both figures present flip bifurcation diagrams. In Figures 1(c) and 1(d), one can see the flip bifurcation on varying the speed parameter, v 2 . In Figure 1(e), the phase plane of period-8 cycle with Nash equilibrium point is depicted. Figure 1(f) shows a one-piece chaotic attractor. In Figures 1(g) and 1(h), a two-piece disconnecting chaotic attractor is given. In these figures the marginal costs are fixed.

The Invariant Manifold
Looking at the map (7), one can see that the point (0, 0) is trapping the map. at is to say, if q 1 (t) � 0 or q 2 (t) � 0, one gets q 1 (t + 1) � 0 or q 2 (t + 1) � 0, and hence both coordinate axes become invariant axes for the map. erefore, the invariant set made of the axes q 1 and q 2 forms an invariant manifold for T. Accordingly, the dynamics of the map (7) on the invariant axes can be described by the one-dimensional map given by It is easy to see that the map (10) is topologically equivalent to through the transformation given below. and

Dynamic Analysis.
For equation (11), suppose the following function.
Function (13) shows the logistic equation, and μ is the parameter responsible for its dynamics. It has a first derivative given by σ � (u) � μ(1 − 2u). Its maximum value occurs at μ/4 at the point u � 1/2. In addition, σ(0) � 0 and σ(1) � 0, which means that σ(u) ∈ [0, 1) on the open interval 0 < μ < 4 and for all u ∈ [0, 1]. Moreover, it is easy to see that (11) has only two fixed points, that is, u � 0 and In addition, one can see that the map (11) has two fixed points, u � 0 and and then u � 0 is stable if μ ∈ (0, 1) and unstable otherwise. e other fixed point e simulation in Figure 2(a) presents a 1D bifurcation diagram on varying μ that confirms that the point u � 1 − 1/μ is stable in μ ∈ (1, 3). In Figure 2(b) a cobweb diagram is given at μ � 2.99 in the stability interval and shows a stable fixed point, while in Figure 2(c), an unstable behavior is plotted at μ � 3.75, and hence the instability of the fixed point is confirmed. Proof.
e proof is straightforward as μ i � 1 + (a − c)] i , and our discussion given above confirms the topological equivalence between the map (7) and (11) on which the trajectories of (11) diverge when μ i > 3, i � 1, 2.
and r 1 � r 2 � r is studied. As one can see, the map T defined in (7) is symmetric. at is to say, its structure does not change if the decision variables q 1 and q 2 are swapped; i.e., T ∘ S � S ∘ T, where S: is means that the diagonal defined by Δ � (q 1 , q 2 ): q 1 � q 2 forms an invariant manifold. It implies that the dynamics of any trajectories starting on the diagonal (q 1 (0) � q 2 (0)) will get back to the diagonal for every t. In this case, the dynamic behavior of map (7) under   Discrete Dynamics in Nature and Society the restriction imposed on Δ will be described by the onedimensional map given below.
at is to say, any synchronized trajectories (i.e., q 1 (t) � q 2 (t) for every time step t) will be governed by T Δ : Δ ⟶ Δ, and the following proposition is raised.

Proposition 7.
e map (14) is unimodal and concave, When the model becomes unimodal, one gets the following critical point.
Proof. e proof is straightforward by looking at the function g(q) that has a global maximum without any local maximum, and hence it is unimodal. e sign of the second derivative of g(q) proves that the map is concave and finally g ′ (q cr ) � 0. Now, the map T Δ has only one fixed point, that is, q, and this point is locally stable if ] < 2/a − c while flip bifurcation takes place at v � v flip : � 2/a − c. Figure 2(d) shows a 1D bifurcation diagram where flip bifurcation takes place at v ≈ 3.33.

Basin of Attraction
Here, some properties of the map (7) regarding its attractive basin are given. According to [29], it is possible to see that the coordinate axes and their preimages of any rank construct boundaries for initial conditions of any nondiverging trajectories. It can be seen that the dynamics of map (7) can be restricted to the invariant axes and can be rewritten in the one-dimensional map as follows.
which possesses a positive fixed point, It is easy to see that g i (q i ) is concave and increasing function provided that q i < 1 + ] i (a − c i )/4v i . In addition, for any points in the form is concave and unimodal. In the last case, one can get the nonnegative points q bound is means that trajectories commencing close to these invariant axes with positive initial states are repelled by those axes. In Figure 3(b), the line segments ω i , i � 1, 2, and their inverses ω − 1 i , i � 1, 2, are depicted at the set of parameter values a � 1, b � 0.35, c 1 � 0.4, c 2 � 0.3, r 1 � 0.7, r 2 � 0.5, v 1 � 4.5696, and v 2 � 3.852 where the point O − 1 3 is the fourth preimage point of the origin and is given by It is also noted that this point is the intersection of both inverses ω − 1 1 and ω − 1 2 .

Critical Curves.
Besides the fact that the map (7) is nonlinear, it is also noninvertible map, which means that the rank-1 preimages of a point (q � 1 , q � 2 ) may not exist or may be multivalued. As one can see, the origin point possesses four rank-1 preimages, and hence the phase plane of map (7) includes Z 4 zone, the zone whose points have four rank-1 preimages. It may also include Z 2 , the zone whose points possess two rank-1 preimages, and Z 0 , the zone with no rank-1 preimages. To investigate these zones, let us calculate the rank-1 preimages for the points along the invariant axes and take the forms (p, 0) and (0, q). Substituting these points into the map (7), one gets the following algebraic systems.
Solving the second equation of (18a) gives the following: , and then one gets two case. Case 1. Substituting q 2 � 0 in the first equation of (18a) gives the following quadratic equation.
is means that the point in the form (p, 0) may have no preimages or two distinct real-1 preimages.
, one gets the algebraic system Solving (19b), we get the following quadratic equation.
Discrete Dynamics in Nature and Society is means that in this case, the point (p, 0) may have no preimages or two distinct real rank-1 preimages. Combining the two cases, one can conclude that such points in the invariant axes may have four, two, or no preimages. For the interior points in the form (p, q) ≠ (0, 0), it is complicated to categorize such points as belonging to one of the zones Z 0 , Z 2 , or Z 4 due to the difficulty of solving the algebraic system (18a) at this point. For this reason, numerical simulation is used to calculate critical curves. e critical curves are responsible for dividing the map's phase plane into these zones. e critical curve of rank-1 is denoted by LC and forms a locus of coincident points of rank-1 preimages located on a set denoted by LC − 1 . Since the map (7) is continuously differentiable, LC − 1 can be defined as the locus of points where the associated Jacobian determinant is vanished; i.e., LC − 1 ⊆ (q 1 , q 2 ) ∈ R 2 : det(J(q 1 , q 2 )) � 0 , where J(q 1 , q 2 ) is defined by (A.1) in the appendix. erefore, LC represents the rank-1 image of LC − 1 under the map T given in (7); i.e., LC � T(LC − 1 ). e expression for det(J(q 1 , q 2 )) � 0 is cumbersome and is omitted here, and we instead do numerical simulation for presenting both LC − 1 and LC. Figure 3(c) shows that LC − 1 consists of two branches LC a − 1 and LC b − 1 while Figure 3(d) presents the critical curve that is also formed by two branches, LC a and LC b . It is also clear that the map (7) is not invertible as its phase plane is divided into three zones, Z 0 , Z 2 , and Z 4 .

Global Analysis
In economic models, global analysis is required to investigate the topological structure of basin of attraction that is not recognized by local analysis and whose qualitative changes occur due to change in parameter values. e long-term behavior of variables of these economic models given initial   conditions is indeed necessary in order to show the influences that may result if these initials are taken far away from a fixed point of an attracting set. e numerical analysis shows two important aspects for the map (7), the asymmetric and symmetric cases. Let us first begin analyzing the asymmetric case by adopting the following set of parameter values: a � 1, b � 0.35, c 1 � 0.4, c 2 � 0.3, r 1 � 0.7, r 2 � 0.5, v 1 � 3.3285, and v 2 � 3.852. is set is given in Figure 4(a) which represents a period-8 cycle with Nash equilibrium point and its attractive basins. ese basins are colored orange and green while the escaping (or infeasible points) domain is plotted in gray color. It is also clear that the escaping domain forms a connecting set. By increasing the speed parameter of the first firm to v 1 � 4.03 while keeping the other values fixed, the dynamic of map (7) turns to a chaotic attractor that consists of two disconnecting areas as shown in Figure 4(b). It can be observed that the branch LC b of the critical curve LC is moving upward. With further increase in this parameter to v 1 � 4.4675, a one-piece chaotic attractor emerges in Figure 4(c). Its attracting domain is displayed in white color and is separated from the infeasible domain by the segments ω − 1 1 and ω − 1 2 . At v 1 � 4.647, with the other parameter values being fixed, the chaotic attractor continues to appear, the critical curve LC b becomes tangent to the boundary ω − 1 2 , and hence contact bifurcation is raised as given in Figure 4(d). is contact takes place only when one branch LC b of the critical curve which separates Z 2 zone from Z 0 zone becomes tangent to ω − 1 2 . is contact with the boundary of infeasible domain occurs from one side due to asymmetry between the two branches of map (7). is contact bifurcation becomes more visible in Figure 4(e) at the parameters set ] 1 � 4.63 and ] 2 � 3.95 while the other parameters are fixed. As one can see, the contact bifurcation causes disconnection in the infeasible domain because an area h 0 from Z 0 zone moves into Z 2 zone; i.e., is kind of contact bifurcation is well-known (see [29]). Just after it, the area h 0 from the infeasible domain enters inside the region Z 2 . As a result, the points belonging to this region have now two distinct preimages. Since this region belongs to Z 0 , their preimages will belong to the same region and form holes lying in the attraction domain that is bounded by the quadrilateral shape Indeed, only at ] i � 2/(a − c i ), i � 1, both h 0 and its preimage h − 1 0 (that is formed by two areas connected by the branch LC b − 1 ) have a contact with coordinate axes. It is also clear that both h 0 and h − 1 0 are entirely located in the region Z 2 . erefore, each of them has further preimages that form smaller holes of the infeasible domain located in the attracting domain and are bounded by rank-2 preimages of the segments ω 1 and ω 2 .
In contrast, the symmetric case is now investigated. is case happens at v 1 � v 2 � v, c 1 � c 2 � c, and r 1 � r 2 � r. is means the map (7) can be rewritten in the following form.
e map (20) is symmetric on the diagonal Δ, and studying its dynamics through local analysis around Δ is related to synchronization (intermittency). In this regard, one can say that a dynamic behavior of the map is synonymized if there exists a step time � t such that q 1 (t) � q 2 (t) for any t > � t. erefore, in this case, the dynamics can be studied by recalling the Jacobian of map (20) given by where, ℓ(q) en, the eigenvalues can be given as follows.

Proposition 8.
e point e * is locally stable provided that v < 2/a − c. At v � v f � 2/a − c, the point e * undergoes flip bifurcation.
Proof. Substituting e * into (22), we get Simple calculations make Combining the two conditions on ] gives v < 2/a − c, and this completes the proof.
From the above proposition, one can deduce that the parameter v plays a key role in the stability of e * . is means that high reactivity of each firm with respect to a marginal change in its average profit should exist. Furthermore, the importance of the above proposition comes from its generalization when an attracting cycle for the map T * on Δ exists. To see that, suppose that a k-cycle given by  (q 1 , q 1 ), . . . , (q m , q m ) of the map T * is embedded in the invariant diagonal Δ which corresponds to the cycle q 1 , . . . , q m . en, the multipliers for this k-cycle are given by whose eigenvectors are (1, 1) and (1, − 1). As shown, the stability of e * is governed by |λ ⊥ | < 1 and is confirmed by |λ ‖ | < 1. In contrast, a complex chaotic attractor H for the map T * says that H is asymptotically stable if and only if all trajectories belonging to H are transversely attracting. is asymptotic stability condition of H is given in terms of transverse Lyapunov exponent as follows.
where q(0) ∈ H and q(i) is the trajectories generated by the map e following definitions are given [29]. A spectrum of Lyapunov exponents can be defined according to the initial conditions as follows [30].
where Λ nat ⊥ denotes a Lyapunov exponent evaluated at a generic trajectory in the chaotic attractor H. If Λ max ⊥ is nonpositive, then a set is called Lyapunov attractor. While Λ max ⊥ > 0 and Λ nat ⊥ < 0, the set becomes no longer Lyapunov stable and becomes Milnor attractor. It is defined as follows.  (20) is represented by period-4 cycle given in Figure 4(f ). Its attractive basin and the basin of Nash equilibrium point are colored green and orange, respectively, and are bounded by ω 1 and ω 2 and their inverses. e infeasible domain of this cycle is given in gray color. By keeping the same set of values and increasing v to 4.2199558, four closed invariant rings emerge around the diagonal as depicted in Figure 4(g). All these rings are symmetric around the invariant diagonal Δ and lie within the region Z 2 only when two of them are tangent to the branch of critical curve LC b . By further increasing ] to 4.2818322, the four rings convert into a period-20 cycle as shown in Figure 4(h) whose basin of attraction is colored cyan while the basin of Nash point is denoted by orange color. At ] � 4.414216, while the other parameter values are fixed, a Milnor attractor is raised on the invariant diagonal as presented in Figure 5(a). is means the dynamics of map (20) are captured by a two-piece chaotic attractor on the invariant diagonal Δ that is in Milnor sense. is indicates that any cycle embedded in the diagonal Δ is transversely unstable, and hence all trajectories approaching that diagonal are exploded away by transversely unstable cycles. By further increasing ] to 4.763885, a one-piece chaotic attractor emerges as shown in Figure 5 is chaotic attractor is symmetric on the diagonal because synchronization still occurs. e symmetry of this chaotic attractor is  given in Figure 5(c) where its boundary at the same set of parameter values is given by the critical curves LC a , LC b , LC a 1 , LC b 1 , LC a 2 , and LC b 2 . In Figure 5(d), the transverse Lyapunov exponent is given. It gives better understanding on the map's dynamics in the long term. As one can see, the bifurcation diagram corresponding to the transverse Lyapunov exponent confirms the existence of chaotic attractor on the invariant diagonal in the interval ] ∈ [4.4, 4.7] except for windows of some periodic cycles. It is also noted that in this interval of the speed parameter ], the Lyapunov exponent changes its sign many times until it eventually becomes positive for large values of the speed parameter ]. From economic perspective, this implies that the long-term dynamics of the map (20) can be very sharp due to small changes in this parameter. In other words, firms commencing with the same initial states may be different in the long term.

Conclusion
In this paper, some dynamic characteristics such as multistability, synchronization phenomena, and complex structure of some attractive basins are investigated for a nonlinear two-dimensional noninvertible map describing a Cournot duopoly game with an invariant one-dimensional manifold where synchronization takes place. It has been proved that, with identical firms, the game's map becomes symmetric and both coordinate axes and diagonal q 1 � q 2 form invariant manifolds. e main results in this paper concern the global analysis in the asymmetric and symmetric case where synchronization has taken place. Such global analysis has focused on the properties of critical curves of the noninvertible map. e phase plane of game's map has been divided into three zones of preimages, that is, Z 0 , Z 2 , and Z 4 . We have shown that for a set of parameter values, the one-   1 − τ + δ � 4ABv 1 v 2 + b 2 AB 1 + r 1 1 + r 2 , (A.6a) (A.6c) Since a > c 1 , c 2 and a > b, then the multiplication AB is always nonnegative, and then condition (A.6a) is always positive. Suppose that conditions (A.6b) and (A.6c) are nonnegative; then, combining them gives 0 < 2Av 1 + 2Bv 2 − 4ABv 1 v 2 + b 2 (1 + r 1 )(1 + r 2 )AB < 4, and hence Nash point is stable. e proof is completed. □ Proof of Proposition 5. Jacobian (A.1) at Nash point e * has two eigenvalues taking the following form.