%I #16 Sep 08 2022 08:45:29
%S 1,0,-4,48,-624,9600,-175680,3790080,-95235840,2752081920,
%T -90328089600,3328103116800,-136191650918400,6131573025177600,
%U -301213549769932800,16030999766605824000,-918678402394841088000,56387623092958789632000,-3690023220507773140992000
%N Expansion of e.g.f.: sqrt(1+4*x)/(1+2*x).
%C A row of an array that is under investigation.
%H G. C. Greubel, <a href="/A126967/b126967.txt">Table of n, a(n) for n = 0..365</a>
%F D-finite with recurrence: a(n) +6*(n-1)*a(n-1) +4*(n-1)*(2*n-3)*a(n-2)=0. - _R. J. Mathar_, Jan 23 2020
%p seq(coeff(series( sqrt(1+4*x)/(1+2*x), x, n+1)*n!, x, n), n = 0..20); # _G. C. Greubel_, Jan 29 2020
%t nmax=20; CoefficientList[Series[Sqrt[1 + 4 x] / (1 + 2 x), {x, 0, nmax}], x] Range[0, nmax]! (* _Vincenzo Librandi_, Jan 24 2020 *)
%o (Magma) m:=20; R<x>:=PowerSeriesRing(Rationals(), m); b:=Coefficients(R!(Sqrt(1+4*x)/(1+2*x))); [Factorial(n-1)*b[n]: n in [1..m]]; // _Vincenzo Librandi_, Jan 24 2020
%o (PARI) my(x='x+O('x^30)); Vec(serlaplace( sqrt(1+4*x)/(1+2*x) )) \\ _G. C. Greubel_, Jan 29 2020
%o (Sage) [factorial(n)*( sqrt(1+4*x)/(1+2*x) ).series(x,n+1).list()[n] for n in (0..30)] # _G. C. Greubel_, Jan 29 2020
%Y Cf. A126966.
%K sign
%O 0,3
%A _N. J. A. Sloane_, Mar 22 2007