Next Article in Journal
Transient Behavior of the MAP/M/1/N Queuing System
Next Article in Special Issue
Deep Learning Models for Predicting Monthly TAIEX to Support Making Decisions in Index Futures Trading
Previous Article in Journal
Sustainability, Big Data and Mathematical Techniques: A Bibliometric Review
Previous Article in Special Issue
Revisiting the Valuable Roles of Global Financial Assets for International Stock Markets: Quantile Coherence and Causality-in-Quantiles Approaches
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Complex Investigations of a Piecewise-Smooth Remanufacturing Bertrand Duopoly Game

Department of Statistics and Operations Research, College of Science, King Saud University, P.O. Box 2455, Riyadh 11451, Saudi Arabia
Mathematics 2021, 9(20), 2558; https://doi.org/10.3390/math9202558
Submission received: 5 September 2021 / Revised: 1 October 2021 / Accepted: 5 October 2021 / Published: 13 October 2021
(This article belongs to the Special Issue Mathematics and Financial Economics)

Abstract

:
This paper considers a Bertrand competition between two firms whose decision variables are derived from a quadratic utility function. The first firm produces new products with their own prices while the second firm re-manufactures returned products and sells them in the market at prices that may be less than or equal to the price of the first firm. Dynamically, this competition is constructed on which boundedly rational firms apply a gradient adjustment mechanism to update their prices in each period. According to this mechanism and the nature of the competition, a two-dimensional piecewise smooth discrete dynamic map was constructed in order to study the complex dynamic characteristics of the game. The phase plane of the map was divided into two different regions, separated by border curve. The equilibrium points of the map, in each region on where they are defined, were calculated, and their stability conditions were investigated. Furthermore, we conducted a global analysis to investigate the complex structure of the map, such as closed invariant curves, periodic cycles, and chaotic attractors and their basins, which cause qualitative changes as some parameters are allowed to vary.

1. Introduction

In the economy, there are different kinds of competitions that attract both scientific researchers and economists (who study the complex characteristics of these competitions). Studies show that some of these economic models may be chaotic. Here, we must cite the individual who implemented the first trial of such studies, essentially opening the way for more studies about such models. Augustin Cournot introduced the first mathematical model of a duopoly game [1]. A duopoly game is a special ’case’ of the oligopoly game, in which only two firms dominate the market. It is more complex than a monopoly, where only one firm dominates the market. Duopoly games assume that firms sell homogeneous goods in the same market, but with independent prices. Consequently, the optimal decisions of these firms are affected by their opponents. Competition in such games is divided into static or dynamic; in this manuscript, we focus on the dynamic. The Cournot and Bertrand models are popular models that describe dynamic competition [2]. In Cournot models, the two firms produce homogeneous commodities; the decision variables here are the quantities they produce. On the other hand, when the decision variables are prices instead of quantities, the model is known as a Bertrand. The Stackelberg game is also interesting [3,4]—it involves two firms, one of which is more powerful than the other. In such a game, the weaker firm does not take action until the stronger firm takes action. The firms are known as the leader and the follower.
Studying the dynamic features of the above-mentioned models shows whether we can make predictions about these models in the long run. Studying the stability of the equilibrium points of such models occurs after modeling; this requires following certain mechanisms, such as bounded rationality and Puu’s approach. The latter is not applied much in the literature because it has some disadvantages when studying the stability of the equilibrium points [5,6,7]. The bounded rationality mechanism has been intensively applied in many economic games investigated in the literature. More information about this mechanism and its application exists [8,9,10]. Studying the stability of the equilibrium reveals that such models open the gate toward discovering other complex dynamic characteristics of systems modeling such a game. Characteristics include the coexistence of different types of bifurcations that make the equilibrium points lose their stability and behave chaotically. There are several studies in the literature [11,12,13,14,15] that have reported that those systems lose their stability due to flip and Neimark–Sacker bifurcations; then routes to chaos coexist.
Remanufacturing in the economy is important and requires more analytical and numerical studies to show its benefits. Remanufacturing is the process of disassembling commodities, which are then repaired and collected, and sent back to the line of production to be sold as new products. The firm, known as a third-party company, handles this part of production. Examples of such companies include Xerox, as well as Mercedes-Benz and Ford for car industries [13]. Studies show that most of the current remanufacturing competitions highlight the decisions taken by all competitors involved in the economy. Only a few studies have focused on the complex dynamic features behind this kind of competition. For example, in [13], a remanufacturing Cournot duopoly game was introduced. It shows that the equilibrium point of the game loses its stability due to the coexistence of two types of bifurcations—period doubling and Neimark–Sacker. Furthermore, in [16], a remanufacturing duopoly game was introduced and its characteristics were investigated. The current paper extends the work introduced in [13]. Here, we study the case of Bertrand adopting the inverse demand functions used in [13]. We assume that the preferences adopted by the consumers follow a quadratic utility function that is concave, by which the inverse demand functions of both firms are derived. However the prices derived from this utility are linear, but the analysis given in this paper is rich in extra dynamic features. Our analysis, either analytic or numeric, differs from those given in [13,16]. In the current paper, we introduce a piecewise smooth Bertrand duopoly game, which is entirely different than smooth Bertrand games in the literature, and gives rise to our motivation. The piecewise Bertrand game is more complicated than classic ones and the bifurcation taking place in this game belongs to border-collision bifurcation. In the literature, several works on tackling piecewise smooth maps were discussed. In [17], a two-dimensional discontinuous piecewise linear map was considered. This map described the dynamics of a discrete-time version of a continuous-time fashion cycle model. The bifurcation discussed in this map belongs to the class of border-collision bifurcation. A three-region footloose-entrepreneur economic geography model was introduced in [18]. The map described this geography model as a piecewise map, defined in three regions. Its local and global dynamics and the corresponding bifurcation were discussed in [18]. In [19], a four-region economic geography model was introduced and analyzed as a piecewise map that was defined in four regions, in which regions are differentiated in size and geographical position. The effects of the constraints were discussed in an adaptive segregation model in [20]. In that study, the model was described by a discrete two-dimensional piecewise smooth dynamic map to investigate the dynamics of the entry and exit of two populations into the model described by such a map. The asymmetric case of the map given in [20] was studied in [21]. In [22], the border collision bifurcations of a piecewise adaptive oligopoly model whose demand function was isoelastic was investigated. Discriminatory individual choices may lead to two groups of people of opposite kinds as a result of segregation. People may be separated for different reasons, such as age, sex, language, income, nationality, skin color, etc. This segregation has led to investigations of such piecewise maps. This piecewise characteristic, which is one of the peculiarities of such models, gives rise to a particular kind of bifurcation due to the change of map definition. This kind of bifurcation is related to closed invariant sets (such as attractors, frontiers, manifolds), having contact with the border region on where the map is defined. Such maps are more general and complicated than the traditional smooth maps and only few studies on those maps have been reported in the literature and its applications to economic models are still undeveloped areas.
In brief, this paper is organized as follows. In Section 2, the piecewise smooth Bertrand map describing the proposed game is formulated. In Section 3, the stability conditions for the equilibrium points in each region where the map is defined are discussed. In Section 4, local and global analyses based on numerical experiments are investigated. Finally, we summarize our obtained results in this paper in the conclusion section.

