D <- function(f,delta=.000001){ function(x){ (f(x+delta) - f(x-delta))/(2*delta)} }