cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A026638 a(n) = A026637(2*n, n).

This page as a plain text file.
%I A026638 #20 Jul 03 2024 01:19:24
%S A026638 1,2,8,26,92,332,1220,4538,17036,64412,244928,935684,3588392,13806704,
%T A026638 53271548,206040506,798600332,3101109164,12062148368,46986821516,
%U A026638 183276382472,715748620424,2798274135368,10951009023716,42895901012792,168167959150232,659793819847040
%N A026638 a(n) = A026637(2*n, n).
%H A026638 Vincenzo Librandi, <a href="/A026638/b026638.txt">Table of n, a(n) for n = 0..1000</a>
%F A026638 From _Vaclav Kotesovec_, Oct 21 2012: (Start)
%F A026638 G.f.: (3 - (x+1)*sqrt(1-4*x))/((x+2)*sqrt(1-4*x)).
%F A026638 Recurrence: 2*n*a(n) = (7*n-4)*a(n-1) + 2*(2*n-1)*a(n-2).
%F A026638 a(n) ~ 2^(2*n+2)/(3*sqrt(Pi*n)) (End)
%t A026638 CoefficientList[Series[1/(2+x)+3/((2+x)*Sqrt[1-4*x])-1,{x,0,20}],x] (* _Vaclav Kotesovec_, Oct 21 2012 *)
%o A026638 (PARI) my(x='x+O('x^66)); Vec( 1/(2+x)+3/((2+x)*sqrt(1-4*x))-1 ) \\ _Joerg Arndt_, May 04 2013
%o A026638 (Magma) [1] cat [n le 2 select 2^(2*n-1) else ((7*n-4)*Self(n-1) + 2*(2*n-1)*Self(n-2))/(2*n): n in [1..40]]; // _G. C. Greubel_, Jul 01 2024
%o A026638 (SageMath)
%o A026638 @CachedFunction
%o A026638 def a(n): # a = A026638
%o A026638     if n<3: return 2^(n*(n+1)/2)
%o A026638     else: return ((7*n-4)*a(n-1) + 2*(2*n-1)*a(n-2))/(2*n)
%o A026638 [a(n) for n in range(41)] # _G. C. Greubel_, Jul 01 2024
%Y A026638 Cf. A026637, A026639, A026640, A026642, A026643, A026644, A026645.
%Y A026638 Cf. A026646, A026647, A026966, A026967, A026968, A026969, A026970.
%K A026638 nonn
%O A026638 0,2
%A A026638 _Clark Kimberling_