Discrete Dynamics in Nature and Society

Volume 2014 (2014), Article ID 908093, 13 pages

http://dx.doi.org/10.1155/2014/908093

Research Article

Chaotic Characteristics and Application of Cooperative Game and Evolutionary Game

Yujing Yang, Junhai Ma, and Hongliang Tu

College of Management Economic, Tianjin University, Tianjin, China

Received 2 October 2013; Accepted 26 December 2013; Published 19 March 2014

Academic Editor: Zbigniew Leśniak

Copyright © 2014 Yujing Yang 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.

Abstract.

According to a dynamical multiteam Cournot game in exploitation of a renewable resource, a new dynamic Cournot duopoly game model with team players in exploitation of a renewable resource is built up in this paper. Based on the theory of bifurcations of dynamical systems, the stability of the system is studied and the local stable region of Nash equilibrium point is obtained. The effect of the output adjustment speed parameters and the weight parameter of the system on the dynamic characteristics of the system are researched. The complexity of the system is described via the bifurcation diagrams, the Lyapunov exponents, the phase portrait, the time history diagram, and the fractal dimension. Furthermore, the chaos control of the system is realized by the parameter adjustment method. At last, an evolutionary game as a special dynamic system is constructed and analyzed which is more useful and helpful in application. The derived results have very important theoretical and practical values for the renewable resource market and companies.

1. Introduction

Chaos has become a hot topic in the competition of oligarchs. Research on the complexity of the oligopoly game model has been paid attention to by researchers and scholars recently. Reference [1] studied different strategies which are the Cournot model, the Stackelberg case, and the dynamic system of a duopoly game and investigated stable equilibrium point, cycles, bifurcation, and chaos of the systems. Reference [2] obtained the explicit stability zones for Cournot game with 3 and 4 competitors, and [3] analyzed the stability, bifurcation, chaos, and chaos control of a Kopel model. Reference [4] studied the chaotic dynamics in nonlinear duopoly game with bounded and naive players. Reference [5] constructed a stochastic optimal control model of pollution abatement and studied the control of the model. Reference [6] studied the dynamical behaviors of a duopoly game with delayed bounded rationality and obtained some practical and theoretical significance in the practice. Reference [7] studied the complex dynamics in nonlinear triopoly market with different expectations and obtained some useful conclusions. Reference [8, 9] made the analysis on the complexity of a Cournot-Bertrand duopoly game model with limited information and introduced a 4D Hyperchaotic System, making the numerical simulation and achieving a second control. Reference [10] introduced a social behavior selection model for evolutionary dynamics of behaviors in social networks which exhibits a rich set of emergent behaviors of evolution and studied the chaotic dynamics of the social behavior selection networks in crowd simulation. Reference [11] studied the dynamics and adaptive control of a duopoly advertising model based on heterogeneous expectations and gave the scope of the convergent condition and control intensity. And the result demonstrated that the adaptive controlled method could lead the chaos track to low cycle track. In particular, cooperation and multiteam games had become popular in the fields of chaos in recent years. Reference [12] made analysis on whether the migration cost influences cooperation among success-driven individuals. Reference [13] studied cooperation in an evolutionary prisoner’s dilemma game with probabilistic strategies. Reference [14] made analysis on dynamical multiteam and signaling games. References [15, 16] formulated the standard and multiteam Bertrand game which is based on Puu’s incomplete information and studied quantum team games, respectively. Reference [17] proposed two different versions of the multiteam model: one is that the firms in each team can help each other, and the other is that they competed with another team. Finally, the conditions of their local asymptotic stability were studied. Reference [18] studied an incomplete information dynamical system and a modification of the system was suggested and applied to the standard Cournot game and formulated and studied the multiteam dynamic Cournot game. Reference [19] mainly established a dynamical multiteam Cournot game in exploitation of a renewable resource and analyzed the asymptotic stability of the equilibrium solution of the game. It is worthy to further analyze and extend it in this area. Reference [20] analyzed the purely competitive evolutionary dynamics for games and proposed that perhaps the system will enter a chaotic state for some parameters.

Reference [19–24] have analyzed the dynamic characters in the economic fields, and drew some valuable conclusions.

