04/21/2015, 08:23 AM
(This post was last modified: 04/21/2015, 08:26 AM by sheldonison.)
(04/21/2015, 03:20 AM)marraco Wrote: ...The kneser.gp program will only generate useful results for you for real bases>exp(1/e). It will take a couple of weeks of work and experimenting to figure out how to rewrite the algorithms to work for bases between 0 and 1, which is actually a different function (in many ways) than tetration for real bases>eta. I haven't done that yet; and things are very busy now
I get crashes for almost all bases
[code]
(23:15) gp > init(0.5)
Try "init(exp(1));"
or init(1.5); /* runs a little slower for bases closer to exp(1/e) */
or init(2);
or init(3);
or init(4);
As I remember, the program will generate results for bases between 1 and eta as well, but what it does for bases between 1 and eta is not what you might expect it to be; so I would recommend avoiding those bases as well. Kneser.gp accurately generates Kneser's merged tetration solution from both complex repelling fixed points for most real bases>eta, for much of the complex plane, unless the iterated function overflows. pari-gp's internal representation has limitations, and sexp(10) is a humongous number. The largest number pari-gp can represent is about \( \text{sexp}_e\left(4.08\right)\approx 7.24676\cdot 10^{131475196{ \)
- Sheldon