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.

A026639 a(n) = A026637(2*n, n-1).

This page as a plain text file.
%I A026639 #9 Jul 02 2024 02:15:13
%S A026639 1,5,20,74,278,1049,3980,15170,58052,222914,858512,3314960,12829070,
%T A026639 49748705,193259660,751954250,2929965020,11431262390,44651369720,
%U A026639 174597927740,683388447260,2677230376490,10496941482680,41188078562324
%N A026639 a(n) = A026637(2*n, n-1).
%H A026639 G. C. Greubel, <a href="/A026639/b026639.txt">Table of n, a(n) for n = 1..1000</a>
%F A026639 a(n) = ((7*n^2 - 4*n + 1)*a(n-1) + 2*n*(2*n-1)*a(n-2))/(2*(n^2-1)), with a(0) = 1, a(1) = 5, a(2) = 20. - _G. C. Greubel_, Jul 01 2024
%t A026639 a[n_]:= a[n]= If[n<4, (5*4^(n-1) -Boole[n==1])/4, ((7*n^2-4*n+1)*a[n- 1] +2*n*(2*n-1)*a[n-2])/(2*(n^2-1))];
%t A026639 Table[a[n], {n,40}] (* _G. C. Greubel_, Jul 01 2024 *)
%o A026639 (Magma)
%o A026639 [1] cat [n le 2 select 5*(3*n-2) else ((7*n^2+10*n+4)*Self(n-1) + 2*(2*n+1)*(n+1)*Self(n-2))/(2*n*(n+2)): n in [1..40]]; // _G. C. Greubel_, Jul 01 2024
%o A026639 (SageMath)
%o A026639 @CachedFunction
%o A026639 def a(n): # a = A026639
%o A026639     if n<4: return (5*4^(n-1) - 0^(n-1))/4
%o A026639     else: return ((7*n^2 - 4*n + 1)*a(n-1) + 2*n*(2*n-1)*a(n-2))/(2*(n^2-1))
%o A026639 [a(n) for n in range(1,41)] # _G. C. Greubel_, Jul 01 2024
%Y A026639 Cf. A026637, A026638, A026640, A026641, A026642, A026643, A026644.
%Y A026639 Cf. A026645, A026646, A026647, A026966, A026967, A026968, A026969.
%Y A026639 Cf. A026970.
%K A026639 nonn
%O A026639 1,2
%A A026639 _Clark Kimberling_