models/Primer/vanderpol.lua

   1 --%  ft| TRAJECTORY_T0_V0_A1_O1 sn| Convergence~to~a~fixed~point n| #0 d| 0.5 n| #1 d| 0.1 n| #2 d| 5 n| #3 d| -0.1 n| #4 d| 1.5 n| #5 d| 0.01 n| #6 d| 0 n| #7 d| 10000 n| #8 d| 10000 n| #9 d| x n| #10 d| y
   2 --@@
   3 name = "VanderPol"
   4 description = "See Model refs in user's guide"
   5 type = "C"
   6 parameters = {"k", "mu", "b"}
   7 variables = {"x", "y"}
   8
   9 function f(k, mu, b, x, y)
  10
  11         x1= k * y + mu * x * (b - y^2)
  12         y1= - x + mu
  13
  14         return x1, y1
  15
  16 end
  17
  18
  19 function Jf(k, mu, b, x, y)
  20
  21 return
  22
  23 mu * (b - y^2), k - 2 * mu * x * y,
  24 -1,                     0
  25
  26 end