octave:>function y=f(x) >y = x / (1 + sin(x)^2); >endfunction octave:>[V, IER, NFUN, ERR]=quad("f", 0, 2*pi)