home *** CD-ROM | disk | FTP | other *** search
/ OS/2 Shareware BBS: 10 Tools / 10-Tools.zip / octa21eb.zip / octave / DLFCN.ZIP / dlfcn / examples / oregonator.m < prev   
Text File  |  1998-11-10  |  465b  |  19 lines
  1. ## The `oregonator'.
  2. ##
  3. ## Reference:
  4. ##
  5. ##   Oscillations in chemical systems.  IV.  Limit cycle behavior in a
  6. ##   model of a real chemical reaction. Richard J. Field and Richard
  7. ##   M. Noyes, The Journal of Chemical Physics, Volume 60 Number 5,
  8. ##   March 1974.
  9.  
  10. function dx = oregonator (x, t)
  11.  
  12.   dx = zeros (3, 1);
  13.  
  14.   dx(1) = 77.27*(x(2) - x(1)*x(2) + x(1) - 8.375e-06*x(1)^2);
  15.   dx(2) = (x(3) - x(1)*x(2) - x(2)) / 77.27;
  16.   dx(3) = 0.161*(x(1) - x(3));
  17.  
  18. end
  19.