OFFSET
1,2
COMMENTS
See A228929 for explanation.
LINKS
G. C. Greubel, Table of n, a(n) for n = 1..500
FORMULA
Given a positive real number x, let z(0) = x - floor(x) and z(k+1) = abs(z(k))*round(1/abs(z(k))) - 1; then a(n) = sign(z(n))*round(1/abs(z(n))) for n>0.
EXAMPLE
e = 2 + 1*(1 - 1/4*(1 + 1/8*(1 + 1/67*(1 + 1/266*(1 + 1/9757*(1 + ...)))))).
MAPLE
ArticoExp := proc (n, q::posint)::list; local L, i, z; Digits := 50000; L := []; z := frac(evalf(n)); for i to q+1 do if z = 0 then break end if; L := [op(L), round(1/abs(z))*sign(z)]; z := abs(z)*round(1/abs(z))-1 end do; return L end proc
# List the first 20 terms of the expansion of exp(1)-2
ArticoExp(exp(1), 20)
MATHEMATICA
ArticoExp[x_, n_] := Round[1/#] & /@ NestList[Round[1/Abs[#]]*Abs[#] - 1 &, FractionalPart[x], n]; Block[{$MaxExtraPrecision = 50000}, ArticoExp[Exp[1] - 2, 20]] (* G. C. Greubel, Dec 26 2016 *)
CROSSREFS
KEYWORD
sign,cofr
AUTHOR
Giovanni Artico, Sep 09 2013
STATUS
approved