Abstract

Based on a nonlinear demand function and a market-clearing price, a cobweb model is introduced in this paper. A gradient mechanism that depends on the marginal profit is adopted to form the 1D discrete dynamic cobweb map. Analytical studies show that the map possesses four fixed points and only one attains the profit maximization. The stability/instability conditions for this fixed point are calculated and numerically studied. The numerical studies provide some insights about the cobweb map and confirm that this fixed point can be destabilized due to period-doubling bifurcation. The second part of the paper discusses the memory factor on the stabilization of the map’s equilibrium point. A gradient mechanism that depends on the marginal profit in the past two time steps is adopted to incorporate memory in the model. Hence, a 2D discrete dynamic map is constructed. Through theoretical and numerical investigations, we show that the equilibrium point of the 2D map becomes unstable due to two types of bifurcations that are Neimark–Sacker and flip bifurcations. Furthermore, the influence of the speed of adjustment parameter on the map’s equilibrium is analyzed via numerical experiments.

1. Introduction

Different disciplines such as biology, engineering, and economy are characterized by real-life models which possess complex dynamic behaviors and multistability criteria. Economists and mathematicians have reported these behaviors in economic models such as monopoly, duopoly, and oligopoly models. In that direction, many studies have reported that the equilibrium points of these models can be destabilized due to different types of bifurcations such as flip and Neimark–Sacker. Our contribution in this paper focuses on further investigating on complex dynamic characteristics of the cobweb model constructed based on the nonlinear demand function and gradient mechanism. We follow in our paper the nice discussions given in the seminal articles of Ezekiel [1] and Naimzada [2]. In [1], the effects of prices in the fluctuations of certain markets have been explained. Some important aspects on supply and market demand have been investigated in [1]. Ezekiel in [3] has claimed that the quantity of production must be set based on time so that firms’ producers can review prices and hence a time lag has to be made between supply and demand.

The current paper follows up the hypothesis of Ezekiel and the discussion made by Naimzada in order to set up and analyze the suggested cobweb model. Building this model does not require assuming that the firms (or producers) have no full knowledge about market information but they have to have only knowledge about the demand. We adopt here an irrational function that is used to represent the demand. Using this kind of knowledge and the clearing price condition, the producer can update its production due to the variations that occurred in the marginal profit. Recalling the bounded rationality mechanism defined elsewhere [2], the producer may decrease or increase its outputs depending on the decrease or increase taking place in the marginal profit. In order to brief the outcomes of this paper, we divide the contribution into two important parts: the first part deals with the studies and investigations on a one-dimensional nonlinear discrete-time map describing the change in price. It has four nonzero real fixed points and only one of them attains the maximum value of the profit. The stability conditions of this point are calculated and some insights about it are provided based on analytical and numerical analysis. In the second part, we introduce the memory and convert the 1D map into a two-dimensional map. The 2D map contains two important parameters; the parameter of memory and the other for the adjustment speed. The memory parameter is represented by some weights that their influences on the stability of the equilibrium point are discussed. Through some analytical and numerical experiments, we show that low and high values for the memory parameter affect the stabilization of the equilibrium point due to the coexistence of Neimark–Sacker and period-doubling bifurcations. Based on the global analysis, some dynamic behaviors of the map such as multistability and chaotic attractors are reported. Some of those behaviors possess attractive basin with lobes from divergent points. Furthermore, we study the influence of the speed of adjustment by assuming equal memory weights. Some dynamic characteristics are obtained and discussed through pushing the global analysis towards the speed of adjustment parameter.

