cp's OEIS Frontend

A081566 Second binomial transform of expansion of exp(3cosh(x)).

%I A081566 #14 Sep 08 2022 08:45:09
%S A081566 1,2,7,26,118,572,3127,18146,114793,765602,5463982,40870436,323326813,
%T A081566 2667777842,23092966267,207651618746,1947316349278,18906249136892,
%U A081566 190564801592107,1982986181092226,21345005629846213,236628248493001202
%N A081566 Second binomial transform of expansion of exp(3cosh(x)).
%C A081566 Binomial transform of A081565.
%H A081566 Vincenzo Librandi, <a href="/A081566/b081566.txt">Table of n, a(n) for n = 0..200</a>
%F A081566 E.g.f.: exp(2*x) * exp(3*cosh(x))/e^3 = exp(3*cosh(x)+2*x-3).
%p A081566 seq(coeff(series(exp(3*cosh(x)+2*x-3), x, n+1)*factorial(n), x, n), n = 0 .. 30); # _G. C. Greubel_, Aug 13 2019
%t A081566 With[{nn = 30}, CoefficientList[Series[Exp[3 Cosh[x] + 2 x - 3], {x, 0, nn}], x] Range[0, nn]!] (* _Vincenzo Librandi_, Aug 08 2013 *)
%o A081566 (PARI) my(x='x+O('x^30)); Vec(serlaplace( exp(3*cosh(x)+2*x-3) )) \\ _G. C. Greubel_, Aug 13 2019
%o A081566 (Magma) m:=30; R<x>:=PowerSeriesRing(Rationals(), m); b:=Coefficients(R!( Exp(3*Cosh(x)+2*x-3) )); [Factorial(n-1)*b[n]: n in [1..m]]; // _G. C. Greubel_, Aug 13 2019
%o A081566 (Sage) [factorial(n)*( exp(3*cosh(x)+2*x-3) ).series(x,n+1).list()[n] for n in (0..30)] # _G. C. Greubel_, Aug 13 2019
%Y A081566 Cf. A081564, A081565.
%K A081566 easy,nonn
%O A081566 0,2
%A A081566 _Paul Barry_, Mar 22 2003