Thread Rating:
  • 0 Vote(s) - 0 Average
  • 1
  • 2
  • 3
  • 4
  • 5
Expansion of base-e pentation
#12
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;



Attached Files
.gp   pentation.gp (Size: 27.26 KB / Downloads: 639)
Reply


Messages In This Thread
Expansion of base-e pentation - by andydude - 10/24/2009, 04:30 AM
RE: Expansion of base-e pentation - by bo198214 - 10/24/2009, 10:13 AM
RE: Expansion of base-e pentation - by andydude - 10/24/2009, 10:49 AM
RE: Expansion of base-e pentation - by andydude - 10/25/2009, 07:14 AM
RE: Expansion of base-e pentation - by andydude - 10/29/2009, 11:36 AM
RE: Expansion of base-e pentation - by andydude - 10/31/2009, 09:47 AM
RE: Expansion of base-e pentation - by sheldonison - 11/02/2010, 03:37 PM

Possibly Related Threads...
Thread Author Replies Views Last Post
  Precision check on [pentation.gp] SOLVED Cherrina_Pixie 7 17,763 07/02/2011, 01:39 AM
Last Post: Cherrina_Pixie



Users browsing this thread: 1 Guest(s)