Specify the correct jac_alg

Author: ErikQQY

docs/src/tutorials/bvp_example.md

@@ -132,7 +132,7 @@ end
 u0 = [1.0, 1.0, 1.0]
 tspan = (0.0, 1.0)
 prob = SecondOrderBVProblem(f!, bc!, u0, tspan)
-sol = solve(prob, MIRKN4(), dt = 0.01)
+sol = solve(prob, MIRKN4(; jac_alg = BVPJacobianAlgorithm(diffmode = AutoForwardDiff())), dt = 0.01)
 ```
 
 ## Example 3: Semi-Explicit Boundary Value Differential-Algebraic Equations
@@ -173,5 +173,5 @@ u0 = [0.0, 0.0, 0.0, 0.0]
 tspan = (0.0, 1.0)
 fun = BVPFunction(f!, bc!, mass_matrix = [1 0 0 0; 0 1 0 0; 0 0 1 0; 0 0 0 0])
 prob = BVProblem(fun, u0, tspan)
-sol = solve(prob, Ascher4(zeta = [0.0, 0.0, 1.0]), dt = 0.01)
+sol = solve(prob, Ascher4(zeta = [0.0, 0.0, 1.0], jac_alg = BVPJacobianAlgorithm(diffmode = AutoForwardDiff())), dt = 0.01)
 ```