Improve the documentation for `pow`

[schillic]

julia> using IntervalArithmetic

julia> x = interval(-4, -2)
[-4.0, -2.0]_com

julia> pow(x, 2)
∅_trv

The problem seems to be in these lines:

IntervalArithmetic.jl/src/intervals/arithmetic/power.jl

Lines 157 to 159 in 544c000

	domain = _unsafe_bareinterval(T, zero(T), typemax(T))
	x = intersect_interval(x, domain)
	isempty_interval(x) && return emptyinterval(BareInterval{R})

It is unclear to me why domain should be nonnegative. Clearly you can take a power of a negative number.

[OlivierHnt]

Thanks for opening this issue. I agree this may be odd looking; for the reasons behind this behavior refer to Table 9.1 from the IEEE standard.

The documentation can always be improved, and in particular that of pow. So I changed the title this issue accordingly.

Note that pow is not an alias for ^,

julia> interval(-4, -2)^2
[4.0, 16.0]_com

julia> pown(interval(-4, -2), 2)
[4.0, 16.0]_com

julia> pow(interval(-4, -2), 2)
∅_trv