erfinv - once more

[G. B|rger]

Hi octave,

I ran into another strange behavior of the erfinv function. Using the
values given below the algorithm enters an infinit loop (using
octave-2.1.34).

This only happens on my office rs2k aix machine. At home using
linux/i386 the error is gone.

           Gerd


------- script --------
format long;
n = 360; k=84;
ww = (k-1)/(n-1);
y = ww*(k-1)/n + (1-ww)*k/n;
a = erfinv(2*y-1);
------- script --------

------- output (with z_new from erfinv) --------
z_new = -0.516197146962978
z_new = -0.516197146962979
z_new = -0.516197146962978
z_new = -0.516197146962979
...
warning: erfinv: iteration limit exceeded
------- output (with z_new from erfinv) --------



-------------------------------------------------------------
Octave is freely available under the terms of the GNU GPL.

Octave's home on the web:  http://www.octave.org
How to fund new projects:  http://www.octave.org/funding.html
Subscription information:  http://www.octave.org/archive.html
-------------------------------------------------------------