%
% 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
% 
d_x1 = diffgen('c_4_p_2_fn', 1, 1, 0.001 );
d_x2 = diffgen('c_4_p_2_fn', 1, 2, 0.001 );
d_x3 = diffgen('c_4_p_2_fn', 1, 3, 0.001 );

disp( sprintf( 'D(x=1)= %10.8f; D(x=2)= %10.8f; D(x=3)= %10.8f', d_x1, d_x2, d_x3 ) );