clear all


f=inline('2./(1+x.^2)', 'x'); % f=@(x) 2./(1+x.^2);
a=-1; b=1; tol=1e-6;


[Q,k]=mid_comp(f,a,b,tol);
fprintf('mid:           Q=%19.15f    fcount=%5d\n', Q,k)

[Q,k]=mid_compextr(f,a,b,tol);
fprintf('mid (extr):    Q=%19.15f    fcount=%5d\n', Q,k)

[Q,k]=trap_comp(f,a,b,tol);
fprintf('trap:          Q=%19.15f    fcount=%5d\n', Q,k)

[Q,k]=trap_compextr(f,a,b,tol);
fprintf('trap (extr):   Q=%19.15f    fcount=%5d\n', Q,k)

[Q,k]=simp_comp(f,a,b,tol);
fprintf('simp:          Q=%19.15f    fcount=%5d\n', Q,k)

[Q,k]=simp_compextr(f,a,b,tol);
fprintf('simp (extr):   Q=%19.15f    fcount=%5d\n', Q,k)

[Q,k]=quad(f,a,b,tol);
fprintf('quad:          Q=%19.15f    fcount=%5d\n', Q,k)
