% f1(x) = 1 + O(h^2) function value = f1(x) value = exp(x) - x; end