|
|
| View previous topic :: View next topic |
| Author |
Message |
Gottfried Helms Guest
|
Posted: Mon Jun 30, 2008 10:06 am Post subject: -- Tetration: x - f(x) + f(f(x)) - f(f(f(x))) + ... - ... = |
|
|
A curious result from study of the tetra-series.
I considered the "reverse" of the tetra-series problem.
Instead of asking for the a_lternating s_um of powertowers of
increasing p_ositive heights (asp)
asp(x,dxp) = dxp°0(x) - dxp°1(x) + dxp°2(x) - ... + ...
where dxp(x) = exp(x)-1 and dxp°h(x) is the h'th integer-iterate.
I asked for a function tf(x) where
asp(x,tf) = (e^x-1)/2 = tf°0(x) - tf°1(x) + tf°2(x) - ... +...
so I ask for the function tf(x) whose tetraseries defines (e^x - 1)/2
Using the matrix-operator-approach I got the result
tf(x) = x - x^2 + 2/3*x^3 - 3/4*x^4 + 11/15*x^5 - 59/72*x^6 + 379/420*x^7 - 331/320*x^8
+ 1805/1512*x^9 - 282379/201600*x^10 + 3307019/1995840*x^11 - 6152789/3110400*x^12
+ 616774003/259459200*x^13 - 3212381993/1117670400*x^14 + 54372093481/15567552000*x^15
- 594671543783/139485265920*x^16 + 58070127447587/11115232128000*x^17
- 1209735800444267/188305108992000*x^18 + 26776614379573099/3379030566912000*x^19
- 209181772596680209/21341245685760000*x^20 + 1034961114326994557/85151570286182400*x^21
- 80852235077445729119/5352384417988608000*x^22
+ 2210690796475549862239/117509166994931712000*x^23
- 18624665294361841906483/793412278431252480000*x^24
+ 379264261780067802109819/12926008369442488320000*x^25
- 6584114267874407529534167/179240649389602504704000*x^26
+ 5046681464320089079803469/109576837040799744000000*x^27
- 326480035696597942691643978259/5646080455772478898176000000*x^28
+ 327920863401689931801359966641/4511103058030460180889600000*x^29
- 418419411682443365665393881223739/4573325169175707907522560000000*x^30
+ 15798888070625404329026746075454779/137047310902965380295426048000000*x^31
+ O(x^32)
which can be determined to arbitrary many coefficients by a
recursive process on rational numbers.
The float-representation is
tf(x) = 1.00000000000*x - 1.00000000000*x^2 + 0.666666666667*x^3 - 0.750000000000*x^4
+ 0.733333333333*x^5 - 0.819444444444*x^6 + 0.902380952381*x^7 - 1.03437500000*x^8
+ 1.19378306878*x^9 - 1.40068948413*x^10 + 1.65695596841*x^11 - 1.97813432356*x^12
+ 2.37715218038*x^13 - 2.87417649515*x^14 + 3.49265533084*x^15 - 4.26332874559*x^16
+ 5.22437379435*x^17 - 6.42433870711*x^18 + 7.92434807834*x^19 - 9.80176020073*x^20
+ 12.1543397362*x^21 - 15.1058348510*x^22 + 18.8129220299*x^23 - 23.4741329327*x^24
+ 29.3411740841*x^25 - 36.7333765543*x^26 + 46.0560972611*x^27 - 57.8241911808*x^28
+ 72.6919467774*x^29 - 91.4912883306*x^30 + 115.280540468*x^31 + O(x^32)
so I assume, that this series tf(x) has radius of convergence limited to
about |x|<0.7
Moreover, the iterates of this function seem always to be of a similar
form, so the alternating sum of the found coefficients is divergent for
each coefficient (but may be Euler-summed). So this result must be considered
in more detail next, since I had inconsistency of the matrix-method
with serial summation either for increasing positive or for increasing
negative heights.
However, for x=1/2 or x=1/3 or smaller we can accelerate convergence
by Euler-summation such that I get good (?) approximation to the
six'th digit for x=1/3 using the truncated series with 31 terms only.
The process for the generation of these coefficients is a bit tedious
yet; so I don't have -for instance - the function, whose iterations
must be non-alternating summed to get the exp(x)-1 value
(or (exp(x)-1)/2 or some other scalar multiple) which -as I guess-
could have better range of convergence.
I'll post the result, if I got it.
Gottfried |
|
| Back to top |
|
 |
