login
This site is supported by donations to The OEIS Foundation.

 

Logo

Annual Appeal: Please make a donation to keep the OEIS running. In 2018 we replaced the server with a faster one, added 20000 new sequences, and reached 7000 citations (often saying "discovered thanks to the OEIS").
Other ways to donate

Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A121351 a(n] = (3*n+1)*a(n-1) - a(n-2), starting a(0)=0, a(1)=1. 3
0, 1, 7, 69, 890, 14171, 268359, 5889727, 146974816, 4109405121, 127244583935, 4322206448669, 159794394016818, 6387453554224051, 274500708437617375, 12620645134576175199, 618137110885794967376, 32130509120926762128353 (list; graph; refs; listen; history; text; internal format)
OFFSET

0,3

LINKS

Table of n, a(n) for n=0..17.

FORMULA

3*a(n)= Pi*BesselJ_{4/3 + n}(2/3)* BesselY_{4/3}(2/3) - Pi*BesselJ_{4/3}(2/3) * BesselY_{4/3 + n}(2/3).

a(n) = sum {k = 0..floor((n-1)/2)} (-1)^k*3^(n-2*k-1)*(n-2*k-1)!*binomial(n-k-1,k)*binomial(n-k+1/3,k+4/3), cf. A058798. - Peter Bala, Aug 01 2013

a(n) ~ n! * BesselJ(4/3, 2/3) * 3^(n+1/3) * n^(1/3). - Vaclav Kotesovec, Jul 31 2014

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]

RecurrenceTable[{a[0]==0, a[1]==1, a[n]==(3n+1)*a[n-1]-a[n-2]}, a, {n, 20}]  (* Vaclav Kotesovec, Jul 31 2014 *)

CROSSREFS

Cf. A053984, A058798, A121323, A121353.

Sequence in context: A243668 A265033 A226270 * A302353 A059321 A217400

Adjacent sequences:  A121348 A121349 A121350 * A121352 A121353 A121354

KEYWORD

nonn,easy

AUTHOR

Roger L. Bagula and Bob Hanlon (hanlonr(AT)cox.net), Sep 05 2006

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 December 18 18:07 EST 2018. Contains 318243 sequences. (Running on oeis4.)