2. The Model

In an economic market, individuals have their own preferences, which are used to measure their satisfaction about certain goods. Those preferences are translated into decisions of consumers, which are represented according to the demand theory by a utility function. Under a certain budget condition, this utility must be maximized in order to achieve a consumer’s perspective [23]. In this paper, we assume the following utility function,
U ( q 1 , q 2 ) = q 1 + θ q 2 1 2 ( q 1 2 + 2 θ q 1 q 2 + θ q 2 2 ) ; q 1 , q 2 > 0
Which was defined by Singh and Vives [24]. It is easy to see that this function is strictly concave for any values of the decision variables q 1 and q 2 and the fraction parameter θ where 0 < θ < 1 . The meaning of this fraction is defined later in this section. Now, assume that we have two competing firms. The first firm is called the manufacturer and it produces new goods that are denoted by the variable q 1 , while the second firm receives returned goods and repairs them, and then sells them again in the market. It reproduces the quantity q 2 . Indeed, we restrict the decision variables space to S = ( q 1 , q 2 ) : q 1 , q 2 > 0 . As the first firm is responsible for providing the market with the q 1 at period time t, the second firm should provide the market with returned products at time t + 1 . Therefore, the consumers will have to differentiate between those products and are willing to pay a fraction θ in order to buy new remanufactured products. Moreover, the second firm’s marginal cost must be less than this fraction θ otherwise it loses the market. So, using a budget constraint i = 1 2 p i q i = 1 with p i = U q i , i = 1 , 2 and with (1), we obtain the inverse demand functions as follows,
p 1 = 1 q 1 θ q 2 , p 2 = θ 1 q 1 q 2
where p i , i = 1 , 2 denotes the price for q i , i = 1 , 2 . θ = 1 indicates customers’ bias to pay same price for both new and remanufactured goods, which may not be accepted in the economy. θ = 0 means remanufactured goods are not approved by customers. Hence, this variety of preferences gives restrictions of the parameter θ ( 0 , 1 ) . Now, the demand functions of (2) are given by,
q 1 = 1 1 1 θ p 1 + 1 1 θ p 2 , q 2 = 1 1 θ p 1 + 1 θ ( 1 θ ) p 2
Using linear cost functions, C i ( q i ) = c i q i , i = 1 , 2 where c i , i = 1 , 2 are the marginal costs. This gives the following profits,
π 1 = p 1 c 1 1 1 1 θ p 1 + 1 1 θ p 2 , π 2 = p 2 c 2 1 1 θ p 1 + 1 θ ( 1 θ ) p 2
Using (4), the marginal profits become:
π 1 p 1 = 1 + c 1 1 θ 2 1 θ p 1 + 1 1 θ p 2 , π 2 p 2 = c 2 θ ( 1 θ ) 2 θ ( 1 θ ) p 2 + 1 1 θ p 1
By setting π i p i = 0 , i = 1 , 2 the equilibrium point E = ( p ¯ 1 , p ¯ 2 ) for the firms is unique and takes the following form:
p ¯ 1 = 2 ( 1 + c 1 θ ) + c 2 4 θ , p ¯ 2 = θ ( 1 + c 1 θ ) + 2 c 2 4 θ
It depends on the parameters c 1 , c 2 and θ . Let us first show the influence of the parameter θ on the equilibrium point. It is easy to get,
p ¯ 1 θ = 2 c 1 + c 2 6 4 θ 2 , p ¯ 2 θ = θ 2 8 θ + 4 + 2 ( 2 c 1 + c 2 ) 4 θ 2
This means that p ¯ 1 θ < 0 under the condition 2 c 1 + c 2 < 6 and p ¯ 2 θ > 0 is always satisfied for θ 2 8 θ + 4 + 2 ( 2 c 1 + c 2 ) > 0 and 0 < θ < 1 . This means that by increasing the parameter θ , the price of the remanufactured goods will be high and the price of the new goods will decrease. However, in the economy, consumers prefer to buy new goods instead of buying remanufactured goods; therefore, it may be preferable to keep the fraction θ low for both firms. Now, we assume a repeated discrete dynamic games between the two firms according to their prices (the so-called Bertrand game). To do that, we recall the bounded rationality mechanism [25,26,27]. This mechanism depends on the change occurring in the marginal profit, Φ i ( p 1 , p 2 ) = π i p i , i = 1 , 2 , whether it is increased or decreased. It is expressed as follows,
p i ( t + 1 ) = p i ( t ) + k i ( p i ) Φ i ( p 1 , p 2 ) , i = 1 , 2
where, k i ( p i ) refers to an adjustment parameter. This parameter is important in such games and has influences on the stability of the equilibrium point of the game. Using (5) and (8) and assuming that α i ( q i ) = k i p i , i = 1 , 2 the following discrete dynamic system is obtained.
p 1 , t + 1 = p 1 , t + k 1 p 1 , t 1 + c 1 1 θ 2 1 θ p 1 , t + 1 1 θ p 2 , t , p 2 , t + 1 = p 2 , t + k 2 p 2 , t c 2 θ ( 1 θ ) 2 θ ( 1 θ ) p 2 , t + 1 1 θ p 1 , t
where k i , i = 1 , 2 is positive constant. Here, we should mention that the price of the second firm at time t + 1 should be less than or equal those of the first firm at time t, i.e., p 2 ( t + 1 ) p 1 ( t ) . Therefore, the system defined by (9) should be modified to: p 2 ( t + 1 ) = min Ω ( t ) , p 1 ( t ) where Ω = p 2 ( t ) + k 2 p 2 c 2 θ ( 1 θ ) 2 θ ( 1 θ ) p 2 + 1 1 θ p 1 then the system (9) is rewritten to,
p 1 , t + 1 = p 1 , t + k 1 p 1 , t 1 + c 1 1 θ 2 1 θ p 1 , t + 1 1 θ p 2 , t , p 2 , t + 1 = min Ω ( t ) , p 1 ( t )
Therefore, the map defined in (10) is considered as a piecewise map. In order to study the stability of this map, one should consider the second equation in the above map. It is governed by the function Ω . This function is responsible for converting the map (9) to a piecewise one. This means that if p 1 , p 2 Ω and Ω < p 1 then we have only one part from the map (10) to be applied, otherwise if Ω > p 1 we have the other part of the map to be used. Here comes the case when Ω = p 1 , which is called the border curve, in which the map (10) is continuous but not differentiable. This border curve divides the phase plane into different regions where the map is defined by different expressions. It has the following form,
p 1 = p 2 1 + k 2 θ ( 1 θ ) c 2 2 p 2 1 k 2 1 θ p 2 : = g ( p 2 )
that makes the map (10) becomes,
p 1 , t + 1 = p 1 , t + k 1 p 1 , t 1 + c 1 1 θ 2 1 θ p 1 , t + 1 1 θ p 2 , t , p 2 , t + 1 = p 2 , t + k 2 p 2 , t c 2 θ ( 1 θ ) 2 θ ( 1 θ ) p 2 , t + 1 1 θ p 1 , t if p 1 , t g ( p 2 ) p 1 , t if p 1 , t g ( p 2 )
or it is equivalent to,
p 1 , t + 1 , p 2 , t + 1 = p 1 , t + k 1 p 1 , t 1 + c 1 1 θ 2 1 θ p 1 , t + 1 1 θ p 2 , t , p 2 , t + k 2 p 2 , t c 2 θ ( 1 θ ) 2 θ ( 1 θ ) p 2 , t + 1 1 θ p 1 , t if p 1 , t g ( p 2 ) ( region 1 ) p 1 , t + 1 , p 2 , t + 1 = p 1 , t + k 1 p 1 , t 1 + c 1 1 θ 2 1 θ p 1 , t + 1 1 θ p 2 , t , p 1 , t if p 1 , t g ( p 2 ) ( region 2 )
It is important to note that the stability of the above map should be studied according to the border curve. This means that the map may have fixed points that are located in different regions of the phase plane. The destabilization of those fixed points means that their dynamics are usually never related to one region and, hence, the dynamic of the map must be studied as a piecewise map.

