lsode-3.m

来自「GNU Octave is a high-level language, pri」· M 代码 · 共 14 行

function xdot = f (x, t)
  xdot = x;
endfunction

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

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

y = [1; e];

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