//SLUGens released under the GNU GPL as extensions for SuperCollider 3, by Nick Collins, http://composerprogrammer.com/index.html
Euler ODE solver implementation of the Brusselator equations (http://en.wikipedia.org/wiki/Brusselator).
x' = x^2*y - (mu+1)*x + gamma
y' = -x^2*y + mu*x
All inputs can have .kr rate UGens plugged in.
Nonlinear oscillators can blow up, treat with caution. This one is relatively stable however, converging to a fixed point, or a limit cycle, in the upper positive quadrant. for some reasonable values. Just be careful if mu gets much bigger than gamma (though making it larger is necessary to get some chaotic oscillation behaviour); you can retrigger to get back to normal, and keep the rate lower to avoid Euler integration blow-ups. You may just need to scale and push down around zero to avoid a DC offset. Fixed point is at (x,y) = (gamma, mu/gamma)
If > 0.0, restart with new initial conditions sampled from initx, inity
update rate for a sample step
equation constant. Set mu > (gamma**2 + 1.0) for the more fun limit cycle regions
reset value for x
reset value for y