<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Archiving and Interchange DTD v1.0 20120330//EN" "JATS-archivearticle1.dtd">
<article xmlns:xlink="http://www.w3.org/1999/xlink">
  <front>
    <journal-meta>
      <journal-title-group>
        <journal-title>Kyiv, Ukraine, June</journal-title>
      </journal-title-group>
    </journal-meta>
    <article-meta>
      <title-group>
        <article-title>Nonlinear Dynamic Model of a Microeconomic System with Different Reciprocity and Expectations Types of Firms: Stability and Bifurcations</article-title>
      </title-group>
      <contrib-group>
        <contrib contrib-type="author">
          <string-name>Vitaliy Kobets</string-name>
          <email>vkobets@kse.org.ua</email>
          <xref ref-type="aff" rid="aff0">0</xref>
          <xref ref-type="aff" rid="aff1">1</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Alexander Weissblut</string-name>
          <xref ref-type="aff" rid="aff0">0</xref>
          <xref ref-type="aff" rid="aff1">1</xref>
        </contrib>
        <aff id="aff0">
          <label>0</label>
          <institution>Kherson State University</institution>
          ,
          <addr-line>27, 40 rokiv Zhovtnya st., Kherson, 73000</addr-line>
          <country country="UA">Ukraine</country>
        </aff>
        <aff id="aff1">
          <label>1</label>
          <institution>ey Terms: Expectations</institution>
        </aff>
      </contrib-group>
      <pub-date>
        <year>2016</year>
      </pub-date>
      <volume>2</volume>
      <fpage>1</fpage>
      <lpage>24</lpage>
      <abstract>
        <p>This paper analyzes the dynamic interaction between selfish and reciprocity firms in the market of homogeneous product. The decisions of both types of firms in respect of their output strategies are investigated under naive, adaptive and generalized expectations. The standard postulate for competitive firms' model has been extended by the assumption that there is a share of reciprocity firms which, unlike selfish firms, maximize both private and social benefits as consumer surplus. It has been proved that the unique Nash equilibrium is stable for all affordable values of parameters in the model with adaptive expectations, and is unstable for the model with naive expectations at sufficiently large number of firms in the market. A special desktop application has been created for animation of model trajectory and demonstration of stable quantity trajectories and bifurcation diagrams of firms' output. Naive expectations of two-thirds of firms result in a state of dynamic chaos in the market leading to degeneration of the existing competition model between the two types of firms. The crucial factor which ensures the stable equilibrium in the market and the ability to predict firms' output is the adaptive approach which takes into account the adaptive expectations of firms planning their product quantity.</p>
      </abstract>
      <kwd-group>
        <kwd>microeconomic system</kwd>
        <kwd>reciprocity</kwd>
        <kwd>naïve expectations</kwd>
        <kwd>adaptive expectations</kwd>
        <kwd>consumer surplus</kwd>
        <kwd>stability</kwd>
        <kwd>bifurcation</kwd>
      </kwd-group>
    </article-meta>
  </front>
  <body>
    <sec id="sec-1">
      <title>DynamicSystem,</title>
    </sec>
    <sec id="sec-2">
      <title>DesktopApplication,</title>
    </sec>
    <sec id="sec-3">
      <title>NashEquilibrium,</title>
      <p>In recent years the researchers are renouncing the assumption of perfect rationality as
unconditional basis of economic agents’ behavior. The neoclassical ‘rational man’
does not exist in reality; economic agents act according to established rules, without
being fully informed and maximizing their own utility [1].</p>
      <p>- 503</p>
      <p>Karl Polanyi identified the alternative economic organization where social norms
are not generated by economic self-interest of the individual. This network of
reciprocal relations is based on mutual economic cooperation of efforts and resources
between the members of non-economic network, dominated by cultural norms rather
than market laws. Under reciprocity relations the exchange donor and recipient can be
transposed. So this is a symmetrical relationship of gifts exchange between members
of horizontal social networks [2].</p>
      <p>This relationship is not regulated by formal institutions but based on informal
commitments giving moral right to mutual help and reciprocal exchange on
sustainable basis in the long run period. But this is a relationship with minimal risk
for participants and the penalty of loss of social capital (reputation and trust) and
social isolation. Society supports the stability of the exchange to ensure their survival
during crises and wars. Reciprocity is not altruism which does not create reciprocal
obligations in quantitative, qualitative or time respects, just vague commitment (e.g.,
you give, if you can).</p>
      <p>Actually, reciprocity relations, commodity exchange and hierarchical
