diff --git a/docs/make.jl b/docs/make.jl
index 2185529..3cca0f6 100644
--- a/docs/make.jl
+++ b/docs/make.jl
@@ -4,4 +4,4 @@ using Documenter, Looping
 DocMeta.setdocmeta!(Looping.Utilities, :DocTestSetup, :(using Looping.Utilities); recursive=true)
 DocMeta.setdocmeta!(Looping.FloquetUtils, :DocTestSetup, :(using Looping.FloquetUtils); recursive=true)
 
-makedocs(sitename="Looping Documentation", modules=map(eval, modules), draft=false, strict=:doctest)
+makedocs(sitename="Looping Documentation", modules=[Looping.Utilities, Looping.FloquetUtils], draft=false, strict=:doctest)
diff --git a/docs/src/index.md b/docs/src/index.md
index bc0c775..210ded9 100644
--- a/docs/src/index.md
+++ b/docs/src/index.md
@@ -2,10 +2,3 @@
 
 Some numerical tools and experiments for the fiber loop reservoir
 engineering project.
-
-
-
-```@contents
-Pages = ["modules/Utilities.md" "modules/FloquetUtilities.md"]
-Depth = 3
-```
diff --git a/docs/src/modules/FloquetUtils.md b/docs/src/modules/FloquetUtils.md
index c3d52a4..3207b79 100644
--- a/docs/src/modules/FloquetUtils.md
+++ b/docs/src/modules/FloquetUtils.md
@@ -1,19 +1,5 @@
 # Floquet Utilities
 
-```@docs
-Looping.FloquetUtils
-```
-
-
-```@meta
-CurrentModule = Looping.FloquetUtils
-```
-
-
-
-```@docs
-time_evolution_op
-floquet_hamiltonian
-trivial_floquet_hamiltonian
-KickOperator
+```@autodocs
+Modules = [Looping.FloquetUtils]
 ```
diff --git a/docs/src/modules/Utilities.md b/docs/src/modules/Utilities.md
index f579146..282f507 100644
--- a/docs/src/modules/Utilities.md
+++ b/docs/src/modules/Utilities.md
@@ -1,17 +1,5 @@
 # Utilities
 
-```@docs
-Looping.Utilities
-```
-
-
-```@meta
-CurrentModule = Looping.Utilities
-```
-
-
-
-```@docs
-periodic_distance
-restrict_to_range
+```@autodocs
+Modules = [Looping.Utilities]
 ```
diff --git a/src/FloquetUtils.jl b/src/FloquetUtils.jl
index 286c534..a84790a 100644
--- a/src/FloquetUtils.jl
+++ b/src/FloquetUtils.jl
@@ -19,7 +19,7 @@ using ..Utilities
 
 Calculate the time evolution operator for a Hamiltonian `H` up to a
 total time `T`. The rest arguments are passed on to
-`DifferentialEquations.solve`.
+[`de.solve`](@ref).
 """
 function time_evolution_op(H::Function, T::Real, args...; kwargs...)
     function u_step!(dU, U, _, t)
@@ -52,7 +52,7 @@ end
 
 Returns the Floquet Hamiltonian given a Hamiltionian `H` and a time
 `T`. The rest arguments are passed on to
-[`DifferentialEquations.solve`](@ref).
+[`de.solve`](@ref).
 """
 function floquet_hamiltonian(H::Function, T::Real, args...; kwargs...)
     U = time_evolution_op(H, T, args...; kwargs...)
@@ -93,12 +93,12 @@ struct KickOperator
     H_F::Matrix
 
 
-    """
-        KickOperator(U(t), H_F)
+"""
+    KickOperator(U(t), H_F)
 
-    Return the Kick operator Given the time evolution operator `U(t)` and the Floquet Hamiltonian
-    `H_F`.
-    """
+Return the Kick operator Given the time evolution operator `U(t)` and the Floquet Hamiltonian
+`H_F`.
+"""
     function KickOperator(U, H_F::Matrix)
         if size(U(0)) != size(H_F)
             throw(DimensionMismatch("`U` and `H_F` should have the same dimension."))
@@ -107,7 +107,11 @@ struct KickOperator
     end
 end
 
+"""
+    KickOperator(H(t), T)
 
+Return the Kick operator Given the Hamiltonian `H(t)` and the total time ``T``.
+"""
 function KickOperator(H::Function, T::Real)
     U = time_evolution_op(H, T)
     H_F = floquet_hamiltonian(U(T), T)
diff --git a/src/Looping.jl b/src/Looping.jl
index 0a7677c..2a57df1 100644
--- a/src/Looping.jl
+++ b/src/Looping.jl
@@ -1,8 +1,4 @@
 module Looping
-export hi
-
-hi() = print("hi")
-
 include("Utilities.jl")
 include("FloquetUtils.jl")
 end