CompatHelper: bump compat for NonlinearSolve to 4 for package docs, (keep existing compat)

Author: CompatHelper Julia

docs/Project.toml

@@ -21,8 +21,8 @@ FiniteDiff = "6a86dc24-6348-571c-b903-95158fe2bd41"
 ForwardDiff = "f6369f11-7733-5829-9624-2563aa707210"
 IncompleteLU = "40713840-3770-5561-ab4c-a76e7d0d7895"
 JLD2 = "033835bb-8acc-5ee8-8aae-3f567f8a3819"
-NonlinearSolve = "8913a72c-1f9b-4ce2-8d82-65094dcecaec"
 ModelingToolkit = "961ee093-0014-501f-94e3-6117800e7a78"
+NonlinearSolve = "8913a72c-1f9b-4ce2-8d82-65094dcecaec"
 ODEProblemLibrary = "fdc4e326-1af4-4b90-96e7-779fcce2daa5"
 Optimization = "7f7a1694-90dd-40f0-9382-eb1efda571ba"
 OptimizationNLopt = "4e6fcdb7-1186-4e1f-a706-475e75c168bb"
@@ -60,8 +60,8 @@ FiniteDiff = "2"
 ForwardDiff = "0.10"
 IncompleteLU = "0.2"
 JLD2 = "0.4, 0.5.1"
-NonlinearSolve = "3.15"
 ModelingToolkit = "9"
+NonlinearSolve = "3.15, 4"
 ODEProblemLibrary = "0.1"
 Optimization = "3, 4"
 OptimizationNLopt = "0.2, 0.3"