Mathematica program for tetration based on the series with q-binomial coefficients
#1
A Mathematica program to compute tetration \( {^z a} \) for \( 1<a<e^{1/e} \) and \( \Re(z)>-2 \) based on the formula with q-binomial coefficients from http://mathoverflow.net/q/259278/9550.

Increase Terms and/or WorkingPrecision option values to get more digits. Examples of values computed using this program are http://goo.gl/MLfEMC and http://goo.gl/LVRkk7.

Code:
ClearAll[Tetration];
Options[Tetration] = {Terms -> 40, WorkingPrecision -> 100};
Tetration[_, -1, OptionsPattern[]] = 0;
Tetration[a_, n_Integer /; n >= 0, OptionsPattern[]] := Power @@ Table[a, {n}];
Tetration[a_ /; 1 < a < E^(1/E), z_ /; Re[z] > -2, OptionsPattern[]] :=
  With[{q = -ProductLog[-Log[a]]},
  With[{limit = SequenceLimit[N[Accumulate[Table[Sum[(-1)^(n + k) q^Binomial[n - k, 2] QBinomial[z, n, q] QBinomial[n, k, q] Tetration[a, k], {k, 0, n}], {n, 0, OptionValue[Terms]}]], OptionValue[WorkingPrecision]]]},
   SetPrecision[limit, 3/4 Precision[limit]]]];
   
(* Examples:
   Tetration[Log[3], 1/2]
   Tetration[Log[3], 1/2, Terms -> 50, WorkingPrecision -> 150] *)


Possibly Related Threads…
Thread Author Replies Views Last Post
  Slick Mathematica "trick" Daniel 0 573 09/29/2022, 10:07 AM
Last Post: Daniel
  GitHub repo for Mathematica packlet FractionalIteration Daniel 2 1,205 08/26/2022, 03:34 PM
Last Post: Daniel
Wink new fatou.gp program sheldonison 33 65,611 07/10/2022, 02:23 AM
Last Post: Catullus
  The beta method program JmsNxn 0 1,227 02/25/2022, 03:05 AM
Last Post: JmsNxn
  My new ABEL_L.gp program JmsNxn 13 10,310 10/06/2021, 07:18 PM
Last Post: Ember Edison
  Natural complex tetration program + video MorgothV8 1 6,587 04/27/2018, 07:54 PM
Last Post: MorgothV8
  complex base tetration program sheldonison 23 85,137 10/26/2016, 10:02 AM
Last Post: Gottfried
  C++ program for generatin complex map in EPS format MorgothV8 0 5,576 09/17/2014, 04:14 PM
Last Post: MorgothV8
  Single-exp series computation code mike3 0 5,322 04/20/2010, 08:59 PM
Last Post: mike3
  Computations with the double-exp series mike3 0 4,636 04/20/2010, 07:32 PM
Last Post: mike3



Users browsing this thread: 1 Guest(s)