# HG changeset patch # User John W. Eaton # Date 1396731795 14400 # Node ID 6fdd3ab55b7823b4291e5f21f5f3dad21229cf43 # Parent c644ed73c6cec0af602641cd53f8f8e70887db06 undo unintended changes to fplot in changeset 37c300acfcfd diff -r c644ed73c6ce -r 6fdd3ab55b78 scripts/plot/draw/fplot.m --- a/scripts/plot/draw/fplot.m Sat Apr 05 15:30:30 2014 -0400 +++ b/scripts/plot/draw/fplot.m Sat Apr 05 17:03:15 2014 -0400 @@ -146,10 +146,10 @@ ## problems with the current solution. while (n < 2^18) # Something is wrong if we need more than 250K points - yi = interp1 (x0, y0, x, "linear") + yi = interp1 (x0, y0, x, "linear"); ## relative error calculation using average of [yi,y] as reference ## since neither estimate is known a priori to be better than the other. - err = 0.5 * max (abs ((yi - y) ./ (yi + y))(:)) + err = 0.5 * max (abs ((yi - y) ./ (yi + y))(:)); if (err < tol || abs (err - err0) < tol/2) ## Either relative tolerance has been met OR ## algorithm has stopped making any reasonable progress per iteration.