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

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A195287 a(n) = (A091137(n)/n!) * Integral_{u=-1..1} u*(u+1)*...*(u+n-1) du. 3
2, 0, 4, 8, 232, 448, 18224, 35424, 1036064, 2025472, 130960832, 257072000, 689908475264, 1358275350528, 8031885897472, 15847920983552, 7981032500085248, 15774370258485248, 12448755354530366464 (list; graph; refs; listen; history; text; internal format)
OFFSET

0,1

COMMENTS

Numerators of the second row of an array based on Adams numerical integration. Take q!*s(m,q) = Integral_{-m-1..1} u*(u+1)*...*(u+q-1) du. a(n) is in the second row (case m=0) numerators of s(m,q) in the comments.

The unreduced array s(m,q), (m=-1,0,1,...,   columns q=0,1,2,...) is

1,   1/2,   5/12,   9/24,   251/720, 475/1440,  = A002657(n)/A091137(n),

2,     0,   4/12,   8/24,   232/720, 448/1440,  = a(n)/A091137(n),

3,  -3/2,   9/12,   9/24,   243/720, 459/1440,

4,  -8/2,  32/12,      0,   224/720, 448/1440,

5, -15/2,  85/12, -55/24,   475/720, 475/1440,

6, -24/2, 180/12, -216/24, 2376/720,        0.

Column numerators: A000027, -A067998(n), A152064(n), A157371(n), A165281(n).

Page 56 of the reference.

(*) 2/2 = 1,

    2/2 + 0 = 1,

    2/3 + 0 + 1/3 = 1,

    2/4 + 0 + 1/6 + 1/3 = 1. Reduced.

REFERENCES

P. Curtz, Intégration numérique des systèmes differentiels à conditions initiales, Centre de Calcul Scientifique de l'Armement, Arcueil, 1969.

LINKS

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

FORMULA

b(n) = a(n)/A091137(n).

b(0)/2 = 1,

b(0)/2 + b(1) = 1,

b(0)/3 + b(1)/2 + b(2) = 1,

b(0)/4 + b(1)/3 + b(2)/2 + b(3) = 1.

First vertical denominators: A028310(n) + 1. See A104661.

Values in (*).

MAPLE

A195287 := proc(n)

        mul(u+i, i=0..n-1) ;

        int(%, u=-1..1) ;

        %/n!*A091137(n) ;

end proc:

seq(A195287(n), n=0..20) ; # R. J. Mathar, Oct 02 2011

MATHEMATICA

(* a7 = A091137 *) a7[n_] := a7[n] = Product[d, {d, Select[Divisors[n] + 1, PrimeQ]}]*a7[n-1]; a7[0]=1; a[n_] := a7[n]/n!*Integrate[ Pochhammer[u, n], {u, -1, 1}]; Table[a[n], {n, 0, 18}] (* Jean-François Alcover, Aug 13 2012 *)

CROSSREFS

Sequence in context: A101189 A295321 A001443 * A070015 A021492 A342201

Adjacent sequences:  A195284 A195285 A195286 * A195288 A195289 A195290

KEYWORD

nonn

AUTHOR

Paul Curtz, Sep 20 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 | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified July 24 05:05 EDT 2021. Contains 346273 sequences. (Running on oeis4.)