clear, clc

a=0;
b=pi/4;
n=8;
%f = @(x) 0.2+25*x-200*x^2+675*x^3-900*x^4+400*x^5;
%f = @(x) (1 + 1/x)^2;
f= @(x) log(1+tan(x));
I = e4_Simpson(f,a,b,n);

fprintf('\n\n I = %0.9f unit square\n\n',I)

fplot(f, [a,b])
grid on