Is there a good way to find fixed points in gp?
I found some code (simplified version below) which finds them by direct iteration, but I have to wonder if there's a more elegant way, either in general or in this particular case:
li(z)=-eint1(-log(z));
Eps_= 4.0*10.0^(-default(realprecision));
z=1+I; zlast=0;
for(k=1, 1e6, z=li(z); if(abs(z-zlast)<Eps_, break); zlast=z)
Charles Greathouse
Case Western Reserve University