19 Jan
2007
19 Jan
'07
1:29 p.m.
On Fri, 19 Jan 2007, R. William Gosper wrote:
(c186) subst(%,x,d184)
What does d184 do?
(interesting that it got imagpart exact. Macsyma doesn't, e.g., force 3.1... + .69b-105*%i to real.) But note that the convergence is *cubic* if we just subst realparts: (c190) subst(realpart(d185),x,d184)
Does this just throw away the imaginary part at each iteration?
x 1 (d245) (%i - -) (----------- + 1) + x n %i x n (---- + 1) n
Am I reading it right? (i-x/n) * (1/ (i*x/n +1)^n +1) + x -J