| |
Ads |
Advertising
Sponsor
|
|
Guest
|
Posted: Mon Jun 30, 2008 10:08 am Post subject: Re: -- Tetration: x - f(x) + f(f(x)) - f(f(f(x))) + ... - .. |
|
|
On 30 juin, 07:06, Gottfried Helms <he...@uni-kassel.de> wrote:
| Quote: |
A curious result from study of the tetra-series.
I considered the "reverse" of the tetra-series problem.
Instead of asking for the a_lternating s_um of powertowers of
increasing p_ositive heights (asp)
asp(x,dxp) = dxp°0(x) - dxp°1(x) + dxp°2(x) - ... + ...
where dxp(x) = exp(x)-1 and dxp°h(x) is the h'th integer-iterate.
I asked for a function tf(x) where
asp(x,tf) = (e^x-1)/2 = tf°0(x) - tf°1(x) + tf°2(x) - ... +....
so I ask for the function tf(x) whose tetraseries defines (e^x - 1)/2
Using the matrix-operator-approach I got the result
tf(x) = x - x^2 + 2/3*x^3 - 3/4*x^4 + 11/15*x^5 - 59/72*x^6 + 379/420*x^7 - 331/320*x^8
+ 1805/1512*x^9 - 282379/201600*x^10 + 3307019/1995840*x^11 - 6152789/3110400*x^12
+ 616774003/259459200*x^13 - 3212381993/1117670400*x^14 + 54372093481/15567552000*x^15
- 594671543783/139485265920*x^16 + 58070127447587/11115232128000*x^17
- 1209735800444267/188305108992000*x^18 + 26776614379573099/3379030566912000*x^19
- 209181772596680209/21341245685760000*x^20 + 1034961114326994557/85151570286182400*x^21
- 80852235077445729119/5352384417988608000*x^22
+ 2210690796475549862239/117509166994931712000*x^23
- 18624665294361841906483/793412278431252480000*x^24
+ 379264261780067802109819/12926008369442488320000*x^25
- 6584114267874407529534167/179240649389602504704000*x^26
+ 5046681464320089079803469/109576837040799744000000*x^27
- 326480035696597942691643978259/5646080455772478898176000000*x^28
+ 327920863401689931801359966641/4511103058030460180889600000*x^29
- 418419411682443365665393881223739/4573325169175707907522560000000*x^30
+ 15798888070625404329026746075454779/137047310902965380295426048000000*x^31
+ O(x^32)
which can be determined to arbitrary many coefficients by a
recursive process on rational numbers.
The float-representation is
tf(x) = 1.00000000000*x - 1.00000000000*x^2 + 0.666666666667*x^3 - 0.750000000000*x^4
+ 0.733333333333*x^5 - 0.819444444444*x^6 + 0.902380952381*x^7 - 1.03437500000*x^8
+ 1.19378306878*x^9 - 1.40068948413*x^10 + 1.65695596841*x^11 - 1.97813432356*x^12
+ 2.37715218038*x^13 - 2.87417649515*x^14 + 3.49265533084*x^15 - 4.26332874559*x^16
+ 5.22437379435*x^17 - 6.42433870711*x^18 + 7.92434807834*x^19 - 9.80176020073*x^20
+ 12.1543397362*x^21 - 15.1058348510*x^22 + 18.8129220299*x^23 - 23.4741329327*x^24
+ 29.3411740841*x^25 - 36.7333765543*x^26 + 46.0560972611*x^27 - 57.8241911808*x^28
+ 72.6919467774*x^29 - 91.4912883306*x^30 + 115.280540468*x^31 + O(x^32)
so I assume, that this series tf(x) has radius of convergence limited to
about |x|<0.7
Moreover, the iterates of this function seem always to be of a similar
form, so the alternating sum of the found coefficients is divergent for
each coefficient (but may be Euler-summed). So this result must be considered
in more detail next, since I had inconsistency of the matrix-method
with serial summation either for increasing positive or for increasing
negative heights.
However, for x=1/2 or x=1/3 or smaller we can accelerate convergence
by Euler-summation such that I get good (?) approximation to the
six'th digit for x=1/3 using the truncated series with 31 terms only.
The process for the generation of these coefficients is a bit tedious
yet; so I don't have -for instance - the function, whose iterations
must be non-alternating summed to get the exp(x)-1 value
(or (exp(x)-1)/2 or some other scalar multiple) which -as I guess-
could have better range of convergence.
I'll post the result, if I got it.
Gottfried
|
Dear Gottfried,
the function f(x) verifies exp(f(x)) + exp(x) = 2*x +2
So f(x)=ln( -exp(x) +2*x+2)
and series near 0, x -x^2 +2/3x^3-11/15x^5 ...
Alain |
|
| Back to top |
|
 |
