% Test of integration using table data
clear
load intdata.dat
x = intdata(:,1);
y = intdata(:,2);
h = x(2) - x(1);
exact = 1.1064840;

% Trap rule
t = trap_fun(h,y);
rt = (t-exact)/exact

% Simpsons rule
s = simp_fun(h,y);
rs = (s-exact)/exact