diff --git a/src/Flowpipes/solutions.jl b/src/Flowpipes/solutions.jl
index cb1c5b6fd1..8d21f9b2a9 100644
--- a/src/Flowpipes/solutions.jl
+++ b/src/Flowpipes/solutions.jl
@@ -100,6 +100,7 @@ Base.getindex(sol::ReachSolution, I::AbstractVector) = getindex(sol.F, I)
 
 # indexing: fp[j, i] returning the j-th reach-set of the i-th flowpipe
 Base.getindex(sol::ReachSolution, I::Int...) = getindex(sol.F, I...)
+Base.eachindex(sol::ReachSolution) = eachindex(sol.F)
 
 # evaluation interface
 Base.getindex(sol::ReachSolution, t::Float64) = getindex(sol.F, t)