subordination exist at the same time. But it is reciprocity that underlies most
decentralized corrections of diverse shortcomings and failures of markets and firms.
These relations are long-run factors of economic efficiency; they set most social
obligations of firms towards individuals without government intervention. No society
can exist without reciprocal relationship [3].</p>
      <p>Reciprocity or social responsibility implies that the firms not only pursue their
selfish goal of increasing profits, but are also ready to sacrifice some of their own
profits for the benefit of consumers without direct compensation for it by the state [1].
Such targets can be stipulated by the firms’ desire to get stable profits in the long run
rather than maximal short-run profits. Such forward-thinking firms-reciprocators are
considered in the model of this paper. Their objective function is a weighted average
of the profits and consumer surplus of their market segment.</p>
      <p>The real economic processes make a clear demonstration that neoclassical "rational
man" is not their subject. In real economy "optimal imperfect decisions" are taken by
simple and non-expensive calculations, well adapted to frequent repetitions, to
evolution; it is more efficient for perfectly rational firm to perform multiple
experiments with quantity to estimate the demand function rather than search for
nonrecurrent, instantaneous achieving of equilibrium. New paradigm of nonlinear
economics is a mix of qualitative theory of nonlinear dynamical system, optimal
control theory, game theory, and theory of stochastic processes [4, 5].</p>
      <p>The evolutionary approach and analysis of the dynamics allow to explain why one
type of firm ousts another from the market, why sometimes the economic system is
stable, but in other cases is unstable [6, 7]. If the system has multiple equilibria, the
dynamics and evolution is the selection mechanism of best equilibrium according to
certain criteria [8]. The evolutionary process is analogous to social learning. An
example of its application is the pricing mechanisms for auctions that occur in agents
social networks, e-commerce and trade through the Internet [9, 10].</p>
      <p>The study of the evolution of the markets with the strategic interaction usually uses
the following assumptions: (a) two firms or two types of similar firms operate in the
market; (b) firms produce homogeneous goods in quantities of x1(t ) and x2 (t ) ; (c)
no firm knows the rivals’ quantities; (d) the firms seek to predict the output of the
competitors using the adaptive scheme.</p>
      <p>Planning of quantity for the following period firms resolve optimization problem:
MaxПi (xi ; xej (t 1)) , where Пi is the objective function of firm i, xej (t 1) - expected
quantity of a competitor j ( i, j  1, 2 ).</p>
      <p>Examples of bounded rationality of firms are: ignoring the impact of competitors'
actions on their own output (local monopoly approximation LMA), naive expectation
(assumption of unchangeable behavior of competitors for a long time and using
x j (t) instead xej (t 1) ). [8] Of course both decision making approaches (adaptive and
naive, bounded rational) coexist in the market with a certain probability.</p>
      <p>Analysis of nonlinear oligopoly with heterogeneous players reveals that a higher
degree of product differentiation may destabilize the Cournot-Nash equilibrium.
Authors showed that a cascade of flip bifurcation may lead to periodic cycles and
chaotic motions [11]. Stability conditions of Nash equilibrium and complex dynamics
are also studied for heterogeneous duopoly with isoelastic demand function. For such
heterogeneous players a cascade of flip bifurcation leads to periodic cycles and chaos
and the Neimark-Sacker bifurcation generates attractive invariant closed curve [12].</p>
      <p>Such scheme serves as the basis for mathematical model of this paper, which
distinguishes from the other models in that: (a) firms use more than one way of
decision-making, and combine different ones; (b) except their own selfish interests,
firms take into account social ones.</p>
      <p>The paper goal is to consider the impact of naive, adaptive and generalized
expectations of egoist and reciprocator firms on stability of equilibrium and the
conditions of transition to dynamic chaos in the numerical experiment using a
specially designed desktop application.</p>
      <p>The paper is organized as follows: part 2 describes two-dimensional market model
with naïve and adaptive expectations; part 3 is devoted to dynamics model in general
case; part 4 demonstrates C#-application model for numerical investigation; part 5
concludes.
2</p>
      <sec id="sec-3-1">
        <title>Two-Dimensional Market Model</title>
        <p>We consider the market of homogeneous product, which consists of n firms, including
k identical firms-reciprocators, each of them producing x units of product and n – k
identical selfish firms, each of them producing y units of product. Thus, the industry
quantity of the two types of firms is Q  k  x  (n  k )  y . Product price P in the
market is given by the inverse market demand function P  P(Q) 
( b  0 ).</p>
        <p>The objective function of a firm-egoist is profit  Y  (P  v)  y , where v is the
firm’s costs per unit in the market. Firm-reciprocator maximizes both its own profit
 X  (P  v)  x and consumer surplus CS of its own market segment:
