cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-2 of 2 results.

A062138 Coefficient triangle of generalized Laguerre polynomials n!*L(n,5,x)(rising powers of x).

Original entry on oeis.org

1, 6, -1, 42, -14, 1, 336, -168, 24, -1, 3024, -2016, 432, -36, 1, 30240, -25200, 7200, -900, 50, -1, 332640, -332640, 118800, -19800, 1650, -66, 1, 3991680, -4656960, 1995840, -415800, 46200, -2772, 84, -1, 51891840, -69189120
Offset: 0

Views

Author

Wolfdieter Lang, Jun 19 2001

Keywords

Comments

The row polynomials s(n,x) := n!*L(n,5,x)= sum(a(n,m)*x^m,m=0..n) have e.g.f. exp(-z*x/(1-z))/(1-z)^6. They are Sheffer polynomials satisfying the binomial convolution identity s(n,x+y) = sum(binomial(n,k)*s(k,x)*p(n-k,y),k=0..n), with polynomials sum(|A008297(n,m)|*(-x)^m, m=1..n), n >= 1 and p(0,x)=1 (for Sheffer polynomials see A048854 for S. Roman reference).
These polynomials appear in the radial part of the l=2 (d-wave) eigen functions for the discrete energy levels of the H-atom. See Messiah reference.
For m=0..5 the (unsigned) column sequences (without leading zeros) are: A001725(n+5), A062148-A062152. Row sums (signed) give A062191; row sums (unsigned) give A062192.
The unsigned version of this triangle is the triangle of unsigned 3-Lah numbers A143498. - Peter Bala, Aug 25 2008

Examples

			Triangle begins:
  {1};
  {6, -1};
  {42, -14, 1};
  {336, -168, 24, -1};
  ...
2!*L(2, 5, x) = 42-14*x+x^2.
		

References

  • A. Messiah, Quantum mechanics, vol. 1, p. 419, eq.(XI.18a), North Holland, 1969.

Crossrefs

For m=0..5 the (unsigned) column sequences (without leading zeros) are: A001725(n+5), A062148, A062149, A062150, A062151, A062152.
Row sums (signed) give A062191, row sums (unsigned) give A062192.
Cf. A143498.

Programs

  • Mathematica
    Flatten[Table[((-1)^m)*n!*Binomial[n+5,n-m]/m!,{n,0,8},{m,0,n}]] (* Indranil Ghosh, Feb 24 2017 *)
  • PARI
    tabl(nn) = {for (n=0, nn, for (m=0, n, print1(((-1)^m)*n!*binomial(n+5, n-m)/m!, ", "); ); print(); ); } \\ Indranil Ghosh, Feb 24 2017
    
  • PARI
    row(n) = Vecrev(n!*pollaguerre(n, 5)); \\ Michel Marcus, Feb 06 2021
    
  • Python
    import math
    f=math.factorial
    def C(n, r):return f(n)//f(r)//f(n-r)
    i=-1
    for n in range(26):
        for m in range(n+1):
            i += 1
            print(str(i)+" "+str(((-1)**m)*f(n)*C(n+5, n-m)//f(m))) # Indranil Ghosh, Feb 24 2017

Formula

T(n, m) = ((-1)^m)*n!*binomial(n+5, n-m)/m!.
E.g.f. for m-th column: ((-x/(1-x))^m)/(m!*(1-x)^6), m >= 0.

A062152 Sixth (unsigned) column of triangle A062138 (generalized a=5 Laguerre).

Original entry on oeis.org

1, 66, 2772, 96096, 3027024, 90810720, 2663781120, 77630192640, 2270683134720, 67111301537280, 2013339046118400, 61498356317798400, 1916698771904716800, 61039483966811750400, 1988143192061868441600
Offset: 0

Views

Author

Wolfdieter Lang, Jun 19 2001

Keywords

Examples

			a(2) = (2+5)! * binomial(2+10,10) / 5! = (5040 * 66) / 120 = 2772. - _Indranil Ghosh_, Feb 24 2017
		

Crossrefs

Programs

  • Magma
    [Factorial(n+5)*Binomial(n+10,10)/Factorial(5): n in [0..20]]; // G. C. Greubel, May 11 2018
  • Mathematica
    Table[(n+5)!*Binomial[n+10,10]/5!,{n,0,14}] (* Indranil Ghosh, Feb 24 2017 *)
  • PARI
    a(n) = (n+5)!*binomial(n+10,10)/5! \\ Indranil Ghosh, Feb 24 2017
    
  • Python
    import math
    f=math.factorial
    def C(n, r):return f(n)/f(r)/f(n-r)
    def A062152(n): return f(n+5)*C(n+10, 10)/f(5) # Indranil Ghosh, Feb 24 2017
    

Formula

a(n) = A062138(n+5, 5).
a(n) = (n+5)!*binomial(n+10, 10)/5!.
E.g.f.: N(5;5, x)/(1-x)^16 with N(5;5, x) := Sum_{k=0..5} A062190(5, k)* x^k = 1 + 50*x + 450*x^2 + 1200*x^3 + 1050*x^4 + 252*x^5.
If we define f(n,i,x) = Sum_{k=i..n} Sum_{j=i..k} binomial(k,j)* Stirling1(n,k)*Stirling2(j,i)*x^(k-j) then a(n-10) = (-1)^n*f(n,10,-6), (n>=10). - Milan Janjic, Mar 01 2009
Showing 1-2 of 2 results.