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