cp's OEIS Frontend

A126958 Define an array by d(m, 0) = 1, d(m, 1) = m; d(m, k) = (m - k + 1) d(m+1, k-1) - (k-1) (m+1) d(m+2, k-2). Sequence gives d(n,4).

%I A126958 #9 Oct 21 2022 21:13:53
%S A126958 36,24,-60,-216,-420,-624,-756,-720,-396,360,1716,3864,7020,11424,
%T A126958 17340,25056,34884,47160,62244,80520,102396,128304,158700,194064,
%U A126958 234900,281736,335124,395640,463884,540480,626076,721344,826980,943704,1072260,1213416,1367964
%N A126958 Define an array by d(m, 0) = 1, d(m, 1) = m; d(m, k) = (m - k + 1) d(m+1, k-1) - (k-1) (m+1) d(m+2, k-2). Sequence gives d(n,4).
%D A126958 V. van der Noort and N. J. A. Sloane, Paper in preparation, 2007.
%H A126958 G. C. Greubel, <a href="/A126958/b126958.txt">Table of n, a(n) for n = 0..1000</a>
%H A126958 <a href="/index/Rec#order_05">Index entries for linear recurrences with constant coefficients</a>, signature (5,-10,10,-5,1).
%F A126958 a(n) = (n+3)*(n+1)*(n^2 -10*n +12).
%F A126958 From _G. C. Greubel_, Jan 29 2020: (Start)
%F A126958 G.f.: 12*(3 -13*x +15*x^2 -3*x^3)/(1-x)^5.
%F A126958 E.g.f.: (36 -12*x -36*x^2 +x^4)*exp(x). (End)
%p A126958 seq( (n+3)*(n+1)*(n^2 -10*n +12), n=0..40); # _G. C. Greubel_, Jan 29 2020
%t A126958 Table[(n+3)*(n+1)*(n^2 -10*n +12), {n,0,40}] (* _G. C. Greubel_, Jan 29 2020 *)
%o A126958 (PARI) vector(41, n, my(m=n-1); (m+3)*(m+1)*(m^2 -10*m +12)) \\ _G. C. Greubel_, Jan 29 2020
%o A126958 (Magma) [(n+3)*(n+1)*(n^2 -10*n +12): n in [0..40]]; // _G. C. Greubel_, Jan 29 2020
%o A126958 (Sage) [(n+3)*(n+1)*(n^2 -10*n +12) for n in (0..40)] # _G. C. Greubel_, Jan 29 2020
%o A126958 (GAP) List([0..40], n-> (n+3)*(n+1)*(n^2 -10*n +12)); # _G. C. Greubel_, Jan 29 2020
%Y A126958 A row of A105937.
%K A126958 sign,easy
%O A126958 0,1
%A A126958 Vincent v.d. Noort, Mar 21 2007