try fix 1.6 again

Author: TorkelE

test/miscellaneous_tests/nonlinear_solve.jl

@@ -1,7 +1,7 @@
 ### Fetch Packages ###
 
 # Fetch packages.
-using Catalyst, NonlinearSolve, OrdinaryDiffEq, SteadyStateDiffEq, DifferentialEquations
+using Catalyst, NonlinearSolve, OrdinaryDiffEq, SteadyStateDiffEq
 using Random, Test
 
 # Sets rnd number.
@@ -25,7 +25,7 @@ let
     nl_prob = NonlinearProblem(steady_state_network_1, u0, p)
 
     # Solves it using standard algorithm and simulation based algorithm.
-    sol1 = solve(nl_prob; abstol=1e-12, reltol=1e-12).u
+    sol1 = solve(nl_prob, LevenbergMarquardt(); abstol=1e-12, reltol=1e-12).u
     sol2 = solve(nl_prob, DynamicSS(Rosenbrock23(); abstol=1e-12, reltol=1e-12); abstol=1e-12, reltol=1e-12).u
 
     # Tests solutions are correct.
@@ -54,7 +54,7 @@ let
     nl_prob = NonlinearProblem(steady_state_network_2, u0, p)
 
     # Solves it using standard algorithm and simulation based algorithm.
-    sol1 = solve(nl_prob; abstol=1e-12, reltol=1e-12).u
+    sol1 = solve(nl_prob, LevenbergMarquardt(); abstol=1e-12, reltol=1e-12).u
     sol2 = solve(nl_prob, DynamicSS(Rosenbrock23(); abstol=1e-12, reltol=1e-12); abstol=1e-12, reltol=1e-12).u
 
     # Computes NonlinearFunction (manually and automatically).
@@ -90,7 +90,7 @@ let
     nl_prob_2 = NonlinearProblem(steady_state_network_3, u0, p)
 
     # Solves it using standard algorithm and simulation based algorithm.
-    sol1 = solve(nl_prob_1; abstol=1e-12, reltol=1e-12)
+    sol1 = solve(nl_prob_1, LevenbergMarquardt(); abstol=1e-12, reltol=1e-12)
     sol2 = solve(nl_prob_2, DynamicSS(Rosenbrock23(); abstol=1e-12, reltol=1e-12); abstol=1e-12, reltol=1e-12)
 
     # Checks output using NonlinearFunction.