Logo

Greetings from The On-Line Encyclopedia of Integer Sequences!

Hints

Search: id:A136668
Displaying 1-1 of 1 results found. page 1
     Format: long | short | internal | text      Sort: relevance | references | number      Highlight: on | off
A136668 Triangle of Coefficients of a Bessel polynomial recursion: P(x, n) = 2*(n-1)*P(x, n - 1)/x - n*P(x, n - 2); with substitution: x->1/y. +0
1
1, 0, 1, -2, 0, 2, 0, -11, 0, 8, 8, 0, -74, 0, 48, 0, 119, 0, -632, 0, 384, -48, 0, 1634, 0, -6608, 0, 3840, 0, -1409, 0, 24032, 0, -81984, 0, 46080, 384, 0, -32798, 0, 389312, 0, -1178496, 0, 645120, 0, 18825, 0, -741056, 0, 6966848, 0, -19270656, 0, 10321920, -3840, 0, 666830, 0, -17232128, 0, 137188224, 0 (list; table; graph; listen)
OFFSET

1,4

COMMENT

Row sums:

{1, 1, 0, -3, -18, -129, -1182, -13281, -176478, -2704119, -46909362}

The relationship between the Bessel polynomials and the Hermite

polynomials is better shown by this set of functions.

It took me a while to figure how to do this right: this result is like a toral inverse (1/x) of the Bessel polynomials.

Reversing the coefficient vectors and doing a Gram-Schmidt on them:

f[n_]:=Table[0,{i,1,n}];

b=Table[Join[a[[n]],f[Length[a]-n]],{n,1,Length[a]}];

Table[Length[b[[n]]],{n,1,Length[b]}];

<<LinearAlgebra`Orthogonalization`;

c=GramSchmidt[b];

Gives that these polynomials are at best alternating orthogonal.

REFERENCES

M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1972, 10th edition, (and various reprintings), p. 631.

FORMULA

p(x,0)=1;p(x,1)=1/x; P(x, n) = 2*(n-1)*P(x, n - 1)/x - n*P(x, n - 2); with substitution of 1/y for x.

EXAMPLE

{1},

{0, 1},

{-2, 0, 2},

{0, -11, 0, 8},

{8, 0, -74, 0, 48},

{0, 119, 0, -632, 0, 384},

{-48, 0,1634, 0, -6608, 0, 3840},

{0, -1409, 0, 24032, 0, -81984, 0, 46080},

MATHEMATICA

P[x, 0] = 1; P[x, 1] = 1/x; P[x_, n_] := P[x, n] = 2*(n-1)*P[x, n - 1]/x - n*P[x, n - 2]; Table[ExpandAll[P[x, n] /. x -> 1/y], {n, 0, 10}]; a = Table[CoefficientList[P[x, n] /. x -> 1/y, y], {n, 0, 10}]; Flatten[a]

CROSSREFS

Cf. A106174, A123956.

Adjacent sequences: A136665 A136666 A136667 this_sequence A136669 A136670 A136671

Sequence in context: A099554 A107729 A113400 this_sequence A057498 A137949 A019214

KEYWORD

uned,tabl,sign

AUTHOR

Roger L. Bagula (rlbagulatftn(AT)yahoo.com), Apr 03 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 January 7 17:35 EST 2009. Contains 152824 sequences.


AT&T Labs Research