• 0 Vote(s) - 0 Average
• 1
• 2
• 3
• 4
• 5
 closed form for regular superfunction expressed as a periodic function mike3 Long Time Fellow Posts: 368 Threads: 44 Joined: Sep 2009 08/31/2010, 08:04 AM (This post was last modified: 08/31/2010, 08:32 AM by mike3.) (08/31/2010, 06:51 AM)mike3 Wrote: The product of exp is the exponential of a formal power series. This can be expressed using the Bell polynomials: $\prod_{n=0}^{\infty} \exp(a_n y^n) = \exp\left(\sum_{n=0}^{\infty} a_n y^n\right) = \exp(a_0) \exp\left(\sum_{n=1}^{\infty} a_n y^n\right)$. This then becomes $\exp(a_0) \exp\left(\sum_{n=1}^{\infty} a_n y^n\right) = \exp(a_0) \sum_{n=1}^{\infty} \frac{\sum_{k=1}^n B_{n,k}(1! a_1, ..., (n-k+1)! a_{n-k+1})}{n!} y^n = \exp(a_0) \sum_{n=1}^{\infty} \frac{B_n(1! a_1, ..., n! a_n)}{n!} y^n$. Thus the equations to solve are $a_n L^n = \exp(a_0) \frac{B_n(1! a_1, ..., n! a_n)}{n!}$. Since $a_0 = L$ and $a_1 = 1$, this is $a_n L^n = L \frac{B_n(1, 2! a_2, ..., n! a_n)}{n!}$. This is derived from Faà di Bruno's formula, see http://en.wikipedia.org/wiki/Fa%C3%A0_di...7s_formula for details. Doing some tests, it appears that $B_n(1, 2! a_2, ..., n! a_n)$ has only one occurrence of $n! a_n$, and no higher powers of it, and it never seems to be multiplied by any sort of n-dependent coefficient. This means that $B_n(1, a_2, ..., a_n) - n! a_n = B_n(1, 2! a_2, ..., (n-1)! a_{n-1}, 0)$. I don't have a proof at this point, but we can then solve this: $a_n n! L^n = L B_n(1, 2! a_2, ..., n! a_n)$ $a_n n! L^{n-1} = B_n(1, 2! a_2, ..., n! a_n)$ $a_n n! L^{n-1} - n! a_n = B_n(1, 2! a_2, 3! a_3, ..., (n-1)! a_{n-1}, 0)$ $n! (L^{n-1} - 1) a_n = B_n(1, 2! a_2, 3! a_3, ..., (n-1)! a_{n-1}, 0)$ $a_n = \frac{B_n(1, 2! a_2, 3! a_3, ..., (n-1)! a_{n-1}, 0)}{n! (L^{n-1} - 1)}$. And this is the recurrent formula for the general coefficients. Together with $a_0 = L$ and $a_1 = 1$, this completes the as-close-to-explicit-as-possible-so-far formula for the coefficients of the Fourier series for the regular iteration. EDIT: posts corrected to include factorials on terms $a_n$ in Bell polynomials sheldonison Long Time Fellow Posts: 626 Threads: 22 Joined: Oct 2008 08/31/2010, 02:37 PM (08/31/2010, 07:08 AM)Gottfried Wrote: Hi Sheldon - I recognize your coefficients....Gottfried, thanks for your reply. I'm trying to make sense of the coefficients, and would probably need to program the matrix into pari-gp to verify it. Earlier, Gottfried wrote: Quote:Hi Sheldon - just to allow me to follow (think I can't involve much) - I don't have a clue from where this is coming, what, for instance, is L at all? I think you've explained it elsewhere before but don't see it at the moment... Would you mind to reexplain in short or to provide the link? Gottfried L is the fixed point of base(e). If you iterate the natural logarithm function hundreds of times (say starting with z=0.5), you get a very good approximation of L. Then the equations describe the regular superfunction for base e, which is complex valued at the real axis, but is analytic and entire. The rest of what I did basically, it comes down to expressing the regular superfunction as an analytic Fourier series, where all of the terms of the series decay to zero at i*infinity. Such a series is guaranteed to be analytic. Here is an example of a generalized view of such a Fourier series, not related to the superfunction, with a period of 2Pi. $f(z) = \sum_{n=0}^{\infty}a_ne^{i*nz}$ A full Fourier series would also have the negative coefficients as well, and is often not analytic. $f(z) = \sum_{n=-\infty}^{\infty}a_ne^{i*nz}$ Both of these can be wrapped around the unit circle, with the substitution $y=e^{iz}$. In the analytic case, we have an analytic Taylor series. $f(z) = \sum_{n=0}^{\infty}a_ny^n$. In the general case, with terms from -infinty to +infinty, we have a Laurent series, with singularities inside the unit circle, and an annular ring of convergence. For the full Fourier series, often, the Laurent series only converges on the edge of the unit circle. Hope that helps. - Sheldon sheldonison Long Time Fellow Posts: 626 Threads: 22 Joined: Oct 2008 08/31/2010, 02:41 PM Mike, Thanks a lot for your links and your post! I'll have to program it into pari-gp, and verify that it matches the discreet Fourier series terms. Next goal; rewrite the equations for other bases (except $\eta$), especially for basespc dFac(1,5)*Mat(VE(pc,5) )    \\ show the coefficients scaled by factorials in window "pc" \\ result: \\                                         1 \\                                        a1 \\                                 a1^2+2*a2 \\                         a1^3+6*a2*a1+6*a3 \\  a1^4+12*a2*a1^2+24*a3*a1+(12*a2^2+24*a4) \\ now begin to substitute a1 = 1 pc1=subst(VE(pc,6),'a1,a1) %box >pc dFac(1,6)*Mat(pc1 )    \\ show the coefficients scaled by factorials \\                                                   1 \\                                                   1 \\                                              2*a2+1 \\                                       6*a2+(6*a3+1) \\                     12*a2^2+12*a2+(24*a3+(24*a4+1)) a2 = 1/2/(L-1) pc2=subst(VE(pc1,6),'a2,a2) %box >pc dFac(1,4)*Mat(VE(pc2,4) ) \\ show only 4 rows                                                     1                                                     1                                               L/(L-1)                          ((6*a3+1)*L+(-6*a3+2))/(L-1) a3 = (1/6 + a2)/(L^2-1) pc3=subst(VE(pc2,6),'a3,a3) %box >pc dFac(1,4)*Mat(VE(pc3,4) ) a4=(1/24 + a3 + a2/2 + a2^2/2)/(L^3-1) pc4=subst(VE(pc3,6),'a4, a4) %box >pc dFac(1,5)*Mat(VE(pc4,5) ) a5=(1/5! + a4+a3/2+a3*a2+a2^2/2+a2/3!)/(L^4-1) pc5=subst(VE(pc4,6),'a5, a5) %box >pc dFac(1,6)*Mat(VE(pc5,6) )                                                                              1 \\                                                                             1 \\                                                                       L/(L-1) \\                                                     (L^3+2*L^2)/(L^3-L^2-L+1) \\                                 (L^6+5*L^5+6*L^4+6*L^3)/(L^6-L^5-L^4+L^2+L-1) \\  (L^10+9*L^9+24*L^8+40*L^7+46*L^6+36*L^5+24*L^4)/(L^10-L^9-L^8+2*L^5-L^2-L+1) \\ see the coefficients at powers of L in the numerators equal that matrix-entries, for instance last row(revert the order and begin at L^4: [24  36  46  40  24  9  1] \\ denominators can be factored by (L-1),(L^2-1) and so on` Thanks also for the other explanations. It surely helps to have an entry for my own understanding of the procedure once I've got the whole idea behind... So far - Gottfried Gottfried Helms, Kassel mike3 Long Time Fellow Posts: 368 Threads: 44 Joined: Sep 2009 08/31/2010, 07:26 PM (This post was last modified: 08/31/2010, 07:51 PM by mike3.) Formulas for the Bell polynomials can be found here: http://en.wikipedia.org/wiki/Bell_polynomials and here: http://mathworld.wolfram.com/BellPolynomial.html There's a method involving taking a matrix determinant, a convolution formula, and a sum over partitions (we are only interested in the "complete" Bell polynomials $B_n$ here). E.g. $(x \diam y)_n = \sum_{j=1}^{n-1} {n \choose j} x_j y_{n-j}$ then $B_{n,k}(x_1, ..., x_{n-k+1}) = \frac{x_n^{k \diam}}{k!}$ (where $x^{k \diam}$ is like "exponentiating" the sequence with the "diamond" operator) and $B_n(x_1, ..., x_n) = \sum_{k=1}^n B_{n,k}(x_1, ..., x_{n-k+1})$ mike3 Long Time Fellow Posts: 368 Threads: 44 Joined: Sep 2009 08/31/2010, 07:34 PM (This post was last modified: 08/31/2010, 07:39 PM by mike3.) For the variable-base case, $a_n (L \log(B))^n = L \frac{B_n(\log(B), \log(B) 2! a_2, ..., \log(B) n! a_n)}{n!}$ thus $a_n L^{n-1} \log(B)^n n! = B_n(\log(B), \log(B) 2! a_2, ..., \log(B) n! a_n)$ $a_n = \frac{B_n(\log(B), \log(B) 2! a_2, ..., \log(B) (n-1)! a_{n-1}, 0)}{(L^{n-1} \log(B)^n - \log(B)) n!}$. mike3 Long Time Fellow Posts: 368 Threads: 44 Joined: Sep 2009 08/31/2010, 08:23 PM Now all we need is some way to express the Riemann mapping sheldonison Long Time Fellow Posts: 626 Threads: 22 Joined: Oct 2008 08/31/2010, 08:33 PM (08/31/2010, 07:34 PM)mike3 Wrote: For the variable-base case, $a_n (L \log(B))^n = L \frac{B_n(\log(B), \log(B) 2! a_2, ..., \log(B) n! a_n)}{n!}$ thus $a_n L^{n-1} \log(B)^n n! = B_n(\log(B), \log(B) 2! a_2, ..., \log(B) n! a_n)$ $a_n = \frac{B_n(\log(B), \log(B) 2! a_2, ..., \log(B) (n-1)! a_{n-1}, 0)}{(L^{n-1} \log(B)^n - \log(B)) n!}$.Wow! Thanks Mike! Only matrices and statistics, "n choose k" formulas is one of my weak areas, so it might be a few days before I catch up. I'm kind of intrigued at actually having a closed form for a real valued superfunction, like sqrt(2). I'm guessing that the number of terms required for convergence gets extremely large as z increases. From the terms I computed a couple of days ago, it looks like the terms are decreasing exponentially, which means the series acts like it has a singularity. However, since the regular superfunction is entire, we must have convergence to infinity, so the terms must eventually decrease faster than exponentially. - Sheldon Gottfried Ultimate Fellow Posts: 757 Threads: 116 Joined: Aug 2007 08/31/2010, 08:34 PM (This post was last modified: 08/31/2010, 08:40 PM by Gottfried.) Please excuse - this post contains some s But it's worth to notice how the bell-polynomials and my -for most fellows here: still cryptic - matrix-method are related. That I didn't see this earlier ----------------- In the mathworld-wolfram-link I find a short description of the Bell-polynomials first kind, fortunately with some example. $\sum_{k=0}^{\infty} \frac{B_k(x)}{k!}t^k = e^{(e^t-1)x}$ This expressed with my matrix-formulae is V(t)~ * dF(-1)*S2 * V(x) = exp( (exp(t)-1)*x) or using my standard-matrix for x->exp(x)-1 plus one intermediate step V(t) ~ * fS2F = V(exp(t)-1)~ write "y" for "exp(t)-1" V(y)~ * dF(-1) * V(x) = F(-1)~* dV(y)*V(x) = F(-1)~ * V( y*x) = exp( y*x) = exp( (exp(t)-1)*x) Or in one expression: V(t) ~ * (fS2F * dF(-1)) * V(x) = exp( (exp(t)-1)*x) (which I can recognize immediately to be correct because I'm extremely used to that notation) or the even simpler definition of the vector of Bell-polynomials: B(x) = S2*V(x) So now I see at least, how the Bell-polynomials of the *first kind* are related to my matrix-lingo. (So I should rename fS2F into "Bell" perhaps...) ---------------------------------- But then follows the Bell-polynomials of *second kind*. And there I'm lost again... No example, no redundancy... as if the reader could not do some error in parsing a complex formula... Is there possibly meant the iteration of the Bell-polynomials (=matrix-power of S2)? That would be simple then... Or, wait - looking at the *subscripted* "x" I think, that they are now coefficients for some arbitrary function developed as a powerseries with reciprocal factorials, ... and then... Well, this seems to be just the definition of what I found out myself at the very beginning of my fiddling with this subject and called it a "matrixoperator" for some function. With the additional property, that it is (lower) triangular (because of the missing x0 in the formula in mathworld), so for instance all my U_t-matrices for the decremented exponentiation (or "U-tetration to base t" in my early speak) were such collections of Bell-polynomials of second kind (but also the schröder-matrices, just all my lower-triangular matrix-operators which I worked with the last years and have the reciprocal factorial scaling ) So now - if we talk about the Bell-polynomials (first or second kind) and the symbolic representation in terms of the log(L) (or log(t)) then this is just what I solved in my earlier posted link (good to know)! Here it is again: http://go.helms-net.de/math/tetdocs/APT.htm . And Faa di Bruno/Bell-polynomials and "matrixoperators" is the same and one needs only read about one side of these notations... Amen - Gottfried Gottfried Helms, Kassel « Next Oldest | Next Newest »

 Possibly Related Threads... Thread Author Replies Views Last Post Half-iterates and periodic stuff , my mod method [2019] tommy1729 0 84 09/09/2019, 10:55 PM Last Post: tommy1729 Is there a function space for tetration? Chenjesu 0 194 06/23/2019, 08:24 PM Last Post: Chenjesu Degamma function Xorter 0 672 10/22/2018, 11:29 AM Last Post: Xorter Periodic analytic iterations by Riemann mapping tommy1729 1 1,973 03/05/2016, 10:07 PM Last Post: tommy1729 Should tetration be a multivalued function? marraco 17 14,929 01/14/2016, 04:24 AM Last Post: marraco Introducing new special function : Lambert_t(z,r) tommy1729 2 3,240 01/10/2016, 06:14 PM Last Post: tommy1729 Natural cyclic superfunction tommy1729 3 2,688 12/08/2015, 12:09 AM Last Post: tommy1729 Tommy-Mandelbrot function tommy1729 0 1,778 04/21/2015, 01:02 PM Last Post: tommy1729 Can sexp(z) be periodic ?? tommy1729 2 3,520 01/14/2015, 01:19 PM Last Post: tommy1729 [2014] Uniqueness of periodic superfunction tommy1729 0 1,760 11/09/2014, 10:20 PM Last Post: tommy1729

Users browsing this thread: 1 Guest(s)