3. The Stability

It is clear that the map (13) is defined in two different regions, so there may be different fixed points in each region. We denote the fixed points in these regions by E r 1 and E r 2 where E r 1 = 2 ( 1 + c 1 θ ) + c 2 4 θ , θ ( 1 + c 1 θ ) + 2 c 2 4 θ and E r 2 = 1 + c 1 θ , 1 + c 1 θ . The Jacobians at those points are,
J r 1 = 1 2 k 1 ( 1 + c 1 ) [ 2 ( c 2 + θ ) θ 2 ] 4 ( 1 + c 1 ) ( 1 + c 2 ) θ 2 θ k 1 [ 2 ( c 2 + θ ) θ 2 ] 4 ( 1 + c 1 ) ( 1 + c 2 ) θ 2 θ 2 k 2 ( 1 + 2 c 1 ) 4 ( 1 + c 1 ) ( 1 + c 2 ) θ 2 1 2 θ k 2 ( 1 + 2 c 1 ) ( c 2 + θ ) 4 ( 1 + c 1 ) ( 1 + c 2 ) θ 2 , J r 2 = 1 2 k 1 ( 1 + c 1 ) θ + 2 ( 1 + c 1 ) θ k 1 θ + 2 ( 1 + c 1 ) 1 0
Studying the stability of those fixed points depends on the eigenvalues that are calculated by λ 1 , 2 = 1 2 τ ± τ 2 4 δ . Where τ and δ refer to the trace and determinant respectively. They are depending on the region on where the fixed points lie in. There may be τ r 1 or τ r 2 if the fixed points lie in region 1 or region 2. The same observation is for δ ( δ r 1 or δ r 2 ). They are given by,
τ r 1 = 2 1 k 1 c 2 + 2 ( 1 + c 1 θ ) 4 θ 1 θ k 2 θ 2 c 2 + θ ( 1 + c 1 θ ) 4 θ 1 θ , τ r 2 = 1 2 ( 1 + c 1 θ ) 1 θ k 1 , δ r 1 = 1 2 k 1 c 2 + 2 ( 1 + c 1 θ ) 4 θ 1 θ 2 k 2 θ 2 c 2 + θ ( 1 + c 1 θ ) 4 θ 1 θ k 1 k 2 1 4 θ 2 c 2 + 2 ( 1 + c 1 θ ) 4 θ 1 θ 2 c 2 + θ ( 1 + c 1 θ ) 4 θ 1 θ , δ r 2 = ( 1 + c 1 θ ) 1 θ k 1
The eigenvalues in both regions are,
λ 1 r 1 , 2 r 1 = 1 k 1 c 2 + 2 ( 1 + c 1 θ ) 4 θ 1 θ k 2 θ 2 c 2 + θ ( 1 + c 1 θ ) 4 θ 1 θ ± ± 1 θ θ 2 k 1 2 c 2 + 2 ( 1 + c 1 θ ) 4 θ 1 θ 2 θ 1 θ k 1 k 2 c 2 + 2 ( 1 + c 1 θ ) 4 θ 1 θ 2 c 2 + θ ( 1 + c 1 θ ) 4 θ 1 θ + k 2 2 2 c 2 + θ ( 1 + c 1 θ ) 4 θ 1 θ 2 , λ 1 r 2 , 2 r 2 = 1 2 k 1 1 + c 1 θ 1 θ ± 1 2 1 + 4 k 1 2 1 + c 1 θ 1 θ 2
Or it can be rewritten in the form,
λ 1 r 1 , 2 r 1 = 1 2 τ r 1 ± τ r 1 2 4 δ r 1 λ 1 r 2 , 2 r 2 = 1 2 τ r 2 ± τ r 2 2 4 δ r 2
It should be noted that the map (13) is continuous and differentiable in both regions separately. However, on the border curve Ω = p 1 , the map is continuous, but not differentiable. Furthermore, there may be no fixed point in one region of the phase plane. This means that the fixed point E r 1 (or E r 2 ) may turn out to be in region 2 (or region 1). This also means that the character of a virtual fixed point exists. In addition, the stability of those fixed points will be governed by the following triangles,
S r 1 = ( τ r 1 , δ r 1 ) : 1 + τ r 1 + δ r 1 > 0 , 1 τ r 1 + δ r 1 > 0 , 1 δ r 1 > 0 S r 2 = ( τ r 2 , δ r 2 ) : 1 + τ r 2 + δ r 2 > 0 , 1 τ r 2 + δ r 2 > 0 , 1 δ r 2 > 0
where,
1 + τ r 1 + δ r 1 = c 2 + 2 1 + c 1 θ 2 c 2 + θ 1 + c 1 θ θ 4 θ 1 θ 2 k 1 k 2 , 1 τ r 1 + δ r 1 = 4 4 c 2 + 2 ( 1 + c 1 θ ) 4 θ 1 θ k 1 4 θ 2 c 2 + θ ( 1 + c 1 θ ) 4 θ 1 θ k 2 1 4 θ c 2 + 2 ( 1 + c 1 θ ) 4 θ 1 θ 2 c 2 + θ ( 1 + c 1 θ ) 4 θ 1 θ k 1 k 2 , 1 δ r 1 = 2 k 1 c 2 + 2 ( 1 + c 1 θ ) 4 θ 1 θ + 4 θ 2 c 2 + θ ( 1 + c 1 θ ) 4 θ 1 θ + 1 8 θ 2 c 2 + 2 ( 1 + c 1 θ ) 4 θ 1 θ 2 c 2 + θ ( 1 + c 1 θ ) 4 θ 1 θ k 1 k 2 , 1 + τ r 2 + δ r 2 = ( 1 + c 1 θ ) 1 θ k 1 , 1 τ r 2 + δ r 2 = 2 3 ( 1 + c 1 θ ) 1 θ k 1 , 1 δ r 2 = 1 + ( 1 + c 1 θ ) 1 θ k 1
Now, the next section carries out some numerical simulation in order to investigate the routes and the types of bifurcation by which the fixed points of the above system lose their stability.