b</p>
        <p>Q
  Q 
CS  k   P(q)dq  PQ  , where  is the parameter defining the segment of the
market, which the reciprocator firm believes its own and optimizes it ( 0    k ), 
is the minimal technologically possible quantity of product. Then</p>
        <p>CS  k  b  ln  Q   Qb  Q   bk  ln  Q  1  bk ln  Q  ,
where    e . The specific choice of ε does not affect the dynamics of the model
because the objective function, as any potential, is set up to an arbitrary constant
accuracy, so further we will write ε instead of ˆ . Then the objective function of
firmreciprocator is:
Consider the dynamics of this two-dimensional model in discrete time t = 0, 1,...;
where xt, yt are the outputs at time t of reciprocator and egoist firm, respectively. On
the basis of these values at time t each firm finds the optimal value for its own
production quantity in the next moment t 1 , maximizing its objective function. It
distinguishes this model, in which the firm responds to changes in output of both their
and other types of firms from traditional competition models, where one type of firm
responds to changes in other types only. So each selfish firm is looking for such value
of yt1 at which it maximizes its own profits, suggesting that SR firms and the other
n  k 1 PI firms leave their quantities unchanged:</p>
        <p>
           b 
 Y    v   yt1 . (
          <xref ref-type="bibr" rid="ref1">1</xref>
          )
        </p>
        <p>
           yt1  k  xt  (n  k 1)  yt 
Obviously, the maximum point for yt1 is found from the condition
 0 , whence:
v( yt1  kxt  (n  k 1) yt )2  b(kxt  (n  k 1) yt ) .
(
          <xref ref-type="bibr" rid="ref2">2</xref>
          )
