diff --git a/Project.toml b/Project.toml
index 3eaaa1e0..86603445 100644
--- a/Project.toml
+++ b/Project.toml
@@ -1,6 +1,6 @@
 name = "BoundaryValueDiffEq"
 uuid = "764a87c0-6b3e-53db-9096-fe964310641d"
-version = "5.6.0"
+version = "5.6.1"
 
 [deps]
 ADTypes = "47edcb42-4c32-4615-8424-f2b9edc5f35b"
@@ -42,7 +42,7 @@ Aqua = "0.7"
 ArrayInterface = "7"
 BandedMatrices = "1"
 ConcreteStructs = "0.2"
-DiffEqBase = "6.138"
+DiffEqBase = "6.138 - 6.143"
 FastAlmostBandedMatrices = "0.1"
 ForwardDiff = "0.10"
 LinearAlgebra = "1.9"
diff --git a/src/sparse_jacobians.jl b/src/sparse_jacobians.jl
index a5225d73..1983d010 100644
--- a/src/sparse_jacobians.jl
+++ b/src/sparse_jacobians.jl
@@ -57,7 +57,7 @@ function __generate_sparse_jacobian_prototype(::MIRKCache, ::StandardBVProblem,
         N)
     fast_scalar_indexing(ya) ||
         error("Sparse Jacobians are only supported for Fast Scalar Index-able Arrays")
-    J_c = BandedMatrix(Ones{eltype(ya)}(M * (N - 1), M * N), (1, 2M))
+    J_c = BandedMatrix(Ones{eltype(ya)}(M * (N - 1), M * N), (1, 2M - 1))
     return ColoredMatrix(J_c, matrix_colors(J_c'), matrix_colors(J_c))
 end