%I M4620 N1973 #33 Oct 15 2023 01:42:25
%S 1,1,9,39,1141,12721,804309,17113719,1886573641,65373260641,
%T 11127809595009,570506317184199,138730500808639741,9867549661639871761
%N From the expansion of sinh(x) / cos(x): a(n) = odd part of A002084(n).
%D N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
%D N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
%H T. D. Noe, <a href="/A002085/b002085.txt">Table of n, a(n) for n = 0..50</a>
%H Eric Duchene, Aviezri S. Fraenkel, Vladimir Gurvich, Nhan Bao Ho, Clark Kimberling, and Urban Larsson, <a href="http://www.wisdom.weizmann.ac.il/~fraenkel/Papers/WythoffWisdomJune62016.pdf">Wythoff Wisdom</a>, 43 pages, no date, unpublished.
%H Eric Duchene, Aviezri S. Fraenkel, Vladimir Gurvich, Nhan Bao Ho, Clark Kimberling, and Urban Larsson, <a href="/A001950/a001950.pdf">Wythoff Wisdom</a>, unpublished, no date. [Cached copy, with permission]
%H J. M. Gandhi, <a href="http://dx.doi.org/10.4153/CMB-1970-059-9">The coefficients of sinh x/ cos x</a>. Canad. Math. Bull. 13 1970 305-310.
%F Divide A002084(n) by 4^(floor((n+1)/4)).
%t a[ n_] := If[ n < 0, 0, SeriesCoefficient[ Sinh[x] / Cos[x], {x, 0, 2 n + 1}] (2 n + 1)! / 2^(n + Mod[n, 2])] (* _Michael Somos_, Apr 10 2011 *)
%o (PARI) {a(n) = local(A, m); if( n<0, 0, m = 2*n + 1; A = x * O(x^m) ; 2^(n + 1 - n%2) * m! * polcoeff( sinh( x/2 + A) / cos( x/2 + A), m))} /* _Michael Somos_, Apr 10 2011 */
%Y Cf. A002084.
%K nonn,easy
%O 0,3
%A _N. J. A. Sloane_