;; same as above (i.e. central) (bind-func ndifferentiate "numeric differentiation - central diff" (lambda (f:[double,double]* x) (/ (- (f (+ x (/ 1e-5 2.0))) (f (- x (/ 1e-5 2.0)))) 1e-5)))