%I #19 Jan 22 2018 02:39:10
%S 0,1,-1,4,-16,91,-616,4894,-44416,453541,-5145856,64222984,-874402816,
%T 12897183391,-204862698496,3486535686274,-63292764553216,
%U 1220795601029641,-24931903045697536,537463528325234764,-12196043454064623616,290588368323926512291,-7253387436695186046976
%N Expansion of e.g.f.: sinh(log(1+sinh(x))).
%H G. C. Greubel, <a href="/A009568/b009568.txt">Table of n, a(n) for n = 0..438</a>
%F a(2*n) = A006154(2*n)/2, n>0. - _Ralf Stephan_, Apr 29 2004
%t With[{nn=20},CoefficientList[Series[Sinh[Log[1+Sinh[x]]],{x,0,nn}],x] Range[0,nn]!] (* _Harvey P. Dale_, Dec 08 2017 *)
%o (PARI) x='x+O('x^30); concat([0], Vec(serlaplace( sinh(log(1+sinh(x)))))) \\ _G. C. Greubel_, Jan 21 2018
%K sign,easy
%O 0,4
%A _R. H. Hardin_
%E Extended with signs by _Olivier GĂ©rard_, Mar 15 1997
%E Prior Mathematica program replaced and definition clarified by _Harvey P. Dale_, Dec 08 2017
%E Terms a(20) to a(22) added by _G. C. Greubel_, Jan 21 2018