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
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.
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.
FORMULA
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
KEYWORD
nonn
AUTHOR
Paul Curtz, Sep 20 2011
STATUS
approved