diff --git a/docs/src/visualizations.md b/docs/src/visualizations.md
index 10393d7a..244d5cf9 100644
--- a/docs/src/visualizations.md
+++ b/docs/src/visualizations.md
@@ -14,80 +14,17 @@ The following GUI is obtained with the function [`interactive_trajectory_timeser
 </video>
 ```
 
-```julia
-using DynamicalSystems, GLMakie, ModelingToolkit
-
-# Define the variables and parameters in symbolic format
-@variables t
-D = Differential(t)
-@parameters begin
-    a = 0.29
-    b = 0.14
-    c = 4.52
+```@setup MAIN
+using DynamicalSystems
+struct ShowFile
+    file::String
 end
-@variables begin
-    x(t) = 10.0
-    y(t) = 0.0
-    z(t) = 1.0
-    nlt(t) # nonlinear term
+function Base.show(io::IO, ::MIME"text/plain", f::ShowFile)
+    write(io, read(f.file))
 end
-
-# Create the equations of the model
-eqs = [
-    D(x) ~ -y - z,
-    D(y) ~ x + a*y,
-    D(z) ~ b + nlt - z*c,
-    nlt ~ z*x, # observed variable
-]
-
-# Create the model via ModelingToolkit
-@named roessler = ODESystem(eqs)
-# Explicitly disable parameter spliting, so that integer
-# indices can also be used (not recommended in general!)
-model = structural_simplify(roessler; split = false)
-# Cast it into an `ODEProblem` and then into a `DynamicalSystem`
-prob = ODEProblem(model)
-ds = CoupledODEs(prob)
-# If you have "lost" the model, use:
-model = referrenced_sciml_model(ds)
-
-# Define which parameters will be interactive during the simulation
-parameter_sliders = Dict(
-    # can use integer indexing
-    1 => 0:0.01:1,
-    # the global scope symbol
-    b => 0:0.01:1,
-    # or the symbol obtained from the MTK model
-    model.c => 0:0.01:10,
-)
-
-# Define what variables will be visualized as timeseries
-norm(u) = sqrt(u[1]*u[1] + u[2]*u[2])
-observables = [
-    1,         # can use integer indexing,
-    z,         # MTK state variable
-    model.nlt, # MTK observed variable
-    norm,      # or arbitrary function of the state
-]
-
-# Define what variables will be visualized as state space trajectory
-# same as above, any indexing works, but ensure to make the vector `Any`
-# so that integers are not converted to symbolic variables
-idxs = Any[1, y, 3]
-
-u0s = [
-    # no fancy indexing here yet; numbers must correspond to state variables
-    [-4.0, -4, 0.1],
-    [4.0, 4, 0.2],
-    [5.72, 0.28, 0.21],
-    [-5.72, 0.0, 0.0],
-]
-
-fig, dsobs = interactive_trajectory_timeseries(ds, observables, u0s;
-    parameter_sliders, idxs,
-)
-
-display(fig)
+```
+```@example MAIN
+ShowFile(joinpath(dirname(pathof(DynamicalSystems)), "../test/mtk_gui.jl")) # hide
 ```
 
 ```@docs