Abstract

The comparison between two nonlinear duopoly models constructed based on symmetric utility function that is derived from Cobb–Douglas is investigated in this paper. The first model consists of two firms which update their outputs using gradient-based mechanism called bounded rationality. The second model contains a bounded rational firm that is competing with a firm whose outputs depend on a trade-off between market share maximization and profit maximization. For the two models, the fixed points are calculated and their conditions of stability are analyzed. The obtained results show that the second model is more stabilizing provided that the second firm adopts low weights of trade-offs. We show that the two models can be destabilized via flip bifurcation only. Furthermore, the noninvertibility of the two models that can give rise to several stable attractors is discussed.

1. Introduction

Literature has reported several works that have studied the complex dynamic characteristics of Cournot duopoly games. Such games have given rise to interesting results about the models describing them. These results have included investigations of the types of bifurcation responsible for instability of games’ equilibrium points and chaos routes. Some of those games have been formed based on popular utility function such as Cobb–Douglas [1], constant elasticity of substitution (CES) [2], and Singh and Vives utility function [3]. There are other utility production functions that have been adopted to model Cournot games like, for example, the isoelastic utility derived from Cobb–Douglas (see [4]). The current paper adopts a symmetric utility function that is derived from the well-known Cobb–Douglas one. Although the adopted symmetric utility in the current paper is simple, the two models formed based on it are complicated and provide interesting results about their dynamics as will be shown later.

In order to make the current paper self-contained, we report in the introduction some important works that have been cited in this direction of research. For instance, we cite from literature the reported works ([510]) which have discussed and investigated the routes of chaos in Cournot games whose models have been constructed based on naive expectations. Such works have analyzed the chaotic behaviors on those models via detecting the types of bifurcations appeared on them. In [11], a gradient-based mechanism called bounded rationality has been adopted to globally study the competition of Cournot game. In this study, a rich analysis has shown that the game’s equilibrium points can be destabilized due to coexistence of periodic cycles and chaotic behaviors. The adoption of bounded rationality approach has opened the gate to several studies on its influences on duopoly games and it has been applied to other games whose players may be homogeneous or heterogeneous. In [1215], the authors have applied the bounded rationality and other approaches such as Puu’s approach in order to analyze Cournot games and the cooperation that might be occurred between their players. These studies have confirmed what has been reported in literature on the destabilization routes of chaos that affected the stability of the equilibrium points in those games. Other useful works have been reported in [16, 17], in which Agiza et al. have investigated oligopoly games with both linear and nonlinear inverse demand function that have been used to study those games under limited rationality. In [18], chaos has been controlled in a duopoly game of master and slave players after the game’s dynamics have been analyzed and routes to chaos have been identified. The complex dynamic characteristics of both triopoly Cournot and Bertrand games whose players produce differentiated commodities have been investigated in [19]. Here, we cannot ignore to cite two important works of Bischi et al. ([20, 21]) which introduced an intensive analysis of bifurcation and chaos occurred in the dynamics of oligopoly games. In [22], a Cournot duopoly game whose players are homogeneous and their inverse demand functions are isoelastic has been introduced and its dynamic characteristics have been discussed. For more works in this direction of research, interesting readers are advised to see the literature [2330].

The current paper belongs to both homogeneous and heterogeneous duopoly games in which games’ players adopt similar and different decisional mechanisms. The main results of this paper concern with the comparison between two nonlinear duopoly games which are formed based on a symmetric utility function derived from Cobb–Douglas well-known production utility. The comparison includes investigation of routes where the equilibrium points can be destabilized, analyzing the global behaviors of the two models describing the games that include investigation of the structure of attractive basins. Our obtained results show that the equilibrium points of both models cannot preserve their stabilities due to flip bifurcation only. Moreover, they show that the model with heterogeneous players seems to be more stable if the second firm uses low trade-off weights between market share and profit maximization. The structure of basins for both models may be quite complicated as we deal with noninvertible maps.

After this introduction, we can summarize the parts of the paper as follows. Section 2 introduces the two maps describing the duopoly games in this paper. In Section 3, we analyze the properties of the first map and this includes discussing the stability conditions of its equilibrium points and proves it is a noninvertible map. As in Section 3, Section 4 discusses the second map and analyzes its complex characteristics. Some comparisons between the two maps are given in Section 5. In Section 6 ,we end the paper with some conclusions and future works.