The structures of this paper are as follows. In Section 2, the dynamic Cournot duopoly game model is established. In Section 3, the existence and the local stability of the Nash equilibrium point are studied, and the stable region of the Nash equilibrium point is derived. In Section 4, dynamical characteristics of the system are analyzed. The effect of the output adjustment speed parameters and the weight parameter on the dynamics behaviors of the system is discussed. Numerical simulation method is used to show complex dynamics of the system by means of the bifurcation diagram, the Lyapunov exponents, the phase portrait, the time history diagram, and the fractal dimension. In Section 5, bifurcation and chaos control of the model are used by the parameter adjustment method. In Section 6, a model of evolutionary game is analyzed by numerical simulation. We obtain the conclusion that evolutionary game is a special dynamic by analyzing the change of Lyapunov exponents. Finally, some results are summarized.

2. The Model

Suppose that there are two mainly oligopoly companies , in exploitation of a renewable resource. The company    makes the optimal output decision, and suppose that the -output is   . At each period , the price is determined by the total output .

According to [19], we also propose the renewable resource market with the linear inverse demand function:

and the cost function of the company    is as follows:

where    is the fixed cost and    is positive parameter.

We can get the profit of the company    as follows:

Here, it is assumed that the renewable resource firms , establish a team, and the profit of the team is

where the is the weight parameter of the profit of the renewable resource company in the team. We propose that the company    takes bounded rational strategy under the maximum of team profit. Since the game between the companies is a continuous and long-term repeated dynamical process, the dynamical adjustment of the player in the duopoly game can be expressed as follows:

where    is the output adjustment speed parameter.

Combining (4) and (5), a new dynamic Cournot duopoly game with team players in exploitation of a renewable resource is obtained. This map has the following form:

3. The Stability of the System

In system (6),    is taken as bifurcation parameter, and the other parameters are as follows: , , , , , and .

Letting the marginal profits be equal to 0, we can get the Nash equilibrium point. The fixed points of system (6) satisfy the following equations:

The equations [21] are solved and three meaningful fixed points , , and are obtained. Here, we only consider the stability of the Nash equilibrium point and denote the total output at the Nash equilibrium point as .

We can calculate the Jacobian matrix of system (6) at the Nash equilibrium point :

where

Furthermore, we can get the characteristic polynomial of system (6) at :

where

According to the Jury test [21], the necessary and sufficient condition of the local stability of Nash equilibrium should satisfy the following three conditions:

(i),(ii),(iii).

The local stable region of Nash equilibrium point can be obtained by solving the above equations. The stable region of Nash equilibrium point about the output adjustment speed parameter    is shown in Figure 1. The Nash equilibrium is stable for the values of and inside the stable region. The meaning of the stable region is that both the two renewable resource companies will eventually maintain at Nash equilibrium output after finite games whatever initial output is chosen in the local stable region. It is valuable to study that the renewable resource companies increase the output adjustment speed in order to get more profit. Though output adjustment speed parameters are unrelated to the Nash equilibrium point, the system will become unstable and even fall into chaos if one renewable resource company that adjusts output speed is too fast and pushes the output adjustment speed parameter    out of the stable region. Numerical simulation method is used to analyze the dynamic characteristics of system (6) with the change of   . Numerical results such as the bifurcation diagrams, the strange attractors, the Lyapunov exponents, the sensitive dependence on initial conditions, and the fractal structure will be considered in the following section.

Figure 1: The stable region of Nash equilibrium point about and .

4. Dynamic Characteristics of the System

In this game, the renewable resource companies make the optimal output decision to get the maximum profit and adjust their output based on the marginal profit of last period. In the oligopoly market, the players have the driving force to increase their output in the hope of achieving more profits. The renewable resource firms can adjust their output speed to increase their output. So, the output adjustment speed parameter    affects the game results very much. The output adjustment speed parameter    affects the dynamic behaviors of system (6), which will be analyzed in the following section.

4.1. The Effect of Output Adjustment Speed on the System

The stability of Nash equilibrium point will change if company accelerates the output adjustment speed and pushes out of the stable region. When , with the output adjustment speed increasing, the output evolution of the duopoly starts with equilibrium state, through period doubling bifurcation, and ends with chaotic state. From the bifurcation diagram in Figure 2 and the corresponding Lyapunov exponents diagram in Figure 3, we can see that is the domain of system (6) in equilibrium state as shown in Figure 4, is the domain of 2 period orbits of system (6), and is the domain of 4 period orbits of system (6), is the domain of 8 period orbits of system (6), for , system (6) is in a chaotic state. Calculation of the Lyapunov exponent is used to analyze the quantitative characteristics of the dynamic system. The system is in chaotic state if the largest Lyapunov exponent is positive. In addition, the larger the positive Lyapunov exponent is, the more obvious the chaotic state of the system is. Figure 5 shows the representative chaotic attractor at initial point and , and equals  .

