%
% Subroutine that calls "fdstencilnonuniform" to obtain an approximation 
% of the value of the kth derivative of a function "f" at a point "xb" 
% using a finite difference approximation based on "n" given points (n>k).

%It also computes the exact derivative value at xb, the
% error between the approximation and the exact value of the derivative, and the
% estimated error using the expression for the local truncation error. For
% these calculations, the kth derivative, the nth and nth+1 derivatives are
% also provided.


function FDnonuniformNew(f,fkp,fnp,fn1p,k,xb,xpts)