login
The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A237622 Interpolation polynomial through n points (0,1), (1,1), ..., (n-2,1) and (n-1,n) evaluated at 2n, a(0)=1. 2
1, 1, 5, 31, 169, 841, 3961, 18019, 80081, 350065, 1511641, 6466461, 27457585, 115892401, 486748081, 2035917451, 8485840801, 35263382881, 146157442201, 604404010981, 2494365759601, 10275832148401, 42264944401681, 173588164506901, 712027089322849 (list; graph; refs; listen; history; text; internal format)
OFFSET

0,3

LINKS

Alois P. Heinz, Table of n, a(n) for n = 0..1000

FORMULA

E.g.f.: exp(x)+2*exp(2*x)*(BesselI(1,2*x)*(x-1)+x*BesselI(0,2*x)).

a(n) ~ sqrt(n)*4^n/sqrt(Pi). - Vaclav Kotesovec, Feb 14 2014

MAPLE

a:= proc(n) option remember; `if`(n<3, [1, 1, 5][n+1],

       (n*(15*n^3-44*n^2+43*n-18) *a(n-1)

        -2*(n-1)*(2*n-3)*(3*n^2-n+2) *a(n-2))/

        ((n-2)*(n+1)*(3*n^2-7*n+6)))

    end:

seq(a(n), n=0..30);

MATHEMATICA

a[n_] := Module[{m}, If[n == 0, 1, InterpolatingPolynomial[Table[{k, If[k == n-1, n, 1]}, {k, 0, n-1}], m] /. m -> 2n]];

a /@ Range[0, 30] (* Jean-Fran├žois Alcover, Dec 22 2020 *)

CROSSREFS

Cf. A002061 (evaluated at n), A158842 (at n+1), A237664 (n+1 points).

Sequence in context: A045904 A034353 A294722 * A239334 A180635 A078526

Adjacent sequences:  A237619 A237620 A237621 * A237623 A237624 A237625

KEYWORD

nonn

AUTHOR

Alois P. Heinz, Feb 10 2014

STATUS

approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified September 21 07:28 EDT 2021. Contains 347596 sequences. (Running on oeis4.)