Figure 2: Bifurcation diagram of system (6) with and .

Figure 3: The corresponding Lyapunov exponents diagram with and .

Figure 4: For and , the wave plots of system (6) with the change of time.

Figure 5: Chaos attractor of system (6) for and .

Furthermore, the bifurcation diagram in Figure 6 and the corresponding Lyapunov exponents diagram in Figure 7 show a one-parameter bifurcation diagram with respect to when is . We can see that is the domain of system (6) in equilibrium state. is the domain of 2 period orbits of system (6). is the domain of 4 period orbits of system (6). is the domain of 8 period orbits of system (6). For , system (6) is in a chaotic state. Figure 8 shows the representative chaotic attractor at initial point and , and equals .

Figure 6: Bifurcation diagram of system (6) with and .

Figure 7: The corresponding Lyapunov exponents diagram with and .

Figure 8: Chaos attractors of system (6) for and .

From the analysis above, we can find that the stability of system (6) will be changed and even the complex dynamic behaviors occur with the increase of the output adjustment speed   .

4.2. The Weight Parameter Effects on the System

In this section, we study the weight parameter which affects the dynamic behaviors of system (6).

Likewise, the bifurcation diagram in Figure 9 and the corresponding Lyapunov exponents diagram in Figure 10 shows a one-parameter bifurcation diagram with respect to when and . is the domain which has economic meaning. One can see that system (6) is stable at different Nash equilibrium points when ; that is to say, the duopoly renewable resource companies are in different equilibrium states. is the domain of 2 period orbits of system (6) as shown in Figure 11. In both sides of the stable domain, system (6) undergoes period doubling bifurcations to chaos with the weight parameter decreasing and increasing, which only has dynamic meaning, and they do not have economic meaning. For , the output , but the output is in chaotic state. For , the output , but the output has 8-cycle output orbits as shown in Figure 13. For , the output , but the output has 4-cycle output orbits as shown in Figure 12. For , the output , but the output has different equilibrium states. For , the output , but the output has 2-cycle output orbits as shown in Figure 14. For , the output , but the output have 4-cycle output orbits as shown in Figure 15. For , the output , but the output has 8-cycle output orbits as shown in Figure 16. For , the output , but the output is in chaotic state as shown in Figure 17.

Figure 9: Bifurcation diagram of system (6) with and and .

Figure 10: The corresponding Lyapunov exponents diagram with and and .

Figure 11: The wave plots of the output with the change of time when , and .

Figure 12: The wave plots of the output with the change of time when , and .

Figure 13: The wave plots of the output with the change of time when and , .

Figure 14: The wave plots of the output with the change of time when , and .

Figure 15: The wave plots of the output with the change of time when , and .

Figure 16: The wave plots of the output with the change of time when , and .

Figure 17: The wave plots of the output with the change of time when , and .

From the analysis above, we can see that the weight parameter not only affects the stability of system (6), but also changes the Nash equilibrium point. In both sides of the stable domain, system (6) undergoes period bifurcations to chaos with the weight parameter decreasing or increasing, which has no economic meaning.

4.3. Sensitive Dependence on Initial Conditions

One of the most important characteristics of the chaos is extremely sensitive dependence on initial conditions. Figures 18 and 19 and Figures 20 and 21 show the output history diagram and difference with the change of time when system (6) has different initial conditions. We can see that there is almost no distinction between them in the beginning, but the difference becomes more and more between them with the number of games increasing. It implies that a slight difference between initial values can lead to a great effect on the game results. For , and and , and , it further affirms that system (6) is in a chaotic state. When the system is in a chaotic state, the market will be destroyed and it is difficult for the renewable resource companies to make long-term plan. So, every action from the renewable resource firms can result in a great loss.

Figure 18: For , and , and the initial conditions which are (0.56, 0.92) and (0.561, 0.92), the history diagram with the change of time about the output .

Figure 19: For , and , and the initial conditions which are (0.56, 0.92) and (0.561, 0.92), the difference with the change of time about the output .

