f = @(x)sin(x) ;
s = simpson(f,0,pi,10)
error1 = s-2

x = linspace(0,pi,11);
y = sin(x);
s2 = simpson(y,0,pi)
error2 = s2-2
