Mercurial > octave-nkf
view examples/oregonator.m @ 14025:9867be070ee1
use pcre/pcre.h if it is present
* configure.ac: Check for pcre/pcre.h.
* liboctave/regexp.cc: Include either pcre.h or pcre/pcre.h depending
on which is available.
author | John W. Eaton <jwe@octave.org> |
---|---|
date | Sun, 11 Dec 2011 22:35:13 -0500 |
parents | 4295d634797d |
children |
line wrap: on
line source
## The `oregonator'. ## ## Reference: ## ## Oscillations in chemical systems. IV. Limit cycle behavior in a ## model of a real chemical reaction. Richard J. Field and Richard ## M. Noyes, The Journal of Chemical Physics, Volume 60 Number 5, ## March 1974. function dx = oregonator (x, t) dx = zeros (3, 1); dx(1) = 77.27*(x(2) - x(1)*x(2) + x(1) - 8.375e-06*x(1)^2); dx(2) = (x(3) - x(1)*x(2) - x(2)) / 77.27; dx(3) = 0.161*(x(1) - x(3)); end