Mercurial > hg > octave-lyh

diff scripts/optimization/fzero.m @ 10394:bc475cd49147

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
avoid div by zero in fzero
author Jaroslav Hajek <highegg@gmail.com>
date Fri, 05 Mar 2010 07:31:19 +0100
parents 7c97da90fc8f
children 9500a66118dc
line wrap: on
line diff
--- a/scripts/optimization/fzero.m
+++ b/scripts/optimization/fzero.m
@@ -287,7 +287,7 @@
 
   ## Check solution for a singularity by examining slope
   if (info == 1)
-    if ((b - a) != 0 && abs ((fb - fa)/(b - a) / slope0) > max (1e6, 0.5/tolx))
+    if ((b - a) != 0 && abs ((fb - fa)/(b - a) / slope0) > max (1e6, 0.5/(eps+tolx)))
       info = -5;
     endif
   endif