Somehow I lost the original post.
For some reason Laguerre's method doesn't converge for this polynomial.
There seem to be large regions of non-convergence. Changes of up to 1e_5
in any of the constants still fail to converge.
The point seems to oscillate in a wide range but doesn't head toward any
solution. In the implementation the initial guess is always 0j0.
I don't see why this would fail. If anyone on this list can help, or
knows someone who can help, I'd appreciate suggestions.
Henry Rich
----------------------------------------------------------------------
For information about J forums see http://www.jsoftware.com/forums.htm