File: lsode-2.m

package info (click to toggle)
octave 2.0.13.95-1
  • links: PTS
  • area: main
  • in suites: slink
  • size: 26,052 kB
  • ctags: 16,420
  • sloc: cpp: 67,184; fortran: 41,514; ansic: 26,607; sh: 7,291; makefile: 4,089; lex: 1,961; yacc: 1,852; perl: 1,676; lisp: 1,664; exp: 123
file content (13 lines) | stat: -rw-r--r-- 221 bytes parent folder | download | duplicates (5) 
1
2
3
4
5
6
7
8
9
10
11
12
13
 
function xdotdot = f (x, t)
  xdotdot = [x(2); -x(1)];
endfunction

x0 = [1; 0];
t = [0; 2*pi];
tol = 100 * dassl_options ("absolute tolerance");

x = lsode ("f", x0, t);

y = [1, 0; 1, 0];

all (all (abs (x - y) < tol))