Research Article  Open Access
NeimarkSacker Bifurcation in DemandInventory Model with StockLevelDependent Demand
Abstract
An analysis of dynamics of demandinventory model with stockleveldependent demand formulated as a threedimensional system of difference equations with four parameters is considered. By reducing the model to the planar system with five parameters, an analysis of oneparameter bifurcation of equilibrium points is presented. By the analytical method, we prove that nondegeneracy conditions for the existence of NeimarkSacker bifurcation for the planar system are fulfilled. To check the sign of the first Lyapunov coefficient of NeimarkSacker bifurcation, we use numerical simulations. We give phase portraits of the planar system to confirm the previous analytical results and show new interesting complex dynamical behaviours emerging in it. Finally, the economical interpretation of the system is given.
1. Introduction
Great economic development after the second world war has released a need of development of mathematical methods to support optimization of economic and business processes. Material flow, production, and inventory are aspects of a business, which, to make it profitable, need to be optimized. Therefore, many models of supply chain were created in the mid20th century. To mention the most noticeable, we ought to list works of Wagner and Whitin [1], Brown [2], and Holt et al. (HMMS model) [3]. They have laid the foundation for supply chain modelling. Those models, despite being relatively simple, have become an inspiration for contemporary researchers: to redefine models in order to fit to the current challenges and to analyse them using available computation power, for example, [4–6].
Prediction of future demand and inventory is an important aspect of running and managing manufacturing or trade company. Methods supporting those tasks have been developed by economists already in the mid20th century; nonetheless, they are still being improved, as economy still changes and creates new challenges. One of those methods is modelling of economic phenomena using mathematical formulations. The topic of demand and inventory needs a contextual approach, since many factors can influence it and different views may be needed. Therefore, the models are created with the usage of different mathematical tools. We can mention here recent works related to demand and inventory that investigate and describe specific economic cases: Chen and Hu in [7] consider an inventory and pricing model and develop polynomial time algorithms to maximize the total profit; Mondal et al. in [8] consider generic algorithm for the case of inventory of a deteriorating item; Qi et al. in [4] consider supply chain that experiences a disruption in demand during the planning horizon.
Ma and Feng in [6] proposed the dynamical model of demand and inventory with mechanism of demand stimulation and inventory limitation. Nowakowska in [9] showed application of such a model to electrical energy market. Hachuła et al. analysed the stability of equilibrium points in [10, 11].
The model describes demand and inventory of a product at one echelon of supply chain at the retailer. The considered supply chain consists of three echelons: manufacturer, retailer, and customers. The following rules are applied to the model: customers buy a good from a retailer; a retailer orders a product in the forecasted amount and forecast is prepared using single exponential smoothing Brown model [12]; a manufacturer produces and delivers exactly the ordered amount and the production capacity is unlimited; customers’ demand depends on a retail price, which can be changed by a discount; price cannot be arbitrarily changed but the retailer can offer a discount depending on stock volume.
The model takes a form of the following system of difference equations:where indicates instance of time; is a stock volume at ; is a real number; interpretation of nonnegative values is obvious and the negative value of stock informs what amount is missing to fulfill demand and needs to be overcompensated in the next periods; is a demand volume at , which assumes constant demand elasticity (see [13, p. 280]), ; is a forecast of demand at and order placed by a retailer at a manufacturer; moreover, by assumption of unlimited capacity, it is also the delivered quantity at , ; is a parameter for discount steering, ; is a parameter for defining the target stock, ; is a price elasticity coefficient that regulates dependence between price, discount, and demand, ; is a forecast smoothing coefficient of Brown model, .
The first equation describes the influence of relation between current stock and target stock on future demand , in a way that when stock is high, the retailer offers a discount to encourage customers to buy a product. The second equation lets us calculate future stock as a sum of current stock and delivery of amount of forecasted demand diminished by actual sales . The third equation describes a method of demand forecasting. Future demand is forecasted using single exponential smoothing by Brown [2]. Smoothing coefficient split importance between actual demand (sales) and forecasted demand .
We analyse properties of a given model (1) on the ground of theory of difference equations. Threedimensional difference systems were studied by many authors, for example, in [14–17]. Under some assumptions, such systems can be rewritten as thirdorder difference equations. Asymptotic properties of these types of equations were investigated in [18–23]. For the background of difference equations theory, see monographs [24, 25]. Firstly, we construct a positive invariant set of the threedimensional system. We notice that the given system has an invariant plane and can be reduced to the family of the planar systems, which we construct positive invariant sets and provide the equilibrium points for. We recall the properties of equilibrium points obtained in [11]. The aim of this paper is to analyse the dynamics of planar systems on the ground of bifurcation theory (see [26]) and numerical simulations. It was shown that nondegeneracy conditions for NeimarkSacker bifurcation are fulfilled and existence of transcritical bifurcation is suggested. The analytical results are confirmed by numerical simulations. Finally, the economical interpretation of the system is given.
2. Main Results
For the sake of convenience, let us write the model as a system of three firstorder difference equations with , , and replacing , , and , respectively, in system (1),with , , or as with a mapping:It is easy to see that for and for are equilibrium points of the system .
Let us shortly elaborate on the meaning of those equilibrium points. The goal of the retailer is to reach target stock . If additionally actual demand and forecasted demand are equal and greater than zero, the retailer’s stock level does not change in the next period, which means that system (2) is in the equilibrium point . On the other hand, when , the retailer’s stock level cannot be changed. This means that the retailer remains in equilibrium point .
We recall that is the positive invariant set of the system if and .
Proposition 1 (see [11]). Here, is the positive invariant set of the system , where is given by (3).
One can observe that the second and the third equation of system (2) are linear ones. Therefore, we treat them with some linear transformations in order to obtain a more convenient form for analysis: both sides of the second equation multiplied by and the second and the third equation added. As an output, we getwhich enables us to state the following.
Corollary 2. System (2) has the invariant plane for any . Hence, system (2) can be reduced to a family of planar systems:with , where , .
Proof. From (5), we get the firstorder recursion of the form which implies that . Constraint of the expression is a direct consequence of the positive invariant set construction.
Let us later write our systems (6) in the form of a family of mappings:for . Taking into consideration economic interpretation, , where .
Proposition 3. Let , , and . The set is the positive invariant set of the system , where is given by (7).
Proof. Let . For any , , we have and
Mapping (7) may have at most two equilibrium points:(i)(ii) iff
Stability of the equilibrium points of (2) has been examined in [11]. Based on the findings of [11], one can formulate the conclusions about the existence and stability of equilibrium points of planar systems (6), which are presented in Table 1. Now, using these results, we analyse the dynamics of the planar systems using bifurcation theory. We are especially interested in proving the existence of NeimarkSacker bifurcation.

