%I #12 May 08 2018 15:11:55
%S 2,3,2,-1,4,-3,10,-691,280,-10851,438670,-1222277,3418052,-1181820455,
%T 1077690978,-23749461029,137853460416080,-84802531453387,
%U 541314450257070,-26315271553053477373,761798417598805340
%N Numerator of (2n+1)(2n+2) B_{2n}, where B_n are the Bernoulli numbers. Also numerators of the asymptotic expansion of the polygamma function psi'''(z).
%H Vincenzo Librandi, <a href="/A076549/b076549.txt">Table of n, a(n) for n = 0..100</a>
%H M. Abramowitz and I. A. Stegun, eds., <a href="http://www.convertit.com/Go/ConvertIt/Reference/AMS55.ASP">Handbook of Mathematical Functions</a>, National Bureau of Standards, Applied Math. Series 55, Tenth Printing, 1972 [alternative scanned copy].
%H Abramowitz and Stegun, Handbook of mathematical functions, <a href="http://members.fortunecity.com/aands/page_260.jpg">p. 260</a> (6.4.14)
%H Abramowitz and Stegun, <a href="http://www.convertit.com/Go/ConvertIt/Reference/AMS55.ASP">Handbook of mathematical functions</a>, <a href="http://members.fortunecity.com/aands/page_260.jpg">p. 260</a> (6.4.14)
%t Join[{2, 3}, Table[Numerator[(4 n^2 + 6 n + 2) BernoulliB[2 n]], {n, 1, 30}]] (* _Vincenzo Librandi_, Aug 02 2013 *)
%Y Denominators are in A006956. Cf. A002427/A006955.
%K easy,frac,sign
%O 0,1
%A _Ralf Stephan_, Oct 19 2002