%I #19 Sep 08 2022 08:45:44
%S 0,4,19,54,124,250,459,784,1264,1944,2875,4114,5724,7774,10339,13500,
%T 17344,21964,27459,33934,41500,50274,60379,71944,85104,100000,116779,
%U 135594,156604,179974,205875,234484,265984,300564,338419,379750,424764
%N a(n) = n^2*(n^2 + 15)/4.
%H Vincenzo Librandi, <a href="/A159833/b159833.txt">Table of n, a(n) for n = 0..1000</a>
%H <a href="/index/Rec#order_05">Index entries for linear recurrences with constant coefficients</a>, signature (5,-10,10,-5,1).
%F a(n) = A008488(n+1)-2 = 4 - 15*A000292(n+1) + 6*A000332(n+4) + 20*A000217(n+1) - 15*(n+1).
%F a(n) = 5*a(n-1) - 10*a(n-2) + 10*a(n-3) - 5*a(n-4) + a(n-5).
%F G.f.: -x*(1+x)*(4*x^2-5*x+4)/(x-1)^5.
%F E.g.f.: x*(16 +22*x +6*x^2 +x^3)*exp(x)/4. - _G. C. Greubel_, May 19 2018
%p seq(n^2*(n^2+15)/4,n=0..80)
%t CoefficientList[Series[-x*(1 + x)*(4*x^2 - 5*x + 4)/(x-1)^5, {x, 0, 40}], x] (* _Vincenzo Librandi_, Dec 18 2012 *)
%t LinearRecurrence[{5,-10,10,-5,1},{0,4,19,54,124},40] (* _Harvey P. Dale_, May 30 2016 *)
%o (Magma) [n^2 * (n^2 + 15)/4: n in [0..40]]; // _Vincenzo Librandi_, Dec 18 2012
%o (PARI) for(n=0, 30, print1(n^2*(n^2 +15)/4, ", ")) \\ _G. C. Greubel_, May 19 2018
%Y Cf. A008488, A000292, A000332, A000217.
%K nonn,easy
%O 0,2
%A _R. J. Mathar_, Apr 23 2009
|