2. The Game of Competition

We consider a game of two firms (or players) whose decision variables in the market are quantities. We take and as the quantities produced by each firm. We recall the Cobb–Douglas function of two variables appearing in many economic models as follows:

This function is popular as Cobb–Douglas (this function is a production function used to organize production possibilities for firms from different and various branches [31]. In the application to production, as a firm’s output is something measurable, the sum or makes sense as watershed between increasing and decreasing returns to scale. In our assumption in this paper, maximizing utility returns or that means the total revenue for both firms is equal and this is due to the assumed symmetry as we take .) function which should properly be called Wicksell function since the Swedish economist Knut Wicksell was the first to introduce this production function [32]. It represents the agent’s preferences. Both and are constants. It has the following properties:(i)It is a monotonic function under the transformation . This means that and and hence the marginal utility of each quantity is positive.(ii)It is a concave utility function. This means that and and hence the marginal utility of each quantity is decreasing.(iii)It is strongly additive. This means that and hence the marginal of the quantity is independent of the quantity . So, the quantities cannot represent substitutes or complements.(iv)It is a homothetic function. This means that if , , and , then .

We assume the symmetric case, , with the budget constraint where is the price of . Now, we have the following maximization problem:Solving (2) gives the following prices:

At the economic market, each competitor wants to maximize its profit given bywhere is the cost function of each quantity. It is a quadratic function and is often met in applications. For example, in the modeling of renewable resources exploitation, such as fisheries [33], in [34], it has been considered in order to investigate the role of convexity. So, the marginal cost is not constant, , and is a constant parameter. The marginal profits take the following form:

Here, we discuss two situations. The first situation assumes that both firms use a gradient-based mechanism such as the bounded rationality defined in the literature [4, 3537]. Using this mechanism, we get the following discrete dynamical system:where is a speed of adjustment parameter. The second situation is a heterogeneous situation. We assume that the second firm shares the market with certain profit. It trades off between market share maximization and profit maximization. The market share means it seeks the optimum output by putting and then we get the optimum, , while the profit maximization gives . Since it is traded off between those, we get the following:where . When the second firm shares market maximization only, we take , while means it seeks profit maximization only. Now, the second firm updates its output according to the following:where is a constant weight. So, we have another model with heterogeneous players given by

It is obvious that if , then the second firm in (9) is naive, while at means it updates its production based on market share maximization. Economically, profit maximization only is mainly a short-term goal and is primarily restricted to the accounting analysis of the financial situation. Furthermore, it is primarily concerned as to how the company will survive and grow in the existing competitive business environment [38]. On the other hand, market share can allow firms to improve profitability either by lowering prices, using advertising, or introducing new or different products. It can also grow the size of its market share by appealing to other audiences or demographics. In our paper, we introduce a weighted average decision-making mechanism between those two mechanisms discussed above. This includes the main contribution in this paper that is to compare between the maps (6) and (9) in order to see which model is more stable and gives large region of stability [39].

3. Properties of Map (6)

3.1. Analytical Results

The map (6) can make negative or unbounded trajectories if the initial condition is chosen far from . It is obvious that if , negative values for can be obtained.

Proposition 1. The map (6) has two fixed points that are and .

Proof. Setting and in (6) and solving algebraically, we get the two fixed points.

Proposition 2. The fixed point is an unstable point.

Proof: The Jacobian matrix of map (6) at becomeswhich contains indefinite form and then the point is unstable and the proof is completed.

Proposition 3. The fixed point is locally asymptotically stable provided that

Proof. At , the Jacobian matrix becomeswhose trace and determinant are given bySubstituting (13) in Jury conditions [40] givesSince and are positive speeds of adjustment parameters then is always satisfied. and giveCombining the above inequalities completes the proof.

Proposition 4. The fixed point can be destabilized due to flip bifurcation only.

Proof. The Jacobian (12) has the following eigenvalues:It is clear that is always positive, and therefore, we have only two real eigenvalues. So the point becomes unstable because of flip bifurcation.

3.2. Numerical Results

