Mercurial > octave-nkf

annotate examples/oregonator.m @ 18479:85d31344fb23 stable rc-3-8-1-2

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
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 4295d634797d
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
2689
8c7955a8d49f [project @ 1997-02-18 09:06:10 by jwe]
jwe
parents:
diff changeset
1 ## The `oregonator'.
3162
7c96e85c76db [project @ 1998-04-08 18:19:35 by jwe]
jwe
parents: 2689
diff changeset
2 ##
7c96e85c76db [project @ 1998-04-08 18:19:35 by jwe]
jwe
parents: 2689
diff changeset
3 ## Reference:
7c96e85c76db [project @ 1998-04-08 18:19:35 by jwe]
jwe
parents: 2689
diff changeset
4 ##
7c96e85c76db [project @ 1998-04-08 18:19:35 by jwe]
jwe
parents: 2689
diff changeset
5 ## Oscillations in chemical systems. IV. Limit cycle behavior in a
7c96e85c76db [project @ 1998-04-08 18:19:35 by jwe]
jwe
parents: 2689
diff changeset
6 ## model of a real chemical reaction. Richard J. Field and Richard
7c96e85c76db [project @ 1998-04-08 18:19:35 by jwe]
jwe
parents: 2689
diff changeset
7 ## M. Noyes, The Journal of Chemical Physics, Volume 60 Number 5,
7c96e85c76db [project @ 1998-04-08 18:19:35 by jwe]
jwe
parents: 2689
diff changeset
8 ## March 1974.
2689
8c7955a8d49f [project @ 1997-02-18 09:06:10 by jwe]
jwe
parents:
diff changeset
9
8c7955a8d49f [project @ 1997-02-18 09:06:10 by jwe]
jwe
parents:
diff changeset
10 function dx = oregonator (x, t)
8c7955a8d49f [project @ 1997-02-18 09:06:10 by jwe]
jwe
parents:
diff changeset
11
8c7955a8d49f [project @ 1997-02-18 09:06:10 by jwe]
jwe
parents:
diff changeset
12 dx = zeros (3, 1);
8c7955a8d49f [project @ 1997-02-18 09:06:10 by jwe]
jwe
parents:
diff changeset
13
8c7955a8d49f [project @ 1997-02-18 09:06:10 by jwe]
jwe
parents:
diff changeset
14 dx(1) = 77.27*(x(2) - x(1)*x(2) + x(1) - 8.375e-06*x(1)^2);
8c7955a8d49f [project @ 1997-02-18 09:06:10 by jwe]
jwe
parents:
diff changeset
15 dx(2) = (x(3) - x(1)*x(2) - x(2)) / 77.27;
8c7955a8d49f [project @ 1997-02-18 09:06:10 by jwe]
jwe
parents:
diff changeset
16 dx(3) = 0.161*(x(1) - x(3));
8c7955a8d49f [project @ 1997-02-18 09:06:10 by jwe]
jwe
parents:
diff changeset
17
8c7955a8d49f [project @ 1997-02-18 09:06:10 by jwe]
jwe
parents:
diff changeset
18 end