models/Economic/cournot.lua

   1 --%  ft| BIFURCATION_1 sn| Period~-~doubling~route~to~chaos n| #0 d| 0.01 n| #1 d| 0.01 n| #2 d| 1 n| #3 d| b n| #4 d| 5.77 n| #5 d| 6.25 n| #6 d| -0.01 n| #7 d| 0.17 n| #8 d| 3000 n| #9 d| 200 n| #10 d| x
   2 --%  ft| TRAJECTORY_T0_V0_A1_O0 sn| Chaotic~attractor n| #0 d| 0.01 n| #1 d| 0.01 n| #2 d| 1 n| #3 d| 6.2 n| #4 d| 100 n| #5 d| 5000 n| #6 d| 5000 n| #7 d| x n| #8 d| y
   3 --@@
   4 name = "Cournot Olygopoly"
   5 description = "See Model refs in user's guide"
   6 type = "D"
   7 parameters = {"a", "b"}
   8 variables ={"x", "y"}
   9
  10 function f(a, b,  x, y)
  11
  12         if (a*y <= 1) then
  13
  14         x1= math.sqrt(y/a)-y
  15
  16         else
  17
  18         x1=0
  19
  20         end
  21
  22         if (b*x <= 1) then
  23
  24         y1= math.sqrt(x/b)-x
  25
  26         else
  27
  28         y1=0
  29
  30         end
  31
  32
  33 return x1, y1
  34
  35 end
  36
  37
  38
  39 function Jf(a,b,x,y)
  40
  41
  42         if (a*y <= 1) then
  43
  44         dxdy =math.sqrt(1/(4*a*y))-1
  45
  46         else
  47
  48         dxdy = 0
  49
  50         end
  51
  52
  53         if (b*x <= 1) then
  54
  55         dydx = math.sqrt(1/(4*b*x))-1
  56
  57         else
  58
  59         dydx = 0
  60
  61         end
  62
  63
  64 return  0,      dxdy,
  65                 dydx,   0
  66
  67 end
  68
  69
  70
  71
  72
  73
  74
  75
  76
  77