4. Numerical Simulation

By looking at the map (13) one can see that it contains five important parameters, k 1 , k 2 , c 1 , c 2 and θ . We start our numerical analysis by investigating the influences of k 1 and k 2 on the stability of the map. Figure 1a,b present 2D-bifurcation diagrams in the k 1 , k 2 -plane at two different sets of parameters’ values. These sets are ( θ , c 1 , c 2 ) = 0.9721801 , 0.41 , 0.45 and θ , c 1 , c 2 = 0.74 , 0.4 , 0.8316 for the same datum p 0 , 1 , p 0 , 2 = 0.12 , 0.11 . Figure 1a shows different colored regions on where different periodic cycles, chaotic attractors, and closed invariant curves coexist due to both Neimark–Sacker and flip bifurcation. Figure 1b presents closed invariant sets due to Neimark–Sacker bifurcation followed by flip bifurcation. Therefore, we study, in the next section, the global influence of the parameters k 1 and k 2 separately in the dynamic of map (13). Now, we study the local impact of these parameters by assuming ( θ , c 1 , c 2 , k 2 ) = 0.9721801 , 0.41 , 0.45 , 0.006 . Figure 1c presents a situation of stable equilibrium point for any value of the parameter k 1 until this parameter reaches the value 0.061 on where there are two routes to two different types of bifurcations, Neimark–Sacker comes first then followed by flip bifurcation. The same observation is presented in Figure 1d for the variable q 2 . At the parameters set ( θ , c 1 , c 2 , k 1 ) = 0.9721801 , 0.41 , 0.45 , 0.006 , Figure 1e,f show the bifurcation diagrams for the map’s variables when varying the parameter k 2 . Another mixed scenario arises when we assume the following set of parameter values, θ , c 1 , c 2 = 0.74 , 0.4 , 0.8316 . At this set, the equilibrium point becomes stable when varying either k 1 or k 2 until closed invariant curves are born because of Neimark–Sacker bifurcation. Afterward, higher periodic cycles and chaotic attractors coexist due to flip bifurcation, as shown in Figure 2a,b. This discussion shows that the speeds of the adjustment parameters k 1 and k 2 have a great impact on the map’s equilibrium points. There is more evidence by which we confirm that the map enters chaotic regions. The evidence is confirmed by plotting the largest Lyapunov exponent (LLE). It is known that when LLE takes positive values, it means that the map’s dynamic behavior enters chaotic regions and, hence, closed invariant curves or higher period cycles, and then routes to chaos coexist. We show in Figure 2c,d the corresponding LLE for the two cases plotted in Figure 2a,b with respect to the parameters k 1 and k 2 .

