function y = func2(x)
    y = exp(x) - (1 + x + x.^2 / 2);
end