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.

A126963 Numerators of sequence defined by f(0)=1, f(1)=5/4; f(n) = ( (6*n-1)*f(n-1) - (2*n-1)*f(n-2) )/(4n).

This page as a plain text file.
%I A126963 #24 Sep 08 2022 08:45:29
%S A126963 1,5,43,177,2867,11531,92479,370345,11857475,47442055,379582629,
%T A126963 1518418695,24295375159,97182800711,777467420263,3109879375897,
%U A126963 199032580597603,796130905791967,6369049515119561,25476202478636219,407619274119811709,1630477163761481141
%N A126963 Numerators of sequence defined by f(0)=1, f(1)=5/4; f(n) = ( (6*n-1)*f(n-1) - (2*n-1)*f(n-2) )/(4n).
%H A126963 G. C. Greubel, <a href="/A126963/b126963.txt">Table of n, a(n) for n = 0..500</a>
%H A126963 D. Doster, <a href="http://www.jstor.org/stable/2691073">Problem 1318, Three Term Recurrence</a>, Math. Magazine, 63 (1990), 127-128.
%F A126963 f(n) = Sum_{k=0..n} binomial(-1/2,k)*(-1/2)^k.
%F A126963 f(n) -> sqrt(2) as n -> oo.
%F A126963 G.f.: (sqrt(-x)*arccsc(1-x)/sqrt(2)-(Pi*i*sqrt(x))/sqrt(2)^3)/x. - _Vladimir Kruchinin_, Oct 10 2012
%F A126963 a(n) = numerator( Sum_{k=0..n} binomial(2*k, k)/8^k ). - _G. C. Greubel_, Jan 29 2020
%p A126963 seq( numer( add(binomial(2*k, k)/8^k, k=0..n) ), n=0..25); # _G. C. Greubel_, Jan 29 2020
%t A126963 a[n_] := Sqrt[2](1-(Gamma[1/2+n] Hypergeometric2F1[n,1/2+n,1+n,-1])/(Sqrt[Pi] Gamma[1+n])); Table[Numerator[FullSimplify[a[n]]], {n,20}] (* _Gerry Martens_, Aug 09 2015 *)
%t A126963 f[n_]:= If[n==0, 1, If[n==1, 5/4, ((6*n-1)*f[n-1]-(2*n-1)*f[n-2])/(4*n)]];
%t A126963 Table[Numerator[f[n]], {n, 0, 25}] (* _G. C. Greubel_, Jan 29 2020 *)
%o A126963 (PARI) A126963(n)=numerator(sum(k=0,n,binomial(-1/2,k)/(-2)^k)) \\ f(n)=if(n>1,((6*n-1)*f(n-1)-(2*n-1)*f(n-2))/(4*n),(5/4)^n) yields the same results. - _M. F. Hasler_, Aug 11 2015
%o A126963 (Magma) [Numerator( &+[Binomial(2*k, k)/8^k: k in [0..n]] ): n in [0..25]]; // _G. C. Greubel_, Jan 29 2020
%o A126963 (Sage) [numerator( sum(binomial(2*k, k)/8^k for k in (0..n)) ) for n in (0..25)] # _G. C. Greubel_, Jan 29 2020
%o A126963 (GAP) List([0..25], n-> NumeratorRat( Sum([0..n], k-> Binomial(2*k,k)/8^k) ));  # _G. C. Greubel_, Jan 29 2020
%Y A126963 Denominators are in A088802.
%K A126963 nonn,frac
%O A126963 0,2
%A A126963 _N. J. A. Sloane_, Mar 20 2007