diff --git a/docs/src/manual/solving_harmonics.md b/docs/src/manual/solving_harmonics.md
index 0f2712a5..65057e61 100644
--- a/docs/src/manual/solving_harmonics.md
+++ b/docs/src/manual/solving_harmonics.md
@@ -8,7 +8,7 @@ having called `get_harmonic_equations`, we need to set all time-derivatives to z
 Once defined, a `Problem` can be solved for a set of input parameters using `get_steady_states` to obtain `Result`.
 
 ```@docs
-Problem
+HarmonicBalance.Problem
 get_steady_states
 HarmonicBalance.Result
 ```
diff --git a/docs/src/tutorials/classification.md b/docs/src/tutorials/classification.md
index 8996b79a..b79178ec 100644
--- a/docs/src/tutorials/classification.md
+++ b/docs/src/tutorials/classification.md
@@ -20,7 +20,7 @@ We performe a 2d sweep in the driving frequency $\omega$ and driving strength $\
 fixed = (ω₀ => 1.0, γ => 0.002, α => 1.0)
 varied = (ω => range(0.99, 1.01, 100), λ => range(1e-6, 0.03, 100))
 
-result_2D = get_steady_states(harmonic_eq, varied, fixed, threading=true)
+result_2D = get_steady_states(harmonic_eq, varied, fixed)
 ```
 By default the steady states of the system are classified by four different catogaries:
 * `physical`: Solutions that are physical, i.e., all variables are purely real.