models/ikeda.lua

   1 --%  ft| BIFURCATION_1 sn| Period-doubling~route~to~chaos n| #0 d| 1 n| #1 d| 1 n| #2 d| 1 n| #3 d| 0.4 n| #4 d| 6 n| #5 d| c n| #6 d| 0.35 n| #7 d| 0.88 n| #8 d| -0.25 n| #9 d| 1.6 n| #10 d| 500 n| #11 d| 200 n| #12 d| x
   2 --%  ft| TRAJECTORY_T0_V0_A1_O0 sn| Chaotic~attractor n| #0 d| 0 n| #1 d| 1 n| #2 d| 1 n| #3 d| 0.4 n| #4 d| 0.9 n| #5 d| 6 n| #6 d| 100 n| #7 d| 20000 n| #8 d| 20000 n| #9 d| x n| #10 d| y
   3 --@@
   4
   5 name = "Ikeda"
   6 description = "See Model refs in user's guide"
   7 type = "D"
   8
   9 parameters = {"a", "b", "c", "d"}
  10
  11 variables = {"x", "y"}
  12
  13
  14 function f(a, b, c, d, x, y)
  15
  16     r = b - d / (1.0 + x^2 + y^2)
  17
  18     x1 = a + c * (x * math.cos(r) - y * math.sin(r))
  19     y1 = c * (x * math.sin(r) + y * math.cos(r))
  20
  21     return x1, y1
  22
  23 end