diff --git a/test/Project.toml b/test/Project.toml
index 50b045e2..754ee9f0 100644
--- a/test/Project.toml
+++ b/test/Project.toml
@@ -1,5 +1,6 @@
 [deps]
 Aqua = "4c88cf16-eb10-579e-8560-4a9242c79595"
+ExplicitImports = "7d51a73a-1435-4ff3-83d9-f097790105c7"
 LinearAlgebra = "37e2e46d-f89d-539d-b4ee-838fcccc9c8e"
 Pkg = "44cfe95a-1eb2-52ea-b672-e2afdf69b78f"
 StaticArrays = "90137ffa-7385-5640-81b9-e52037218182"
@@ -11,6 +12,7 @@ Test = "8dfed614-e22c-5e08-85e1-65c5234f0b40"
 
 [compat]
 Aqua = "0.8"
+ExplicitImports = "1"
 StaticArrays = "1"
 SymEngine = "0.8, 0.9.1, 0.10, 0.11, 0.12, 0.13"
 SymPy = "2"
diff --git a/test/runtests.jl b/test/runtests.jl
index 3c558343..ee15d96b 100644
--- a/test/runtests.jl
+++ b/test/runtests.jl
@@ -18,6 +18,7 @@ Pkg.build("SymPy")
 using SymPy: SymPy
 
 using Aqua: Aqua
+using ExplicitImports: check_no_implicit_imports, check_no_stale_explicit_imports
 
 @testset "BSeries" begin
     @testset "lazy representation of exact ODE solution" begin
@@ -2678,4 +2679,10 @@ using Aqua: Aqua
             Aqua.test_ambiguities([BSeries])
         end
     end
+
+    @testset "ExplicitImports" begin
+        @test check_no_implicit_imports(BSeries) === nothing
+
+        @test check_no_stale_explicit_imports(BSeries) === nothing
+    end
 end # @testset "BSeries"