login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A191499 E.g.f. sqrt(1+tan(2*x)). 1
1, 1, -1, 11, -47, 601, -5521, 86771, -1296287, 25482481, -527699041, 12800059931, -335639304527, 9794548687561, -308817517422961, 10573293809103491, -388317397661640767, 15275057087004591841, -639584224876056953281, 28426125263460829489451, -1335823888802587475761007 (list; graph; refs; listen; history; text; internal format)
OFFSET

0,4

LINKS

Robert Israel, Table of n, a(n) for n = 0..390

FORMULA

a(n) = Sum_{m=0,..,(n-1)/2} ( binomial(2*n-4*m-2,n-2*m-1)*Sum_{j=0,..,2*m} ( binomial(j+n-2*m-1,n-2*m-1)*(j+n-2*m)!*2^(6*m-n-j+1)*(-1)^(m+n+j+1)*stirling2(n,j+n-2*m)))/(n-2*m) ) ), n>0, a(0)=1.

MAPLE

S:= series(sqrt(1+tan(2*x)), x, 31):

seq(coeff(S, x, j)*j!, j=0..30); # Robert Israel, Feb 28 2017

MATHEMATICA

With[{nn = 50}, CoefficientList[Series[Sqrt[1 + Tan[2*x]], {x, 0, nn}], x]*Range[0, nn]!] (* G. C. Greubel, Feb 28 2017 *)

PROG

(Maxima)

a(n):=sum((binomial(2*n-4*m-2, n-2*m-1)*sum(binomial(j+n-2*m-1, n-2*m-1)*(j+n-2*m)!*2^(6*m-n-j+1)*(-1)^(m+n+j+1)*stirling2(n, j+n-2*m), j, 0, 2*m))/(n-2*m), m, 0, (n-1)/2);

(PARI) x='x + O('x^50); Vec(serlaplace(sqrt(1 + tan(2*x)))) \\ G. C. Greubel, Feb 28 2017

CROSSREFS

Sequence in context: A354590 A067355 A138362 * A072372 A230982 A024530

Adjacent sequences: A191496 A191497 A191498 * A191500 A191501 A191502

KEYWORD

sign

AUTHOR

Vladimir Kruchinin, Jun 03 2011

STATUS

approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified January 31 06:19 EST 2023. Contains 359947 sequences. (Running on oeis4.)