%I #25 Jan 02 2024 08:58:54
%S 0,0,0,1,4,10,22,49,112,256,580,1309,2956,6682,15106,34141,77152,
%T 174352,394024,890473,2012404,4547866,10277806,23227033,52491280,
%U 118626160,268085740,605852581,1369179004,3094236490,6992730202,15803018149
%N a(n) = 4*a(n-1) - 6*a(n-2) + 6*a(n-3) - 3*a(n-4), with initial terms 0, 0, 0, 1.
%C Essentially the partial sums of A052103. - _R. J. Mathar_, Apr 01 2008
%H G. C. Greubel, <a href="/A137247/b137247.txt">Table of n, a(n) for n = 0..1000</a>
%H <a href="/index/Rec#order_04">Index entries for linear recurrences with constant coefficients</a>, signature (4,-6,6,-3).
%F From _R. J. Mathar_, Apr 01 2008: (Start)
%F O.g.f.: x^3/((1-x)*(1-3*x+3*x^2-3*x^3)).
%F A052103(n) = a(n+2) - a(n+1). (End)
%p a[0]:=0: a[1]:=0: a[2]:=0: a[3]:=1: for n from 4 to 30 do a[n]:=4*a[n-1]-6*a[n-2]+6*a[n-3]-3*a[n-4] end do: seq(a[n],n=0..30); # _Emeric Deutsch_, Mar 17 2008
%t LinearRecurrence[{4,-6,6,-3}, {0,0,0,1}, 41] (* _G. C. Greubel_, Apr 15 2021 *)
%o (Magma) I:=[0,0,0,1]; [n le 4 select I[n] else 4*Self(n-1) -6*Self(n-2) +6*Self(n-3) -3*Self(n-4): n in [1..41]]; // _G. C. Greubel_, Apr 15 2021
%o (Sage)
%o def A137247_list(prec):
%o P.<x> = PowerSeriesRing(ZZ, prec)
%o return P( x^3/((1-x)*(1-3*x+3*x^2-3*x^3)) ).list()
%o A137247_list(40) # _G. C. Greubel_, Apr 15 2021
%Y Cf. A052103.
%K nonn,easy
%O 0,5
%A _Paul Curtz_, Mar 10 2008
%E More terms from _R. J. Mathar_, _Rolf Pleisch_ and _Emeric Deutsch_, Apr 01 2008
%E Name edited by _Michel Marcus_, Jan 29 2019