models/msmith.lua

   1 --%  ft| BIFURCATION_1 sn| Route~to~chaos~by~a~Neimark-Sacker~bifurcation n| #0 d| 0.2 n| #1 d| 0.1 n| #2 d| 0.1 n| #3 d| epsilon n| #4 d| 1.27 n| #5 d| 1.48 n| #6 d| -0.2 n| #7 d| 1.05 n| #8 d| 500 n| #9 d| 100 n| #10 d| x
   2 --%  ft| TRAJECTORY_T0_V0_A1_O0 sn| Chaotic~attractor n| #0 d| 0.2 n| #1 d| 0.1 n| #2 d| -1.375 n| #3 d| 1.8 n| #4 d| 100 n| #5 d| 10000 n| #6 d| 10000 n| #7 d| x n| #8 d| y
   3 --@@
   4 name = "MS"
   5 description = "See Model refs in user's guide"
   6 type = "D"
   7 parameters = {"epsilon", "mu"}
   8 variables = {"x", "y"}
   9
  10 function f (epsilon, mu, x, y)
  11
  12         x1 = epsilon * x + mu - y^2
  13         y1 = x
  14
  15 return x1, y1
  16
  17 end
  18
  19 function Jf (epsilon, mu, x, y)
  20
  21         return
  22
  23         epsilon,        -2 * y,
  24         1,               0
  25
  26 end
  27
  28
  29