models/tinkerbell.lua

   1 --%  ft| TRAJECTORY_T0_V0_A1_O0 sn| Chaotic~Tinkerbell~attractor n| #0 d| .1 n| #1 d| .1 n| #2 d| 0.9 n| #3 d| -0.6014 n| #4 d| 2 n| #5 d| 0.5 n| #6 d| 1000 n| #7 d| 100000 n| #8 d| 100000 n| #9 d| x n| #10 d| y
   2 --%  ft| BIFURCATION_1 sn| Route~to~chaos~by~Neimark-Sacker~bifurcations n| #0 d| .1 n| #1 d| .1 n| #2 d| 0.9 n| #3 d| -0.6014 n| #4 d| 2 n| #5 d| c4 n| #6 d| -0.27 n| #7 d| 0.52 n| #8 d| -1.25 n| #9 d| 0.37 n| #10 d| 500 n| #11 d| 100 n| #12 d| x
   3 --@@
   4 name = "Tinkerbell Map"
   5 description = "See Model refs in user's guide"
   6 type = "D"
   7 parameters = {"c1", "c2", "c3", "c4"}
   8 variables = {"x", "y"}
   9
  10 function f(c1, c2, c3, c4, x, y)
  11
  12         x1 = x^2 - y^2 + c1 * x + c2 * y
  13         y1 = 2 * x * y + c3 * x + c4 * y
  14
  15         return x1, y1
  16
  17 end
  18
  19 function Jf(c1, c2, c3, c4, x, y)
  20
  21         return
  22
  23         2 * x + c1,             -2 * y + c2,
  24       2 * y + c3,               2 * x + c4
  25
  26 end
  27