• 0 Vote(s) - 0 Average
• 1
• 2
• 3
• 4
• 5
 Naturally solving the sexp bo198214 Administrator Posts: 1,389 Threads: 90 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: 765 Threads: 119 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 Sexp redefined ? Exp^[a]( - 00 ). + question ( TPID 19 ??) tommy1729 0 1,790 09/06/2016, 04:23 PM Last Post: tommy1729 Solving tetration using differintegrals and super-roots JmsNxn 0 2,121 08/22/2016, 10:07 PM Last Post: JmsNxn Can sexp(z) be periodic ?? tommy1729 2 4,179 01/14/2015, 01:19 PM Last Post: tommy1729 pseudo2periodic sexp. tommy1729 0 1,975 06/27/2014, 10:45 PM Last Post: tommy1729 [2014] tommy's theorem sexp ' (z) =/= 0 ? tommy1729 1 3,086 06/17/2014, 01:25 PM Last Post: sheldonison Multiple exp^[1/2](z) by same sexp ? tommy1729 12 15,055 05/06/2014, 10:55 PM Last Post: tommy1729 entire function close to sexp ?? tommy1729 8 10,319 04/30/2014, 03:49 PM Last Post: JmsNxn Is sexp(z) pseudounivalent for Re(z) > 0 ? tommy1729 3 3,760 03/26/2014, 01:24 PM Last Post: tommy1729 Vincent's theorem and sin(sexp) ? tommy1729 0 1,876 03/22/2014, 11:46 PM Last Post: tommy1729 sexp for base (1/e)^e ~= 0.0660? sheldonison 10 12,129 11/22/2013, 11:20 PM Last Post: mike3

Users browsing this thread: 1 Guest(s)