Naturally solving the sexp bo198214 Administrator Posts: 1,616 Threads: 102 Joined: Aug 2007 08/21/2007, 06:41 PM (This post was last modified: 08/21/2007, 06:45 PM by bo198214.) The Abel equation for the slog was $\text{slog}_b(b^x)=\text{slog}_b(x)+1$. Our original recurring for the super exponential (which I will call similarly) sexp is $\text{sexp}_b(x+1)=b^{\text{sexp}_b(x)}$. So if we develop it at 0, say $\text{sexp}_b(x)=\sum_{n=0}^\infty \rho_n x^n$ satisfying the inverted Abel equation: $\sum_{n=0}^\infty \rho_n (x+1)^n = \exp\left(\log(b)\sum_{n=0}^\infty \rho_n x^n\right)$ we get $\sum_{n=0}^\infty \rho_n \sum_{k=0}^n \left(n\\k\right) x^k = 1 + \sum_{n=0}^\infty B(\log(b)\rho_1,\dots,\log(b)\rho_n) \frac{x^n}{n!}$ where $B(x_1,\dots,x_n)$ is the complete Bell polynomial. The left side develops to $\sum_{k=0}^\infty x^k \sum_{n=k}^\infty \rho_n \left(n\\k\right)$ and so we have the infinite equation system $\sum_{i=k}^\infty \rho_i \left(i\\k\right) = B(\log(b)\rho_1,\dots,\log(b)\rho_k) \frac{1}{k!}$ And I wonder if we solve it the natural way whether we get exactly the inverse of the slog (which I assume). Unfortunately there is no complete Bell polynomial in Maple (at least I didnt find it) and I am too lazy in the moment to program it myself And yes, it is not a linear equation system. Perhaps it is despite solvable, who knows ... jaydfox Long Time Fellow Posts: 440 Threads: 31 Joined: Aug 2007 08/22/2007, 03:15 AM (This post was last modified: 08/22/2007, 03:29 AM by jaydfox.) I've tried doing something similar. I used Faa di Bruno's formula, which, if I understand correctly, is similar to using a Bell polynomial (I can even see the reference in the article you linked to). However, the system is of non-linear equations (as you point out, and Andrew points this out in his paper). The sytem can be solved iteratively by plugging approximations back into the system, and it converges exponentially for the first 11 terms. When I tried to solve the 12th term (with 11 unknowns, the first term being known), the system was unstable. I was able to manually pick the 11th unknown and solve the 10 unknown system, but this required two sets of iterations. Solving beyond this would appear to require manually tweaking all the terms beyond the 10th and then iterating, and this becomes infeasible beyond about three or four additional terms, due to the non-linearity. If someone has a spiffier method that come overcome the non-linearity, it might be doable, but I don't see it being possible to solve for, say, a 100x100 system. To be honest, I was blown away that the slog function was a simple linear system, given how chaotic the sexp function is. Hats off to Andrew for his approach. ~ Jay Daniel Fox Gottfried Ultimate Fellow Posts: 898 Threads: 130 Joined: Aug 2007 08/22/2007, 08:58 AM bo198214 Wrote:The Abel equation for the slog was $\text{slog}_b(b^x)=\text{slog}_b(x)+1$. Our original recurring for the super exponential (which I will call similarly) sexp is $\text{sexp}_b(x+1)=b^{\text{sexp}_b(x)}$. So if we develop it at 0, say $\text{sexp}_b(x)=\sum_{n=0}^\infty \rho_n x^n$ satisfying the inverted Abel equation: $\sum_{n=0}^\infty \rho_n (x+1)^n = \exp\left(\log(b)\sum_{n=0}^\infty \rho_n x^n\right)$ we get $\sum_{n=0}^\infty \rho_n \sum_{k=0}^n \left(n\\k\right) x^k = 1 + \sum_{n=0}^\infty B(\log(b)\rho_1,\dots,\log(b)\rho_n) \frac{x^n}{n!}$ where $B(x_1,\dots,x_n)$ is the complete Bell polynomial. The left side develops to $\sum_{k=0}^\infty x^k \sum_{n=k}^\infty \rho_n \left(n\\k\right)$ and so we have the infinite equation system $\sum_{i=k}^\infty \rho_i \left(i\\k\right) = B(\log(b)\rho_1,\dots,\log(b)\rho_k) \frac{1}{k!}$ And I wonder if we solve it the natural way whether we get exactly the inverse of the slog (which I assume). Unfortunately there is no complete Bell polynomial in Maple (at least I didnt find it) and I am too lazy in the moment to program it myself And yes, it is not a linear equation system. Perhaps it is despite solvable, who knows ... :-) Now I have a name for my matrix-computations. Expliciting the matrix-formulae into sum-notation for the multiplications I arrived exactly at these "complete Bell-polynomials". Thanks! I'll have to read and adapt my naming conventions. Gottfried Gottfried Helms, Kassel « Next Oldest | Next Newest »

 Possibly Related Threads… Thread Author Replies Views Last Post When Does $$\displaystyle\int_{-1}^0\text{sexp}(x)\,\mathrm{d}x$$ Equal $$\frac12$$? Catullus 0 278 10/31/2022, 11:47 PM Last Post: Catullus solving f(g(x)) = f(x) converging to f(exp(x)) = f(x) tommy1729 2 821 05/26/2022, 11:07 PM Last Post: JmsNxn Revitalizing an old idea : estimated fake sexp'(x) = F3(x) tommy1729 0 663 02/27/2022, 10:17 PM Last Post: tommy1729 Sexp redefined ? Exp^[a]( - 00 ). + question ( TPID 19 ??) tommy1729 0 3,919 09/06/2016, 04:23 PM Last Post: tommy1729 Solving tetration using differintegrals and super-roots JmsNxn 0 4,145 08/22/2016, 10:07 PM Last Post: JmsNxn Can sexp(z) be periodic ?? tommy1729 2 8,359 01/14/2015, 01:19 PM Last Post: tommy1729 pseudo2periodic sexp. tommy1729 0 3,678 06/27/2014, 10:45 PM Last Post: tommy1729 [2014] tommy's theorem sexp ' (z) =/= 0 ? tommy1729 1 5,999 06/17/2014, 01:25 PM Last Post: sheldonison Multiple exp^[1/2](z) by same sexp ? tommy1729 12 29,175 05/06/2014, 10:55 PM Last Post: tommy1729 entire function close to sexp ?? tommy1729 8 19,655 04/30/2014, 03:49 PM Last Post: JmsNxn

Users browsing this thread: 1 Guest(s)