octave:>function y=f(x) 
>y = x / (1 + sin(x)^2);
>endfunction 

octave:>[V, IER, NFUN, ERR]=quad("f", 0, 2*pi)