Non-monotonic behavior of exp_remainder

[schillic]

Around order p = 20 the result is the smallest and after that it grows again.

julia> m = IntervalMatrix([-1.1..0.9 -4.1.. -3.9; 3.9..4.1 -1.1..0.9]);

julia> for p in [10, 15, 20, 25, 30]
           println(IntervalMatrices._exp_remainder(m, 1.0, p))
       end
Interval{Float64}[[-3.20833, 3.20833] [-3.20833, 3.20833]; [-3.20833, 3.20833] [-3.20833, 3.20833]]
Interval{Float64}[[-0.0194887, 0.0194887] [-0.0194887, 0.0194887]; [-0.0194887, 0.0194887] [-0.0194887, 0.0194887]]
Interval{Float64}[[-2.77417e-05, 2.77417e-05] [-2.77417e-05, 2.77417e-05]; [-2.77417e-05, 2.77417e-05] [-2.77417e-05, 2.77417e-05]]
Interval{Float64}[[-0.0922537, 0.0922537] [-0.0922537, 0.0922537]; [-0.0922537, 0.0922537] [-0.0922537, 0.0922537]]
Interval{Float64}[[-452.145, 452.145] [-452.145, 452.145]; [-452.145, 452.145] [-452.145, 452.145]]