% batch_quadint
f = @(x) (sqrt(16 - x .^2)/8); % true value of the integral: pi

v = quadint(f, -4, 4, 100000)