%I #26 Apr 16 2022 11:35:52
%S 0,64,1600,9216,30976,78400,166464,313600,541696,876096,1345600,
%T 1982464,2822400,3904576,5271616,6969600,9048064,11560000,14561856,
%U 18113536,22278400,27123264,32718400,39137536,46457856,54760000
%N Squares of even octagonal numbers.
%H <a href="/index/Rec#order_05">Index entries for linear recurrences with constant coefficients</a>, signature (5,-10,10,-5,1).
%F G.f.: 64*x*(1+20*x+29*x^2+4*x^3)/(1-x)^5. - _Colin Barker_, Jan 06 2012
%F a(n) = 16n^2*(3n-1)^2. - _Vincenzo Librandi_, Jan 07 2012
%F E.g.f.: 16*exp(x)*x*(4 + 46*x + 48*x^2 + 9*x^3). - _Stefano Spezia_, Apr 16 2022
%t Table[16*n^2*(3*n-1)^2,{n,0,30}] (* _Vincenzo Librandi_, Jan 07 2012 *)
%t LinearRecurrence[{5,-10,10,-5,1},{0,64,1600,9216,30976},30] (* _Harvey P. Dale_, Nov 27 2015 *)
%o (Magma) [16*n^2*(3*n-1)^2: n in [1..50]]; // _Vincenzo Librandi_, Jan 07 2012
%o (PARI) a(n) = 16*n^2*(3*n-1)^2 \\ _Vincenzo Librandi_, Jan 07 2012
%Y Cf. A000567, A014641, A014642, A014793.
%K nonn,easy
%O 0,2
%A _Mohammad K. Azarian_
%E More terms from _Patrick De Geest_, Aug 2000
%E a(8) corrected by _Vincenzo Librandi_, Jan 07 2012