diff --git a/docs/src/tutorials/bvp_example.md b/docs/src/tutorials/bvp_example.md
index 5402b9733..7c2881136 100644
--- a/docs/src/tutorials/bvp_example.md
+++ b/docs/src/tutorials/bvp_example.md
@@ -145,10 +145,10 @@ Consider a semi-explicit boundary value differential-algebraic equation formulat
 
 ```math
 \begin{cases}
-x_1'=(ε+x_2-p_2(t))y+p_1'(t) \\
-x_2'=p_2'(t) \\
-x_3'=y \\
-0=(x_1-p_1(t))(y-e^t)
+x_1' = \left(ε + x_2 - \sin(t)\right) y + \cos(t) \\
+x_2' = \cos(t) \\
+x_3' = y \\
+0 = \left(x_1 - \sin(t)\right) \left(y - e^t\right)
 \end{cases}
 ```
 
@@ -180,7 +180,6 @@ u0 = [0.0, 0.0, 0.0, 0.0]
 tspan = (0.0, 1.0)
 fun = BVP.BVPFunction(f!, bc!, mass_matrix = [1 0 0 0; 0 1 0 0; 0 0 1 0; 0 0 0 0])
 prob = BVP.BVProblem(fun, u0, tspan)
-sol = BVP.solve(prob,
-    BVP.Ascher4(; zeta = [0.0, 0.0, 1.0], jac_alg = BVP.BVPJacobianAlgorithm(BVP.AutoForwardDiff()));
-    dt = 0.01)
+solver = BVP.Ascher4(; zeta = [0.0, 0.0, 1.0], jac_alg = BVP.BVPJacobianAlgorithm(BVP.AutoForwardDiff()))
+sol = BVP.solve(prob, solver; dt = 0.01)
 ```