cp's OEIS Frontend

A126935 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,3).

%I A126935 #11 Oct 18 2022 15:30:54
%S A126935 0,-12,-24,-30,-24,0,48,126,240,396,600,858,1176,1560,2016,2550,3168,
%T A126935 3876,4680,5586,6600,7728,8976,10350,11856,13500,15288,17226,19320,
%U A126935 21576,24000,26598,29376,32340,35496,38850,42408,46176,50160,54366,58800,63468,68376,73530
%N A126935 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,3).
%D A126935 V. van der Noort and N. J. A. Sloane, Paper in preparation, 2007.
%H A126935 G. C. Greubel, <a href="/A126935/b126935.txt">Table of n, a(n) for n = 0..1000</a>
%H A126935 <a href="/index/Rec#order_04">Index entries for linear recurrences with constant coefficients</a>, signature (4,-6,4,-1).
%F A126935 a(n) = n*(n+2)*(n-5).
%F A126935 From _G. C. Greubel_, Jan 29 2020: (Start)
%F A126935 G.f.: (-6)*x*(2 - 4*x + x^2)/(1-x)^4.
%F A126935 E.g.f.: x*(-12 + x^2)*exp(x). (End)
%p A126935 seq( n*(n+2)*(n-5), n=0..50); # _G. C. Greubel_, Jan 29 2020
%t A126935 T[n_, k_]:= T[n, k]= If[k==0, 1, If[k==1, n, (n-k+1)*T[n+1, k-1] - (k-1)*(n+1)* T[n+2, k-2]]]; Table[T[n, 3], {n,0,50}] (* _G. C. Greubel_, Jan 29 2020 *)
%o A126935 (PARI) vector(50, n, my(m=n-1); m*(m+2)*(m-5) ) \\ _G. C. Greubel_, Jan 29 2020
%o A126935 (Magma) [n*(n+2)*(n-5): n in [0..50]]; // _G. C. Greubel_, Jan 29 2020
%o A126935 (Sage) [n*(n+2)*(n-5) for n in (0..50)] # _G. C. Greubel_, Jan 29 2020
%Y A126935 A row of A105937.
%K A126935 sign,easy
%O A126935 0,2
%A A126935 Vincent v.d. Noort, Mar 21 2007