Google has failed me. Any responses are greatly appreciated.
I don't have it with me right now but I think it reads 10^-15 as zero, so the easiest way would be approximate it with the interval x,x+10^-14 where x is the point. It is optional in the nDeriv function to state an interval but if you don't input one I can only assume it takes the smallest number it has.
You would lose a lot of significant digits if you do it that way. I think that if derivatives have to be evaluated numerically, then you would be better off using some extraplation orseries technique. E.g. you can write the Taylor expansion of a function formally as:
f(x+t) = exp(t d/dx) f(x)
The symmetric difference with step t is thus given by:
Delta_t f(x) = [exp(t d/dx) - exp(-td/dx)]/2 f(x) = sinh(t d/dx) f(x)
So, this means that formally we can express the derivative operator in terms of the finite symmetric difference operator as:
d/dx = 1/t arcsinh(Delta_t) = 1/t [Delta_t - 1/6 Delta_t^3 + ...]
So, to comnpute the derivative at a point, all you need to do is to repeatedly apply the symmteric finite difference operator with some stepsize t. The smaler you take t, the faster te series converges, but then you lose significant digits. So, you should take t not too small and a few terms of the series.
Most graphing calculators compute a derivative by taking the symmetric difference quotient with the value of the difference being a small number close to zero such as .001.
Ref: Calculus: Graphical, Numerical, Algerbraic, by Ross Finney et. al. p. 111.
Separate names with a comma.