| |
Ads |
Advertising
Sponsor
|
|
Gottfried Helms Guest
|
Posted: Wed Jul 02, 2008 11:02 am Post subject: Re: -- Tetration: x - f(x) + f(f(x)) - f(f(f(x))) + ... - .. |
|
|
Am 30.06.2008 17:00 schrieb alainverghote@gmail.com:
| Quote: |
Gottfried
Well, inverse function of f(x) using lambertW()
is -Lambert(-1/2exp(1/2e^x - 1)) +1/2e^x - 1
series near 0 , x +x^2 +4/3x^3 +29/12x^4 +51/10x^5 ...
Yepp, I got the same series - good! |
But now - iterations and especially the sum of iterations
of this functions in the Lambert-representation should be
intractable, so I don't assume this can be helpful to get
more insight in the source of the inconsisteny-problem;
remember:
by the formal application of the matrix-approach
asn(x)+asp(x)-x = 0 // expected
or
asn(x,f(x)) = x - asp(x,f(x)) // expected
or
asp(x,f°(-1)(x)) = x - asp(x,f(x)) // expected
is not true, at least for the function f(x) = exp(x)-1
Well - it was a try...
--------------------------------------
Meanwhile I refined my computation-process, so I've
now even the function fz(x) with the condition
e^x - 1 = fz(x) + fz(fz(x)) + fz(fz(fz(x))) + ..
= sum{h=1..inf} fz°h(x)
I got
fz(x) = 2*(x/4)/1! + 6*(x/4)^2/2! + 10*(x/4)^3/3! - 46*(x/4)^4/4! - 554*(x/4)^5/5!
- 1690*(x/4)^6/6! + 27882*(x/4)^7/7! + 505986*(x/4)^8/8! + 2529590*(x/4)^9/9!
- 61918794*(x/4)^10/10! - 1726391798*(x/4)^11/11! - 14268435022*(x/4)^12/12!
+ 352044609814*(x/4)^13/13! + O(x^14)
where the integer parts of the coefficients are
[2]
[6]
[10]
[-46]
[-554]
[-1690]
[27882]
[505986]
[2529590]
[-61918794]
[-1726391798]
[-14268435022]
[352044609814]
which have to be divided py powers of 4 and by factorials to give the
coefficients of the function.
The float representation of this function is
fz(x) = 0.500000000000*x^1 + 0.187500000000*x^2 + 0.0260416666667*x^3 - 0.00748697916667*x^4
- 0.00450846354167*x^5 - 0.000573052300347*x^6 + 0.000337655203683*x^7 + 0.000191486449469*x^8
+ 0.0000265917660278*x^9 - 0.0000162726971838*x^10 - 0.0000103115449999*x^11
- 0.00000177549511523*x^12 + 0.000000842437121499*x^13 + 0.000000632647393830*x^14
+ O(x^15)
Can we give a range for x where this converges ?
The quotients of subsequent coefficients give the following sequence
0.375000000000,0.138888888889,-0.287500000000,
0.602173913043,0.127105896510,-0.589222316145,
0.567106466538,0.138870223462,-0.611944959460,
0.633671534805,0.172185168687,-0.474480112208,
0.750972835462,0.212965249501,-0.333234112499,
0.928624638217,0.256835516535,-0.218432174828,
1.22327362079,0.303250445581,-0.126531683586,
1.82969444504,0.353026162019,-0.0509774937397,
3.94499797418,0.407780103462,0.0133041813614,-13.1183083394,
0.470026310618,0.0698970628404,-2.15481130004,
0.543690120454,...
-----------------------------------------------------------
Using 32 coefficients for the function and 60 iterates for the sum
I could approximate e^1 -1 relatively well. I got
sum(h=1,60,fz°h(1.0)) - ( exp(1)-1 ) = -3.24385306514 E-13
where the quality of approximation increased when terms of the
function and numbers of iterates are increased.
Fun... :-)
Gottfried Helms |
|
| Back to top |
|
 |
