OFFSET
1,1
COMMENTS
The function f defined at A229350 is the continued fraction transform; to define f(x), start with x > 0: let p(i)/q(i), for i >=0, be the convergents to x; then f(x) is the number [p(0)/q(0), p(1)/q(1), p(2)/q(2), ... ].
Thus, f(pi) = 3.291191..., f(f(pi)) = 3.276718..., f(f(f(pi))) =3.276503 ...; let L(x) = lim(f(n,x)), where f(0,x) = x, f(1,x) = f(x), and f(n,x) = f(f(n-1,x)). Then L(pi) = 3.276503 ..., as in A228993.
EXAMPLE
f(pi) = 3.29119170409194159105202823875110454700455275620391586...
MATHEMATICA
$MaxExtraPrecision = Infinity;
z = 600; x[0] = Pi; c[0] = Convergents[x[0], z]; x[n_] := N[FromContinuedFraction[c[n - 1]], 80]; c[n_] := Convergents[x[n]]; Table[x[n], {n, 1, 20}] (* A228492, f(pi), f(f(pi)), ... *)
t1 = RealDigits[x[1]] (* A228493 *)
t2 = Numerator[c[1]] (* A228992 *)
t3 = Denominator[c[1]] (* A228993 *)
CROSSREFS
KEYWORD
nonn,cons
AUTHOR
Clark Kimberling, Oct 01 2013
STATUS
approved