1:
2:
3:
4:
5:

f = @(x)cos(x)

resabs = sum_simpson(abs(f),a,b,tol);