Mercurial > octave-nkf
annotate examples/oregonator.cc @ 18479:85d31344fb23 stable rc-3-8-1-2
3.8.1 release candidate 2
* configure.ac (OCTAVE_VERSION): Bump to 3.8.1-rc2.
(OCTAVE_MINOR_VERSION): Bump to 1-rc2.
(OCTAVE_RELEASE_DATE): Set to 2014-02-14.
author | John W. Eaton <jwe@octave.org> |
---|---|
date | Mon, 17 Feb 2014 02:20:39 -0500 |
parents | 224e76250443 |
children |
rev | line source |
---|---|
2689 | 1 #include <octave/oct.h> |
2 | |
3 DEFUN_DLD (oregonator, args, , | |
17791
224e76250443
Use GNU style coding conventions for code in examples/
Rik <rik@octave.org>
parents:
9053
diff
changeset
|
4 "The `oregonator'.\n\ |
3162 | 5 \n\ |
6 Reference:\n\ | |
7 \n\ | |
8 Oscillations in chemical systems. IV. Limit cycle behavior in a\n\ | |
9 model of a real chemical reaction. Richard J. Field and Richard\n\ | |
10 M. Noyes, The Journal of Chemical Physics, Volume 60 Number 5,\n\ | |
11 March 1974.") | |
2689 | 12 { |
13 ColumnVector dx (3); | |
14 | |
3681 | 15 ColumnVector x (args(0).vector_value ()); |
2689 | 16 |
4205 | 17 dx(0) = 77.27 * (x(1) - x(0)*x(1) + x(0) - 8.375e-06*pow (x(0), 2.0)); |
2689 | 18 dx(1) = (x(2) - x(0)*x(1) - x(1)) / 77.27; |
19 dx(2) = 0.161*(x(0) - x(2)); | |
20 | |
21 return octave_value (dx); | |
22 } |