x = [0.1:0.1:2];
h = 0.01
f = @sin
df = cos(x);
d2f = -sin(x)

mean(abs(df - d1_forward(f, x, h)))
mean(abs(df- d1_forward2(f, x, h)))
mean(abs(df- d1_backward(f, x, h)))
mean(abs(df- d1_backward2(f, x, h)))
mean(abs(df - d1_central(f, x, h)))
mean(abs(df - d1_central4(f, x, h)))
mean(abs(d2f - d2_central(f, x, h)))
mean(abs(d2f - d2_central4(f, x, h)))
