diff --git a/Project.toml b/Project.toml
index d86a3434..8ae47c1a 100644
--- a/Project.toml
+++ b/Project.toml
@@ -28,7 +28,7 @@ Requires = "1.0"
 StaticArraysCore = "1.1"
 SymbolicIndexingInterface = "0.1, 0.2"
 Tables = "1"
-Zygote = "< 0.6.56"
+Zygote = "0.6.56"
 julia = "1.6"
 
 [extensions]
@@ -56,4 +56,4 @@ test = ["SafeTestsets", "Aqua", "ForwardDiff", "LabelledArrays", "NLsolve", "Ord
 
 [weakdeps]
 Tracker = "9f7883ad-71c0-57eb-9f7f-b5c9e6d3789c"
-Zygote = "e88e6eb3-aa80-5325-afca-941959d7151f"
\ No newline at end of file
+Zygote = "e88e6eb3-aa80-5325-afca-941959d7151f"
diff --git a/src/vector_of_array.jl b/src/vector_of_array.jl
index e734125c..88b88c6e 100644
--- a/src/vector_of_array.jl
+++ b/src/vector_of_array.jl
@@ -353,6 +353,11 @@ function Base.copy(VA::AbstractDiffEqArray)
 end
 Base.copy(VA::AbstractVectorOfArray) = typeof(VA)(copy(VA.u))
 Base.sizehint!(VA::AbstractVectorOfArray{T, N}, i) where {T, N} = sizehint!(VA.u, i)
+
+Base.reverse!(VA::AbstractVectorOfArray) = reverse!(VA.u)
+Base.reverse(VA::VectorOfArray) = VectorOfArray(reverse(VA.u))
+Base.reverse(VA::DiffEqArray) = DiffEqArray(reverse(VA.u), VA.t, VA.sc, VA.observed, VA.p)
+
 function Base.push!(VA::AbstractVectorOfArray{T, N}, new_item::AbstractArray) where {T, N}
     push!(VA.u, new_item)
 end