//SLUGens released under the GNU GPL as extensions for SuperCollider 3, by Nick Collins, http://composerprogrammer.com/index.html
NOTE: experimental, not working well
Naive Euler ODE solver implementation of a neural oscillator model by Terman and Wang (see Brown, G.J., and Wang, D. 2006. "Neural and Perceptual Modeling." In Wang, D., and Brown, G. J., eds. Computational Auditory Scene Analysis: Principles, Algorithms, and Applications. Hoboken, NJ: John Wiley and Sons/IEEE Press, p.355).
Dx = ratex * ( 3*x - x^3 (+ 2) - y + input )
Dy = ratey * (eta ( alpha* (1 + tanh(x/beta)) - y ))
All parameters can have .kr rate UGens plugged in. The +2 is left out and can be added back in as part of the input. eta can also be left out as part of yrate.
audio rate input
update rate for x
update rate for y
reset value for u
reset value for w