Figure 20: For , and , and the initial conditions which are (0.76, 0.92) and (0.76, 0.921), the history diagram with the change of time about the output .

Figure 21: For , and , and the initial conditions which are (0.76, 0.92) and (0.76, 0.921), the difference with the change of time about the output .

4.4. Fractal Dimension

Fractal dimension can be used as another criterion to judge whether the system is in a chaotic state or not. There are many specific definitions of fractal dimension, but none of them can be taken as the universal one. According to [22], the following definition of fractal dimension is adopted in this paper:

where are the Lyapunov exponents and is the largest positive integer for which satisfies and . If , , then . If , , then .

The Lyapunov exponents of system (6) are and for and . System (6) is in a chaotic state because the largest Lyapunov exponent is positive. Fractal dimension demonstrates that the chaotic motion has self-similar structure, which is an important difference between the chaotic motion and the stochastic motion. The fractal dimension of system (6) is . The fractal dimension also reflects the space density of the chaotic attractor. The larger the dimension of the chaotic attractor is, the larger the occupied space is. Thus, the structure of the chaotic attractor is more compact, and the system is more complex and vice versa. The fractal dimension of the 2D discrete system (6) is more than 1.5, so the occupied space is broad, and the structure is tight, which can be seen in Figure 5.

5. Chaos Control

We can see that system (6) becomes unstable and eventually falls into chaos if the output adjustment speed parameter    exceeds a critical value. Both the two renewable resource companies will be harmed and the renewable resource market will become irregular when chaos occurs. Therefore, no one is able to make good strategies and decide reasonable output. To avert risk, it is a good idea for the duopoly renewable resource firms to stay at Nash equilibrium state.

In this part, the parameter adjustment method is used to control the chaos of system (6). The system which is under controlled is as follows:

where is the adjustment parameter and other parameters are the same as above.

It can be seen from the bifurcation diagram in Figure 22 and the corresponding Lyapunov exponents path in Figure 23 that, for and , controlled system (13) stabilized at Nash equilibrium point for . It demonstrates that chaos control of the system (6) can be realized by a small parameter adjustment.

Figure 22: Bifurcation diagram of system (13) with , and .

Figure 23: The corresponding Lyapunov exponents with , and .

When and , the stable region of expands from the original 2-cycle bifurcation point in Figure 2 to in Figure 24, which implies that once system is under control, chaos is delayed or eliminated completely.

Figure 24: Bifurcation diagram of system (13) with , and and .

The purpose of system control is to take effective measures to regulate market behaviors for avoiding the occurrence of chaos. The system can regain equilibrium state by the parameter adjustment method when the system (6) is in a chaotic state, and the bifurcation and chaos behaviors of the system can be delayed or eliminated. The parameter adjustment method can ensure that the renewable resource market develop orderly, and the renewable resource firms compete rationally and healthily.

6. Model of Evolutionary Game

Assume that all the firms in the market take the strategy of team combination. Every combination is based on the last stage, which is similar to the view of biological evolution [20]. Thus, we analyze the game process using the team combination strategy under the theory of evolutionary game. The types of free combination which can take are , , and ; the probability of free combination for is ; the probability of free combination for is . The types of free combination which can take are , , and ; the probability of free combination for is ; the probability of free combination for is . The types of free combination which can take are , , and ; the probability of free combination for is ; the probability of free combination for is . Consider

We make , , and ; thus, , , and , and the model can be written:

Solve the equation; then we can obtain the equilibrium points:

The dynamic equation can be written as follows:

where () is the adjustment speed parameter.

6.1. Lyapunov Exponents and Chaotic Attractors

