Marcelo Forets

The main PyPlot - specific option used seems to be the step function, may be generalized with sth like ```julia pyplot() j = (0:(s.n-1))/s.n; Plots.plot([s.u[:];s.u[s.n];s.d[s.n]], [j;1;1], linetype=:steppre) i = (1:(s.n))/s.n;...

An easy alternative (if you are not interested in supporting other plotting backends in this package) might be to add PyPlot as an optional dependency through Requires.jl, and rename the...

Taking products of interval matrices is not even associative, so this problem is (unfortunately) not uncommon. A general rule is to avoid computing the same terms, some articles call this...

I also find the [definition](https://github.com/JuliaIntervals/IntervalArithmetic.jl/blob/master/src/multidim/intervalbox.jl#L72) confusing; from a set perspective: *all* intervals should be included and *at least* one of them should be strictly included.

FWIW, enclosure between interval boxes is also used in Picard iteration (see [TaylorModels.iscontractive](https://github.com/JuliaIntervals/TaylorModels.jl/blob/7441eb2f30f841c7f8b88cceda4227302adc5d3e/src/validatedODEs.jl#L67)). That said, I think that the (less strict) idea in my [previous comment](https://github.com/JuliaIntervals/IntervalArithmetic.jl/issues/490#issuecomment-921848476) also applies.

Bump :smile: Is this "only" needing to add testing and docs?

I think it's fair to say IntervalArithmetic.jl is used in most if not all the publications listed [here](https://github.com/JuliaReach/LazySets.jl#blue_book-publications).

Is there a chance of getting this function in `master` but with another name, eg. `fast_pow(a, n)`?

I thought this branch gives a faster `pow` than master but it hasn't been resolved yet because it sacrifices tightness for speed. On `master` i measure ~1us for `Interval(0, 1)^2`,...

.. so if that's the case i was just proposing to export the "fast" version, possibly with a different name, keeping the slow one as default.