L2cfrac := proc(L) local a, i; a := op(-1, L) ; for i from 2 to nops(L) do a := op(-i, L)+1/a ; od: RETURN(a) ; end:
A128294 := proc() local b, n, bnxt; b := [1] ; for n from 2 to 10 do bnxt := op(-1, b)+L2cfrac(b) ; b := [op(b), bnxt] ; od: [seq( denom(b[i]), i=1..nops(b))] ; end:
A128294() ; #
R. J. Mathar, Sep 24 2007