Both local and global analyses regarding bifurcation analysis and chaotic attractors within this paper have been carried out based on similar investigations and discussions. For instance, Shabbir et al. in [4] have introduced a new discrete-time system based on cannibalism in the prey population with the addition of Allee effects. They have calculated the system’s fixed points and have discussed their stability that becomes unstable due to period-doubling and Neimark–Sacker bifurcation. In [5], the dynamic characteristics such as topological classification, Lyapunov exponents, and manifold theory have been investigated for the Goodwin model. Din and Haider in [6] have studied the complex dynamic characteristics of the discrete-time version of the Schnakenberg model. In [7], a deterministic oligopoly model whose players are heterogeneous has been introduced and its dynamic behaviors that are formed based on best responders and imitators have been discussed. Agliari et al. in [8] have given an intensive dynamic analysis for a Cournot duopoly game with differentiated goods whose players are homogenous and have adopted a gradient-based mechanism in order to update their outputs. Other interesting studies that have adopted discrete dynamic systems and their complex dynamic characteristics in different applications such as engineering have been reported in Bao et al. [9], Li et al. [10], and Bao et al. [11].

Briefly, the paper consists of six sections. The first section introduces the main contribution of the paper. In Section 2, the literature review is reported. In Section 3, the 1D cobweb map is given and studied. The influence of memory is investigated in Section 4. The influence of the speed of adjustment is analyzed in Section 5. Finally, we conclude the obtained results in Section 6.

2. Literature Review

Several economic contexts have reported in literature different studies and investigations on the cobweb model. We highlight in this section some of those studies. For example, different sectors in the economic market have reported the cobweb model where there are academia [2], real estate [12], nurses [13], potatoes [14], and bioenergy crops [15]. Based on nonlinear demand and supply functions the equilibrium points of cobweb mode have been calculated and their stability conditions have been investigated in [16]. In [17], the adaptive expectation was used to form complex cobweb models. In [18], the traditional cobweb model has been constructed using a nonlinear supply function and it has been reported that period-doubling cycles can be found when firms adopt the adaptive expectation mechanism. In [19], the authors have detected chaotic behaviors for a standard cobweb model whose players adopted adaptive expectations using monotonic demand and supply functions. Based on linear demand and nonlinear supply functions, the dynamics of the cobweb map have been investigated under adaptive expectations in [20].

The studies mentioned above have proposed certain assumptions about the producers (or firms). They have assumed that those producers possessed somehow the cost function and hence they can calculate the profit which depends on produced quantity and its price. Indeed, this gives producers some knowledge on the determination of quantities sent to the economic market as price-based functions. Such hypotheses have given rise to many mechanisms such as the bounded rationality mechanism to model the dynamic economic behaviors behind these hypotheses. For example, a cobweb model has been constructed and described by a nonlinear discrete dynamical system using a gradient mechanism in [2]. The authors in this study assumed that producers have no complete knowledge of demand and they instead perform empirical estimations on marginal profit. The bounded rationality approach is an effective mechanism adopted when dealing with such economic systems and information about it has been reported elsewhere in the literature ([3, 2130]). Furthermore, the bounded rationality has been adopted in many studies in the literature that require knowledge and computational experiments. For instance, in [31], a cobweb model whose supply and demand are nonlinear functions has been investigated. In that study, the producers have recalled the mechanism of backward expectation to make prices forecasting in the future evolution. This mechanism adopted in [31] depended on prices given in the last two periods. In [32], another mechanism that is called naive expectation has been used by producers with more general demand and cost functions in a cobweb model. Other studies in the literature have suggested rational producers in [33], heterogeneous producers in [34], and replicator dynamics in [35].

3. The Cobweb Model

Let be the consumer demand, where refers to price. This demand function may be linear or nonlinear depending on the preferences of the consumer. In this paper, we consider throughout this paper the following nonlinear demand:where represents the time. It is clear that means the market is not supported by quantities and then the consumer is willing to buy a unit of production with a maximum price that is . From (1) we getthat represents the inverse of demand. In order to provide the market with produced quantities at time , producers make the estimation of some factors raised at time . Of these factors are the production volume and gained profit. They receive lack of market information and then perform market experiments so that they can decide the market state. Based on these factors, producers decide to supply the market with quantities demanded. They increase the amount of production in the next period of time if the profit increases and they reduce the amount if the profit is decreased. The profit of quantity supplied to the market may be given bywhere represents supplied quantity at time and denotes the total cost. This total cost is taken as a linear function as follows:where is the marginal cost . The cobweb theory states that attaining market equilibrium can be obtained by assuming at any time . Therefore, (3) can be rewritten as follows:

