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 _______________________________________________ Bug-gsl mailing list [email protected] http://lists.gnu.org/mailman/listinfo/bug-gsl
