function A=simp(f,a,b,n)
i1=trap(f,a,b,n);
i2=trap(f,a,b,n/2);
A=1/3*(4*i1-i2);