%I #32 Feb 19 2023 03:33:24
%S 5,77,221,437,725,1085,1517,2021,2597,3245,3965,4757,5621,6557,7565,
%T 8645,9797,11021,12317,13685,15125,16637,18221,19877,21605,23405,
%U 25277,27221,29237,31325,33485,35717,38021,40397,42845,45365,47957,50621,53357,56165,59045
%N a(n) = (6*n+1)*(6*n+5).
%H T. D. Noe, <a href="/A001513/b001513.txt">Table of n, a(n) for n = 0..1000</a>
%H <a href="/index/Rec#order_03">Index entries for linear recurrences with constant coefficients</a>, signature (3,-3,1).
%F Sum_{k>=0} 1/a(k) = Pi/(8*sqrt(3)) = 0.22672492... - _Jaume Oliver Lafont_, May 30 2010
%F a(n) = 72*n + a(n-1) with a(0)=5. - _Vincenzo Librandi_, Nov 12 2010
%F G.f.: (-5 - 62*x - 5*x^2) / (x-1)^3. - _R. J. Mathar_, Jan 19 2013
%F From _Amiram Eldar_, Feb 19 2023: (Start)
%F a(n) = A016921(n)*A016969(n).
%F Sum_{n>=0} (-1)^n/a(n) = log(2+sqrt(3))/(4*sqrt(3)).
%F Product_{n>=0} (1 - 1/a(n)) = 2*cos(sqrt(5)*Pi/6).
%F Product_{n>=0} (1 + 1/a(n)) = 2*cos(sqrt(3)*Pi/6). (End)
%t a[n_] := (6*n + 1)*(6*n + 5); Array[a, 40, 0] (* _Amiram Eldar_, Feb 19 2023 *)
%o (PARI) a(n)=(6*n+1)*(6*n+5) \\ _Charles R Greathouse IV_, Jun 17 2017
%Y Cf. A016921, A016969.
%K nonn,easy
%O 0,1
%A _N. J. A. Sloane_