- About this Journal ·
- Abstracting and Indexing ·
- Aims and Scope ·
- Annual Issues ·
- Article Processing Charges ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Recently Accepted Articles ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Subscription Information ·
- Table of Contents
Abstract and Applied Analysis
Volume 2013 (2013), Article ID 323290, 11 pages
Nonlinear Dynamics in a Cournot Duopoly with Different Attitudes towards Strategic Uncertainty
1Department of Economics and Management, University of Pisa, Via Cosimo Ridolfi 10, 56124 Pisa, Italy
2Department of Law, University of Genoa, Via Balbi 30/19, 16126 Genoa, Italy
Received 2 September 2013; Accepted 12 November 2013
Academic Editor: Luca Guerrini
Copyright © 2013 Luciano Fanti et al. 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.
This paper analyses the dynamics of a duopoly with quantity-setting firms and different attitudes towards strategic uncertainty. By following the recent literature on decision making under uncertainty, where the Choquet expected utility theory is adopted to allow firms to plan their strategies, we investigate the effects of the interaction between pessimistic and optimistic firms on economic dynamics described by a two-dimensional map. In particular, the study of the local and global behaviour of the map is performed under three assumptions: (1) both firms have complete information on the market demand and adjust production over time depending on past behaviours (static expectations—“best reply” dynamics); (2) both firms have incomplete information and production is adjusted over time by following a mechanism based on marginal profits; and (3) one firm has incomplete information on the market demand and production decisions are based on the behaviour of marginal profits, and the rival has complete information on the market demand and static expectations. In cases 2 and 3 it is shown that complex dynamics and coexistence of attractors may arise. The analysis is carried forward through numerical simulations and the critical lines technique.
In this paper, we analyse the dynamics of a Cournot duopoly under strategic uncertainty with pessimistic and optimistic firms within the framework of a nonlinear dynamic oligopoly as those developed by a recent burgeoning literature (see  and the papers cited therein).
The issue of decision making under uncertainty as distinct from risk has recently been revisited, amongst others, by [2–6]. In these papers, strategic uncertainty is represented by means of the Choquet expected utility (CEU) theory  where agents exhibit different attitudes towards uncertainty, that is, pessimism or optimism, overweighing less or more uncertain events. This theory has also been adopted by  to represent strategic behaviour à la Cournot with different firms’ attitude towards uncertainty (CEU theory has been applied to other economic contexts, where optimism and pessimism can explain the paradox of people buying insurance and gambling, the equity premium puzzle, and the small stock puzzle ).
In a strategic context such as a duopoly game, it is crucial to forecast the behaviour of the competitor in order to make a decision and to specify the information set available to each player. In the literature on nonlinear oligopolies, two distinct assumptions with regard to available information are usually made: players have a complete knowledge of the market demand and use some form of expectations about the rival’s strategic variable decision (e.g., naive, rational, or adaptive expectations or, alternatively, some weighted sum of previous rules) to set the price or the quantity in the future period (e.g., [9, 10]); players have limited information about the market demand and use some forms of estimation of their own current marginal profits (e.g., [11–14]) or other adjustment mechanisms such as the local monopolistic approximation to determine the price or quantity in the future period . This is because, under the hypothesis of limited information, players are unable to solve the optimisation problem by accounting for expectations about the value of the strategic variable that the competitor will choose for the next period, but they are able to get either a correct estimate of the slope of their own profit function in the current period, that is, the partial derivative of the profit function computed at the current state of production, or use a linear approximation of the demand function by market experiments without any guess about the influence of the competitors (i.e., monopolistic approximation).
In this paper we study local and global dynamics of a Cournot duopoly model with strategic uncertainty as in  by considering different information sets of players. In particular, both firms have complete information on the market demand and use the “best reply” adjustment mechanism (static expectations) to vary production period by period; both firms have incomplete information on the market demand and adjust production by following a mechanism based on marginal profits; one firm has incomplete information and production decisions are based on the behaviour of marginal profits, and the rival has complete information and static expectations. These assumptions make the topological structure of the map different and then comparing the local and global properties of the three dynamic systems is relevant. In particular, in cases 2 and 3 we find that complex dynamics as well as coexistence of attractors may occur. These phenomena depend on the relative value of the parameter that weights the strategic uncertainty of firms. The analysis is performed by applying the critical lines techniques as well as through numerical simulations.
The rest of the paper proceeds as follows. Section 2 develops the Cournot model with strategic uncertainty. Section 3 studies the dynamics of the model under complete information and static expectations of both firms. Section 4 introduces the adjustment mechanism of production based on marginal profits for both firms and performs the local and global dynamics. Section 5 analyses the mixed case. Section 6 outlines the conclusions.
2. The Model
We consider a Cournot duopoly for a single homogenous product with a liner negatively sloped inverse demand given by where , denotes the consumers’ marginal willingness to pay for product , and (resp., ) is the output produced by firm (resp., firm ). The average and marginal costs for each single firm to provide one additional unit of output in the market are equal and given by .
Profits of firm , , can be written as follows: We assume that firms are risk neutral and their single productive capacity is able to satisfy the whole market. With regard to strategic uncertainty, each firm does not know whether the rival behaves as a quantity setter à la Cournot. However, each firm has some beliefs about rival’s behaviour and such beliefs are different because each firm has a specific attitude, that is, it is either pessimistic or optimistic, towards uncertainty. In particular, we assume that the pessimistic firm assigns a positive probability of being in the worst case (which is realised when the market price equals zero) (the firm considers it possible that the supply of the other firm is large enough to get the price to zero), while the optimistic firm assigns a positive probability of being in the best case, (which is realised when it behaves as a monopolist in the market).
By using CEU theory, it is assumed that each firm maximises its own CEU function which is given by a weighted average (with the parameter ) of its expected profits and the profits in the worst case (resp., best case) for the pessimistic (resp., optimistic) firm, where (assumed to be the same for both firms) represents the degree of confidence that each firm has about the Cournot behavioural assumption [5, 8]. An economic interpretation of is the following : values of close to zero may represent a relatively new market characterised by a high level of strategic uncertainty (in which pessimistic firms are more likely to decide to do not produce because the level of uncertainty may be too high for them); values of close to one may represent mature markets characterised by a relatively stable environment, where firms are relatively confident about their conjecture.
Let firm (resp., firm ) be the pessimistic (resp., optimistic ) firm. They have to solve the following problem of maximisation of expected utility, respectively: From (3) we see that if the pessimistic firm will expect , while the optimistic firm will expect monopolistic profits given by . When the model boils down to the standard Cournot model without strategic uncertainty. Maximisation of (3) gives the following marginal CEU for firm and firm , respectively:
From (4) and (5), if we get the following reaction functions of the pessimistic and optimistic firms; that is: The interior Nash equilibrium of the game therefore is obtained as follows: From (4) and (5), we note that if no interior Nash equilibrium exists. Thus, for economic reasons we impose the following.
Assumption 1 (). In addition, from (8) the condition to guarantee that the quantities are positive is as follows: Then, we introduce the following.
Assumption 2 (). We note that the values of price and profits corresponding to Nash equilibrium (8) are given by which are positive without imposing any other conditions than (9). In fact, for (11) we note that the condition implies (being for all ). This is the solution of the static game. In the following sections, we study the Cournot duopoly model with strategic uncertainty from a dynamic point of view. In particular, we introduce dynamic adjustment mechanisms to look at whether the Nash equilibrium represents the long-term stable allocation of the market or nonconvergent dynamics exist. With this regard, Section 3 introduces the case in which players have complete information and static expectations as in . Section 4 analyses the case in which both firms have limited information about the market demand and use an adjustment mechanism of production based on an estimate of their own marginal CEU period by period [17, 18]. Section 5 considers the case in which one firm has limited information and the rival has complete information about the market demand with static expectations , that is, the “best-reply” dynamics. This is because under heterogeneous adjustment mechanisms the map has a topological structure different than that when both firms adjust production period by period by using an estimate of their own marginal CEU.
3. Dynamics under Complete Information (“Best Reply” Dynamics)
One of the first dynamic adjustment mechanisms studied in the literature on nonlinear oligopolies is the one proposed by , which is based on firms’ reaction functions. In this case, players play their best replies by assuming that the rival does not modify production with respect to the previous period. We now therefore introduce time, which is discrete and indexed by and assume that players have complete information about the market demand and use static expectations to adjust production period by period. Then, by using (6) and (7) the two-dimensional map describing the dynamics of the economy is as follows: where is the unit-time advancement operator; that is, if the right-hand side variables are defined at time , the left-hand side ones are defined at time . In order to make the reading easier, we recall that , , and .
Since the market demand is linear and average (and marginal) costs are constant, map (12) is piecewise linear and Nash equilibrium (8) is the unique interior fixed point of the map. We note that is not a fixed point of map (12). In fact, given the subsequent iterate leads to a positive value of . In addition, corner fixed points and are avoided by Assumption 2 and the max operator in (12), respectively. In what follows, we focus on dynamics starting from the set that is, dynamics that start with a positive price. We note that considering initial conditions that lie on set does guarantee that the dynamics remain on for every iteration. In fact, we have that and , and thus also the inequality holds. In particular, if the dynamics lie on for any , where is the number of iterations, then dynamics are described by where are fixed for a given initial condition . However, depending on initial conditions it is possible that an iterate lies on the border of . This implies that at least one of the two firms does not produce. In any case, this is a temporary result since the dynamics definitely lie on and they are captured by Nash equilibrium (8) which is the global attractor of map (12) on .
Proof. First of all, we focus on
which is obtained by (12) by removing the nonnegative constraints. It means that if for then for . We note that according to the signs of parameters of the model we have that map is a contraction; that is there exists such that for any . Thus, fixed point (8) is the global attractor of . We now prove that also is a contraction; that is, there exists such that for any . Let us consider . Then, we have and . The following cases are possible.
If , then and , and the property for follows from the property stated for .
If ,, then The result follows by assuming . A similar argument may be applied for the case ,.
If , then The result follows by assuming . A similar argument may be applied for the case ,.
If , then The result follows by assuming . A similar argument may be applied for the case .
Then, we can conclude that fixed point (8) is the global attractor for on . From the inequality for any , we have the result.
In addition, we note that in contrast with the standard Cournot game (), under strategic uncertainty () it is more likely that the pessimistic (resp., optimistic) firm decides not to produce (resp., to produce) a positive quantity at a certain date (as shown in Figure 1).
4. Dynamics under Limited Information
This section studies the dynamics of the Cournot model with strategic uncertainty by using an adjustment mechanism of production introduced by  in a model with continuous time and used by  in a model with discrete time. With regard to the information set of players, in this section we assume that both firms have limited information (no knowledge of the market demand) as in [11–14]. In order to overcome this informational lacuna, we assume that at any time each player uses an adjustment mechanism based on local estimates of its own marginal CEU at time () to determine production at time . The adjustment mechanism is as follows: where is a coefficient that captures the speed of adjustment of each firm’s quantity with respect to a marginal change in its marginal CEU and is the intensity of the reaction of each player. Therefore, in this case the pessimistic firm and the optimistic firm increase or decrease their production at time depending on whether and are positive or negative, respectively. This type of adjustment mechanism implies that although players have incomplete information about demand and cost functions, they are able to get a correct estimate of their marginal CEU in the current period.
By using (4), (5), and (6), the two-dimensional system that characterises the dynamics of a Cournot duopoly, under strategic uncertainty and incomplete information is as follows: From a mathematical point of view, map (20) is defined on the whole plane but for economic reasons only the dynamics that lie on set are meaningful (i.e., ), where is the number of iterations and the last inequality guarantees a positive value of the market price along the trajectory. We note that starting from initial conditions that lie on set does not guarantee that the dynamics remain on for every iteration (i.e., ). To avoid this problem,  and  introduce nonnegativity constraints on the mechanism of adjustment. This allows having well-defined trajectories by starting on whatever feasible initial condition from an economic point of view. By adapting this idea to the model under scrutiny we have that the map defined on set becomes as follows: Although in this paper we concentrate on the study of map (20), this does not represent a loss of generality, since it is possible to show that trajectories that exit from coincide with trajectories of map (21) for which there exists a such that for any at least one of the two variables is zero; that is, one of the two firm exits from the market (see ).
4.1. Local Analysis
The fixed points of map (20) are obtained as nonnegative solutions of the algebraic system: Thus, the stationary equilibria are the following: , , , and .
The study of local stability of equilibrium solutions is based on the study of the Jacobian matrix: evaluated at the fixed point.
Proposition 4. The fixed points , , and are unstable.
Proof. We have
The eigenvalues associated with (24) are and . Then, is a source: The eigenvalues associated with (25) are , where the sign of the inequality follows from Assumption 2, and . Then, is a source or a saddle depending on the value of : The eigenvalues associated with (26) are , where the sign of the inequality follows from Assumption 2, and . Then, is a source or a saddle depending on the value of .
Proposition 5. The fixed point is locally stable if and only if , where is defined in proof. For , undergoes a supercritical flip bifurcation. No other local bifurcations may arise around .
Proof. The Jacobian matrix of map (20) evaluated at is the following: The trace and determinant associated with (27) are therefore given by
The local stability of Nash equilibrium is given by using Jury’s conditions; that is,
which for our system become
From condition we have that second inequality in (30) is always verified; with regard to the third condition, we note that is a second-degree polynomial with respect to . In addition, since (in fact, ) it follows that admits two positive roots and , if the discriminant associated with the equation results to be positive. We now verify the sign of the discriminant. By simple calculations, we have that
This expression is quadratic in and by a simple study we find that is always positive for any . Then, defines a concave function.
Now, since , and , then for any feasible value of .
To sum up, for any such that , where and their explicit expressions are given by
With regard to the first condition in (30), we note that the sign of is given by the sign of where is the first derivative of . Then, for a value , which corresponds to the minimum point of . By considering the sign of the coefficients and inequalities involved, we have the result.
4.2. Critical Curves
An important feature of map (20) is that it is a noninvertible endomorphism. In fact, for a given the rank- preimage (that is the backward iterate defined as ) may not exist or may be multivalued. In the specific case, if we want to compute in terms of in (20) we have to solve a fourth-degree algebraic system that may have four, two, or no solutions. In a natural way, we are led to subdivide the plane in regions , , and according to the number of such preimages (where the subscripts in indicate this number). A direct consequence of this fact is that if we let vary in the plane , , the number of the rank- preimages changes as the point crosses the boundary that separates these regions. Such boundaries are generally characterised by the existence of two coincident preimages. In this regard, by following  we introduce the definition of the critical curves. The critical curve of rank , denoted by is defined as the locus of points that have two (or more) coincident rank- preimages located on a set called . It is quite intuitive to interpret the set as the two-dimensional generalisation of the notion of critical value, local minimum or maximum, of a one-dimensional map, and as the generalization of the notion of critical point (local extremum point). Arcs of separate the regions of the plane characterised by a different number of real preimages.
Since is a continuously differentiable map, belongs to the locus of points where the Jacobian determinant of vanishes (i.e., the points where is not locally invertible). In our case, and is the rank- image of under ; that is, .
From direct computation, we have that if and only if It is easy to check that (35) is the equation of a hyperbola in the plane . Thus, is formed by two branches, denoted by and . This implies that also and subsequent iterations of the critical curves may be seen as the union of two different branches. In Figure 2 we can read and . Each branch of the critical curve separates the phase plane of into regions whose points have the same number of distinct rank- preimages. Specifically, separates region from region , and separates region from region . This allows us to study some global properties of the map and the evolution of basins of attraction and their qualitative changes (or bifurcations) as some parameters are varied.
4.3. Basins of Attraction
In this section, we describe the properties of the basins of attraction of map (20). We begin by studying the projection of the map on the Cartesian axes. For axes and we have, respectively, We note that these equations are both conjugated to the logistic map through the transformation with for (36) and with for (37).
It follows that the dynamics on axes and can be obtained from the well-known behaviour of the standard logistic map by a homeomorphism (see ). In particular, (a) if (resp., ), then we can deduce that bounded trajectories along (resp., ) are generated if the initial conditions are taken inside the segment line (resp., ; (b) from the computation of the eigenvalues of the cycles belonging to one of axes, we have that the direction transverse to the coordinate axes is always repelling; (c) initial conditions with or generate a divergent trajectory. From (a), (b), and (c) it follows that and and their preimages belong to the boundary of . In addition, under the conditions in (a), when the preimages belong to , we can show that is given by the region outside the quadrilateral with , , and where and so that the quadrilateral represents the region of trajectories that converge to a finite distance attractor that does not lie on the axes. Finally, from an economic point of view we note that, in order to preserve a positive value of the price inside the quadrilateral , we have to impose the condition ; that is, .
4.4. Global Analysis and Numerical Simulations
In this section we study the dynamic system for the following parameter values: , , , and and let vary. Figures 3(a) and 3(b) show the existence of a period doubling cascade that (starting from the flip bifurcation value ) generates cyclic attractors of higher period until a global attractor is born.
By increasing the value of , we have an important topological change in the structure of the basins of attraction. In particular, when we have a tangency between and the upper side of the quadrilateral (the grey region in the figures). For higher values of portions of the basin of attraction of the attractors on the axes enter region (at least one of the firms exit the market). After the bifurcation (tangency), one main lake lies inside . Hence, the lake has further preimages which form smaller lakes within the grey region. These ones lie inside the quadrilateral in the region complement to . When increases further, continues to move upwards, the portion grows up and then the holes become larger (see Figures 4(a) and 4(b)). At this stage it is really difficult to predict the long-term dynamics of the economic model, because slight changes in the initial conditions may lead to very different long-term outcomes (it is possible that only one firm produces or both firms produce but with erratic patterns). We note that the attractors located on the axes are not locally attracting. However, they attract a set of initial conditions with positive Lebesgue measure; that is, they are attractors à la Milnor (we recall that a closed invariant set is a Milnor attractor if its stable set has positive Lebesgue measure (see  for details)). Finally, for sufficiently large a contact bifurcation destroys the interior attractor.
5. Dynamics under Both Limited and Complete Information
In this section, we assume that firms adjust production period by period by using different mechanisms. In particular, the pessimistic firm has limited information and modifies production depending on the value of its own marginal profits, while the optimistic firm has complete information and static expectations (“best reply” dynamics). Given this type of heterogeneity and using (4) (pessimistic firm) and (7) (optimistic firm), the two-dimensional system that characterises the dynamics of the economy is now given by which is defined on .
Map (40) admits two fixed points: defined by (8) and . To guarantee that the unique Nash equilibrium of the game has economic meaning, we assume that the inequality continues to be fulfilled. This condition also implies that the equilibrium values of profits and price are positive.
5.1. Local Analysis
The study of local stability of equilibrium solutions is based on the study of the Jacobian matrix of the dynamic system. The Jacobian matrix of map computed in a generic point has the following form:
The following propositions hold.
Proposition 6. The fixed point is unstable.
Proposition 7. The fixed point is locally stable if and only if . For , undergoes a supercritical flip bifurcation. No other local bifurcations may arise around .
Proof. The Jacobian matrix of map (40) evaluated at is as follows: The trace and determinant associated with (43) are therefore given by The proof follows by using Jury’s conditions in (29), which are now given by
5.2. Global Dynamics
This section develops the global analysis of map (40). First of all, we note that if we relax the nonnegative constraints on variables and , map (40) results to be invertible (i.e., given there exists one and only one ). However, map (40) is noninvertible with nonnegative constraints on and . This can be ascertained by looking at Figure 5. From a global perspective, it is possible to identify regions on set corresponding to which at least one of the best replies of the two firms is zero (Figure 5). It is important to note that the introduction of heterogeneity in the mechanism of adjustment between the two firms leads to different behaviours of the map along the Cartesian axes with respect to the model developed in Section 4. In particular, the axis is a trapping subset. In fact, when the subsequent iteration results to be . In contrast, the axis does not have this property: the points that start on axis can lead to different fates (Figure 6(b)). Figure 5 shows that in the black region the subsequent iteration leads to a point , with , while in the yellow region the subsequent iterate leads to a point with and . In the red region the subsequent iteration leads to and finally the successive iterates will be located on the segment line .
We note, however, that by starting from points that lie in the yellow region it is possible that the subsequent iteration continues to lie in the yellow region or, alternatively, it leads to either the red region or black region and the dynamics will definitely end up on the point . In this case the side equilibrium is not locally attracting. However, it attracts a set of initial conditions with positive Lebesgue measure. From a mathematical point of view, these portions of the phase plane can be identified through the union of the preimages of any rank of the points in black and red regions.
Analogously with map , complex dynamics can be observed also for map . Figures 6(a) and 6(b) show the bifurcation diagram and the corresponding Lyapunov exponent () for by using the following parameter values: , , , and . When , the interior fixed point is stable. Starting from , there exists a cascade of period-doubling bifurcations until the occurrence of a chaotic attractor, which is shown in Figure 6(b) that depicts the basin of attraction and the corresponding chaotic attractor for .
In addition, in the model with heterogeneous adjustment mechanisms we note that coexistence of interior attractors can occur by slightly reducing the value of (this result is not usual in the nonlinear oligopoly literature with heterogeneous adjustment mechanisms [9, 10]). With this regard, Figure 7 shows that a cycle of period four (black points in the light-grey region) coexists with a cycle of period six (yellow points in the dark-grey region) when .
This paper developed a nonlinear Cournot duopoly to study the role of strategic uncertainty on the dynamics of the model economy. We characterised the local and global properties of a discrete two-dimensional map by considering that both firms have complete information on the market demand and adjust production over time depending on past behaviours (static expectations, “best reply” dynamics); both firms have incomplete information and production is adjusted over time by following a mechanism based on marginal profits; and one firm has incomplete information and production decisions depend on marginal profits, and the rival has complete information with static expectations (mixed case). In cases 2 and 3 we showed the existence of complex phenomena such as chaotic attractors and coexistence of attractors. The different behaviour of the model depending on whether the adjustment mechanisms are homogeneous or heterogeneous is interesting from an economic point of view, because of the different long-term dynamics that can be observed. In particular, in the cases in which (a) both firms have incomplete information and the adjustment mechanism of production is homogeneous and (b) the information set available to firms and the adjustment mechanisms of production are mixed, there exists the possibility that at least one firm decides to exit from the market (i.e., to do not produce). This result does not hold when both firms have complete information on the market demand and production is based on the behaviour of the previous period.
The authors gratefully acknowledge that this work has been performed within the activity of the PRIN-2009 project “Local interactions and global dynamics in economics and finance: models and tools,” MIUR (Ministry of Education), Italy, and PRIN-2009 project “Structural change and growth,” MIUR, Italy. Numerical simulations have benefited from algorithms that can be found in http://dysess.wikispaces.com/. The authors also thank an anonymous reviewer for valuable comments. The usual disclaimer applies.
- G. I. Bischi, C. Chiarella, M. Kopel, and F. Szidarovszky, Nonlinear Oligopolies. Stability and Bifurcations, Springer, Berlin, Germany, 2010.
- J. Dow and S. R. D. C. Werlang, “Nash equilibrium under Knightian uncertainty: breaking down backward induction,” Journal of Economic Theory, vol. 64, no. 2, pp. 305–324, 1994.
- J. Eichberger and D. Kelsey, “Uncertainty and strategic interaction in economics,” Discussion Papers 99-26, Department of Economics, University of Birmingham, 1999.
- J. Eichberger and D. Kelsey, “Non-additive beliefs and strategic equilibria,” Games and Economic Behavior, vol. 30, no. 2, pp. 183–215, 2000.
- M. Marinacci, “Ambiguous games,” Games and Economic Behavior, vol. 31, no. 2, pp. 191–219, 2000.
- A. Chateauneuf, J. Eichberger, and S. Grant, “Choice under uncertainty with the best and worst in mind: neo-additive capacities,” Journal of Economic Theory, vol. 137, no. 1, pp. 538–567, 2007.
- D. Schmeidler, “Subjective probability and expected utility without additivity,” Econometrica, vol. 57, no. 3, pp. 571–587, 1989.
- F. Fontini, “Cournot oligopoly under strategic uncertainty with optimistic and pessimistic firms,” Metroeconomica, vol. 56, no. 3, pp. 318–333, 2005.
- H. N. Agiza and A. A. Elsadany, “Chaotic dynamics in nonlinear duopoly game with heterogeneous players,” Applied Mathematics and Computation, vol. 149, no. 3, pp. 843–860, 2004.
- F. Tramontana, “Heterogeneous duopoly with isoelastic demand function,” Economic Modelling, vol. 27, no. 1, pp. 350–357, 2010.
- G.-I. Bischi, L. Stefanini, and L. Gardini, “Synchronization, intermittency and critical curves in a duopoly game,” Mathematics and Computers in Simulation, vol. 44, no. 6, pp. 559–585, 1998.
- G.-I. Bischi, M. Gallegati, and A. Naimzada, “Symmetry-breaking bifurcations and representative firm in dynamic duopoly games,” Annals of Operations Research, vol. 89, no. 1–4, pp. 253–272, 1999.
- L. Fanti, L. Gori, and M. Sodini, “Nonlinear dynamics in a Cournot duopoly with relative profit delegation,” Chaos, Solitons & Fractals, vol. 45, no. 12, pp. 1469–1478, 2012.
- L. Fanti, L. Gori, C. Mammana, and E. Michetti, “The dynamics of a Bertrand duopoly with differentiated products: synchronization, intermittency and global dynamics,” Chaos, Solitons & Fractals, vol. 52, pp. 73–86, 2013.
- G. I. Bischi, A. K. Naimzada, and L. Sbragia, “Oligopoly games with Local Monopolistic Approximation,” Journal of Economic Behavior and Organization, vol. 62, no. 3, pp. 371–388, 2007.
- T. Puu, “Chaos in duopoly pricing,” Chaos, Solitons & Fractals, vol. 1, no. 6, pp. 573–581, 1991.
- A. Dixit, “Comparative statics for oligopoly,” International Economic Review, vol. 27, no. 1, pp. 107–122, 1986.
- G. I. Bischi and A. Naimzada, “Global analysis of a dynamic duopoly game with bounded rationality,” in Advances in Dynamics Games and Application, J. A. Filar, V. Gaitsgory, and K. Mizukami, Eds., vol. 5, pp. 361–385, Birkhäuser, Boston, Mass, USA, 2000.
- J. S. Cánovas, T. Puu, and M. Ruíz, “The Cournot-Theocharis problem reconsidered,” Chaos, Solitons and Fractals, vol. 37, no. 4, pp. 1025–1039, 2008.
- C. Mira, L. Gardini, A. Barugola, and J.-C. Cathala, Chaotic Dynamics in Two-Dimensional Noninvertible Maps, World Scientific Publishing, Singapore, 1996.
- R. L. Devaney, An Introduction to Chaotic Dynamical Systems, Addison-Wesley, Redwood City, Calif, USA, 2nd edition, 1989.