Graphs illustrating the cases in Table 1 are presented in Figures 1–4.
The planar system is dependent on five parameters , , , , and . The most interesting case in the bifurcation analysis for the planar system (6) is dependence on parameter with fixed , , which now we will carry on.
From the conditions in Table 1, in order to examine bifurcations on parameter , we have to consider two cases dependent on .
Case 1 (). For , the equilibrium point is the only equilibrium point of the system which is stable. For , the considered system possesses two equilibrium points: unstable and stable . This suggests that is the bifurcation parameter for the transcritical bifurcation.
Case 2 (). For , probably the transcritical bifurcation occurs like in the previous case. For , the equilibrium point is stable and loses its stability for . We prove that, for , the NeimarkSacker bifurcation occurs. We use classical analytical results from [26] and numerical simulations. For the sake of convenience, let us quote from [26] Theorem 4.5 as Lemma 4 and Theorem 4.6 as Lemma 5.
Lemma 4 (Theorem 4.5, [26]). Suppose that a twodimensional discretetime system , , , with smooth , has, for all sufficiently small , the fixed point with multipliers , where , . Let the following conditions be satisfied:
(C.1). (C.2) for . Then, there are smooth invertible coordinate and parameter changes transforming the system into with and , where is given by formula (4.21) in [26].
Lemma 5 ((Theorem 4.6, [26]) (generic NeimarkSacker bifurcation)). For any generic twodimensional oneparameter system , having at the fixed point with complex multipliers , there is a neighbourhood of in which a unique closed invariant curve bifurcates from as passes through zero.
Remark 6 (see [26]). The genericity conditions assumed in the theorem are the transversality condition (C.1) and the nondegeneracy condition (C.2) from Lemma 4 and the additional nondegeneracy condition: (C.3) .
We start with translation mapping to the mapping which possesses fixed point in for sufficiently small . We consider the mapping with parameter :The Jacobi matrix of at has the form with the complex eigenvalues for sufficiently small of the formand modulus Notice that , , which means that condition (C.1) of Lemma 4 is satisfied. Now, we check that condition (C.2) (called no strong resonance) is satisfied. We have to prove that , where and is such that . We get for because are solutions to equation . Fulfilling conditions (C.1) and (C.2) from Lemma 4 means that nondegeneracy conditions are satisfied. To get the existence of NeimarkSacker bifurcation, we have to check the last condition (C.3) of the form , named the first Lyapunov coefficient, given in Lemma 5. If , we get the supercritical NeimarkSacker bifurcation, and if we get the subcritical NeimarkSacker bifurcation. Due to complicated conditions for any parameters , , and , we check the sign of the first Lyapunov coefficient only for specific , , , and in Section 3.
Remark 7. In the analogous way, we can prove the existence of the NeimarkSacker oneparameter bifurcation of the planar system for any other parameters , , , and with the bifurcation parameter which fulfills the equality .
3. Numerical Simulations
In this section, by using numeral simulation, we give the bifurcation diagrams and phase portraits of system (6) to confirm the previous analytical results and show new interesting complex dynamical behaviours emerging for system (6).
Now, let us give illustration to the cases provided in Table 1. Graphs in Figures 1–4 depict behaviour of (6) depending on the parameter . The other parameters remain unchanged at , , and , and, for C1, C2, C4, and C5, ; for (C3), . The graphs clearly confirm the analytical results.
Interesting behaviour can be observed at border point depicted in Figures 12 in Graphs (C2ab). At the first glance, in Graph (C2a), one can assume stability of the equilibrium point , which would be a contradiction of the analytical result. However, looking from a closer perspective in Graph (C2b) one can see that the orbits leave the surroundings of the equilibrium point, herewith confirming the instability proved with analytical methods.
Let us have again a look at the border case of depicted in Figure 5. In each of the graphs, the same parameters are applied in order to observe the behaviour of the system for varying initial conditions of and (, , , , and ). Each time, we set the initial conditions closer and closer to the equilibrium point . Appearance of this kind of circle suggests supercritical bifurcation of NeimarkSacker type. Indeed, negative sign of the first Lyapunov coefficient obtained numerically confirms that suggestion.
(a)
(b)
(c)
(d)
4. Economical Interpretation
System (6) depends on five parameters: , , , , and . Each of the parameters carries the economical meaning described in Section 1. In a business practice, the parameters and are set in advance. Nonetheless, as it is shown in Table 1, changes of parameter cause the change of behaviour of the equilibrium points. Therefore, we present bifurcation diagrams of and as changes in Figures 6 and 7.
Parameter , which represents price elasticity, is defined depending on the types of a product and a market and often is determined experimentally and must be adjusted frequently. Therefore, it is interesting to analyse behaviour of system (6) at varying values of this parameter. For that purpose, let us illustrate changes of the dynamics of the equilibrium points in the bifurcation diagrams (see Figures 8 and 9). Model (6) includes steering mechanism, which allows the user to influence the behaviour of the system. With the parameter , one can decide on the share of influence between the forecast and actuals of sales. The less is, the greater the impact of actuals is, and vice versa. Figure 11 presents bifurcation diagram for changing parameter . Similarly, in Figure 10, we present bifurcation for varying . is a parameter to control the influence of the discrepancy between the actual stock and the target one .
(a)
(b)
(a)
(b)
Business Example. Let us consider a product that follows demand()price() dependence with price elasticity . The target is set on stock volume and assumed demand at least . Actual demand is insufficient , and hence the actual stock exceeds the target . Smoothing coefficient , which means calculation of the amount to be delivered takes into account actuals and forecast with equal wages. Parameter . One can easily calculate . From Table 1, it is observed that the product satisfies condition (C4), and hence the system has two equilibrium points and , unstable and stable, respectively.
Additional Points
Numerical analysis and graph plotting have been performed using Matcontm (according to [27]) and Matlab R2016a. The bifurcation diagrams depict 200 iterations starting from iteration .
Competing Interests
The authors declare that there are no competing interests regarding the publication of this paper.
References
 H. M. Wagner and T. M. Whitin, “Dynamic version of the economic lot size model,” Management Science, vol. 5, no. 1, pp. 89–96, 1958. View at: Google Scholar  MathSciNet
 R. G. Brown, Statistical Forecasting for Inventory Control, McGrawHill Book Co, New York, NY, USA, 1959.
 C. C. Holt, F. Modigliani, J. F. Muth, and H. A. Simon, Planning Production, Inventories, and Work Forces, PrenticeHall, 1960.
 X. Qi, J. F. Bard, and G. Yu, “Supply chain coordination with demand disruptions,” Omega, vol. 32, no. 4, pp. 301–312, 2004. View at: Publisher Site  Google Scholar
 S. R. T. Kumara, P. Ranjan, A. Surana, and V. Narayanan, “Decision making in logistics: a chaos theory based analysis,” CIRP Annals—Manufacturing Technology, vol. 52, no. 1, pp. 381–384, 2003. View at: Publisher Site  Google Scholar
 J. Ma and Y. Feng, “The study of the chaotic behavior in retailer's demand model,” Discrete Dynamics in Nature and Society, vol. 2008, Article ID 792031, 12 pages, 2008. View at: Publisher Site  Google Scholar  MathSciNet
 X. Chen and P. Hu, “Joint pricing and inventory management with deterministic demand and costly price adjustment,” Operations Research Letters, vol. 40, no. 5, pp. 385–389, 2012. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 S. K. Mondal, J. K. Dey, and M. Maiti, “A single period inventory model of a deteriorating item sold from two shops with shortage via genetic algorithm,” Yugoslav Journal of Operations Research, vol. 17, no. 1, pp. 75–94, 2007. View at: Publisher Site  Google Scholar  MathSciNet
 L. Nowakowska, “Dynamic discrete model for electricity price forecasting,” in Proceedings of the 5th International Youth Conference on Energy (IYCE '15), Pisa, Italy, May 2015. View at: Publisher Site  Google Scholar
 P. Hachuła and E. Schmeidel, “Stability analysis of demandinventory model in a certain business case,” in Proceedings of the 11th European Scientific Conference (TRANSCOM '15), University of Zilina, Žilina, Slovakia, June 2015. View at: Google Scholar
 M. NockowskaRosiak, P. Hachuła, and E. Schmeidel, “Stability of equilibrium points of demandinventory model with stockdependent demand,” Journal of Difference Equations and Applications, vol. 22, no. 10, pp. 1490–1500, 2016. View at: Publisher Site  Google Scholar
 E. S. Gardner Jr., “Exponential smoothing: the state of the artpart II,” International Journal of Forecasting, vol. 22, no. 4, pp. 637–666, 2006. View at: Publisher Site  Google Scholar
 H. R. Varian, Intermediate Microeconomics. A Modern Approach, W. W. Norton & Company, New York, NY, USA, 8th edition, 2010.
 M. Migda, E. Schmeidel, and M. Zdanowicz, “Existence of nonoscillatory solutions for system of neutral difference equations,” Applicable Analysis and Discrete Mathematics, vol. 9, no. 2, pp. 271–284, 2015. View at: Publisher Site  Google Scholar  MathSciNet
 E. Thandapani, R. Karunakaran, and I. M. Arockiasamy, “Bounded nonoscillatory solutions of neutral type difference systems,” Electronic Journal of Qualitative Theory of Differential Equations, p. 25, 2009. View at: Google Scholar  MathSciNet
 E. Thandapani and B. Ponnammal, “Oscillatory properties of solutions of threedimensional difference systems,” Mathematical and Computer Modelling, vol. 42, no. 56, pp. 641–650, 2005. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 E. Thandapani, M. Vijaya, and T. Li, “On the oscillation of third order halflinear neutral type difference equations,” Electronic Journal of Qualitative Theory of Differential Equations, vol. 76, pp. 1–13, 2011. View at: Google Scholar
 M. F. Aktaş, A. Tiryaki, and A. Zafer, “Oscillation of thirdorder nonlinear delay difference equations,” Turkish Journal of Mathematics, vol. 36, no. 3, pp. 422–436, 2012. View at: Publisher Site  Google Scholar  MathSciNet
 A. AndruchSobiło and M. Migda, “Bounded solutions of third order nonlinear difference equations,” The Rocky Mountain Journal of Mathematics, vol. 36, no. 1, pp. 23–34, 2006. View at: Publisher Site  Google Scholar  MathSciNet
 Z. Dolá and A. Kobza, “On thirdorder linear difference equations involving quasidifferences,” Advances in Difference Equations, vol. 2006, Article ID 065652, 2006. View at: Publisher Site  Google Scholar
 Z. Liu, M. Jia, S. M. Kang, and Y. C. Kwun, “Bounded positive solutions for a third order discrete equation,” Abstract and Applied Analysis, vol. 2012, Article ID 237036, 12 pages, 2012. View at: Google Scholar  MathSciNet
 J. Popenda and E. Schmeidel, “Nonoscillatory solutions of third order difference equations,” Portugaliae Mathematica, vol. 49, no. 2, pp. 233–239, 1992. View at: Google Scholar  Zentralblatt MATH  MathSciNet
 E. Schmeidel, “Oscillatory and asymptotically zero solutions of third order difference equations with quasidifferences,” Opuscula Mathematica, vol. 26, no. 2, pp. 361–369, 2006. View at: Google Scholar  Zentralblatt MATH  MathSciNet
 R. P. Agarwal, M. Bohner, S. R. Grace, and D. O'Regan, Discrete Oscillation Theory, Hindawi Publishing Corporation, New York, NY, USA, 2005.
 P. E. Kloeden and M. Rasmussen, Nonautonomous Dynamical Systems, vol. 176 of Mathematical Surveys and Monographs, American Mathematical Society, Providence, RI, USA, 2011. View at: Publisher Site  MathSciNet
 Y. A. Kuznetsov, Elements of Applied Bifurcation Theory, vol. 112, Springer, New York, NY, USA, 2nd edition, 1998. View at: MathSciNet
 W. Govaerts, Yu. A. Kuznetsov, R. Khoshsiar Ghaziani, and H. G. E. Meijer, Cl MatContM: A Toolbox for Continuation and Bifurcation of Cycles of Maps, Universiteit Gent and Utrecht University, 2008.
Copyright
Copyright © 2017 Piotr Hachuła 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.