octave:1> function y=f(x) > y = x * x; > endfunction octave:5> [v,iner,nfun,err]=quad("f",0,1) v = 0.33333 iner = 0 nfun = 21 err = 3.7007e-15 octave:6> quad("f", 0, 2) ans = 2.6667 octave:7> quad("f", 0, 1) ans = 0.33333
