function [val] = di(f, a, b, c, d, n)
    gs = newton(f, a, b, n);
    val = double(rect(gs, c, d, n));
end