OFFSET
1,1
EXAMPLE
sin(1) = 0.84147... = A049469 has the Engel expansion 1/2 + 1/(2*2) + 1/(2*2*3) + ...
MATHEMATICA
EngelExp[A_, n_]:=Join[Array[1&, Floor[A]], First@Transpose@NestList[{Ceiling[1/Expand[ #[[1]]#[[2]]-1]], Expand[ #[[1]]#[[2]]-1]}&, {Ceiling[1/(A-Floor[A])], A-Floor[A]}, n-1]]; EngelExp[N[Sin[1], 6! ], 50] (* Vladimir Joseph Stephan Orlovsky, Jun 13 2009 *)
PROG
(PARI) : s=sin(1); for(i=1, 30, s=s*ceil(1/s)-1; print1(ceil(1/s), ", "); );
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Benoit Cloitre, Mar 03 2002
EXTENSIONS
a(1) inserted by Hauke Worpel (hw1(AT)email.com), Jun 01 2003
Edited by N. J. A. Sloane, Nov 01 2008 at the suggestion of R. J. Mathar
STATUS
approved