%
% Call quadtx 
%
%addpath C:/cygwin/home/bai/NCM_Moler/NCM

f = inline('1/x'); 
a = 1; 
b = 2; 
exactval = log(2);  

tol = 1e-6;   % specify error tolerance  
[Q,fcount] = quadtx(f,a,b,tol) 

exactval 
relative_error = abs(Q-exactval)/exactval
