models/sinecircle.lua

   1 name = "sine circle"
   2 description = "See Model refs in user's guide"
   3 type = "D"
   4 parameters = {"omega", "k"}
   5 variables = {"theta"}
   6
   7 function f(omega, k, theta)
   8
   9 two_pi = 2 * math.pi
  10
  11         y1 = theta + omega + (k/two_pi)*math.sin(two_pi*theta)
  12
  13         y2 = math.mod(y1, 1)
  14
  15         return y2
  16
  17 end
  18
  19 function Jf(omega, k, theta)
  20
  21         return 1 + k*cos(two_pi*theta)
  22 end
  23