The OEIS is supported by the many generous donors to the OEIS Foundation.

 Hints (Greetings from The On-Line Encyclopedia of Integer Sequences!)
 A269820 a(n) = 2*(n-1)*a(n-1) - a(n-2) + 2*(n-2)*a(n-3) with a(0)=a(1)=a(2)=1. 3
 1, 1, 1, 5, 33, 265, 2657, 31949, 447809, 7170193, 129126849, 2583427349, 56849678689, 1364649651865, 35486043526177, 993722660727005, 29814403968536193, 954131785638216481, 32442465432874088321, 1167988313532759193637, 44385462193094951063201 (list; graph; refs; listen; history; text; internal format)
 OFFSET 0,4 LINKS G. C. Greubel, Table of n, a(n) for n = 0..250 FORMULA E.g.f. A(x) satisfies 0 = A'(x)*(A'(x) + A'''(x)) - (A(x) + A''(x))*(A(x) + 2*A''(x)). 0 = a(n)*(+2*a(n+1) - a(n+2) + 4*a(n+3) - a(n+4)) + a(n+1)*(+a(n+1) + 2*a(n+3)) + a(n+2)*(-a(n+2) + 2*a(n+3) - a(n+4)) + a(n+3)*(+a(n+3)) for all n=>0. a(n) ~ c * 2^n * (n-1)!, where c = (BesselJ(0, 1/2) - BesselJ(1, 1/2))/2 = 0.348100674782969508922225... . - Vaclav Kotesovec, Mar 07 2016 EXAMPLE G.f. = 1 + x + x^2 + 5*x^3 + 33*x^4 + 265*x^5 + 2657*x^6 + 31949*x^7 + ... MATHEMATICA RecurrenceTable[{a[n] == 2*(n-1)*a[n-1] - a[n-2] + 2*(n-2)*a[n-3], a[0]==a[1]==a[2]==1}, a, {n, 0, 20}] (* Vaclav Kotesovec, Mar 07 2016 *) FullSimplify[Assuming[Element[x, Reals], CoefficientList[Series[Pi/4 *((BesselJ[1, 1/2] - BesselJ[0, 1/2]) * BesselY[0, -1/2 + x] + BesselJ[0, 1/2 - x] * (BesselY[0, -1/2] + BesselY[1, -1/2])), {x, 0, 20}], x] * Range[0, 20]!]] (* Vaclav Kotesovec, Mar 07 2016 *) PROG (PARI) {a(n) = if( n<3, n>=0, 2*(n-1)*a(n-1) - a(n-2) + 2*(n-2)*a(n-3))}; (PARI) {a(n) = my(A); if( n<3, n>=0, A = vector(n, i, 1); A[3]=5; for(i=4, n, A[i] = 2*(i-1)*A[i-1] - A[i-2] + 2*(i-2)*A[i-3]); A[n])}; (Magma) I:=[1, 1, 5]; [1] cat [n le 3 select I[n] else 2*(n-1)*Self(n-1) - Self(n-2) + 2*(n - 2)*Self(n-3): n in [1..50]]; // G. C. Greubel, Jul 30 2018 CROSSREFS Cf. A269450, A269509. Sequence in context: A367946 A087633 A135075 * A302075 A215671 A049377 Adjacent sequences: A269817 A269818 A269819 * A269821 A269822 A269823 KEYWORD nonn AUTHOR Michael Somos, Mar 07 2016 STATUS approved

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

Last modified August 3 05:44 EDT 2024. Contains 374875 sequences. (Running on oeis4.)