close all; clc; clear; 

% Problem EPage 82
% 
for k=[0,4,100]
  approx = filon( @ (x) sin(x), 1, k, 0, pi, 64 );
  truth = 2/(1-k^2);
  disp( sprintf('k= %5d; filon approx= %10.6f, truth= %10.6f', k, approx, truth) );
end