ORIGINAL RESEARCH article
Extended Dynamic Oligopolies with Flexible Workforce and Isoelastic Price Function
- 1Department of Economics, Chuo University, Tokyo, Japan
- 2Department of Psychology, Center for Cognitive Science, University of Torino, Torino, Italy
- 3Department of Applied Mathematics, University of Pécs, Pécs, Hungary
Single-product oligopolies without product differentiation are examined with linear production, production adjustment, flexible workforce and investment costs. The price function is assumed to be hyperbolic which makes the non-linearity of the model much stronger than in the case of linear price function examined earlier in the literature. The best responses of the firms are determined which are not monotonic in contrast to the linear case. The set of all steady states is then characterized and in the case of a duopoly it is illustrated. The asymptotical behavior of the steady states is examined by using simulation. We analyze the effects of the different types of costs on the industry dynamics and compare them to the prediction by the well known model with hyperbolic price function and no product adjustment and investments costs.
Oligopoly theory and its applications became one of the central issues in the literature of mathematical economics since the pioneering work of Cournot . As a non-cooperative game it models the competition of firms producing similar products or offering similar services. Several variants and extensions of the classical model were introduced and examined including models with and without product differentiation, multi-product oligopolies, rent-seeking games and labor managed firms among others. The existence and uniqueness of the equilibrium was the main focus of research in early stages and later the focus of studies turned to the dynamic extensions of these model variants. The earlier results up to the mid 70s were summarized in Okuguchi  including some of his fundamental contributions. With linear price and cost functions the dynamic models with both gradient adjustments and partial adjustments toward best responses were linear, the asymptotic behavior of which were relatively simple since local asymptotical stability implied global stability. The multiproduct extensions of these models were discussed in Okuguchi and Szidarovszky . More recently non-linear models became the main research focus. There are several ways to introduce non-linearities into oligopoly models. Keeping the linearity of the price and production cost functions, production adjustment costs were introduced and their effect on the asymptotic properties of the equilibrium were examined by Howroyd and Rickard , Macleod , Reynolds [6, 7], Szidarovszky and Yen  among others.
A complete equilibrium analysis was offered in Zhao and Szidarovszky  and Matsumoto et al. . Non-linearities were introduced also by considering cartelizing groups and antitrust thresholds in Matsumoto et al. [11–13], by introducing contingent workforce and investment costs in Merlone and Szidarovszky  and Matsumoto et al.  and by adding adjustment constraints in Burr et al. . The introduction of non-linear price and/or cost functions also leads to non-linear dynamics. Hyperbolic price functions result in interesting dynamic properties. Such oligopoly models are equivalent to rent-seeking games [16–18] as well as to market-share attraction models [19, 20]. A comprehensive summary of different versions of non-linear oligopolies and their asymptotic analysis are offered in Bischi et al. .
The flexibility of workforce is well known to be an important aspect in terms of manufacturers competitivity . More recently, the role of flexibility has been examined as it concerns recession and possible recover . In this paper we reconsider the model of Matsumoto et al.  with keeping linear production, flexible workforce and additional adjustment costs but introducing hyperbolic price function which makes the non-linearity of the model much stronger leading to more interesting dynamic properties. Hyperbolic price functions were introduced into duopolies by Puu  based on general Cobb-Douglas type utility functions of the consumers . They also have the interesting property that the consumers always spend a constant sum on the goods, regardless of price. The choice of linear cost function serves mathematical convenience for the possibility of deriving simple analytic expressions in the different segments of the profit function. More complex cost functions would result in implicit forms of the best responses making the analysis much more difficult if not impossible. Nevertheless, by comparing the dynamics of our model to those in the original model presented in Puu  we can better understand what are the consequences of flexible workforce in the industry.
This paper develops as follows. The mathematical model is introduced and the best responses of the firms are determined in Section 2. The set of all steady states is characterized in Section 3, and the asymptotic behavior of the steady states is examined in Section 4 by using simulation. The last Section 5 concludes the paper with future research directions.
2. The Mathematical Model
An n-firm single-product oligopoly without product differentiation is considered with isoelastic price function, where s is the output of the industry. Let xk denote the output of firm k, then . It is assumed that the firms have linear cost functions, Ck(xk) = ck + dkxk with ck, dk > 0. In addition to these production costs we consider the following cost types. Hiring new workers requires their training and possibly higher wages. Layoff of workers costs the company the unemployment insurance and usually severance pays. The decrease in production levels requires layoffs and any increase is possible only by increasing the workforce. It is assumed that the additional cost of production level changes linearly depend on the levels of decrease or increase in production. This can be modeled as
as the additional cost in time period t. Here xk is the production level of the firm in time period t as decision variable and not the actual value. Increasing the capacity limit beyond the already built up level
also has the investment cost:
Therefore, the profit of firm k can be given as follows,
where is the output of the rest of the industry and Lk is the maximum possible capacity limit that cannot be increased further. There is a difference between Xk(t − 1) and Lk. While Lk is the maximum possible production level that firm k is able to produce, Xk(t − 1) is the built up capacity during the previous time periods. In other words, we consider the cost of workforce flexibility (parameters γk and δk), the cost of adding new machinery (parameters αk) and, finally some structural limits which bound the firms capacity. While we assume that workforce flexibility and new machinery costs are piece-wise linear, as it concerns the structural limits they would need more time to be overcome. For the sake of simplicity let φ1, φ2, and φ3 denote these functions.
so for all feasible xk, , and for all l,
implying that Πk is a strictly concave, piece-wise differentiable, continuous function.
In order to find out the shape of the profit function Πk and determine the best responses Rk we have to consider the following cases.
(i) occurs when , which can be rewritten as
where we assume the reasonable condition dk > δk. For the sake of simplified notation let and . Notice that if sk = 0, then in the first terms of Πk, xk cancels out and xk = 0 is the best choice. However, at sk = xk = 0, Πk is undefined, so this is only a fictitious best response. So this case occurs when with best response Rk = 0. At there is no best response.
(ii) and occur when
The second condition is a quadratic inequality in sk:
The discriminant of the left hand side is
If this is non-positive, then (8) holds for all sk, in which case define and . If , then there are two positive roots . We can prove that . This inequality has the form:
which can be rewritten as
This inequality is obviously satisfied. So this case occurs when
and Rk is the stationary point in interval (0, xk (t − 1)):
(iii) and is the case when and
Notice that if , then as well and if , then as well. In the case when , the left hand side of (11) has two positive roots , and since the left hand side of Inequality (11) is larger than that of (8), and . Otherwise Equation (11) holds for all sk and so we can select and . This case occurs when
and the best response is xk (t − 1).
(iv) and occur when and
The discriminant of the left hand side is
and notice that implies that and if then as well. If then we can select and . Otherwise Expression (13) has two positive roots and since the left hand side of (13) is larger than that of (11), and . Clearly this is the case when
and the best response is the stationary point between xk (t − 1) and Xk (t − 1):
(v) and have two conditions:
The discriminant is
and notice that implies that furthermore implies that . If then we select and , and if , then Expression (16) has two positive roots , where and . This case occurs when
and in this case the best response is Xk (t − 1).
(vi) and is the case when and
The discriminant is
and similarly to the other cases implies that , and implies that . In case of we may chose and . Otherwise Expression (18) has two positive roots , where and . This case occurs when
and the best response is the stationary point between Xk (t − 1) and Lk:
(vii) And finally, the case of occurs when and the best response is Rk = Lk.
For easier understanding of the different cases notice that , so if any one of these quantities in non-positive, then the same holds for all others with larger subscripts.
The computation of Rk can be done in the following algorithm:
• Step 1. Set
Compute values for l = 1, 2, 3, 4, 5. If all are non-positive then Rk is determined from Figure 2, otherwise go to next step.
• Step 2. Let and be the smaller and larger ones of the roots
If , then Rk is determined from Figure 3, otherwise go to next step.
• Step 3. Let and be the smaller and larger ones of the roots
If , then Rk is determined from Figure 4, otherwise go to next step.
• Step 4. Let and be the smaller and larger ones of the roots
If , then Rk is determined from Figure 5, otherwise go to next step.
• Step 5. Let and be the smaller and larger ones of the roots
If , then Rk is determined from Figure 6, otherwise go to next step.
• Step 6. Let and be the smaller and larger ones of the roots
Rk is determined from Figure 7.
3. Steady States Analysis
By denoting the best response function of firm k by Rk(sk, xk, Xk), the dynamic model with positive adjustment toward best responses can be written as
for k = 1, 2, …, n. Clearly a vector is a steady state if and only if for all k,
In determining the set of all steady states we have the following possibilities:
(c) is interior if
Notice that in case (b), , so segment of φ3 is eliminated. The condition of case (b) can be rewritten as
and that of case (c) is equivalent to the following:
So feasible solution of Inequality (30) for exists if and only if the right hand side is positive, which occurs if , and if the left hand side is below Lk, which is the case when (29) is violated.
4. Simulation Study
The model introduced and analyzed in the previous sections is a clear generalization of the duopoly model of Puu [24, 26, Chapter 7]. His special model can be obtained by selecting α1 = α2 = γ1 = γ2 = δ1 = δ2 = 0 and sufficiently large values of L1 and L2. In this particular case, the only fixed point, except the origin, is, of course, the Cournot equilibrium point:
In Puu  it is proved that whenever one of the ratios of the marginal costs of the duopolists falls outside the interval , the Cournot point is not stable. Furthermore, the complexity of the dynamics is illustrated considering a bifurcation diagram of firms' output vs. the marginal cost ratio. When we set α1 = α2 = γ1 = γ2 = δ1 = δ2 = 0 and sufficiently large L1 and L2, we obtain the bifurcation diagram reported in Figure 10, which is identical to the one presented in Puu [26, p. 271]. This Figure has been obtained with A = 1, c1 = c2 = 0, d1 = 1 and initial condition x1(0) = x2(0) = 0.1 when d2 varies in the interval [5.75, 6.25].
The model we present here is much more complex. To understand the effects of the different costs and production constraints we introduced, we analyze each of them separately and compare the dynamics to the one originally presented in Puu . These analyses are reported in Figures 11–14. Each of these figures consists of a bifurcation diagram in the parameter plane (d2/d1, α = α1 = α2) in which the regions of different periodicity are represented by different colors. Bifurcation diagrams are also shown in which amplitudes are plotted vs. marginal cost ratio as in the original figure in Puu [26, p. 271] which is here reproduced as Figure 10. The two bifurcation diagrams with marginal cost ratio as bifurcation variable help to understand how the different cost coefficients affect the dynamics. Therefore, we will next examine how the dynamic properties of the dynamics depend on the different values of the different cost coefficients. The comparison of these bifurcation diagrams helps to understand the respective roles of isoelastic demand function, workforce flexibility costs and structural limits.
Figure 11. Investment cost. (A) Bifurcation diagram in the parameter plane (d2/d1, α = α1 = α2); the regions of different periodicity are represented by different colors. The horizontal yellow lines represent the line on which, in (B), the bifurcation diagram is shown with α = α1 = α2 = 0.20 and in (C), the bifurcation diagram is shown with α = α1 = α2 = 0.48.
Let us start our analysis by considering the effects of the investment costs αk. As the investment costs are small (Figure 11B) the bifurcation diagram is identical to the case when investment costs are zero. In this case the trajectories after a transient of 5000 iterations are also identical. By contrast, when investments costs are larger (α = α1 = α2 = 0.48) the trajectory becomes less complex than in the case of smaller costs (α = α1 = α2 = 0.2). This can be explained because larger investment costs dampen firms' reply.
The effects of the hiring costs γk are much more pronounced. In fact even for small values of hiring costs the dynamics is much simpler as it can be seen in Figures 12B,C. For even larger values of the hiring costs we do not have cycles as it can be seen from the top of Figure 12A. In this case also, hiring costs dampen firms' reply.
Figure 12. Hiring cost. (A) Bifurcation diagram in the parameter plane (d2/d1, γ = γ1 = γ2); the regions of different periodicity are represented by different colors. The horizontal yellow line represents the line on which, in (B), the bifurcation diagram is shown with γ = γ1 = γ2 = 0.05, and in (C), the bifurcation diagram is shown with γ = γ1 = γ2 = 0.4.
When considering layoff costs δk, the effects are similar to those of investment costs. For small values of layoff costs the dynamics is similar but not identical to the one with no costs. In fact, although Figures 10, 13C look similar, an inspection of the trajectories after discarding 5000 transient shows some differences. Also in this case, for even larger values of the hiring costs we do not have cycles as it can be seen from the top of Figure 13A. Large enough layoff costs dampen the dynamics.
Figure 13. Layoff cost. (A) Bifurcation diagram in the parameter plane (d2/d1, δ = δ1 = δ2); the regions of different periodicity are represented by different colors. The horizontal yellow lines represent the line on which, in (B), the bifurcation diagram is shown with δ = δ1 = δ2 = 0.1 and in (C), the bifurcation diagram is shown with δ = δ1 = δ2 = 0.2.
Finally, when considering capacity limits Lk we can see that, unless they are influencing firms' response, they have no effect on the dynamics. As a matter of fact, Figures 10, 14C are identical as they are the respective trajectories after a transient of 5000 iterations. In this case large production limits do not have any effect on the dynamics.
Figure 14. Capacity limit. (A) Bifurcation diagram in the parameter plane (d2/d1, L = L1 = L2); the regions of different periodicity are represented by different colors. The horizontal yellow line represents the line on which, in (B), the bifurcation diagram is shown with L = L1 = L2 = 0.15, and in (C), the bifurcation diagram is shown with L = L1 = L2 = 0.2.
Non-linear single product oligopolies without product differentiation were introduced and examined where the production cost was linear, and the piecewise linear production adjustment and investment costs made the model non-linear. The non-linearity of the model became even stronger by assuming hyperbolic price function. These models are equivalent with rent-seeking and market-share attraction games as well. The profit functions of the firms are continuous, piece-wise differentiable and strictly concave implying the uniqueness of the best responses. The best response functions of the firms were then determined which are not monotonic in contrast with the case of linear price functions assumed earlier in the literature. The set of all steady states were characterized and illustrated in the case of a duopoly. The asymptotical properties of the steady states were investigated by using simulation. Comparing the dynamics to the one of Puu's original model, showed the different roles these costs have on the dynamics. Although some of them seem to have little effects on the dynamics, others –such as the hiring costs– have a deep impact. In particular, as recruitment and selection cost can be staggering, see for instance Gusdorf , they should be considered in order to have a more realistic model. In our analysis these costs make the dynamics less complex than the one predicted by the theoretical model. Furthermore, the dynamics in real word seems to be less complex than the one predicted by models which do not consider these costs. A reason for this discrepancy could be that the hiring costs are high and therefore the dynamics is less complex than predicted.
It will be interesting to consider the case of a generic isoelastic function. Also, further non-linearities can be introduced into the models by assuming non-linear cost functions and different types of the price functions. We will elaborate on these ideas in our next research project.
UM: Original Idea and Simulations. FS: Mathematical Derivation. AM: General Comments.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
AM has received research grants from the MEXT-Supported Program for the Strategic Research Foundation at Private Universities 2013–2017 and Chuo University (Joint Research Grant). UM has developed this work in the framework of the research project on “Dynamic Models for behavioural economics” financed by DESP-University of Urbino.
10. Matsumoto A, Merlone U, Szidarovszky F. Dynamic oligopoly models with production adjustment and investment costs. In: Matsumoto A, Szidarovszky F, Asada T, editors. Essays in Economic Dynamics. Berlin: Springer-Verlag (2016). p. 99–109. doi: 10.1007/978-981-10-1521-2_6
23. Nash BJ, Romero J. Flexible workforce: the role of temporary workers in recession and recovery. Region Focus. (2011) First Quarter:21–3, 38. Available online at: https://www.richmondfed.org/~/media/richmondfedorg/publications/research/region_focus/2011/q1/pdf/feature1.pdf
25. Agliari A, Puu T. A Cournot duopoly with bounded inverse demand function. In: Puu T, Sushko I, editors. Oligopoly Dynamics: Models and Tools. Heidelberg: Springer-Verlag (2002). p. 171–94. doi: 10.1007/978-3-540-24792-0_7
Keywords: oligopolies, repeated games, complex dynamics, contingent workforce, investment cost
JEL code: C72,C73.
Citation: Matsumoto A, Merlone U and Szidarovszky F (2016) Extended Dynamic Oligopolies with Flexible Workforce and Isoelastic Price Function. Front. Appl. Math. Stat. 2:19. doi: 10.3389/fams.2016.00019
Received: 13 July 2016; Accepted: 06 October 2016;
Published: 04 November 2016.
Edited by:Laura Gardini, University of Urbino, Italy
Copyright © 2016 Matsumoto, Merlone and Szidarovszky. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Ugo Merlone, firstname.lastname@example.org