f = @(x) exp(x);
d = richardson(f,1.0, 5, 0.5)