Expansion of base-e pentation

[sheldonison]

I updated the pentation.gp code. The first thing the pentation code does is generate the sexp fixed point. The sexp fixed point gets very close to -2 as the base increases, and since there is a singularity at sexp(-2), that was breaking the previous code for bases>3. I decided to patch the slog function, and make a few other changes so that pentation at least works for bases<=10. It also works for all bases>=1.47, bases for which kneser.gp works.

This code includes all of the kneser.gp code, with the penation.gp code in one program, to make things simpler.

I also included the "sexpupfixed" routing, to iteratively generate the sexp base with a parabolic upper fixed point. The "sexpupfixed" routine requires about a dozen iterations of "init(B);loop" to get results accurate to 32 decimal digits, which takes about 15 minutes; I posted results for the value of the fixed point in this thread..
- Sheldon

Code:

\r pentation.gp
init(2);loop;  /* generate kneser.gp sexp series */
genpent;