The following numerical simulation experiments validate the above theoretical results. We start our numerical experiments by assuming the following parameters’ set , . At this set, and the Jacobian becomeswhere we assume and . The eigenvalues then become with and then is locally asymptotically stable. Keeping the costing shift parameters and as previously, we now study the influence of or on the stability of . Figures 1(a) and 1(b) show that is locally asymptotically stable for all values of and ; until we get the cycle of period 2, the point O2 becomes unstable. Indeed, the simulation shows some interesting behavior of the dynamic of map (6) around the point . To present that dynamics, we have to study the effects of or separately. Assuming the following parameters set, , and . As increases to 6.54, a period-2 cycle (represented by squares in Figure 1(c)) is born around the fixed point . It is plotted with its basins of attraction represented by the light green colors in Figure 1(c) while the grey color refers to divergent (or nonconvergent) points in the phase plane. This cycle turns into a cycle of period four as increases to 7.14 with quite complicated basins of attraction as displayed in Figure 1(d). Other basins of attraction for the cycle of period four are formed in Figure 1(e) at where nonconvergent points appeared in white color within the basins. The basins of attraction become more complicated as the period-8 starting to appear at as depicted in Figure 1(f). Increasing further gives rise to four disconnected pieces and two chaotic attractors and then one chaotic attractor. We give in Figures 1(g) and 1(h) examples of two and one chaotic attractors and their basins at and , respectively. The white color in the basins denotes the nonconvergent points in the phase plane. The same discussion can be performed for the other parameter under the same values of costs. This makes us to display the 2D-bifurcation diagram in the on which interesting readers can detect more about the dynamics of map (6) around . Figure 2(a) presents the 2D-bifurcation in the with examples of attractive basins of period-7 (Figure 2(b)), period-9 (Figure 2(c)), and period-10 (Figure 2(d)). Furthermore, the 2D-bifurcation diagram shows that the dynamics of map (6) do not permit the appearance of period-3 and period-5.

Now, we study the case when by assuming the set of parameters’ values, . In this case, we do not want to repeat the above discussions and only we give some simulations represented by Figures 3(a) to 3(f) at this set in order to see the difference in dynamics in this case and the previous one. All these results so far show a peculiar shape of the basins of attraction and require from us to investigate other characteristics of map (6). From the above discussion and the numerical experiments, we highlight that the assumptions and affect the stability region obtained by the second condition of (14) for the interior equilibrium point in the . We have seen that at , we get and while in the case , we have and for the same condition.

3.3. Critical Curves and Phase Plane Zones

The structure of the basins of attraction of map (6) at any that may represent or periodic cycle or any chaotic attractor is constructed by some boundaries. These boundaries are calculated as follows.

Proposition 5. The origin point has four real preimages.

Proof. Setting and in (6) and solving the system algebraically, we getwhere is obtained by solving algebraically the following system:

Proposition 6. Let and be the two line segments of the coordinates and . Also, let and be their preimages. Then, the boundary of the basins of attraction of any is denoted by and is given bywhere refers to the set of all preimages of rank-n.

These boundaries are displayed in Figure 4. Furthermore, the figure shows the critical curve that is calculated by vanishing the determinant of the Jacobian of (6). It is represented by the following curve:

It is represented by red color in Figure 4 as it contains two parts, . It is hard to plot or to get an analytical form for in order to identify the zones in the phase plane; however, it is simple to see that has two real preimages of rank-1 and hence it belongs to . Furthermore, both and have no preimages and they belong to . With the preimages of the origin, we detect that the phase plane of map (6) is divided into three zones, , and . Therefore, map (6) is a noninvertible map.

4. Properties of Map (9)

4.1. Analytical Results

Proposition 7. The map (9) has only two fixed points that are and where

Proof. Setting and in (9) and solving the system algebraically completes the proof.

Proposition 8. The fixed point is unstable.

Proof. The proof is the same as Proposition 2.

Proposition 9. The fixed point is locally asymptotically stable provided that where

Proof. The Jacobian at becomeswhereand the trace and determinant take the following form:and then,It is clear that and since and is always positive. While is attained, if . The proof is completed.
We should highlight here that if , then can be destabilized due to flip bifurcation only.

4.2. Numerical Results