| |
Ads |
Advertising
Sponsor
|
|
Gottfried Helms Guest
|
Posted: Thu Jul 03, 2008 10:34 am Post subject: Re: -- Tetration: x - f(x) + f(f(x)) - f(f(f(x))) + ... - .. |
|
|
Am 02.07.2008 20:00 schrieb Gottfried Helms:
| Quote: |
Now, for the alternating sum of iterates of f(x)
asp(x) = x - f(x) + f°2(x) - f°3(x) + ... - ...
and of g(x)
asn(x) = x - g(x) + g°2(x) - g°3(x) + ... - ...
where for asp(x) we have the closed form
asp(x) = (exp(x)-1)/2
it would be good to have a closed form for asn(x)
as well. This seems out of reach.
I just tried with a "poor man's" lambertW-implementation. |
To have convergence/summability I used the parameter x=-0.05
using the closed-form for alternating sum of f-iterates:
asp(-0.05) = (exp(-0.05)-1)/2 = -0.0243852877496
using alternating sum of g-iterates:
asn(-0.05) = sumalt(k=0,(-1)^k*g(-0.05,k)) = -0.0243852877496
The difference:
asn(-0.05) - asp(-0.05) = 2.11329439961 E-99
.... which is again surprising, and the test of the conjecture
of zero-value
asn(-0.05) + asp(-0.05) -(-0.05) = 0.00122942450071 =/= 0
is negative, so indeed the serial summation of terms seems
seem to give the correct result - fine.
That the difference of asn and asp is near zero remains now
as another surprise...
Gottfried Helms |
|
| Back to top |
|
 |
| |
Ads |
Advertising
Sponsor
|
|
Gottfried Helms Guest
|
Posted: Thu Jul 03, 2008 10:48 am Post subject: Re: -- Tetration: x - f(x) + f(f(x)) - f(f(f(x))) + ... - .. |
|
|
Am 03.07.2008 07:34 schrieb Gottfried Helms:
| Quote: |
I just tried with a "poor man's" lambertW-implementation.
To have convergence/summability I used the parameter x=-0.05
using the closed-form for alternating sum of f-iterates:
asp(-0.05) = (exp(-0.05)-1)/2 = -0.0243852877496
|
I made an error in the computation of the sum, giving crap.
I'll repost corrected results soon.
Gottfried Helms |
|
| Back to top |
|
 |
