quad-2.m

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

function y = f (x)
  y = x .* sin (1 ./ x) .* sqrt (abs (1 - x));
endfunction
[v, ier, nfun, err] = quad ("f", 0, 3);
(ier == 0 || ier == 1) && v - 1.98194122455795 < sqrt (eps) && nfun > 0