We assume the following , with . At this set, and the Jacobian becomes, at taking and ,and the eigenvalues are with and then is locally asymptotically stable. Increasing above the value gives rise to a period-doubling bifurcation, as shown in Figure 5(a). It is clear that the fixed point is stable for all the values of until reaches the value where the first period-2 cycle is born. As discussed in map (6), we give here some complex behaviors of map (9) at different values of the bifurcation parameter . For instance, at , a cycle of period-2 is emerged and is plotted with its basins of attraction in Figure 5(b). At , a period-4 cycle is coexisted with its attractive basins in Figure 5(c). It is clear that they are quite complicated basins. The basins become more complicated when the period-8 cycle appears as given in Figure 5(d) at . After that higher periodic cycles exist with quite complicated basins and then the map enters the chaos region at and where two chaotic attractors and then one chaotic attractor are emerged and displayed in Figures 5(e) and 5(f), respectively.

The numerical simulation shows that increasing which means the second firm gives more interesting to the market share maximization than the profit maximization reduces the stability interval of . The same observation when the firm gives more interesting to the traded-offs between the market share maximization and profit maximization is obtained. It is noticed that low values of the traded-off parameter keep the stability interval of larger. Furthermore, when we take values of and in the form , the dynamics of map (9) do not give any new complex dynamics and hence the dynamic is still governed by the change that may be happened on both parameters and .

4.3. Critical Curves and Phase Plane Zones

which has no real solutions as the expression under the square root is negative. Therefore, belongs to .

Proposition 10. The origin point has only one real preimage that is .

Proof. Setting and in (9), we getThe first equation in (29) gives or . Substituting in the second equation of (29) gives , then we get . Now, we have the following case:Solving (30) gives

Proposition 11. The fixed point has no real preimages.

Proposition 12. Any point in the form has two real preimages.

Proposition 13. Any point in the form has no real preimages.

Now, for map (9) is given by the following curve:

It is plotted in Figure 6 with red color at the parameters’ values, , and . Simulation shows that any point belonging to has two real preimages. However, it is hard to plot but we can conclude that the phase plane of (9) may be divided into only two zones that are and . Therefore, map (9) is a noninvertible map.

5. Comparison between the Maps and

The above analysis and discussion about the two maps show that the model described by map (9) is quite complicated than those of map (6); however, map (9) gives a large stability region for the fixed point under the same set of parameters. This is clear when comparing the 1D bifurcation diagram for the two maps. This means that the firm adopting the trade-offs between the market share maximization and the profit maximization will be more stable in the market competition than those using a gradient-based mechanism such as the bounded rationality. It is also observed from the numerical experiments that when the second firm updates its outputs based on the average of market share and profit maximization, the stability region of the fixed point reduces and becomes almost like those provided by map (6) under the same set of parameters’ values. Figure 7(a) shows the bifurcation diagram at , and . Increasing further to 0.7 and keeping the other parameters’ values fixed, the region of stability region of the fixed point reduces, as shown in Figure 7(b). Furthermore, Figures 7(c) and 7(d) show the influences of the parameters and on the dynamics of map (9). From an economic perspective, one can see that low values of the parameter means the second firm pays more attention for the profit maximization (just as the first firm decision-making mechanism) than market share and hence the interior equilibrium point becomes stable for any value of in the interval . As increases further, the stability interval of the equilibrium point with respect to is reduced, as shown in Figure 7(c). Similar discussion is for the influences of the parameter given in Figure 7(d).

6. Conclusion

In this manuscript, we have studied two nonlinear duopoly models whose players adopted gradient-based mechanism and trade-off mechanism. We have shown that the fixed points in both models become unstable due to flip bifurcation only. The duopoly in the first model has been compared with the duopoly of the second model whose second firm uses trade-off between market share maximization and profit maximization. Our contributions showed that the second model was more stabilizing than the first model. It has been investigated that the trade-off mechanism can give better stability region for the fixed point. We found that when the second firm in the second model adopts the average between market share and profit maximization, the region of stability was reduced. Furthermore, using higher weights of the trade-off makes the stability region to shrink and hence the second firm should adopt low weights to become more stable in the market. We have also investigated that both maps are noninvertible, and for this reason, the basins of attraction of any dynamic behavior of the maps were quite complicated. In future study, we intend to apply the trade-off mechanism to oligopoly models with more than two firms.

Data Availability

The datasets generated during the current study are available from the corresponding author on reasonable request.

Conflicts of Interest

The author declares no conflicts of interest.

Acknowledgments

This research was supported by King Saud University, Riyadh, Saudi Arabia, Project number (RSP-2020/167).