%I #10 Aug 18 2015 00:14:39
%S 1,1,2,10,560,18393200,11307340057302083200,
%T 79095479027242971758816977525848827652668769392000
%N a(n) = numerator of b(n): b(1)=1; b(n+1) = b(n) * [b(1);b(2),...,b(n)], where [...] is a continued fraction of rational terms.
%C a(9) and a(10) have 131 and 343 digits, respectively and are too large to include here. - _R. J. Mathar_, Oct 08 2007
%e a(5) = the numerator of b(5). b(5) = (10/3) * (1 +1/(1 +1/(2 +3/10))) = 560/99.
%p 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: A128295 := 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( numer(b[i]),i=1..nops(b))] ; end: A128295() ; # _R. J. Mathar_, Oct 08 2007
%Y Cf. A128296.
%K frac,nonn
%O 1,3
%A _Leroy Quet_, Feb 25 2007
%E More terms from _R. J. Mathar_, Oct 08 2007
|