login
The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.

 

Logo

Thanks to everyone who made a donation during our annual appeal!
To see the list of donors, or make a donation, see the OEIS Foundation home page.

Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A076729 a(n) = A001147(n+1) * Integral_{x=0..1} (1 + x^2)^n dx. 10
1, 4, 28, 288, 3984, 70080, 1506240, 38384640, 1133072640, 38038533120, 1431213235200, 59645279232000, 2726781752217600, 135661078090137600, 7295806823277772800, 421717409630060544000, 26071235813929033728000 (list; graph; refs; listen; history; text; internal format)
OFFSET

0,2

COMMENTS

Numerator of the integral where denominator is equal to (2n+1)!! = A001147(n+1).

Also numerator of the integral (1-x^2)^-(n+1/2) for x from 0 to sqrt(1/2). Here the sequence starts at n=1; at n=2 the function is 4.

a(n) = Integral_{x=0..log(1+sqrt2)} cosh(x)^(2*n-1) dx where the denominators are b(n)=(2*n)!/(n!*2^n). E.g., a(3)=28 and b(3)=15; both offsets are 1. - Al Hakanson (hawkuu(AT)excite.com), Mar 02 2004

Self-convolution of A001813. - Vladimir Reshetnikov, Oct 11 2016

LINKS

Table of n, a(n) for n=0..16.

FORMULA

a(n) = 2*n*a(n-1) + (2*n)!/n!.

a(n) = 2^n*Sum_{k=0..n} A001147(k)*A001147(n-k).

a(n) = (2*n+1)!*Sum_{k=0..n} k!*(-2)^k/((2*k+1)!*(n-k)!).

a(n) = (2*n+1)!!*hypergeom([1/2, -n], [3/2], -1). - Vladeta Jovovic, Dec 05 2002

E.g.f.: 1/((1-2*x)*sqrt(1-4*x)). - Vladeta Jovovic, May 11 2003

G.f.: hypergeom([1,1/2],[],4*x)^2  - Mark van Hoeij, May 16 2013

a(n) ~ 2^(2*n+3/2)*n^n/exp(n). - Vaclav Kotesovec, Oct 05 2013

a(n) = (2n+1)!!*Sum_{i=0..n} binomial(n,i)/(2i+1). - John M. Campbell, Feb 06 2016

From Vladimir Reshetnikov, Oct 08 2016: (Start)

a(n) = 2^n*A034430(n) = -(2*n+1)!! * Im(Beta(2, n+1, 1/2))/2.

Recurrence: 2*(3*n+2)*a(n) = a(n+1) + 4*n*(2*n+1)*a(n-1). (End)

Expansion of square of continued fraction 1/(1 - 2*x/(1 - 4*x/(1 - 6*x/(1 - 8*x/(1 - 10*x/(1 - ...)))))). - Ilya Gutkovskiy, Apr 19 2017

EXAMPLE

For n=3, (2n+1)!!=105 and the integral is 96/35 = 288/105, so a(3) = 288.

MAPLE

seq((doublefactorial(2*n+1))*sum((binomial(n, i))/(2*i+1), i=0..n), n=0..20) ; # John M. Campbell, Feb 06 2016

MATHEMATICA

a[n_] := (2n + 1)!!*Integrate[(1 + x^2)^n, {x, 0, 1}]; Table[ a[n], {n, 0, 16}] (* Robert G. Wilson v, Feb 27 2004 *)

Round@Table[-(2 n + 1)!! Im[Beta[2, n + 1, 1/2]]/2, {n, 0, 20}] (* Round is equivalent to FullSimplify here, but is much faster - Vladimir Reshetnikov, Oct 08 2016 *)

PROG

(PARI) a(n)=if(n<0, 0, subst(intformal((1+x^2)^n), x, 1)*(2*n+1)!/2^n/n!)

CROSSREFS

Cf. A001147, A077595, A077745, A086891, A034430, A001813.

Sequence in context: A071212 A090353 A201595 * A078634 A091485 A201354

Adjacent sequences:  A076726 A076727 A076728 * A076730 A076731 A076732

KEYWORD

nonn

AUTHOR

Al Hakanson (hawku(AT)hotmail.com), Oct 28 2002

STATUS

approved

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

License Agreements, Terms of Use, Privacy Policy. .

Last modified January 21 13:52 EST 2020. Contains 331113 sequences. (Running on oeis4.)