From equation (
          <xref ref-type="bibr" rid="ref2">2</xref>
          ) we obtain the response function of the PI firm:
        </p>
        <p>
          yt1  bv (kxt  (n  k 1) yt  kxt  (n  k 1) yt . (
          <xref ref-type="bibr" rid="ref3">3</xref>
          )
Similarly, firm-reciprocator finds such value of xt1 at which the maximum value
of its objective function is:
ПX    xt1  (k 1) bxt  (n  k)  yt  xt1  vxt1   (1 )  bk ln  xt1  (k 1)xt  (n  k) yt  (
          <xref ref-type="bibr" rid="ref4">4</xref>
          )

yt1
share of profit as  
        </p>
        <p>. Then
v(xt1  (k 1)xt  (n  k ) yt )2  b((k 1)xt  (n  k ) yt )  1  bk (xt1  (k 1)xt  (n  k ) yt )
.</p>
        <p>
          Assuming z  xt1  (k 1)xt  (n  k) yt , we present (
          <xref ref-type="bibr" rid="ref5">5</xref>
          ) as:
z2  b ((k 1)xt  (n  k ) yt )  1  b z .
        </p>
        <p>
          v  vk
Hence, in view of (
          <xref ref-type="bibr" rid="ref3">3</xref>
          ), we obtain a system of dynamics equations of this model:
 b  1 1 b 2 1 1 b
xt1  v ((k 1)xt  (n  k ) yt )   2  vk   (k 1)xt  (n  k) yt  2  vk

 b
 yt1  v (kxt  (n  k 1) yt  kxt  (n  k 1) yt .
(
          <xref ref-type="bibr" rid="ref5">5</xref>
          )
(
          <xref ref-type="bibr" rid="ref6">6</xref>
          )
2.2
        </p>
        <p>
          Equilibrium Conditions for the Model with Naive Expectations
In the Nash equilibrium point xt+1=xt=x, yt+1 =yt=y for all t = 0, 1, … . Therefore, at
this point, by (
          <xref ref-type="bibr" rid="ref2">2</xref>
          ) and (
          <xref ref-type="bibr" rid="ref5">5</xref>
          ) we obtain:
kx  (n  k ) y 2  b kx  (n  k 1) y   b ((k 1)x  (n  k ) y) 
v v
1  b (kx  (n  k ) y)
 vk
(
          <xref ref-type="bibr" rid="ref7">7</xref>
          )
        </p>
        <p>From the last equation we obtain x  y 
1  n  k )  y , i.e. the response functions of both types of firms are
 k
2 1</p>
        <p>x  (1 

respectively:
x   k  (1  )  (n  k )  y , y  (2 1)  k  x</p>
        <p>
          (2 1)  k  k  (1 )  (n  k)
To calculate the coordinates of a fixed point, we substitute the expression of y
through x in the first equation (
          <xref ref-type="bibr" rid="ref7">7</xref>
          ). Thus the following proposition is proved.
Proposition 1. There is unique Nash equilibrium point in a dynamic system (
          <xref ref-type="bibr" rid="ref6">6</xref>
          ):
(
          <xref ref-type="bibr" rid="ref8">8</xref>
          )
1

(x 
n  k y) , whence it follows that
        </p>
        <p>
          k
 b  2 1   1 n  k 
x*  1   1  ,
 vn   n    k 
 (
          <xref ref-type="bibr" rid="ref9">9</xref>
          )
        </p>
        <p>b 2 1  2 1 
 y*  vn    1   n .</p>
        <p>However, is this point stable?</p>
        <p>
          Proposition 2. For any set b , v  0 and  ( 0    1 ) Nash equilibrium (
          <xref ref-type="bibr" rid="ref9">9</xref>
          ) is
unstable for sufficiently large number of firms n if k   and k  3   for any
n n 4
  0 .
        </p>
        <p>
          The destabilizing role of number of players n is well known for the evolution of
firms’ strategies in oligopoly games [8]. However, in this case, according to
calculations, point (
          <xref ref-type="bibr" rid="ref9">9</xref>
          ) is unstable even at n  5 .
        </p>
        <p>
          Proof. We show that in dynamic system (
          <xref ref-type="bibr" rid="ref6">6</xref>
          ) at Nash equilibrium point (
          <xref ref-type="bibr" rid="ref9">9</xref>
          ) modulus
of Jacobian J is greater than 1: |det J|&gt;1. This implies that at least one eigenvalue of
the Jacobian is greater than 1 in absolute value, which means instability of the fixed
point (
          <xref ref-type="bibr" rid="ref9">9</xref>
          ). Here, the Jacobian of the system (
          <xref ref-type="bibr" rid="ref6">6</xref>
          ):
J xy    xt
But for point (
          <xref ref-type="bibr" rid="ref9">9</xref>
          ) in the denominator
        </p>
        <p>((k 1)x * (n  k 1) y*) 
b
v</p>
        <p>Jxy 
J yy 
k  3   at a certain   0 , which guarantees that the factors
n 4
b
v
b
v
1 and
2 b ((k 1)xt  (n  k) yt  d 2 2 b</p>
        <p>
          v v
zero for all possible n, k, b, v&gt; 0 and α (0≤α≤1), Q.E.D.
((k 1)xt  (n  k 1) yt
1 do not equal
(
          <xref ref-type="bibr" rid="ref10">10</xref>
          )
(
          <xref ref-type="bibr" rid="ref11">11</xref>
          )
(13)
2.3
        </p>
        <p>Dynamic Model Equations with Adaptive Expectations
Since all selfish firms are assumed as identical, it is natural to suggest that they have
the same planning at moment t , so their production quantities yt1 at moment t 1
will be equal too. Given these expectations, each selfish firm is looking for such value
yt1 at which it obtains the highest profit, suggesting that production quantity of SR
firms will remain unchanged:
Obviously, the maximum point for yt1 is found from the condition  Y  0 , which
yt1
gives us:</p>
        <p> b
 Y  
 kxt  (n  k ) yt1</p>
        <p>
 v   yt1 .</p>
        <p>
kxt  (n  k) yt1 2 </p>
        <p>
          kxt .
b
v
b
v
Then kxt  (n  k ) yt1 
kxt , from here response function of PI firms is:
(n  k) yt1  b kxt  kxt . (
          <xref ref-type="bibr" rid="ref12">12</xref>
          )
        </p>
        <p>v</p>
        <p>Similarly, firm-reciprocator naturally expects that the quantity of production of all
these firms at moment t  1 would be the same. Based on this expectation, each
firmreciprocator finds the value of xt1 at which the objective function is maximal,
assuming that the output of PI firms does not change:
 X    b xt1  vxt1   (1  ) b
 kxt1  (n  k ) yt  k</p>
        <p> kxt1  (n  k ) yt  .
ln </p>
        <p>
            
In the Nash equilibrium point xt+1=xt=x, yt+1=yt=y for all t  0,1,... . Therefore, at this
point in view of (
          <xref ref-type="bibr" rid="ref11">11</xref>
          ) and (14) we get:
(kx  (n  k) y)2  b kx  b (n  k ) y  b 1  (kx  (n  k)  y) . (16)
v v v 
From the second equation we get x  1   x  n k k  (1  1   ) y , whence
response functions for selfish and reciprocator firms are, respectively:
y  n k k   ((11 ))  x x  n k k   ((11 ))  y (17)
To calculate the coordinates of the fixed point, we substitute this expression y in
terms of x at first equation (16):
 kx  (n  k ) n k k   ((11 ))  x 2  bv kx kx2    ((11 ))  1  bv kx
Hence, we obtain:
Proposition 3. There is unique Nash equilibrium point in the dynamic system (15)
with adaptive expectations:
x  vbk (  (21 ) )2 ,
 y  v(nb k )  2  ((2(1)2 ) )2 . .
        </p>
        <p>(18a)
As before, without loss of generality, let   1 , otherwise we can override the share
of profit as  </p>
        <p>
J xy    xt
,
obviously, J xx  J yy  0 where trJ  J xx  J yy  0 . Thus, to test Shur conditions it is
sufficient to establish that det J  1 . But at point (18) y* 
and therefore
b
v(n  k)
 2  ((1 ) )2
(2 )2
kJxy 
Consequently, det J  J xx  J yy  J xy  J yx  0 , Q.E.D.
function P  P(Q) </p>
        <p>The price of product P in the market is given by the inverse market demand
b</p>
        <p>( b  0 ), and the price is not less than a cent, i.e. P  0.01 .</p>
        <p>Q
Therefore, the product quantity of each firm-reciprocator is x 
. Similarly, the
Suppose that in planning under the given market model adaptive expectations are
used with probability p , naïve ones - with probability q  1  p . Then the profit
function for a typical (representative) firm-egoist has the form:</p>
        <p> b
 Y  
 yt1  kxt  p(n  k 1) yt1  q(n  k 1) yt

 v   yt1 ,

and the objective function for the representative firm-reciprocator</p>
        <p> b
 X    
 xt1  p(k 1)xt1  q(k 1)xt  (n  k) yt

xt1  vxt1   ,

(1 ) b ln( xt1  p(k 1)xt1  q(k 1)xt  (n  k) yt ) .</p>
        <p>
          k 
Obviously, for p  0 ( q  1 ) objective functions  Y and  X are consistent with the
results of naive model (
          <xref ref-type="bibr" rid="ref1">1</xref>
          ) and (
          <xref ref-type="bibr" rid="ref4">4</xref>
          ), for p  1 ( q  0 ), they are consistent with the
results of the adaptive model (
          <xref ref-type="bibr" rid="ref10">10</xref>
          ) and (13) respectively. Let us assume
z X  yt1  kxt  (n  k 1)( pyt1  qyt ) zX  xt1  (n  k ) yt  (k 1)( pxt1  qxt )
 b   b  b
In this notation  Y   z X  v   yt1 ;  X     zX xt1  vxt1   (1 ) k
Then the point yt1 of maximum profit function  X is found from the condition
ln zX .
 X  0 , here
yt1
whence
condition

(19)
(20)
(21)
z2X  bv
        </p>
        <p>(kxt  q(n  k 1) yt
yt1  (1 p(n  k 1))  b (kxt  q(n  k 1) yt  k  xt1  (n  k 1)  q  yt (22)
v
The maximum point xt1 for the objective function  X is found from the first order</p>
        <p>X  0 . Forth without loss of generality we assume here   1 ,
xt1
otherwise as above we redefine the share of profit as  </p>
        <p>z2X  bv  ((n  k )  yt  (k 1)  qxt )  1 b(1  pvk(k 1)) zX (23)
Thus, in view of (22), we obtain the dynamics model of equations system of in the
general case:
</p>
        <p>y  p n k  q(1 1  n k k )  x (1 p(k 1))  21 ,
where response function in this case
p(n  k)  q(  (1 ) n  k )</p>
        <p>k
x  G  . (26)
y (2 1)(1 p(k 1))
To calculate the coordinates of the fixed point we substitute from (26) expression for
x in the first equation of (25) y2 (kG  (n  k))2  b y(kG  q(n  k 1)) . Hence
y v
Proposition 5. There is unique Nash equilibrium point in a general dynamical system
(24):</p>
        <p>b (kG  q(n  k 1))
y*  v</p>
        <p>x*  Gy*  v
(kG  (n  k))2 (k  (1/ G)(n  k ))2
where the function G  G( p, q, n, k, ) is given in (26).</p>
        <p>
          Proposition 6. For p  0 ( q  1 ) the equilibrium point (x*; y*) coincides with point
(
          <xref ref-type="bibr" rid="ref9">9</xref>
          ) of a dynamic system with naive expectations. When p  1 ( q  0 ) the
equilibrium point coincides with point (18b) of the dynamic system with adaptive
expectations.
        </p>
        <p>(27)
b (k  q(1/ G)(n  k 1))</p>
        <p>C# - Application Model for Numerical Investigation
C# window application Model has been created specifically for the numerical
investigation of the model of this paper, using a graphical interface of C# system
libraries System.Drawing and System.Windows.Forms. Note that all the calculations
associated with the model, are localized in the method calc of the application Model
that makes it easy to modify the equations of the model and use the Model to study
the other two-dimensional dynamical systems. Fig. 1 shows the application window.</p>
        <p>The right side presents 6 kinds of graphs displayed by the application; their
examples are set forth in the paper. Selected switch indicates that here the graph of
trajectory x(t) is selected. On the left side counters allow us to specify the parameters
of the model and the initial values of the trajectory. After their setting the calculation
results of the iterations’ coordinates below and their image in the center of the
window. This displays an animation of a selected path, the number of iterations been
set on the scroll bar above. Pressing the button Model view left displays information
about the model, its equations and parameter information.</p>
        <p>
          With the increasing probability of naive expectations q, that is with decreasing p,
the market becomes unstable, evolving from simple dynamics (15) with a single
stable equilibrium point to the unpredictable behavior of system (
          <xref ref-type="bibr" rid="ref6">6</xref>
          ). From the proof
of Proposition 2 it follows that the market volatility is proportional to the number n of
firms in the market. Therefore, for fixed q market instability increases with increasing
n. Thus, model (24) has two parameters: the number of firms n and the probability of
a naive approach q, whose growth leads to instability. The transition from stability to
chaos is the same in both cases. Consider this transition for parameter q.
Let n n=20, k=6, b=200, v=2, α=0.9, q=0.5. The trajectory of the dynamical system
(24) with the following parameters and the initial point x0=0.1, y0=0.1 is shown in the
following figures 2 and 3. In Fig. 2 on the x-axis of the system are given iterations of
system (24) from m = 1 to m = 100, on the y-axis – corresponding quantity product of
selfish firm ym.
        </p>
        <p>As we can see from the graph, the path quickly converges to the equilibrium value
y*≈2.488. The graph for the trajectory of firm-reciprocator xm on y-axis is similar.
The equilibrium value of x* is about 6.72. Let us consider the graph of the trajectory
for the same parameters except q. Now q = 0.55 (Fig. 3).
It still has stable Nash equilibrium, but 100 iterations does not suffice for
convergence. Further, let q = 0.6 (Fig. 4).</p>
        <p>Fig. 4. Bifurcation quantity trajectory of selfish firm under probability of naive expectations
q  0.6
As we can see, bifurcation occurred, and instead of equilibrium point there was a
steady cycle, where values of ym are approaching the point of y*≈4 for even m and the
point of y*≈1 for odd m. By doubling the lag between iterations only even or only
odd iterations will be considered, and thus either point y*≈4, or y*≈1 respectively
would be the equilibrium steady state.</p>
        <p>Stable cycle has four cycles for q=0.64 (fig. 5). There was a new cycle doubling
bifurcation. Calculations show that with increasing parameter q doubling bifurcation
cycle continues, following Sharkovskii’s scale. According to this scale, when q≈0.675
there is the state of dynamic chaos (fig. 6). Similarly, the graph of product xm on
yaxis by firm-reciprocator looks like trajectory of a selfish firm.</p>
        <p>Fig. 5. Doubling bifurcation cycle of quantity by selfish firm under probability of naive
expectations q  0.64
Fig. 6. The state of dynamic chaos of quantity by selfish firms under probability of naive
expectations q  0.675</p>
        <p>Note that the ratio between the quantity of output by selfish firms and reciprocators
remains almost unchanged. It is demonstrated in the graph of fig. 9, where each
iteration on x-axis shows the value of output by firms-reciprocators xm, and the
vertical axis - the appropriate output of quantity ym of selfish firms (fig. 7.).
In detail the process of loss of stability and transition to chaos of dynamic system (24)
can be presented in the following bifurcation diagram (fig. 8).
probability of naive expectations ( q ) in a general dynamical system</p>
        <p>Here the horizontal axis represents the parameter value of q multiplied by 10. The
ordinate values quantity volumes of selfish firm on stable cycle, multiplied by 0.3.
This rescaling is done for the sake of clarity. The values of the other parameters are
the same as above. The bifurcation diagram, where on vertical axis are placed the
values of output of firms-reciprocators xm looks similar.</p>
        <p>As noted in numerical simulations, the bifurcation may be interpreted as separation
of equilibrium into several ways, one of which is selected by the market due to
evolution of firms’ strategies, such as repeated interactions and adaptations.
Numerical experiments with n firms as the variable parameter are analogous to those
described above.
5</p>
      </sec>
      <sec id="sec-3-2">
        <title>Conclusion</title>
        <p>Thus, we have designed the strategic model of cooperation between the two types of
firms in the market of homogeneous product, where reciprocator and selfish firms
plan their output using the adaptive approach with probability p and naïve (bounded
rationality) one with a probability of q = 1-p, which distinguishes this model from
existing analogues, where each type of firm adheres to one strategy rather than their
combination and maximizes only its own profit rather than social welfare.</p>
        <p>Desktop C# application Model using a graphical interface to animate the model
trajectories has been created specifically for the numerical investigation of the model.</p>
        <p>It has been proved that in the model with adaptive expectations the unique Nash
equilibrium in a dynamic system is stable for all possible values of the parameters.
The trajectories of the dynamical system converge to the fixed point for any possible
initial values. In the model with naive expectations the unique Nash equilibrium is
unstable for sufficiently large values of n for all possible values of other parameters.
According to the calculations, this point is unstable even at n ≥ 5.</p>
        <p>As a result of numerical experiment we have found that bifurcations of cycle
doubling occur with an increase in naive expectations. This bifurcation can be
interpreted as separation of equilibrium state into several ways, one of which is
selected by the market in the evolution of firms’ strategies. If two-thirds of firms use
naive expectation (q≈0.675), then in accordance with the Sharkovskii scale there
appears the state of dynamic chaos in the market, leading to degeneration of the
existing competition model between two types of firms.</p>
        <p>Thus, the crucial factor, which ensures sustainable equilibrium in the market and
the ability to predict the product quantity of firms, is the adaptive approach, i.e. the
one taking into account adaptive expectations of the firms when they plan their
production.</p>
        <p>Similar results are obtained if instead of q we use parameter n - number of firms in
the market, where system also moves from stability to chaos if n increases.</p>
      </sec>
    </sec>
  </body>
  <back>
    <ref-list>
      <ref id="ref1">
        <mixed-citation>
          1.
          <string-name>
            <surname>Barsukova</surname>
          </string-name>
          , S.Y.:
          <article-title>Reciprocal interactions: the nature</article-title>
          , function, specificity.
          <source>Sociological studies 8</source>
          .
          <fpage>20</fpage>
          --
          <lpage>30</lpage>
          (
          <year>2004</year>
          )
        </mixed-citation>
      </ref>
      <ref id="ref2">
        <mixed-citation>
          2.
          <string-name>
            <surname>Polanyi</surname>
            ,
            <given-names>K.</given-names>
          </string-name>
          :
          <article-title>Economics as institutional process</article-title>
          .
          <source>Economic sociology 2</source>
          .
          <fpage>62</fpage>
          --
          <lpage>73</lpage>
          (
          <year>2002</year>
          )
        </mixed-citation>
      </ref>
      <ref id="ref3">
        <mixed-citation>
          3.
          <string-name>
            <surname>Kolm</surname>
            ,
            <given-names>S.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Ythier</surname>
            ,
            <given-names>J.M.</given-names>
          </string-name>
          :
          <year>2006</year>
          .
          <article-title>Handbook of the Economics of Giving, Altruism and Reciprocity</article-title>
          . Elsevier: Worth-Holland (
          <year>2006</year>
          )
        </mixed-citation>
      </ref>
      <ref id="ref4">
        <mixed-citation>
          4.
          <string-name>
            <surname>Bischi</surname>
            ,
            <given-names>G.I.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Lamantia</surname>
            ,
            <given-names>F.</given-names>
          </string-name>
          <article-title>A dynamic model of oligopoly with R&amp;D externalities along networks</article-title>
          .
          <source>Part I. Mathematics and Computers in Simulation</source>
          ,
          <volume>84</volume>
          .
          <fpage>51</fpage>
          --
          <lpage>65</lpage>
          (
          <year>2012</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref5">
        <mixed-citation>
          5.
          <string-name>
            <surname>Puu</surname>
            ,
            <given-names>T.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Panchuk</surname>
            ,
            <given-names>A.</given-names>
          </string-name>
          <article-title>Nonlinear economic dynamics</article-title>
          . New York: Nova Science Publishers (
          <year>2011</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref6">
        <mixed-citation>
          6.
          <string-name>
            <surname>Puu</surname>
            ,
            <given-names>T.</given-names>
          </string-name>
          <article-title>Chaos in duopoly pricing</article-title>
          .
          <source>Chaos, Solitons&amp;Fractals</source>
          <volume>1</volume>
          .
          <fpage>573</fpage>
          --
          <lpage>581</lpage>
          (
          <year>1991</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref7">
        <mixed-citation>
          7.
          <string-name>
            <surname>Szidarovszky</surname>
            ,
            <given-names>F.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Okuguchi</surname>
            ,
            <given-names>K.</given-names>
          </string-name>
          <article-title>An Oligopoly model of commercial fishing</article-title>
          .
          <source>Seoul Journal of Economics</source>
          <volume>11</volume>
          .
          <fpage>321</fpage>
          --
          <lpage>330</lpage>
          (
          <year>1998</year>
          )
        </mixed-citation>
      </ref>
      <ref id="ref8">
        <mixed-citation>
          8.
          <string-name>
            <surname>Bischi</surname>
            ,
            <given-names>G.I.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Naimzada</surname>
            ,
            <given-names>A.K.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Sbragia</surname>
            ,
            <given-names>L.</given-names>
          </string-name>
          :
          <article-title>Oligopoly games with Local Monopolistic Approximation</article-title>
          .
          <source>Journal of Economic Behavior&amp;Organization</source>
          <volume>62</volume>
          .
          <fpage>371</fpage>
          --
          <lpage>388</lpage>
          (
          <year>2007</year>
          )
        </mixed-citation>
      </ref>
      <ref id="ref9">
        <mixed-citation>
          9.
          <string-name>
            <surname>Kobets</surname>
            ,
            <given-names>V.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Poltoratskiy</surname>
            ,
            <given-names>M:</given-names>
          </string-name>
          <article-title>Forming an Evolutionarily Stable Firm Strategy under Cournot Competition Using Social Preferences</article-title>
          . In: Ermolayev,
          <string-name>
            <surname>V.</surname>
          </string-name>
          et al. (eds.) ICT in Education, Research, and
          <string-name>
            <given-names>Industrial</given-names>
            <surname>Applications</surname>
          </string-name>
          .
          <source>Revised Extended Papers of ICTERI</source>
          <year>2014</year>
          , CCIS 469, pp.
          <fpage>343</fpage>
          -
          <lpage>361</lpage>
          , Springer Verlag, Berlin Heidelberg (
          <year>2014</year>
          ) DOI:
          <fpage>10</fpage>
          .1007/978-3-
          <fpage>319</fpage>
          -13206-8_
          <fpage>17</fpage>
        </mixed-citation>
      </ref>
      <ref id="ref10">
        <mixed-citation>
          10.
          <string-name>
            <surname>Kobets</surname>
            ,
            <given-names>V.</given-names>
          </string-name>
          <string-name>
            <surname>Yatsenko</surname>
            ,
            <given-names>V.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Poltoratskiy</surname>
          </string-name>
          .
          <article-title>Dynamic model of double electronic Vickrey auction</article-title>
          . In: Ermolayev,
          <string-name>
            <surname>V.</surname>
          </string-name>
          et al.
          <source>(eds.) Proc. 11-th Int. Conf. ICTERI</source>
          <year>2015</year>
          , Kherson, Lviv, May
          <volume>14</volume>
          -16,
          <year>2015</year>
          , CEUR-WS.org/Vol-
          <volume>1356</volume>
          , ISSN 1613-
          <issue>0073</issue>
          , pp.
          <fpage>236</fpage>
          --
          <lpage>251</lpage>
          , online CEURWS.org/Vol-1356/ICTERI-2015
          <string-name>
            <surname>-CEUR-</surname>
          </string-name>
          WS-Volume.pdf
        </mixed-citation>
      </ref>
      <ref id="ref11">
        <mixed-citation>
          11.
          <string-name>
            <surname>Fanti</surname>
            ,
            <given-names>L.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Gori</surname>
            ,
            <given-names>L.</given-names>
          </string-name>
          <article-title>The dynamics of a differentiated duopoly with quantity competition</article-title>
          .
          <source>Economic Modelling</source>
          ,
          <volume>29</volume>
          .
          <fpage>421</fpage>
          --
          <lpage>427</lpage>
          (
          <year>2012</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref12">
        <mixed-citation>
          12.
          <string-name>
            <surname>Tramontana</surname>
            ,
            <given-names>F.</given-names>
          </string-name>
          <article-title>Heterogeneous duopoly with isoelastic demand function</article-title>
          .
          <source>Economic Modelling</source>
          ,
          <volume>27</volume>
          .
          <fpage>350</fpage>
          --
          <lpage>357</lpage>
          (
          <year>2010</year>
          ).
        </mixed-citation>
      </ref>
    </ref-list>
  </back>
</article>