Terse Schroeder & Abel function code
#1
The following provides the Schroeder and Abel functions in Mathematica using the support for solving recursive equations and Bell polynomials. The code takes about a hour to run for the first ten derivatives.
Code:
order=10;
H[0]=0;
H[1]=f'[0]^t ;
Do[H[max]=First[r[t]/.RSolve[{r[0]==0,r[t]==Sum[Derivative[k][f][0]BellY[max,k,Table[H[j]/.t->t-1,{j,max}]],{k,2,max}]+ f'[0] r[t-1]},r[t],t]],{max,2,order}];
Schroeder=f'[0]^t z+Sum[1/k! H[k]z^k,{k,2,order}]
Abel=Limit[Schroeder,{f'[0]->1}]
Daniel
Reply
#2
The Pari/GP code for Bell polynomials is
Code:
  Bell(k,n=-1)=
  { my(x, v, dv, var = i->eval(Str("X",i)));
 
    v = vector(k, i, if (i==1, 'E, var(i-1)));
    dv = vector(k, i, if (i==1, 'X*var(1)*'E, var(i)));
    x = diffop('E,v,dv,k) / 'E;
    if (n < 0, subst(x,'X,1), polcoef(x,n,'X));
  }

So to Mathematica code at the beginning of this thread can be migrated to Pari/GP.
Daniel
Reply


Possibly Related Threads…
Thread Author Replies Views Last Post
  Code to calculate tetration using my method Shanghai46 10 3,922 10/19/2023, 09:15 PM
Last Post: marracco
  Quickest way to compute the Abel function on the Shell-Thron boundary JmsNxn 1 1,572 10/10/2023, 05:17 AM
Last Post: leon
  FractionalIteration now in Wolfram Function Repository Daniel 0 795 07/10/2023, 07:16 AM
Last Post: Daniel
  Revisting my accelerated slog solution using Abel matrix inversion jaydfox 22 49,682 05/16/2021, 11:51 AM
Last Post: Gottfried
  C++ code for tet, ate and hexp MorgothV8 0 6,266 07/10/2014, 04:24 PM
Last Post: MorgothV8
  "Kneser"/Riemann mapping method code for *complex* bases mike3 2 13,477 08/15/2011, 03:14 PM
Last Post: Gottfried
  An incremental method to compute (Abel) matrix inverses bo198214 3 16,804 07/20/2010, 12:13 PM
Last Post: Gottfried
  Single-exp series computation code mike3 0 5,827 04/20/2010, 08:59 PM
Last Post: mike3
  SAGE code for computing flow matrix for exp(z)-1 jaydfox 4 17,762 08/21/2009, 05:32 PM
Last Post: jaydfox
  Computing Abel function at a given center jaydfox 10 27,953 11/30/2007, 06:44 PM
Last Post: andydude



Users browsing this thread: 1 Guest(s)