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.

A258114 E.g.f.: Sum_{n>=0} x^n * cosh(n*x).

Original entry on oeis.org

1, 1, 2, 9, 72, 665, 6960, 85057, 1199744, 19070865, 336372480, 6522635801, 137996694528, 3163206890857, 78085740701696, 2065239729737745, 58263449436979200, 1746433243580269217, 55428341343200280576, 1856918215298125692073, 65483209810866254643200, 2424691204935999655757241
Offset: 0

Views

Author

Paul D. Hanna, May 20 2015

Keywords

Examples

			E.g.f.: A(x) = 1 + x + 2*x^2/2! + 9*x^3/3! + 72*x^4/4! + 665*x^5/5! +...
where A(x) = 1 + x*cosh(x) + x^2*cosh(2*x) + x^3*cosh(3*x) + x^4*cosh(4*x) +...
		

Crossrefs

Cf. A030178.

Programs

  • Mathematica
    CoefficientList[Series[(1-x*Cosh[x])/(1-2*x*Cosh[x]+x^2), {x, 0, 20}], x] * Range[0, 20]! (* Vaclav Kotesovec, May 21 2015 *)
  • PARI
    {a(n) = sum(k=0,n, n!/k! * ((n-k)^k + (-n+k)^k)/2)}
    for(n=0,30,print1(a(n),", "))
    
  • PARI
    {a(n) = local(A=1); A = sum(m=0,n, x^m*cosh(m*x +x*O(x^n))); n!*polcoeff(A,n)}
    for(n=0,30,print1(a(n),", "))
    
  • PARI
    {a(n) = local(X=x+x*O(x^n),A=1); A = (1 - x*cosh(X)) / (1 - 2*x*cosh(X) + x^2); n!*polcoeff(A,n)}
    for(n=0,30,print1(a(n),", "))

Formula

E.g.f.: (1 - x*cosh(x)) / (1 - 2*x*cosh(x) + x^2).
a(n) = Sum_{k=0..n} n!/k! * ((n-k)^k + (-n+k)^k)/2.
a(n) ~ n! * (1-c*cosh(c)) / (2*(cosh(c)+c*(sinh(c)-1)) * c^(n+1)), where c = A030178 = LambertW(1) = 0.56714329040978387299996866... . - Vaclav Kotesovec, May 21 2015