From (5), one gets a positive profit if . This indicates that the maximum price paid for buying a unit of commodity must be greater than the value . On the other hand, if , the profit becomes negative that has no economic meaning. Differentiating (5) gives the following marginal profit:

Now, the quantity supplied to the market at can be stated as follows:where we take as a positive parameter known as the speed of adjustment. Using (1) in (7), one can simplify the supply in the formand then market-clearing price gives

From (1) and (8), we get

After some simple calculations, the price at time can be expressed by the following one-dimensional discrete dynamic map:

It possesses four fixed points given by

Since the parameters , and are positive, then both and are positive. The other two points are positive under certain conditions. Before we discuss the stability conditions for these fixed points and detect their complex dynamic behaviors, we investigate the relation between the demand and supply depicted in Figures 1(a) and 1(b). One can see that both and intersect in the equilibrium price . It must not exceed the maximum price that is guaranteed based on the condition ; otherwise, a negative price may be raised if . Starting at while and moving vertically until reaching which is represented by the blue curve, we get the quantity in the next period that is denoted by in Figure 1(a). Since our model is in an equilibrium state at each period of time, we have equal demand and supply and consequently, we move horizontally to (represented by the red curve in Figure 1(a)) to get the price that is being used to set in the next period. After setting , we move horizontally to (represented by the red curve in Figure 1(a)). This procedure is repeated till reaching the equilibrium of market price . Figure 1(b) shows this procedure for an initial price . From this discussion, we conclude that the trajectories of cobweb will not approach negative values.

Now, we begin discussing the dynamics of map (11) around the fixed points which depend on the parameters and . The profit represented by (5) can be rewritten as follows:

Differentiating (13), we get but the profit attains its maximum value at the equilibrium price . Furthermore, we have and if . This means that when we increase the price in the interval , the profit variations become positive while they get negative for . This may be explained due to the monotonic or nonmonotonic way of the production level carried out by producers and the way they deal with the change in profit. Figure 2 confirms that the marginal profit intersects the abscissa axis exactly in (the market equilibrium price).

Proposition 1. The stability/instability of market equilibrium price is constrained by the following possibilities:(i)The equilibrium price is locally stable if (ii)The equilibrium price is unstable if where

Proof. Recalling the marginal demand and marginal supply at givesUsing the conditions given in [36], we get(i)The equilibrium price is locally stable if (ii)The equilibrium price is unstable if So, the proof is completed after substituting (15) in the above conditions.
Now, we study the influences of the parameters , and on the dynamics of the map. We perform some numerical experiments to investigate more these dynamics. The influence of the speed of adjustment parameter while fixing the other parameters is given in Figure 3(a). We assume the following parameters’ values . Figure 3(a) shows two situations: the first one presents the period-doubling diagram with respect to at and . It shows that the equilibrium point becomes locally stable; then, it becomes unstable. For this situation, several numerical experiments have been carried out using low marginal cost and we have concluded that such low marginal cost guarantees stable situation for the equilibrium price. On the other hand and as shown in Figure 3(b), the increase in marginal cost causes a reduction in the region of stability for the equilibrium price. It would be a stability region as long as the condition that preserves nonnegative price and stability of the equilibrium market price too. For the parameters and , we fix the speed of adjustment parameter and study the impact of and . As given in Figure 3(b), the parameter affects the stability of the equilibrium point and makes it loses its stability via period-doubling bifurcation for low and high marginal cost. The same discussion is given for the parameter as depicted in Figure 3(c). These numerical experiments confirm the impact of the map’s parameters on the local stability of the equilibrium price. This urges us to investigate more some of the global investigation of the dynamics of the map. We begin with the set of parameters values, . Figure 3(d) shows that there is a monotonic convergence to the equilibrium market price . Further increase in the parameter to the value 0.9525 gives rise to a sequence of period-doubling whose attractive basin is given in Figure 3(e). As increases further and the other parameters are fixed, higher cycles are obtained until where a chaotic behavior of the map is born in Figure 3(f). Other numerical experiments have been carried out for different values of and have shown that the attracting set of any period cycles or chaotic behavior will have an attractive basin that is bounded within the box . Therefore, any initial prices that are chosen out of this box will lead to negative or unbounded trajectories that have no meaning in economy.