| |
Ads |
Advertising
Sponsor
|
|
Gottfried Helms Guest
|
Posted: Thu Jul 03, 2008 11:01 am Post subject: Re: -- Tetration: x - f(x) + f(f(x)) - f(f(f(x))) + ... - .. |
|
|
Am 02.07.2008 20:00 schrieb Gottfried Helms:
| Quote: |
Now, for the alternating sum of iterates of f(x)
asp(x) = x - f(x) + f°2(x) - f°3(x) + ... - ...
and of g(x)
asn(x) = x - g(x) + g°2(x) - g°3(x) + ... - ...
where for asp(x) we have the closed form
asp(x) = (exp(x)-1)/2
it would be good to have a closed form for asn(x)
as well. This seems out of reach.
I just tried with a "poor man's" lambertW-implementation. |
Recall:
f(x) = log( -exp(x) +2*x+2)
g(x) = -Lambert(-1/2*exp(1/2*exp(x) - 1)) +1/2*exp(x) - 1
where g(x) = f°(-1)(x)
I used Pari/GP with 200 digits float-precision
--------------------------------------------
To have convergence/summability I used the parameter x=-0.5
x = -0.5
using the closed-form for alternating sum of f-iterates:
asp(x) = -0.196734670144... // = (exp(x)-1)/2
using explicitely the alternating sum of f-iterates:
asp(x) = -0.196734670144... + 1.42929732695... E-178*I
// =sumalt(h=0,(-1)^h*f°h(x))
using alternating sum of g-iterates:
asn(x) = -0.303265329856... // =sumalt(h=0,(-1)^h*g°h(x))
The difference:
asn(x) - asp(x) = -0.106530659713...
The old conjecture "as(x) = asp(x) + asn(x) - x = 0"
as(x) = asn(x) + asp(x) - x = -6.09612952428... E-111
So, surprisingly, the old conjecture, that as(x) = 0 matches
this result to very good approximation.
It is surprising, since with f(x) = e^x-1 and g(x) = log(1+x)
I had the systematic deviation from zero as(x) =/= 0 in general.
Take another x:
x = -1
asp(x) = -0.316060279414 // closed form
asp(x) = -0.316060279414 + 2.06700192007 E-187*I // sumalt
asn(x) = -0.683939720586 // sumalt
asn(x)-asp(x) = -0.367879441171
as(x) = asp(x) + asn(x) - x = 5.59598064001 E-111
also here the old conjecture holds.
The explanation, which I could guess at the moment, for this
is, that f(x) and g(x) seem to have the same fixpoint = zero
if iterated to infinity, while exp(x)-1 and log(1+x) converge
to different fixpoints (as far as they both converge for a
certain x).
Gottfried Helms |
|
| Back to top |
|
 |
| |
Ads |
Advertising
Sponsor
|
|
amy666 Guest
|
Posted: Fri Jul 04, 2008 7:56 pm Post subject: lim n-> oo (x - f(x) + f(f(x)) - f(f(f(x))) + ... - ... ) / |
|
|
you might wanna shoot me for not starting a new thread, or not, but i post it here since i feel it is related.
recently i had an old idea again which relates to gottfrieds.
let T(x) be an arbitrary analytic function.
as the title says :
lim n-> oo (x - f(x) + f(f(x)) - f(f(f(x))) + ... - ... ) / n = T(x)
where the sum is taken n times of course , thus an anternating average.
since /n implies that else it doesnt converge , that implies the iterations of f do not converge.
if the iterations are bounded then this can be related to a periodic (analytic) function g where the n th iterate of f = g(n).
if g also has an addition formula , we might even get a closed form for f(x).
and g would then probably be ( here we go again :p )
an inverse hypergeo function !
i believe that with a very good understanding of ergodic theory, umbral calculus and the concepts of the (periodic) inverse hypergeo functions we can solve the equations of type
lim n-> oo (x - f(x) + f(f(x)) - f(f(f(x))) + ... - ... ) / n = T(x)
note :
the limit might perhaps in simple cases be equal to the integral over its period ( contour method ! ) of g.
( this makes sence , if it sounds strange stop for a moment and think about it ! )
it is striking how i once again end up with inverse hypergeo and contour integration !!
however what might interest gottfried more is if this also related to his OP.
we can of course rewrite f(f(x)) as g(2) etc but will that help ?
f(x) = g(y)
f(f(x)) = g(y+1)
x - g(y) + g(y+1) - g(y+2) + ... = T(x)
my guess is the following :
the + and - together with the fact that lim g(oo) is probably = oo leads me to the consideration that the left side is an analytic continuation of a similar series ( of g ).
then also T(x) must be an analytic continuation of something or related to a self-similar looking function ( if it has no poles e.g. ).
( see my concept sum over the reals of f(x) , which is somewhat related )
both approaches are hard but very different.
do the 2 methods relate ?
does the second also relate to contour integration and the concept of inverse hypergeometric functions ??
i admit i dont know the answer.
***
regards
tommy1729 |
|
| Back to top |
|
 |
