Merge pull request #122 from SciML/ap/patch

Patch Watson Function
SciML · Dec 25, 2023 · 3d0defd · 3d0defd · avik-pal · Dec 25, 2023
2 parents 1e8c196 + 6052b7f
commit 3d0defd
Show file tree

Hide file tree

Showing 2 changed files with 5 additions and 2 deletions.
diff --git a/lib/NonlinearProblemLibrary/Project.toml b/lib/NonlinearProblemLibrary/Project.toml
@@ -1,12 +1,14 @@
 name = "NonlinearProblemLibrary"
 uuid = "b7050fa9-e91f-4b37-bcee-a89a063da141"
-version = "0.1.1"
+version = "0.1.2"
 
 [deps]
 LinearAlgebra = "37e2e46d-f89d-539d-b4ee-838fcccc9c8e"
 SciMLBase = "0bca4576-84f4-4d90-8ffe-ffa030f20462"
 
 [compat]
+Aqua = "0.8"
+LinearAlgebra = "1.6"
 SciMLBase = "1, 2"
 julia = "1.6"
 

diff --git a/lib/NonlinearProblemLibrary/src/NonlinearProblemLibrary.jl b/lib/NonlinearProblemLibrary/src/NonlinearProblemLibrary.jl
@@ -90,6 +90,7 @@ p5_dict = Dict("n" => n, "start" => x_start, "sol" => x_sol,
 # ------------------------------------- Problem 6 ------------------------------------------
 function p6_f!(out, x, p = nothing)
     n = length(x)
+    out .= 0
     for i in 1:29
         ti = i / 29.0
         sum1 = 0.0
@@ -113,7 +114,7 @@ function p6_f!(out, x, p = nothing)
         end
     end
 
-    out[1] = out[1] + 3.0 * x[1] - 2.0 * x[1] * x[1] + 2.0 * x[1]^3
+    out[1] = out[1] + 3.0 * x[1] - 2.0 * x[1] * x[2] + 2.0 * x[1]^3
     out[2] = out[2] + x[2] - x[2]^2 - 1.0
     nothing
 end