Global Analysis: Multistability and Basin of Attraction

The local analysis discussed previously provides observations about the stability of the map’s equilibrium points around some neighborhoods; it is important to study the global analysis to gain more insight into the complex analysis of the map’s dynamic behavior. Detecting information about the multistability of the map may provide us with some predictions about its future evolution. For that, we perform, in this subsection, a global analysis of the dynamic behavior of the model to obtain projections about the long run of the map’s trajectories. In the previous section, we concluded that the equilibrium point can be destabilized via flip and Neimark–Sacker Bifurcations. Now, we show a variety of different complex dynamic scenarios for the map (13) using numerical simulations. In particular, the chaotic trajectories that can emerge beside the multistability of the map’s behavior where different chaotic attractors may coexist reflect the wide spectrum of the parameter sets. Such complex characteristics are interesting because the demand function adopted by the two firms is linear not nonlinear. Let us start our global analysis by assuming the following parameters set, c 1 = 0.41 , c 2 = 0.45 , θ = 0.9721801 , k 1 = 0.006 and k 2 = 0.006 . At this set, we have two coincided equilibrium points, E r 1 = E r 2 = 0.4378199 , 0.4378199 whose Jacobian matrices (14) become,
J r 1 = 0.81115 0.094426 0.094426 0.80574 , J r 2 = 0.81115 0.094426 1 0
and the eigenvalues in both regions are λ 1 r 1 , 2 r 1 = 0.90291 , 0.71398 and λ 1 r 2 , 2 r 2 = 0.91441 , 0.10326 . They are real eigenvalues and lie within the unit circle λ 1 r 1 , 2 r 1 < 1 and λ 1 r 2 , 2 r 2 < 1 . In addition, we have τ r 1 = 1.61689 , δ r 1 = 0.64466 , τ r 2 = 0.81115 and δ r 2 = 0.094426 . It is clear that δ r 1 < 1 and δ r 2 < 1 that means the map (13) is dissipative in both regions of the border curve Ω = p 1 and the triangles of stability are satisfied. Figure 3a shows the basin of attraction of the coincided equilibrium points born on the border curve. It is plotted at the same set of the parameter values. Now, we show the rise of the multistability of attractors, as the speed of the adjustment parameter k 1 increases, the other parameter values remain fixed. We should note that all of the numerical simulations were carried out based on the initial data of decision variables p 1 , 0 = 0.12 and p 2 , 0 = 0.11 . In the following figures, iterations of the map (13) present different attractors with their basins, using sets of initial conditions at the long run. Increasing the bifurcation parameter k 1 to 0.0659 , the equilibrium point becomes unstable and the dynamics of the map converts into a spiral on one side from the border curve, as shown in the top left in Figure 3b. Increasing that parameter further gives rise to larger spiral points and then they convert into a closed invariant curve due to the Neimark–Sacker bifurcation at k 1 = 0.0664 . This closed invariant curve is tangent to the border curve. This closed invariant curve continues to appear until k 1 = 0.0796 , where the period-2 cycle is born due to the coexistence of flip bifurcation. It is plotted in Figure 3c with its basin of attraction. As one can see, this period is born on region 1, but its basin is distributed in both regions from the border curve. Increasing k 1 further gives rise to the period-4 cycle, which has a complicated basin of attraction distributed in both regions in Figure 3d. Others increasing in this parameter give four chaotic areas, which start gathering into two chaotic areas, then finally become one chaotic attractor, as shown in Figure 3e,f. It is obvious from the 1D bifurcation diagrams presented in Figure 1c,d that the bifurcation parameter k 1 affects the complex dynamic of the map (13). The dynamics start first with a stable equilibrium point, then closed invariant curves in region 2 appear due to Neimark–Sacker. Afterward, higher periodic cycles and chaotic attractors are born in region 1 due to flip bifurcation. This means that the map’s dynamics jumps from Neimark–Sacker to flip bifurcation.
Now, we study another set of parameter values. We assume c 1 = 0.4 , c 2 = 0.8316 , θ = 0.74 , k 1 = 0.1 and k 2 = 0.1 . At this set, we get two coincided stable equilibrium points born on the border curve Ω = p 1 . It is depicted in Figure 4a with its basin of attraction. As previously, we studied the influences of the bifurcation parameter k 1 on the dynamics of the map (13). Increasing k 1 makes the equilibrium point unstable and gives rise to spiral points. Those points are located in region 2, as shown in Figure 4b. They are turned into closed invariant curves in region 2 as k 1 increases further as given in Figure 4c. The closed invariant curves continue to appear because of Neimark–Sacker bifurcation until they turn into a period-11 cycle due to flip bifurcation at k 1 = 0.451 , and keeping the other parameter values as presented in Figure 4d. It is clear that this period cycle is located in region 2 and it does not cross the border curve that means the dynamic of the map (13) is only obtained by the part defined in this region. In Figure 4e, we obtain a period-11 cycle located in region 2 with a chaotic attractor that crosses the border curve and is distributed in both regions at k 1 = 0.452 . It is also clear that the phase plane becomes quite complicated. Increasing k 1 to 0.46 gives a chaotic attractor of the map in both regions and passes through the border curve. This chaotic attractor continues to appear until a period-9 cycle is born at k 1 = 0.4736 with the other parameter values fixed. Figure 5a displays this cycle with its complicated basin of attraction, and it is clear it has points jumping from region 1 to region 2. This indicates that the map (13) must be manipulated and studied as a piecewise smooth map. When k 1 increases further, a one-piece chaotic attractor continues to appear and pass through the border curve. After that, further increasing in the parameter k 1 gives birth to an additional period-10 cycle crossing the border curve. Figure 5b displays this period cycle with its basin of attraction. As k 1 increases, further higher period cycles that are followed by chaotic attractors passing through the border curve arise. After that, the dynamic of the map changes to a period-2 cycle at k 1 = 0.5 , which is followed by a period-4 cycle at k 1 = 0.501 and both cycles are located in region 1. We end our discussion here by presenting the influences of the parameter θ when fixing the map’s parameters. It is obvious that both equilibrium points depend on that parameter, which means we will get two different equilibrium points. In that case, the dynamic of the map will be related to only one region and we will have only one equilibrium point; the other point will be called the virtual equilibrium point. We present, in Figure 5c,d, the bifurcation diagrams when varying the parameter θ . Influences of that parameter on the map’s dynamic are left to the readers.

