|
|
|
|
1, 2, 13, 130, 1807, 32280, 705421, 18237164, 544505521, 18438430990, 698246022001, 29239344782022, 1341545985079903, 66926098621724300, 3606825675219961657, 208826700420103831480, 12926842112341879416001, 851962999949978920707834, 59561112879709434549509941
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
LINKS
|
|
|
FORMULA
|
a(n) = y(n,n), where y(m+2,n) = (m + n)*y(m+1,n) + y(m,n), with y(0,n)=0, y(1,n)=1 for all n. - Benedict W. J. Irwin, Nov 03 2016
a(n) = round(2*BesselI(n-1,2)*BesselK(2*n-1,2)). - Mark van Hoeij, Nov 08 2022
|
|
MAPLE
|
seq(round(2*BesselI(n-1, 2)*BesselK(2*n-1, 2)), n=1..30); # Mark van Hoeij, Nov 08 2022
|
|
MATHEMATICA
|
Table[DifferenceRoot[Function[{y, m}, {y[2+m]==(m+n)y[1+m]+y[m], y[0]==0, y[1]==1}]][n], {n, 1, 20}] (* Benedict W. J. Irwin, Nov 03 2016 *)
|
|
PROG
|
(Haskell)
a247365 n = a102473 (2 * n - 1) n
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|