%I #7 Mar 17 2021 08:02:47
%S 1,3,10,36,133,501,1918,7440,29180,115522,461044,1852938,7492846,
%T 30464306,124461782,510696350,2103708187,8696498477,36066269640,
%U 150015248758,625664295594,2615929689642,10962436020878,46037427169060
%N A lower diagonal of pendular trinomial triangle A119369.
%H G. C. Greubel, <a href="/A119374/b119374.txt">Table of n, a(n) for n = 0..1000</a>
%F G.f.: A(x) = B(x)^3/(1+x - x*B(x)) = B(x)^3*G(x) = B(x)^2*H(x) = B(x)*I(x), where B(x) is g.f. of A119370, G(x) is g.f. of A119371, H(x) is g.f. of A119372 and I(x) is g.f. of A119373.
%F G.f.: 16*(1+x)/( ((1+x^2) +sqrt((1+x^2)^2-4*x*(1+x)))^3*(1+4*x+x^2 +sqrt((1+4*x+x^2)^2 - 4*x*(1+x)*(3+2*x))) ).
%t CoefficientList[Series[16*(1+x)/( ((1+x^2) +Sqrt[(1+x^2)^2 -4*x*(1+x)])^3*(1+4*x +x^2 +Sqrt[(1+4*x+x^2)^2 -4*x*(1+x)*(3+2*x)])), {x,0,30}], x] (* _G. C. Greubel_, Mar 16 2021 *)
%o (PARI) {a(n)=polcoeff(16*(1+x)/((1+x^2)+sqrt((1+x^2)^2-4*x*(1+x)+x*O(x^n)))^3 /(1+4*x+x^2 + sqrt((1+4*x+x^2)^2 - 4*x*(1+x)*(3+2*x)+x*O(x^n))),n)}
%o (Sage)
%o def A119374_list(prec):
%o P.<x> = PowerSeriesRing(QQ, prec)
%o return P( 16*(1+x)/( ((1+x^2) +sqrt((1+x^2)^2-4*x*(1+x)))^3*(1+4*x+x^2 +sqrt((1+4*x+x^2)^2 - 4*x*(1+x)*(3+2*x))) ) ).list()
%o A119374_list(30) # _G. C. Greubel_, Mar 16 2021
%o (Magma)
%o R<x>:=PowerSeriesRing(Rationals(), 30);
%o Coefficients(R!( 16*(1+x)/( ((1+x^2) +Sqrt((1+x^2)^2-4*x*(1+x)))^3*(1+4*x+x^2 +Sqrt((1+4*x+x^2)^2 - 4*x*(1+x)*(3+2*x))) ) )); // _G. C. Greubel_, Mar 16 2021
%Y Cf. A119369, A119370, A119371, A119372, A119373, A119375, A119376.
%K nonn
%O 0,2
%A _Paul D. Hanna_, May 17 2006
|