(bind-func ndiff "numeric differentiation - central diff" (lambda (f:[double,double]* x dx) (/ (- (f (+ x (/ dx 2.0))) (f (- x (/ dx 2.0)))) dx)))