diff --git a/Project.toml b/Project.toml
index ce120d7d..1673a4ec 100644
--- a/Project.toml
+++ b/Project.toml
@@ -9,7 +9,7 @@ Reexport = "189a3867-3050-52da-a836-e630ba90ab69"
 SparseArrays = "2f01184e-e22b-5df5-ae63-d93ebab69eaf"
 
 [compat]
-ReachabilityAnalysis = "0.18.6 - 0.28"
+ReachabilityAnalysis = "0.18.6 - 0.30"
 ReachabilityBase = "0.2.3 - 0.3"
 Reexport = "1"
 SparseArrays = "<0.0.1, 1.6"
diff --git a/docs/Project.toml b/docs/Project.toml
index cf39636f..4e031203 100644
--- a/docs/Project.toml
+++ b/docs/Project.toml
@@ -14,5 +14,5 @@ Literate = "2"
 MAT = "0.10 - 0.11"
 ModelingToolkit = "8 - 11"
 Plots = "1"
-ReachabilityAnalysis = "0.22 - 0.28"
+ReachabilityAnalysis = "0.22 - 0.30"
 ReachabilityBase = "0.2.3 - 0.3"
diff --git a/src/models/brusselator/brusselator_meta.jl b/src/models/brusselator/brusselator_meta.jl
index a2c28501..fc61a177 100644
--- a/src/models/brusselator/brusselator_meta.jl
+++ b/src/models/brusselator/brusselator_meta.jl
@@ -1,6 +1,6 @@
 using ReachabilityAnalysis
 Dict(["opts" => Dict(:T => 18.0, :alg => TMJets(; orderT=6, orderQ=2)),
-      "X0"   => (0.8 .. 1.0) × (0.0 .. 0.2),
+      "X0"   => BallInf([1.0, 0.1], 0.1),
       "info" => Dict("name" => "brusselator",
       "dim" => 2,
       "linear" => false,
diff --git a/src/models/roessler/roessler_meta.jl b/src/models/roessler/roessler_meta.jl
index 2a7e5902..c16dbee4 100644
--- a/src/models/roessler/roessler_meta.jl
+++ b/src/models/roessler/roessler_meta.jl
@@ -4,7 +4,7 @@ x₀ = Interval(-0.2, 0.2)
 y₀ = Interval(-8.6, -8.2)
 z₀ = Interval(-0.2, 0.2)
 
-X₀ = x₀ × y₀ × z₀;
+X₀ = Hyperrectangle([0, -8.4, 0], [0.2, 0.2, 0.2]);
 
 Dict(["opts" => Dict(:T => 6.0, :alg => TMJets(; abstol=1e-10, orderT=6, orderQ=2, maxsteps=1_000)),
       "X0"   => X₀,
diff --git a/src/models/spring_pendulum/spring_pendulum_meta.jl b/src/models/spring_pendulum/spring_pendulum_meta.jl
index bed8ed2e..8284c1fa 100644
--- a/src/models/spring_pendulum/spring_pendulum_meta.jl
+++ b/src/models/spring_pendulum/spring_pendulum_meta.jl
@@ -5,7 +5,7 @@ r₀ = Interval(1.19, 1.21)
 vr₀ = Singleton([0.0])
 vθ₀ = Singleton([0.0])
 
-X₀ = r₀ × θ₀ × vr₀ × vθ₀;
+X₀ = Hyperrectangle([1.2, 0.5, 0, 0], [0.01, 0.01, 0, 0]);
 
 Dict(["opts" => Dict(:T => 5.0, :alg => TMJets(; orderT=8, orderQ=2, maxsteps=1_000)),
       "X0"   => X₀,
diff --git a/src/models/steam_governor/steam_governor_meta.jl b/src/models/steam_governor/steam_governor_meta.jl
index 2f1c63f3..02594ccd 100644
--- a/src/models/steam_governor/steam_governor_meta.jl
+++ b/src/models/steam_governor/steam_governor_meta.jl
@@ -1,10 +1,6 @@
 using ReachabilityAnalysis
 
-x₀ = Interval(0.9, 1.1)
-y₀ = Interval(-0.1, 0.1)
-z₀ = Interval(0.9, 1.1)
-
-X₀ = x₀ × y₀ × z₀;
+X₀ = Hyperrectangle([1.0, 0, 1.0], [0.1, 0.1, 0.1]);
 
 Dict(["opts" => Dict(:T => 10.0, :alg => TMJets(; abstol=1e-10, orderT=6, orderQ=2, maxsteps=1_000)),
       "X0"   => X₀,