| |
Ads |
Advertising
Sponsor
|
|
amy666 Guest
|
Posted: Fri Jul 04, 2008 9:47 pm Post subject: Re: lim n-> oo (x - f(x) + f(f(x)) - f(f(f(x))) + ... - ... |
|
|
| Quote: |
you might wanna shoot me for not starting a new
thread, or not, but i post it here since i feel it is
related.
recently i had an old idea again which relates to
gottfrieds.
let T(x) be an arbitrary analytic function.
as the title says :
lim n-> oo (x - f(x) + f(f(x)) - f(f(f(x))) + ... -
... ) / n = T(x)
where the sum is taken n times of course , thus an
anternating average.
since /n implies that else it doesnt converge , that
implies the iterations of f do not converge.
|
unless T(x) = 0 of course.
| Quote: |
if the iterations are bounded then this can be
related to a periodic (analytic) function g where the
n th iterate of f = g(n).
|
however if T(x) is non-constant -> the iterations of f(x) go to oo , the 'period' is reached by iterates of f NOT at x but for a suitable complex value z.
( if f is itself an iteration )
with ' period ' note that sin(n) is periodic but the countable sequence of sin(n) never repeats.
| Quote: |
if g also has an addition formula , we might even get
a closed form for f(x).
|
depends for instance on how smooth and symmetric the periodic function is.
very helpful :
solutions to differential equation :
A' = B(A)
where B°z = B
for some complex z
if there are two solutions for z that dont have a common factor like e.g. i and pi then
we have a double periodic analytic function with periods i and pi.
and probably a double periodic inverse hypergeometric function. ( DPIH function )
| Quote: |
and g would then probably be ( here we go again :p )
an inverse hypergeo function !
i believe that with a very good understanding of
ergodic theory, umbral calculus and the concepts of
the (periodic) inverse hypergeo functions we can
solve the equations of type
lim n-> oo (x - f(x) + f(f(x)) - f(f(f(x))) + ... -
... ) / n = T(x)
note :
the limit might perhaps in simple cases be equal to
the integral over its period ( contour method ! ) of
g.
( this makes sence , if it sounds strange stop for a
moment and think about it ! )
it is striking how i once again end up with inverse
hypergeo and contour integration !!
|
quick note : is there always a solution to
lim n-> oo (x - f(x) + f(f(x)) - f(f(f(x))) + ... - ... ) / n = T(x)
for every analytic T(x) ?
i guess so.
its a quick note ; ill think about it.
| Quote: |
however what might interest gottfried more is if this
also related to his OP.
|
if lim T(y) = 0 at imput y and x.
then perhaps substitute h = 1/n and x = y + h
lim n-> oo (x - f(x) + f(f(x)) - f(f(f(x))) + ... - ... ) = n T(y - 1/n)
now express the ( hopefully ) finite right side limit in terms of a function of x...
and you get the equation in the original form of the OP.
| Quote: |
we can of course rewrite f(f(x)) as g(2) etc but will
that help ?
f(x) = g(y)
f(f(x)) = g(y+1)
x - g(y) + g(y+1) - g(y+2) + ... = T(x)
my guess is the following :
the + and - together with the fact that lim g(oo) is
probably = oo leads me to the consideration that the
left side is an analytic continuation of a similar
series ( of g ).
then also T(x) must be an analytic continuation of
something or related to a self-similar looking
function ( if it has no poles e.g. ).
( see my concept sum over the reals of f(x) , which
is somewhat related )
both approaches are hard but very different.
do the 2 methods relate ?
|
analytic continuation matters are not always present of course.
as for the computation of f(x) that has already been posted by alain.
and is the best way possible.
perhaps i answered my own question ; does a solution always exist ?
yes , but we might need analytic continuation might be the answer.
| Quote: |
does the second also relate to contour integration
and the concept of inverse hypergeometric functions
??
i admit i dont know the answer.
***
regards
tommy1729
|
well , hope i clarified everything now.
the case where T(x) is a constant, yet the iterations of f(x) do not converge has already been well studied , so i give no attention to it.
( see ergodic theory , umbral calculus , logistic map and dynamics )
regards
tommy1729 |
|
| Back to top |
|
 |
| |
Ads |
Advertising
Sponsor
|
|
|
|
You cannot post new topics in this forum You cannot reply to topics in this forum You cannot edit your posts in this forum You cannot delete your posts in this forum You cannot vote in polls in this forum
|
|