5. Conclusions

In this paper, we introduced and investigated a piecewise smooth Bertrand game. Our results—by the global analysis—showed multiple attractors due to border collision bifurcations. This is because the map’s phase plane was divided into two regions by a border curve, which has a crucial role in the global dynamics. We analyzed the behavior of the equilibrium points and their stabilities. For two coincided equilibrium points born on the border curve, we showed that they become unstable due to two types of bifurcations, which occurred simultaneously. First, the equilibrium points became unstable due to the coexistence of closed invariant curves by Neimark–Sacker, then those closed invariant curves were destroyed, and different periodic cycles and chaotic attractors emerged because of flip bifurcation. This situation was shown by a 1D bifurcation analysis and then globally by presenting other complex dynamic characteristics, which occurred in each region where the map was defined.

Funding

This research was funded by King Saud University.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

Research supporting project number (RSP-2021/167), King Saud University, Riyadh, Saudi Arabia.

Conflicts of Interest

The author declares that there are no conflicts of interest.

References

  1. Cournot, A. Récherches sur les Principes Mathématiques de la Théorie des Richesses, Dunod, Paris. 1838. Available online: https://www.jstor.org/stable/27579580 (accessed on 5 September 2021).
  2. Bertrand, J. Theorie Mathematique de la Richesse Sociale. J. Savants 1883, 48, 499–508. [Google Scholar]
  3. Askar, S.S. On complex dynamics of cournot-bertrand game with asymmetric market information. Appl. Math. Comput. 2021, 393, 125823. [Google Scholar]
  4. Askar, S.S. Duopolistic Stackelberg game: Investigation of complex dynamics and chaos control. Oper. Res. 2018, 20, 1685–1699. [Google Scholar] [CrossRef]
  5. Askar, S.S.; Abouhawwash, M. Quantity and price competition in a differentiated triopoly: Static and dynamic investigations. Nonlinear Dyn. 2018, 91, 1963–1975. [Google Scholar] [CrossRef]
  6. Askar, S.S.; Elettreby, M.F. The impact of cost uncertainty on Cournot oligopoly games. Appl. Math. Comput. 2017, 312, 169–176. [Google Scholar] [CrossRef]
  7. Askar, S.S. Tripoly Stackelberg game model: One leader versus two followers. Appl. Math. Comput. 2018, 328, 301–311. [Google Scholar] [CrossRef]
  8. Li, Y.; Wang, L. Chaos in a duopoly model of technological innovation with bounded rationality based on constant conjectural variation. Chaos Solitons Fractals 2019, 120, 116–126. [Google Scholar] [CrossRef]
  9. Elsadany, A.A.; Awad, A.M. Dynamics and chaos control of a duopolistic Bertrand competitions under environmental taxes. Ann. Oper. Res. 2018, 274, 211–240. [Google Scholar] [CrossRef]
  10. Agiza, H.N.; Elsadany, A.A. Nonlinear dynamics in the Cournot duopoly game with heterogeneous players. Phys. A 2003, 320, 512–524. [Google Scholar] [CrossRef] [Green Version]
  11. Askar, S.S.; Alshamrani, A.M.; Alnowibet, K. Dynamic Cournot duopoly games with nonlinear demand function. Appl. Math. Comput. 2015, 259, 427–437. [Google Scholar] [CrossRef]
  12. Hommes, C. Behavioral Rationality and Heterogeneous Expectations in Complex Economic Systems; Cambridge University Press: Cambridge, UK, 2013. [Google Scholar]
  13. Peng, Y.; Lu, Q.; Xiao, Y.; Wu, X. Complex dynamics analysis for a remanufacturing duopoly model with nonlinear cost. Phys. A 2019, 514, 658–670. [Google Scholar] [CrossRef]
  14. Elsadany, A.A.; Awad, A.M. Dynamical analysis and chaos control in a heterogeneous Kopel duopoly game. Indian J. Pure Appl. Math. 2016, 47, 617–639. [Google Scholar] [CrossRef]
  15. Askar, S.S.; Alshamrani, A.M.; Alnowibet, K. The arising of cooperation in Cournot duopoly games. Appl. Math. Comput. 2016, 273, 535–542. [Google Scholar] [CrossRef]
  16. Shi, L.; Sheng, Z.H.; Xu, F. The dynamics of competition in remanufacturing: A stability analysis. Econ. Model. 2015, 50, 245–253. [Google Scholar] [CrossRef]
  17. Gardini, L.; Sushko, I.; Matsuyama, K. 2D discontinuous piecewise linear map: Emergence of fashion cycles. Chaos 2018, 28, 055917. [Google Scholar] [CrossRef] [Green Version]
  18. Commendator, P.; Kubin, I.; Petraglia, C.; Sushko, I. Regional integration, international liberalisation and the dynamics of industrial agglomeration. J. Econ. Dyn. Control. 2014, 48, 265–287. [Google Scholar] [CrossRef] [Green Version]
  19. Commendator, P.; Kubin, I.; Mossay, P.; Sushko, I. The role of centrality and market size in a four-region asymmetric new economic geography model. J. Evol. Econ. 2017, 27, 1095–1131. [Google Scholar] [CrossRef] [Green Version]
  20. Radi, D.; Gardini, L.; Avrutin, V. The role of constraints in a segregation model: The symmetric case. Chaos Solitons Fractals 2014, 66, 103–119. [Google Scholar] [CrossRef]
  21. Radi, D.; Gardini, L.; Avrutin, V. The role of constraints in a segregation model: The asymmetric case. Discret. Dyn. Nat. Soc. 2014, 2014, 569296. [Google Scholar]
  22. Agliari, A.; Gardini, L.; Puu, T. Global bifurcations in duopoly when the Cournot point is destabilized via a subcritical Neimark bifuraction. Int. Game Theory Rev. 2006, 8, 1–20. [Google Scholar] [CrossRef]
  23. Puu, T. A new approach to modeling Bertrand duopoly. Rev. Behav. Econ. 2017, 4, 51–67. [Google Scholar] [CrossRef]
  24. Singh, N.; Vives, X. Price and quantity competition in a differentiated duopoly. Rand J. Econ. 1984, 15, 546–554. [Google Scholar] [CrossRef]
  25. Askar, S.S. The rise of complex phenomena in Cournot duopoly games due to demand functions without inflection points. Commun. Nonlinear Sci. Numer. Simul. 2014, 19, 1918–1925. [Google Scholar] [CrossRef]
  26. Askar, S.S.; Al-Khedhairi, A. Analysis of nonlinear duopoly games with product differentiation: Stability, global dynamics, and control. Discret. Dyn. Nat. Soc. 2017, 2017, 2585708. [Google Scholar] [CrossRef] [Green Version]
  27. Agiza, H.N.; Elsadany, A.A. Chaotic dynamics in nonlinear duopoly game with heterogeneous players. Appl. Math. Comput. 2004, 149, 843–860. [Google Scholar] [CrossRef]
