clear, clc

a=0;
b=.8;
n=4000;
f = @(x) 0.2+25*x-200*x^2+675*x^3-900*x^4+400*x^5;
%f = @(x) (1+1/x)^2;

I = e2_trap(f,a,b,n);

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