%
% Written by:
% --
% John L. Weatherwax                2009-04-21
%
% email: wax@alum.mit.edu
%
% Please send comments and especially bug reports to the
% above email address.
%
%-----

addpath('../../BookCode');

close all; clc; clear; 

% Problem EPage 81; Example EPage 66
% 
global r;

for r = [0,1,2] 
  approx = galag( 'c_4_p_7_fn', 8 );
  truth = pi / ( 2 * ( 1 + r + r^3 ) );
  disp( sprintf( 'r = %10.6f, truth= %10.6f, approx= %10.6f', r, truth, approx ) );
end