diff --git a/docs/Project.toml b/docs/Project.toml
index 79cd1f4..be72387 100644
--- a/docs/Project.toml
+++ b/docs/Project.toml
@@ -1,5 +1,4 @@
 [deps]
-CTFlows = "1c39547c-7794-42f7-af83-d98194f657c2"
 Documenter = "e30172f5-a6a5-5a46-863b-614d45cd2de4"
 LinearAlgebra = "37e2e46d-f89d-539d-b4ee-838fcccc9c8e"
 NLPModelsIpopt = "f4238b75-b362-5c4c-b852-0801c9a21d71"
@@ -10,12 +9,11 @@ Plots = "91a5bcdd-55d7-5caf-9e0b-520d859cae80"
 Suppressor = "fd094767-a336-5f1f-9728-57cf17d0bbfb"
 
 [compat]
-CTFlows = "0.4.5"
-Documenter = "1.5.0"
-NLPModelsIpopt = "0.10.2"
-NonlinearSolve = "3.13.1"
-OptimalControl = "0.9.4"
-OrdinaryDiffEq = "6.85.0"
-Plots = "1.40.5"
-Suppressor = "0.2.7"
+Documenter = "1.5"
+NLPModelsIpopt = "0.10"
+NonlinearSolve = "3.13"
+OptimalControl = "0.11"
+OrdinaryDiffEq = "6.88"
+Plots = "1.40"
+Suppressor = "0.2"
 julia = "1.10"
diff --git a/docs/src/saturation.md b/docs/src/saturation.md
index 3dce9a2..dd2739f 100644
--- a/docs/src/saturation.md
+++ b/docs/src/saturation.md
@@ -179,14 +179,9 @@ We define a custom plot function to plot the solution inside the Bloch ball.
     <code>spin_plot</code> — <span class="docstring-category">Function</span>
 </header>
 <section style="display: none;"><div><pre><code class="language-julia hljs">using Plots.PlotMeasures
-using CTFlows
 
 function spin_plot(sol; kwargs...)
 
-    if sol isa CTFlows.OptimalControlFlowSolution
-        sol = CTFlows.OptimalControlSolution(sol)
-    end
-
     y2 = cos(asin(zs))
     y1 = -y2
 
@@ -231,13 +226,8 @@ end</code><button class="copy-button fa-solid fa-copy" aria-label="Copy this cod
 
 ```@example main
 using Plots.PlotMeasures # hide
-using CTFlows # hide
 function spin_plot(sol; kwargs...) # hide
 
-    if sol isa CTFlows.OptimalControlFlowSolution # hide
-        sol = CTFlows.OptimalControlSolution(sol) # hide
-    end # hide
-
     y2 = cos(asin(zs)) # hide
     y1 = -y2 # hide
 
@@ -376,10 +366,6 @@ and its derivative along the solution computed by the direct method.
 </header>
 <section style="display: none;"><div><pre><code class="language-julia hljs">function switching_plot(sol, H1, H01; kwargs...)
 
-    if sol isa CTFlows.OptimalControlFlowSolution
-        sol = CTFlows.OptimalControlSolution(sol)
-    end
-
     t  = sol.times
     u  = sol.control
     q  = sol.state
@@ -427,10 +413,6 @@ the first singular arc, the switching function is not always zero.
 ```@example main
 function switching_plot(sol, H1, H01; kwargs...) # hide
 
-    if sol isa CTFlows.OptimalControlFlowSolution # hide
-        sol = CTFlows.OptimalControlSolution(sol) # hide
-    end # hide
-
     t  = sol.times  # hide
     u  = sol.control # hide
     q  = sol.state # hide