Mercurial > hg > octave-nkf

annotate examples/oregonator.cc @ 17791:224e76250443

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Use GNU style coding conventions for code in examples/ * examples/embedded.cc, examples/fortdemo.cc, examples/helloworld.cc, examples/make_int.cc, examples/mycell.c, examples/myfeval.c, examples/myfunc.c, examples/mypow2.c, examples/myprop.c, examples/myset.c, examples/mystring.c, examples/mystruct.c, examples/oct_demo.cc, examples/oregonator.cc: Use GNU style coding conventions for code in examples/
author Rik <rik@octave.org>
date Tue, 29 Oct 2013 10:28:39 -0700
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 #include <octave/oct.h>
8c7955a8d49f [project @ 1997-02-18 09:06:10 by jwe]
jwe
parents:
diff changeset
2
8c7955a8d49f [project @ 1997-02-18 09:06:10 by jwe]
jwe
parents:
diff changeset
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
7c96e85c76db [project @ 1998-04-08 18:19:35 by jwe]
jwe
parents: 2689
diff changeset
5 \n\
7c96e85c76db [project @ 1998-04-08 18:19:35 by jwe]
jwe
parents: 2689
diff changeset
6 Reference:\n\
7c96e85c76db [project @ 1998-04-08 18:19:35 by jwe]
jwe
parents: 2689
diff changeset
7 \n\
7c96e85c76db [project @ 1998-04-08 18:19:35 by jwe]
jwe
parents: 2689
diff changeset
8 Oscillations in chemical systems. IV. Limit cycle behavior in a\n\
7c96e85c76db [project @ 1998-04-08 18:19:35 by jwe]
jwe
parents: 2689
diff changeset
9 model of a real chemical reaction. Richard J. Field and Richard\n\
7c96e85c76db [project @ 1998-04-08 18:19:35 by jwe]
jwe
parents: 2689
diff changeset
10 M. Noyes, The Journal of Chemical Physics, Volume 60 Number 5,\n\
7c96e85c76db [project @ 1998-04-08 18:19:35 by jwe]
jwe
parents: 2689
diff changeset
11 March 1974.")
2689
8c7955a8d49f [project @ 1997-02-18 09:06:10 by jwe]
jwe
parents:
diff changeset
12 {
8c7955a8d49f [project @ 1997-02-18 09:06:10 by jwe]
jwe
parents:
diff changeset
13 ColumnVector dx (3);
8c7955a8d49f [project @ 1997-02-18 09:06:10 by jwe]
jwe
parents:
diff changeset
14
3681
df54d394acc0 [project @ 2000-06-26 17:46:58 by jwe]
jwe
parents: 3659
diff changeset
15 ColumnVector x (args(0).vector_value ());
2689
8c7955a8d49f [project @ 1997-02-18 09:06:10 by jwe]
jwe
parents:
diff changeset
16
4205
6b7bd9e406be [project @ 2002-11-27 22:04:51 by jwe]
jwe
parents: 3681
diff changeset
17 dx(0) = 77.27 * (x(1) - x(0)*x(1) + x(0) - 8.375e-06*pow (x(0), 2.0));
2689
8c7955a8d49f [project @ 1997-02-18 09:06:10 by jwe]
jwe
parents:
diff changeset
18 dx(1) = (x(2) - x(0)*x(1) - x(1)) / 77.27;
8c7955a8d49f [project @ 1997-02-18 09:06:10 by jwe]
jwe
parents:
diff changeset
19 dx(2) = 0.161*(x(0) - x(2));
8c7955a8d49f [project @ 1997-02-18 09:06:10 by jwe]
jwe
parents:
diff changeset
20
8c7955a8d49f [project @ 1997-02-18 09:06:10 by jwe]
jwe
parents:
diff changeset
21 return octave_value (dx);
8c7955a8d49f [project @ 1997-02-18 09:06:10 by jwe]
jwe
parents:
diff changeset
22 }