[Date Index]
[Thread Index]
[Author Index]
an [0,1] iterative based on the Farey tree that isn't b-normal
*To*: mathgroup at smc.vnet.net
*Subject*: [mg52089] an [0,1] iterative based on the Farey tree that isn't b-normal
*From*: Roger Bagula <tftn at earthlink.net>
*Date*: Wed, 10 Nov 2004 04:45:38 -0500 (EST)
*Reply-to*: tftn at earthlink.net
*Sender*: owner-wri-mathgroup at wolfram.com
This is an iteration based on two transcental numbers E and Pi.
Any other pair of rational transcendentals can be used.
The result is limited ( not very well here!) to the interval [0,1]
It isn't actually a very good way to manufacture transcendentals, but it
works most of the time.
The idea is to present a [0,1] iterative
that doesn't depend on the Modulo one fractional part
and isn't b-normal by the Bailey definition.
Clear[x,y,a,b,s,g,a0,f]
(* transcendental Farey tree recursive function*)
t1=Pi;
t2=E;
f[a_,b_]:=((a/b)/(1-a/b))/;0<=a/b<=1/2
f[a_,b_]:=((1-a/b)/(a/b))/;1/2<a/b<=1
g[n_]:=g[n]=If[( g[n-1]<1/2)&&(g[n-1]>0),N[g[n-1]/(1-g[n-1]),20],If[
g[n-1]<1,N[(1-g[n-1])/g[n-1],20],N[(E/Pi)^2,20]]]
g[0]=N[f[t2,t1],20];
a0=Table[If [g[n]<1,g[n],0],{n,0,400}]
ListPlot[a0,PlotJoined->True,PlotRange->All]
a1=Sort[a0]
ListPlot[a1,PlotJoined->True,PlotRange->All]
Respectfully, Roger L. Bagula
tftn at earthlink.net, 11759Waterhill Road, Lakeside,Ca 92040-2905,tel: 619-5610814 :
alternative email: rlbtftn at netscape.net
URL : http://home.earthlink.net/~tftn
Prev by Date:
**Re: Converting Dates to Mathematica Format**
Next by Date:
**Re: ListContourPlot - different levels**
Previous by thread:
**Re: nonlinear programming with differential-algebraic constraints (2)**
Next by thread:
**iterative b -normal examples from one of Bailey's papers**
| |