[bug-gsl] Re: diff.c (numerical differentiation) (2008-04-08)
|
From: |
Brian Gough <bjg@network-theory.co.uk> |
|
Subject: |
Re: diff.c (numerical differentiation) |
|
Date: |
Tue, 08 Apr 2008 21:55:31 +0100 |
|
To: |
spasmous <spasmous@.........> |
|
Cc: |
bug-gsl@gnu.org |
At Sun, 6 Apr 2008 13:29:41 -0700,
spasmous wrote:
> h = sqrt (GSL_SQRT_DBL_EPSILON / (2.0 * a2));
> ...
> h = (x+h)-x; // ensure difference between (x+h) and x is exactly h
> *result = (GSL_FN_EVAL (f, x + h) - GSL_FN_EVAL (f, x)) / h;
Thanks for the suggestion.
--
Brian Gough
