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!)
A136247 Triangle T(n,k) read by rows: coefficient [x^k] of the polynomial h(n,x) with h(0,x)=1, h(1,x)=1-x and recursively h(n,x) = 1 + n -(1-x)*(1-h(n-1,x)) - n*h(n-2,x). 3

%I #10 Jul 22 2015 02:52:57

%S 1,1,-1,1,-1,1,1,2,2,-1,1,6,-4,-3,1,1,-4,-20,6,4,-1,1,-40,8,44,-8,-5,

%T 1,1,-12,188,-6,-80,10,6,-1,1,308,136,-546,-10,130,-12,-7,1,1,416,

%U -1864,-628,1256,50,-196,14,8,-1,1,-2664,-3640,6696,1984,-2506,-126,280,-16,-9,1

%N Triangle T(n,k) read by rows: coefficient [x^k] of the polynomial h(n,x) with h(0,x)=1, h(1,x)=1-x and recursively h(n,x) = 1 + n -(1-x)*(1-h(n-1,x)) - n*h(n-2,x).

%C Row sums are 1, 0, 1, 4, 1, -14, 1, 106, 1, -944, 1, 10396, 1, -135134, 1, 2027026, 1, -34459424, 1, 654729076, 1...

%C [Row sums s(n) appear to obey s(n) -2*s(n-1) +(n+1)*s(n-2) +2*(1-n)*s(n-3) +(n-2)* s(n-4)=0. - _R. J. Mathar_, Dec 04 2011]

%D Harry Hochstadt, The Functions of Mathematical Physics, Dover, New York, 1986, pp. 8, 42-43.

%e 1;

%e 1, -1;

%e 1, -1, 1;

%e 1, 2, 2, -1;

%e 1, 6, -4, -3, 1;

%e 1, -4, -20,6, 4, -1;

%e 1, -40, 8, 44, -8, -5, 1;

%e 1, -12, 188, -6, -80,10, 6, -1;

%e 1, 308, 136, -546, -10, 130, -12, -7, 1;

%e 1, 416, -1864, -628, 1256, 50, -196, 14, 8, -1;

%e 1, -2664, -3640, 6696, 1984, -2506, -126,280, -16, -9, 1;

%p h := proc(n,x)

%p if n = 0 then

%p 1 ;

%p elif n = 1 then

%p 1-x ;

%p else

%p 1+n-(1-x)*(1-procname(n-1,x)) -n*procname(n-2,x) ;

%p expand(%) ;

%p end if;

%p end proc:

%p A136247 := proc(n,k)

%p coeftayl(h(n,x),x=0,k) ;

%p end proc:

%p seq(seq(A136247(n,k),k=0..n),n=0..12) ; # _R. J. Mathar_, Dec 04 2011

%t Clear[h, a, n, x, y, c, d] (*Solve linear Shabat transform for Hermite type recursion*) Solve[c*x0 + d - x*(c*x1 + d) + n*(c*x2 + d) == 0, x0] c = -1; d = 1; Solve[y = c*x + d == 0, x] h[x, 0] = 1; h[x, 1] = 1 - x; h[x_, n_] := h[x, n] = -(-1 - n + (1 - x) - (1 - x)* h[ x, n - 1] + n *h[x, n - 2]); Table[ExpandAll[h[x, n]], {n, 0, 10}]; a = Table[CoefficientList[h[x, n], x], {n, 0, 10}]; Flatten[a] Table[Apply[Plus, CoefficientList[h[x, n], x]], {n, 0, 10}];

%Y Cf. A137286.

%K easy,tabl,sign

%O 0,8

%A _Roger L. Bagula_, Mar 17 2008

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 March 28 18:04 EDT 2024. Contains 371254 sequences. (Running on oeis4.)