%I #21 Sep 08 2022 08:44:37
%S 1,0,-1,6,-39,295,-2589,26026,-295754,3755016,-52729546,811837268,
%T -13601028438,246312375686,-4793646410058,99730695561532,
%U -2207506199891660,51757576802639920,-1280121375188444260
%N Expansion of e.g.f.: cos(log(1+log(1+x))).
%H G. C. Greubel, <a href="/A009018/b009018.txt">Table of n, a(n) for n = 0..250</a>
%p seq(coeff(series(factorial(n)*cos(log(1+log(1+x))), x,n+1),x,n),n=0..20); # _Muniru A Asiru_, Jul 21 2018
%t With[{nmax = 30}, CoefficientList[Series[Cos[Log[1 + Log[1 + x]]], {x, 0, nmax}], x]*Range[0, nmax]!] (* _G. C. Greubel_, Jul 21 2018 *)
%o (PARI) x='x+O('x^30); Vec(serlaplace(cos(log(1+log(1+x))))) \\ _G. C. Greubel_, Jul 21 2018
%o (Magma) m:=30; R<x>:=PowerSeriesRing(Rationals(), m); b:=Coefficients(R!(Cos(Log(1 + Log(1 + x))))); [Factorial(n-1)*b[n]: n in [1..m]]; // _G. C. Greubel_, Jul 21 2018
%K sign,easy
%O 0,4
%A _R. H. Hardin_
%E Extended with signs by _Olivier GĂ©rard_, Mar 15 1997
|