d_x = (b - a) / 2                # Half interval length