%
% 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
% 
d1_h1 = diffgen('c_4_p_1_fn', 1, 1, 0.1 );
d1_h2 = diffgen('c_4_p_1_fn', 1, 1, 0.01 );

d2_h1 = diffgen('c_4_p_1_fn', 2, 1, 0.1 );
d2_h2 = diffgen('c_4_p_1_fn', 2, 1, 0.01 );

disp( sprintf( 'D(h=0.1)= %10.8f; D(h=0.01)= %10.8f', d1_h1, d1_h2 ) );
disp( sprintf( 'D^2(h=0.1)= %10.8f; D^2(h=0.01)= %10.8f', d2_h1, d2_h2 ) );