As it can be seen, Figure 25(a) is the Lyapunov exponents with the varies of , when , , , , , , , and ; Figure 25(b) is the Lyapunov exponents with the varies of , when , , , , , , , and ; Figure 25(c) is the Lyapunov exponents with the varies of , when , , , , , , , and ; Figure 25(d) is the Lyapunov exponents with the varies of , when , , , , , , , and ; Figure 25(e) is the Lyapunov exponents with the varies of , when , , , , , , , and ; Figure 25(f) is the Lyapunov exponents with the varies of , when , , , , , , , and . The evolutionary game model is a special nonlinear dynamic model. In the process of evolution, when the adjust parameter is 2.01, the system enters a chaotic state. According to the numerical simulation, it shows that the state has nothing to do with the initial probability and the coefficient of the dynamic equation. Whatever initial values and parameters change into, it is at the point of 2.01 that the first largest Lyapunov index changes into zero, and the system becomes chaotic, the second and third largest lyapunov index changing into zero. Thus, the 2.01 is a key point, which is a border of a bifurcation point from the steady state to chaos in the process of evolutionary game. In order to prevent the system enters a chaotic state, it requires us to control the parameters. When the adjustment speed parameter changes too fast, the system enters a chaotic or hyperchaotic state. It is known to us that the direction of evolutionary game will be decided by the initial value. Whether the system goes into a chaotic or hyperchaotic state, it has no relation with the initial values and the direction of the evolution, but the rate of evolution. That is to say, there are three sides in the process of evolution. One part will enter a state of chaos when the adjust parameter is 2.01 and the first largest lyapunov exponent is zero; thus the system becomes chaotic; when the second largest lyapunov exponent is zero and the adjust parameter is 2.01, two positive lyapunov exponents exist and the system goes into hyperchaotic state. When the evolutionary rate increases to 2.01 for the third party, the system appears there are three positive lyapunov exponents, and system enters a state of trichaos. Figure 26(a) is the Lyapunov exponents with the varies of , when , , , , , , , and ; Figure 26(b) is the Lyapunov exponents with the varies of , when , , , , , , , and ; we obtain that when the adjust rate parameter changes, the point of the third lyapunov exponent becoming zero is still at 2.01; Figure 26(c) is the projection onto the - plane of the attractors and Figure 26(d) is the projection onto the - plane of the attractors when , , , , , , , , and .

(a) The Lyapunov exponents with the varies of , when , , , , , , , and

(b) The Lyapunov exponents with the varies of , when , , , , , , , and

(c) The Lyapunov exponents with the varies of , when , , , , , , , and

(d) The Lyapunov exponents with the varies of , when , , , , , , , and

(e) The Lyapunov exponents with the varies of , when , , , , , , , and

(f) The Lyapunov exponents with the varies of , when , , , , , , , and

