close all; clc; clear; 

% Problem EPage 82
% 
fn = @ (x) ( exp( -x.^2 ) .* cos( x ) );
approx = gaherm( fn, 16 )
truth = sqrt( pi ) * exp( -1/4 )
truth - approx