use latest IntervalArithmetic version in docs

Author: schillic

docs/Project.toml

@@ -1,9 +1,7 @@
 [deps]
 Documenter = "e30172f5-a6a5-5a46-863b-614d45cd2de4"
 DocumenterCitations = "daee34ce-89f3-4625-b898-19384cb65244"
-IntervalArithmetic = "d1acc4aa-44c8-5952-acd4-ba5d80a2a253"
 
 [compat]
 Documenter = "1"
 DocumenterCitations = "1.3"
-IntervalArithmetic = "0.21, =0.21"  # v0.22 changed the output format

docs/src/index.md

@@ -108,17 +108,17 @@ Then, at ``t = 1.0``,
 ```jldoctest quickstart
 julia> A + 1/2 * A^2
 2×2 IntervalMatrix{Float64, Interval{Float64}, Matrix{Interval{Float64}}}:
-  [1.0, 4.50001]  [-3.0, 2.0]
- [-4.0, 2.50001]  [-1.0, 9.0]
+  [1.0, 4.5]   [-3.25, 2.5]
+ [-4.25, 3.0]  [-2.25, 9.0]
 ```
 However, that result is not tight. The computation can be performed exactly via
 single-use expressions implemented in this library:
 
 ```jldoctest quickstart
 julia> quadratic_expansion(A, 1.0, 0.5)
 2×2 IntervalMatrix{Float64, Interval{Float64}, Matrix{Interval{Float64}}}:
-  [1.0, 4.50001]  [-2.0, 1.0]
- [-3.0, 1.50001]   [1.0, 7.0]
+  [1.0, 4.5]  [-2.0, 1.0]
+ [-3.0, 1.5]   [1.0, 7.0]
 ```
 We now obtain an interval matrix that is strictly included in the one obtained from
 the naive multiplication.

src/matrix.jl

@@ -18,10 +18,10 @@ parameterized in the number field, the interval type, and the matrix type.
 ### Examples
 
 ```jldoctest
-julia> A = IntervalMatrix([-1 .. -0.8 0 .. 0; 0 .. 0 -1 .. -0.8])
+julia> A = IntervalMatrix([interval(-1, -0.8) interval(0); interval(0) interval(-1, -0.8)])
 2×2 IntervalMatrix{Float64, Interval{Float64}, Matrix{Interval{Float64}}}:
- [-1.0, -0.7999999]   [0.0, 0.0]
-  [0.0, 0.0]         [-1.0, -0.7999999]
+ [-1.0, -0.8]   [0.0, 0.0]
+  [0.0, 0.0]   [-1.0, -0.8]
 ```
 
 An interval matrix proportional to the identity matrix can be built using the