Figure 25: The Lyapunov exponents with the varies of parameters () and (when the step length is 0.001.

(a) The Lyapunov exponents with the varies of , when , , , , , , , and

(b) The Lyapunov exponents with the varies of , when , , , , , , , and

(c) Projection onto the - plane of the attractors

(d) Projection onto the - plane of the attractors

Figure 26: The Lyapunov exponents with the varies of parameters () and (when the step length is 0.001.

7. Conclusions

A new dynamics of a nonlinear duopoly game with team players in exploitation of a renewable resource is built up in this paper. The stability of the Nash equilibrium point and the bifurcation and chaotic behaviors of the dynamic system are studied. It is found that bifurcation and chaos occur with the increasing of the output adjustment speed parameters    and the change of weight parameter . The weight parameter not only affects the stability of system (6), but also changes the Nash equilibrium point of the system. The oligopoly market may become unstable and even fall into chaos while the output adjustment speed parameters are out of the stable region. Furthermore, the chaos can be delayed or eliminated by the parameter adjustment method.

A model of evolutionary game is constructed in this paper, and the dynamic behaviors are discussed. We can obtain some useful results indicating that evolutionary game model is a special dynamic system, because it is at the point of 2.01 that every lyapunov index changes into chaos, no matter how the parameters change. Thus, it will be helpful for us to adopt some measures in the process of game.

It has a very theoretical and practical significance to research the complexity of new nonlinear dynamical system. As the traditional nonrenewable energy is being depleted, the development of renewable resource is an inevitable choice. This paper also shows a guidance for the renewable resource firms to make strategies of their output and exploit the renewable resource and is helpful for the government to formulate relevant policies to manage the renewable resource market.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors would like to seize this opportunity to thank the National Nature Science Foundation of China (Grant no. 61273231). This work was supported by Doctoral Fund of Ministry of Education of China (Grant no. 20130032110073).

References

T. Puu, “The chaotic duopolists revisited,” Journal of Economic Behavior and Organization, vol. 33, no. 3-4, pp. 385–394, 1998.

H. N. Agiza, “Explicit stability zones for Cournot game with 3 and 4 competitors,” Chaos, Solitons and Fractals, vol. 9, no. 12, pp. 1955–1966, 1998.

H. N. Agiza, “On the analysis of stability, bifurcation, chaos and chaos control of Kopel map,” Chaos, Solitons and Fractals, vol. 10, no. 11, pp. 1909–1916, 1999.

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.

T. Dubiel-Teleszynski, “Nonlinear dynamics in a heterogeneous duopoly game with adjusting players and diseconomies of scale,” Communications in Nonlinear Science and Numerical Simulation, vol. 16, no. 1, pp. 296–308, 2011.

D. Dragone, L. Lambertini, G. Leitmann, and A. Palestini, “A stochastic optimal control model of pollution abatement,” Nonlinear Dynamics and Systems Theory, vol. 10, no. 2, pp. 117–124, 2010.

A. A. Elsadany, “Dynamics of a delayed duopoly game with bounded rationality,” Mathematical and Computer Modelling, vol. 52, no. 9-10, pp. 1479–1489, 2010.

J. Ma and X. Pu, “Complex dynamics in nonlinear triopoly market with different expectations,” Discrete Dynamics in Nature and Society, vol. 2011, Article ID 902014, 12 pages, 2011.

J. Ma and Y. Yang, “Hyperchaos numerical simulation and control in a 4D hyperchaotic system,” Discrete Dynamics in Nature and Society, vol. 2013, Article ID 980578, 16 pages, 2013.

X. Xu, Z. Chen, G. Si, X. Hu, Y. Jiang, and X. Xu, “The chaotic dynamics of the social behavior selection networks in crowd simulation,” Nonlinear Dynamics, vol. 64, no. 1-2, pp. 117–126, 2011.

J. Ding, Q. Mei, and H. Yao, “Dynamics and adaptive control of a duopoly advertising model based on heterogeneous expectations,” Nonlinear Dynamics, vol. 67, no. 1, pp. 129–138, 2012.

Y. Liu, X. Chen, L. Zhang, F. Tao, and L. Wang, “Does migration cost influence cooperation among success-driven individuals,” Chaos Solitons Fractals, vol. 45, pp. 1301–1308, 2012.

H. Li, Q. Dai, H. Cheng, and J. Yang, “Cooperation in an evolutionary prisoner's dilemma game with probabilistic strategies,” Chaos Solitons Fractals, vol. 45, pp. 1397–1403, 2012.

E. Ahmed and A. S. Hegazi, “On dynamical multi-team and signaling games,” Applied Mathematics and Computation, vol. 172, no. 1, pp. 524–530, 2006.

E. Ahmed, M. F. Elettreby, and A. S. Hegazi, “On Puu's incomplete information formulation for the standard and multi-team Bertrand game,” Chaos, Solitons & Fractals, vol. 30, no. 5, pp. 1180–1184, 2006.

E. Ahmed, M. F. Elettreby, and A. S. Hegazi, “On quantum team games,” International Journal of Theoretical Physics, vol. 45, no. 5, pp. 907–913, 2006.

M. F. Elettreby and S. Z. Hassan, “Dynamical multi-team Cournot game,” Chaos, Solitons & Fractals, vol. 27, no. 3, pp. 666–672, 2006.

M. F. Elettreby and M. Mansour, “On Cournot dynamic multi-team game using incomplete information dynamical system,” Applied Mathematics and Computation, vol. 218, no. 21, pp. 10691–10696, 2012.

S. S. Asker, “On dynamical multi-team Cournot game in exploitation of a renewable resource,” Chaos, Solitons and Fractals, vol. 32, no. 1, pp. 264–268, 2007.

C. Veller and V. Rajpaul, “Purely competitive evolutionary dynamics for games,” Physical Review E, vol. 86, no. 4, Article ID 041907, 2012.

L. Edelstein-Kashet, Mathematical Models in Biology, Random House, New York, NY, USA, 1992.

J. L. Kaplan and J. A. Yorke, “Preturbulence: a regime observed in a fluid flow model of Lorenz,” Communications in Mathematical Physics, vol. 67, no. 2, pp. 93–108, 1979.

H. Wang and J. Ma, “Complexity analysis of a Cournot-Bertrand duopoly game model with limited information,” Discrete Dynamics in Nature and Society, vol. 2013, Article ID 287371, 6 pages, 2013.

J. Ma and H. Tu, “Complexity of a duopoly game in the electricity market with delayed bounded rationality,” Discrete Dynamics in Nature and Society, vol. 2012, Article ID 698270, 13 pages, 2012.