models/Primer/conb.lua

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