From help-octave-request at bevo dot che dot wisc dot edu Fri Jan 11 09:37:42 2002
Subject: erfinv - once more
From: "G. B|rger" <b280002 at niesel dot dkrz dot de>
To: help-octave at bevo dot che dot wisc dot edu
Date: Fri, 11 Jan 2002 09:37:15 -0600


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
-------------------------------------------------------------