Logo

Greetings from The On-Line Encyclopedia of Integer Sequences!

Hints

Search: id:A142983
Displaying 1-1 of 1 results found. page 1
     Format: long | short | internal | text      Sort: relevance | references | number      Highlight: on | off
A142983 a(1) = 1, a(2) = 2, a(n+2) = 2*a(n+1) + (n+1)*(n+2)*a(n). +0
13
1, 2, 10, 44, 288, 1896, 15888, 137952, 1419840, 15255360, 186693120, 2387093760, 33898314240, 502247692800, 8123141376000, 136785729024000, 2483065912320000, 46822564905984000, 942853671825408000 (list; graph; listen)
OFFSET

1,2

COMMENT

This is the case m = 1 of the general recurrence a(1) = 1, a(2) = 2*m, a(n+2) = 2*m*a(n+1) + (n+1)*(n+2)*a(n) (we suppress the dependence of a(n) on m), which arises when accelerating the convergence of the series 1/2 + 1/2*sum {k = 1..inf} (-1)^(k+1)/(k*(k+1)) = log(2). For other cases see A142984 (m=2), A142985 (m=3), A142986 (m=4) and A142987 (m=5). The solution to the general recurrence may be expressed as a sum: a(n) = n!*p_m(n+1)*sum {k = 1..n} (-1)^(k+1)/(p_m(k)*p_m(k+1)), where p_m(x) = sum {k = 1..m} 2^(k-1)*C(m-1,k-1)*C(x,k) is the polynomial that gives the regular polytope numbers for the m-dimensional cross polytope as defined by [Kim](see A142978). The first few values are p_1(x) = x, p_2(x) = x^2, p_3(x) = (2*x^3+x)/3 and p_4(x) = (x^4+2*x^2)/3.

The polynomial p_m(x) is the unique polynomial solution of the difference equation x*(f(x+1)-f(x-1)) = 2*m*f(x), normalised so that f(1) = 1. The o.g.f. for the p_m(x) is 1/2*((1+t)/(1-t))^x = 1/2 + x*t + x^2*t^2 + (2*x^3+x)/3*t^3 + ... . Thus p_m(x) is, apart from a constant factor, the Meixner polynomial of the first kind M_m(x;b,c) at b = 0, c = -1, also known as a Mittag-Leffler poynomial.

The general recurrence in the first paragraph above has a second solution b(n) = n!*p_m(n+1) with b(1) = 2*m, b(2) = m^2+2. Hence the behaviour of a(n) for large n is given by lim n-> infinity a(n)/b(n) = sum {k = 1..inf} (-1)^(k+1)/(p_m(k)*p_m(k+1)) = 1/((2*m)+ 1*2/((2*m)+ 2*3/((2*m)+ 3*4/((2*m)+...+ *n*(n+1)/((2*m)+...))))) = 1 + (-1)^(m+1) * (2*m)*(log(2) - (1-1/2+1/3- ...+ (-1)^(m+1)/m)), where the final equality follows by a result of Ramanujan (see [Berndt, Chapter 12, Entry 32(i)]).

See A142979, A142988 and A142992 for similar results. For corresponding results for Napier's constant e, the constant zeta(2) and Apery's constant zeta(3) refer to A000522, A142995 and A143003 respectively.

REFERENCES

Bruce C. Berndt, Ramanujan's Notebooks Part II, Springer-Verlag.

LINKS

Hyun Kwang Kim, On regular polytope numbers

Weisstein Eric, W. Meixner polynomial of the first kind

Weisstein Eric, W. Mittag-Leffler polynomial

FORMULA

a(n) = n!*p(n+1)*sum {k = 1..n} (-1)^(k+1)/(p(k)*p(k+1)), where p(n) = n. Recurrence: a(1) = 1, a(2) = 2, a(n+2) = 2*a(n+1)+(n+1)*(n+2)*a(n). The sequence b(n):= n!*p(n+1) satisfies the same recurrence with b(1) = 2, b(2) = 6. Hence we obtain the finite continued fraction expansion a(n)/b(n) = 1/(2 +1*2/(2 +2*3/(2 +3*4/(2 +...+(n-1)*n/2)))), for n >=2. The behaviour of a(n) for large n is given by lim n -> infinity a(n)/b(n) = 1/(2 +1*2/(2 +2*3/(2 +3*4/(2 +...+n*(n+1)/(2 +...))))) = sum {k = 1..inf} (-1)^(k+1)/(k*(k+1)) = 2*log(2) - 1;

MAPLE

a := n -> (n+1)!*sum ((-1)^(k+1)/(k*(k+1)), k = 1..n): seq(a(n), n = 1..20);

CROSSREFS

Cf. A000522, A142979, A142984, A142985, A142986, A142987, A142988, A142992.

Sequence in context: A105485 A151313 A144896 this_sequence A065805 A145239 A068561

Adjacent sequences: A142980 A142981 A142982 this_sequence A142984 A142985 A142986

KEYWORD

easy,nonn

AUTHOR

Peter Bala (pbala(AT)toucansurf.com), Jul 17 2008

page 1

Search completed in 0.002 seconds

Lookup | Welcome | Find friends | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Transforms | Puzzles | Hot | Classics
More pages | Superseeker | Maintained by N. J. A. Sloane (njas@research.att.com)

Last modified December 17 13:29 EST 2009. Contains 170826 sequences.


AT&T Labs Research