Quote:Wait... the series above only works for even approximants, I also got:
for odd approximants, and they can't both be right... what am I doing wrong?
Tetra-series
|
Quote:Wait... the series above only works for even approximants, I also got: for odd approximants, and they can't both be right... what am I doing wrong?
10/31/2009, 12:58 PM
(10/31/2009, 10:38 AM)andydude Wrote: This is probably an old result. For the function Hmmm, don't recognize this coefficients. Would you mind to show more of your computations?
Gottfried Helms, Kassel
(10/31/2009, 11:01 AM)andydude Wrote:Quote:Wait... the series above only works for even approximants, I also got: Alternating summing of coefficients requires averaging if the sequence to be summed is not converging fast enough. For example, for the coefficients at the linear term the alternating sum of iterates give 1x-1x+1x-1x... for the coefficient in AS and to get 0.5 x here in the limit one must employ cesaro or Euler-summation. Note, that the same problem strengthens if the sequence of coefficients at some x have also a growthrate (is divergent) Then for each coefficient you need an appropriate order for Cesaro-/Eulersum. If you use powers of a *triangular* Carleman-matrix X for the iterates, then you can try the geometric series for matrices AS = I - X + X^2 - X^3 + ... - ... = (I + X)^-1 in many cases and use the coefficients of the second row/column. If X is not triangular you have a weak chance, that you get a usable approximate, since the eigenvalues of AS may behave nicer than that of X, because for an eigenvalue x>1 the eigenvalue in AS is 1/(1+x) which is smaller than 1, and from a set of increasing eigenvalues (all >=1) as well from one of decreasing eigenvalues (0<all<=1) you get a set of eigenvalues in AS (0<all<1), which makes the associated powerseries for AS(x) nice-behaved.
Gottfried Helms, Kassel
If I Euler-sum the list of coefficients, which I get with the formula at the superroot-msg for heights h=0..63 , then the resulting powerseries seems to begin with
Code: 0 + 1/2x -1/2 x^2 + 1/4 x^3 -1/6 x^4 + 5/12 x^5 + ??? + ??? , ... ] Because this is tidy to avoid the Euler-summation at all we can do a trick: Let S(x,h) be the formal powerseries for the height h for S(x,h) = (1+x)^(1+x)^...^(1+x) - 1 and S(x) the series for the limit when h-> infinity, then by definition, AS(x) = S(x,0) - S(x,1) + S(x,2) - ... + // Euler-sum Since the coefficients of any height converge to that of the S(x)-series I compute the difference D(x,h) = S(x,h) - S(x) and rewrite AS(x) = D(x,0) - D(x,1) + D(x,2) ... + aeta(0)*S(x) where aeta(0) is the alternating zeta-series zeta(0) meaning aeta(0) = 1-1+1-1+1-... = 1/2 Because the coefficients with index k<h in D(x,h) vanish, I get exact rational values in the coefficients of the formal powerseries in AS(x) Code: 0 * x^0 In float-numerical display this is Code: 0 * x^0 Check, for base x=1/2 I get AS(x) = 0.938253002500
Gottfried Helms, Kassel
(10/31/2009, 02:40 PM)Gottfried Wrote: Check, for base x=1/2 I get AS(x) = 0.938253002500 Yes, with the average of the two series, I get AS(0.5) = 0.938253. The way that I got the coefficients is slightly different than your method. I did this: Let then so I thought, if we know G (1, -1, 1, -1, ...), then and when the matrix size is even I get the first series, and when the matrix size is odd, I get the second series. (10/31/2009, 09:37 PM)andydude Wrote: Hmm, I didn't catch the actual computation, but that's possibly not yet required. I guess, you need the technique of divergent-summation ; it may be, that the implicite series, which occur by multiplication of I found a very nice method to get at least an overview, whether such matrix-product suffers from non-convergence (but which is reparable by Cesaro/Euler-summation). I've defined a diagonal-matrix dE(o) of coefficients for Euler-summation of order o, which can simply be included in the matrix-product. Write In my implementation order o=1 means no change; simply dE(1) = I ; dE(1.7)..dE(2.0) can sum 1-1+1-1... and similar, dE(2.5)..dE(3.0) can sum divergence of type 1-2+4-8+...-... and so on (details for toy-implementation in Pari/GP see below, I can send you some example scripts if this is more convenient) Gottfried Code: // Pari/GP
Gottfried Helms, Kassel
(10/31/2009, 09:37 PM)andydude Wrote: The way that I got the coefficients is slightly different than your method. I did this:Ah, now I understand, B is the matrix which transforms the f- into g - coefficients, G is given and F is sought... B is not triangular here: how do you get the correct entries for its inverse, btw? But whatever: I use this idea too, frequently. However in many instances I found in our context of exponentiation and especially iterated exponentiation, that the inverse of some matrix X represents highly divergent series, such that systematically results which are correct using the non-inverted matrix are not correct for the inverse problem using the (naive) inverse of X. This is especially the case for some matrix X, whose triangular LR-factors have the form of a q-binomial matrix. Such LR-factors occur by a square matrix X = x_{r,c} = base^(r*c) or X = x_{r,c} = base^(r*c)/r! or the like, and if X shall be inverted by inversion of its triangular factors. Such matrices X occur for example in the interpolation which I called "exponential polynomial interpolation" for the T-tetration (or sexp)-Bell-matrices. I used that matrix X also in the example for the "false interpolation for logarithm"-discussion. (But I could not yet find a workaround for the occuring inconsistencies with the inverse) Now I don't see the precise characteristics of your B-matrix so far; I've just to actually construct one and to look into it to be able to say more. Let's see... Gottfried
Gottfried Helms, Kassel
11/03/2009, 03:56 AM
11/03/2009, 04:12 AM
O wait... I was wrong
which means sorry.
If F is all zeroes except for a one somewhere, then that represents an
The first coefficient seems to have a pattern in it, but this is just because Oh, and another weird thing: |
« Next Oldest | Next Newest »
|