(bind-func ndiff2 "numeric differentiation (2nd deriv) - central diff" (lambda (f:[double,double]* x) (/ (- (- (f (+ x 1e-5)) (f x)) (- (f x) (f (- x 1e-5)))) (* 1e-5 1e-5))))