4. Memory Effect on the Cobweb Model (2D Map)

This section introduces the memory in the map (11). The memory here means that producers decide whether they increase or decrease the production in the next period of time. It is given by [37]

This memory is estimated by producers based on the marginal profits at two periods, the previous period and the current period with weight in order to provide the market with the supplied quantity at . Using (6), we get

Imposing the market clearing gives

Setting and in (18), a two-dimensional (2D) map is obtained:

For the map (19), we get four fixed points and we focus only on the one in the form where is given by (11). The following propositions are used in discussing the stability/instability conditions of this point.

Proposition 2. Suppose that and are two eigenvalues for the Jacobian matrix of (19) at ; then,(i) is an attracting node and is locally stable if .(ii) is repelling node and unstable if .(iii) is a saddle point and is unstable if and (or and ).(iv) is a nonhyperbolic point if and (or and ).

Proposition 3. The point is local stable if the following conditions are satisfied:

Proof. Proof. We recall Jury conditions [38] that are given bywhere and represent trace and determinant of the Jacobian matrix of the map (19) at and are given byWe should highlight here the types of bifurcations that may have occurred based on these conditions. If while and , then becomes unstable due to transcritical or fold bifurcation. If while and , then becomes unstable due to period-doubling bifurcation. But if and and , then can be destabilized due to Neimark–Sacker bifurcation ([8, 29]). Substituting (22) in (21) completes the proof.
Hence, the conditions in (20) are used to detect the conditions of stability/instability of . Because of the complicated form of the equilibrium point , some numerical simulations are performed in order to get more insights about the conditions given in (20). It is clear that both and contain the parameter that represents the memory. Therefore, we analyze the effects of on the stability of the equilibrium point . We begin with the set , and . It gives and then the Jacobian becomesFor , we get two complex conjugate eigenvalues, with , and hence is locally stable point. Furthermore, we have and and hence means that the map (19) is dissipative. It is easy to see that the conditions in (20) are all satisfied and then is locally asymptotically stable. In Figure 4(a), we depict how the parameter affects the equilibrium point . It shows the bifurcation diagram and that there are two types of bifurcations found between and the variable . The assumed set of parameters’ values shows that the equilibrium point becomes unstable because of the coexisting of Neimark–Sacker bifurcation for values in the interval while, in the interval , the equilibrium point becomes unstable due to flip bifurcation. Now, we perform some numerical experiments in order to investigate more the global dynamic behavior of the map (19) around . Keeping the parameters’ values fixed and changing to , the dynamics of the map around the equilibrium are converted into four closed rings. The phase portrait for those closed rings is given in Figure 4(b). Increasing to 0.33254 gives rise to a period-4 cycle (denoted by squares) with a disconnected chaotic attractor (red) around the unstable equilibrium for the same set used for the bifurcation diagram. It is plotted in Figure 4(c) with its basins of attraction that it consists of two colors: yellow and cyan. A further increase in the memory parameter to makes the map becomes more chaotic and we get a chaotic behavior around as shown in Figure 4(d). At , the attractive basins of period-9 cycle (denoted by squares) together with the basin of a closed chaotic attractor (red) that surrounds the unstable are displayed in Figure 4(d). This period-9 cycle (as shown in Figure 4(e)) disappears when we increase to 0.5. At , we get a period-4 cycle together with a closed invariant curve surrounding the point . This dynamic situation is given in Figure 4(f) where the yellow and cyan colors refer to the attractive basins of the period-4 cycle. For further increase in , the previous dynamic situation is disappeared and we get instead a spiral point. We give in Figure 4(g) different dynamic behaviors for the map at and . In Figure 4(h), we depict the stability region for the point at the parameters set, in the . It is clear that when and , the point becomes unstable via the Neimark–Sacker bifurcation as increases. As becomes close to a uniform distribution, the stability region of is extended as increases. For values of above 0.75, this stability region of is reduced with respect to and then becomes unstable through flip bifurcation. This investigation makes us analyze in the the 2D bifurcation diagram. This bifurcation diagram is given in Figure 5 at the set of parameters values, and . We give at the end of this section a dynamic situation of chaotic attractor that is formed by four disconnected bands around for a value of near to 1 as shown in Figure 6.