Figure 1. (a,b) 2D bifurcation diagram of the system in the ( k 1 , k 2 ) parameter plane. The light gray color refers to the stability region of the fixed point. Other colors are for different types of period cycles. The white color refers to the non-convergent points. (c,d) Bifurcation diagrams with respect to p 1 and p 2 on varying k 1 at: c 1 = 0.41 , c 2 = 0.45 , θ = 0.9721801 and k 2 = 0.006 . (e,f) Bifurcation diagrams with respect to p 1 and p 2 on varying k 2 at: c 1 = 0.41 , c 2 = 0.45 , θ = 0.9721801 and k 1 = 0.006 .
Figure 1. (a,b) 2D bifurcation diagram of the system in the ( k 1 , k 2 ) parameter plane. The light gray color refers to the stability region of the fixed point. Other colors are for different types of period cycles. The white color refers to the non-convergent points. (c,d) Bifurcation diagrams with respect to p 1 and p 2 on varying k 1 at: c 1 = 0.41 , c 2 = 0.45 , θ = 0.9721801 and k 2 = 0.006 . (e,f) Bifurcation diagrams with respect to p 1 and p 2 on varying k 2 at: c 1 = 0.41 , c 2 = 0.45 , θ = 0.9721801 and k 1 = 0.006 .
Mathematics 09 02558 g001
Figure 2. (a) Bifurcation diagrams with respect to p 1 and p 2 on varying k 1 at: c 1 = 0.4 , c 2 = 0.8316 , θ = 0.74 and k 2 = 0.1 . (b) Bifurcation diagrams with respect to p 1 and p 2 on varying k 2 at: c 1 = 0.4 , c 2 = 0.8316 , θ = 0.74 and k 1 = 0.1 . (c) Largest Lyapunov exponent on varying k 1 for c 1 = 0.4 , c 2 = 0.8316 , θ = 0.74 and k 2 = 0.1 (d) Largest Lyapunov exponent on varying k 2 for c 1 = 0.4 , c 2 = 0.8316 , θ = 0.74 and and k 1 = 0.1 .
Figure 2. (a) Bifurcation diagrams with respect to p 1 and p 2 on varying k 1 at: c 1 = 0.4 , c 2 = 0.8316 , θ = 0.74 and k 2 = 0.1 . (b) Bifurcation diagrams with respect to p 1 and p 2 on varying k 2 at: c 1 = 0.4 , c 2 = 0.8316 , θ = 0.74 and k 1 = 0.1 . (c) Largest Lyapunov exponent on varying k 1 for c 1 = 0.4 , c 2 = 0.8316 , θ = 0.74 and k 2 = 0.1 (d) Largest Lyapunov exponent on varying k 2 for c 1 = 0.4 , c 2 = 0.8316 , θ = 0.74 and and k 1 = 0.1 .
Mathematics 09 02558 g002
Figure 3. (a) The basin of attraction of the equilibrium points. (b) Spiral and closed invariant curves at different values of k 1 . (c) Basin of attraction of period-2 cycle at k 1 = 0.0827 . (d) Basin of attraction of period-4 cycle at k 1 = 0.0834 . (e) Phase plane for the four pieces chaotic areas at k 1 = 0.0836 . (f) Phase plane for the one piece chaotic attractor at k 1 = 0.0845 . Other parameters values are: c 1 = 0.41 , c 2 = 0.45 , θ = 0.9721801 and k 2 = 0.006 .
Figure 3. (a) The basin of attraction of the equilibrium points. (b) Spiral and closed invariant curves at different values of k 1 . (c) Basin of attraction of period-2 cycle at k 1 = 0.0827 . (d) Basin of attraction of period-4 cycle at k 1 = 0.0834 . (e) Phase plane for the four pieces chaotic areas at k 1 = 0.0836 . (f) Phase plane for the one piece chaotic attractor at k 1 = 0.0845 . Other parameters values are: c 1 = 0.41 , c 2 = 0.45 , θ = 0.9721801 and k 2 = 0.006 .
Mathematics 09 02558 g003
Figure 4. (a) The stable equilibrium point with its basin of attraction around it. Parameters values are c 1 = 0.4 , c 2 = 0.8316 , θ = 0.74 , k 1 = 0.1 and k 2 = 0.1 . (b) Spiral points at different values of k 1 . (c) Closed invariant curves at different values of k 1 . (d) The basin of attraction for the period-11 cycle at k 1 = 0.451 . The gray color denotes divergent and infeasible trajectories. (e) Basin of attraction of the period-11 cycle with the chaotic attractor. (f) Phase plane of the chaotic attractor.
Figure 4. (a) The stable equilibrium point with its basin of attraction around it. Parameters values are c 1 = 0.4 , c 2 = 0.8316 , θ = 0.74 , k 1 = 0.1 and k 2 = 0.1 . (b) Spiral points at different values of k 1 . (c) Closed invariant curves at different values of k 1 . (d) The basin of attraction for the period-11 cycle at k 1 = 0.451 . The gray color denotes divergent and infeasible trajectories. (e) Basin of attraction of the period-11 cycle with the chaotic attractor. (f) Phase plane of the chaotic attractor.
Mathematics 09 02558 g004
Figure 5. (a) The basin of attraction for the period-9 cycle. (b) The basin of attraction for the period-10 cycle. (c) Bifurcation diagram with respect to p 1 on varying θ at: c 1 = 0.2 , c 2 = 0.1 , k 1 = 0.1 and k 2 = 0.1 . (d) Bifurcation diagram with respect to p 2 on varying θ at: c 1 = 0.2 , c 2 = 0.1 , k 1 = 0.1 and k 2 = 0.1 .
Figure 5. (a) The basin of attraction for the period-9 cycle. (b) The basin of attraction for the period-10 cycle. (c) Bifurcation diagram with respect to p 1 on varying θ at: c 1 = 0.2 , c 2 = 0.1 , k 1 = 0.1 and k 2 = 0.1 . (d) Bifurcation diagram with respect to p 2 on varying θ at: c 1 = 0.2 , c 2 = 0.1 , k 1 = 0.1 and k 2 = 0.1 .
Mathematics 09 02558 g005
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Askar, S. Complex Investigations of a Piecewise-Smooth Remanufacturing Bertrand Duopoly Game. Mathematics 2021, 9, 2558. https://doi.org/10.3390/math9202558

AMA Style

Askar S. Complex Investigations of a Piecewise-Smooth Remanufacturing Bertrand Duopoly Game. Mathematics. 2021; 9(20):2558. https://doi.org/10.3390/math9202558

Chicago/Turabian Style

Askar, Sameh. 2021. "Complex Investigations of a Piecewise-Smooth Remanufacturing Bertrand Duopoly Game" Mathematics 9, no. 20: 2558. https://doi.org/10.3390/math9202558

APA Style

Askar, S. (2021). Complex Investigations of a Piecewise-Smooth Remanufacturing Bertrand Duopoly Game. Mathematics, 9(20), 2558. https://doi.org/10.3390/math9202558

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop