;; numeric differentation (bind-func ndiff "numeric differentiation - central diff" (lambda (f:[double,double]* x) (/ (- (f (+ x (/ 1e-5 2.0))) (f (- x (/ 1e-5 2.0)))) 1e-5)))