|
| |
|
|
A121351
|
|
An A053984-like Bessel-Binet recursion found by Bob Hanlon's new survey program: a[n] = (3*n + 1)*a[n - 1] - a[n - 2].
|
|
0
| |
|
|
0, 1, 7, 69, 890, 14171, 268359, 5889727, 146974816, 4109405121, 127244583935, 4322206448669, 159794394016818, 6387453554224051, 274500708437617375, 12620645134576175199, 618137110885794967376, 32130509120926762128353
(list; graph; refs; listen; history; internal format)
|
|
|
|
OFFSET
| 1,3
|
|
|
COMMENTS
| I had theorized that there were a related kind of Bessel-Binets of the form: a[n]=(a0*n+c0)*a[n-1]+b0*a[n-2] Where the a0,b0 and c0 were Integers. Bob Hanlon found ones with a0=2 in a program he wrote. A053984 was already in OEIS. This one is a0=3
|
|
|
FORMULA
| a(n) = (3*n + 1)*a(n - 1) - a(n - 2) a(n)=(1/3) pi BesselJ[4/3 + n, 2/3] BesselY[4/3, 2/3] - (1/3) pi BesselJ[4/3, 2/3] BesselY[4/3 + n, 2/3]
|
|
|
MATHEMATICA
| f[n_Integer] = Module[{a}, a[n] /. RSolve[{a[n] == (3*n + 1)*a[n - 1] - a[n - 2], a[0] == 0, a[1] == 1}, a[n], n][[1]] // FullSimplify] Rationalize[N[Table[f[n], {n, 0, 25}], 100], 0]
|
|
|
CROSSREFS
| Cf. A053984, A001503.
Sequence in context: A180911 A084774 A025757 * A059321 A077683 A082797
Adjacent sequences: A121348 A121349 A121350 * A121352 A121353 A121354
|
|
|
KEYWORD
| nonn,uned
|
|
|
AUTHOR
| Roger Bagula and Bob Hanlon (rlbagulatftn(AT)yahoo.com), Sep 05 2006
|
| |
|
|
