Close issue #188 (#189)

* rm duplicate link * fix bug in A42 method; close issue #188
JuliaMath · Jul 6, 2020 · 183193b · 183193b · jverzani · Jul 6, 2020
1 parent bb11bf7
commit 183193b
Show file tree

Hide file tree

Showing 4 changed files with 8 additions and 5 deletions.
diff --git a/Project.toml b/Project.toml
@@ -1,6 +1,6 @@
 name = "Roots"
 uuid = "f2b01f46-fcfa-551c-844a-d8ac1e96c665"
-version = "1.0.2"
+version = "1.0.3"
 
 [deps]
 Printf = "de0858da-6303-5e67-8744-51eddeeeb8d7"

diff --git a/README.md b/README.md
@@ -2,9 +2,6 @@
 Linux: [![Build Status](https://travis-ci.org/JuliaMath/Roots.jl.svg?branch=master)](https://travis-ci.org/JuliaMath/Roots.jl)
 Windows: [![Build status](https://ci.appveyor.com/api/projects/status/goteuptn5kypafyl?svg=true)](https://ci.appveyor.com/project/jverzani/roots-jl)
 
-https://juliahub.com/docs/Roots/
-Documentation
-
 # Root finding functions for Julia
 
 

diff --git a/src/bracketing.jl b/src/bracketing.jl
@@ -669,7 +669,7 @@ function update_state(M::A42, f, state::UnivariateZeroState{T,S}, options::Univa
     if abs(cb - u) > 0.5 * (b-a)
         ch = _middle(an, bn)
     end
-    fch::S = f(cb)
+    fch::S = f(ch)
     incfn(state)
     if check_zero(M, state, ch, fch)
         # tighten up bracket

diff --git a/test/test_find_zero.jl b/test/test_find_zero.jl
@@ -286,6 +286,12 @@ end
 
     ## issue #178 passinig through method
     @test fzero(sin, 3, 4, Roots.Brent()) ≈ π
+
+    ## issue #188 with A42
+    f(x) = x*(1-x^2)/((x^2+a^2)*(1+a^2*x^2))
+    a,r = .18, 0.05
+    xs = (r + 1e-12, 1.0)
+    @test find_zero(x -> f(r)-f(x), xs, Roots.A42()) ≈ 0.4715797678171889
 end