5. The Speed of Adjustment and Its Effect

We assume that represents a symmetric case in the map (19) as follows:

It means that producers adopt the average in marginal profits at the time steps and . Now, the map (24) contains the parameters , and and we study in this section the effect of on the dynamics of the map. Assuming and , Figure 7(a) shows the 1D bifurcation diagram when varying the bifurcation parameter . It shows that the equilibrium point becomes stable for the values of until approaches the value where the point becomes unstable due to the Neimark–Sacker bifurcation. Increasing to 1.423, the equilibrium point converts into an unstable spiral point as shown in Figure 7(b). A further increase in to 1.58 gives rise to a period-4 cycle. The basins of attraction of this cycle consist of yellow and cyan colors (as shown in Figure 7(c)) while the white color refers to the nonconvergent points. This cycle occurs as increases until where this cycle is converted into a chaotic attractor consisting of four disconnected bands (see Figure 7(d)). At , a period-5 cycle is born which has attractive basins with yellow and cyan colors. The dark gray refers to divergent points as given in Figure 7(e). After that, the dynamic of the map becomes chaotic for any increase in the parameter . Figure 7(f) gives an example of the chaotic situation when . Regarding the parameters and numerical experiments show that increasing both parameters affect the stability region of with respect to .

6. Conclusion

In this work, we have investigated a cobweb model whose producers do not possess complete knowledge about the market and update their outputs according to the mechanism of bounded rationality. The producers have estimated their marginal profit by observing profit variations that might have occurred at the beginning of production. In the first part of this paper, we have studied the complex and chaotic behaviors of a one-dimensional discrete cobweb model. Our obtained results have shown that any increase in marginal cost leads to an unstable equilibrium price. Furthermore, we have shown that the equilibrium price for the 1D map becomes locally stable when the adjustment speed parameter takes low values. Higher values of that parameter have given rise to unstable price and therefore periodic cycles or chaotic attractors for the dynamics of price have coexisted.

The second part of this paper contains other contributions that are the inclusion of memory in the 1D model and converting it into a 2D discrete dynamic model. Based on a convex combination of marginal profit at the past two periods, the memory parameter has been incorporated ([31, 37]). Our investigations have analyzed the influence of the memory parameter on the equilibrium price and found that it has a qualitative effect on its stability. Using both analytical and numerical analysis we have discussed the stability/instability conditions of the equilibrium price. We have found that the equilibrium price may be unstable through period-doubling and Neimark–Sacker bifurcations when the memory parameter is used as the bifurcation parameter. Through simulation, we have reported multistability situations. Furthermore, the qualitative impact of the adjustment parameter has been analyzed in a symmetric case on which the average change in marginal profits has been considered.

Data Availability

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

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding this article.

Acknowledgments

The authors are grateful to the Research Supporting Project number (RSP-2020/167), King Saud University, Riyadh, Saudi Arabia. This project is funded by Saud University, Riyadh, Saudi Arabia.