clear;
clf;
format long
a=0.1;
b=2;
%Mencari luas f(x) dengan sb.x dalam 
%interval [a,b]
n=10;
for m=1:10
   k(m)=m
   clc
   dx=(b-a)/n;
	for i=1:n
      LuasTrap(i)=(1/2)*dx*(fungsi(a+((i-1)*dx))+fungsi(a+(i*dx)));
   	LuasSimp(i)=(1/6)*dx*(fungsi(a+((i-1)*dx))+(4*(fungsi(a+(i-0.5)*dx)))+fungsi(a+(i*dx)));
	end
	LuasTrapTot(m)=sum(LuasTrap);
	LuasSimpTot(m)=sum(LuasSimp);
   Error(m)=abs(LuasTrapTot(m)-LuasSimpTot(m));
   if (Error(m))<0.00000001
  		break
   else
      n=2*n;
   end
end
Hasil=[k' LuasTrapTot' LuasSimpTot' Error']
plot(k,Error), grid on