annotate examples/oregonator.cc @ 2689:8c7955a8d49f

[project @ 1997-02-18 09:06:10 by jwe]
author jwe
date Tue, 18 Feb 1997 09:09:12 +0000
parents
children 7c96e85c76db
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 #include <iostream.h>
8c7955a8d49f [project @ 1997-02-18 09:06:10 by jwe]
jwe
parents:
diff changeset
4
8c7955a8d49f [project @ 1997-02-18 09:06:10 by jwe]
jwe
parents:
diff changeset
5 DEFUN_DLD (oregonator, args, ,
8c7955a8d49f [project @ 1997-02-18 09:06:10 by jwe]
jwe
parents:
diff changeset
6 "The `oregonator'.")
8c7955a8d49f [project @ 1997-02-18 09:06:10 by jwe]
jwe
parents:
diff changeset
7 {
8c7955a8d49f [project @ 1997-02-18 09:06:10 by jwe]
jwe
parents:
diff changeset
8 ColumnVector dx (3);
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 ColumnVector x = args(0).vector_value ();
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(0) = 77.27 * (x(1) - x(0)*x(1) + x(0) - 8.375e-06*pow (x(0), 2));
8c7955a8d49f [project @ 1997-02-18 09:06:10 by jwe]
jwe
parents:
diff changeset
13 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
14 dx(2) = 0.161*(x(0) - x(2));
8c7955a8d49f [project @ 1997-02-18 09:06:10 by jwe]
jwe
parents:
diff changeset
15
8c7955a8d49f [project @ 1997-02-18 09:06:10 by jwe]
jwe
parents:
diff changeset
16 return octave_value (dx);
8c7955a8d49f [project @ 1997-02-18 09:06:10 by jwe]
jwe
parents:
diff changeset
17 }