cp's OEIS Frontend

A019427 Continued fraction for tan(1/4).

%I A019427 #22 Sep 08 2022 08:44:44
%S A019427 0,3,1,10,1,18,1,26,1,34,1,42,1,50,1,58,1,66,1,74,1,82,1,90,1,98,1,
%T A019427 106,1,114,1,122,1,130,1,138,1,146,1,154,1,162,1,170,1,178,1,186,1,
%U A019427 194,1,202,1,210,1,218,1,226,1,234,1,242,1,250,1,258,1,266,1,274,1,282,1,290,1,298,1
%N A019427 Continued fraction for tan(1/4).
%H A019427 Harry J. Smith, <a href="/A019427/b019427.txt">Table of n, a(n) for n = 0..20000</a>
%H A019427 G. Xiao, <a href="http://wims.unice.fr/~wims/en_tool~number~contfrac.en.html">Contfrac</a>
%H A019427 <a href="/index/Rec#order_04">Index entries for linear recurrences with constant coefficients</a>, signature (0,2,0,-1).
%F A019427 From _Colin Barker_, Sep 08 2013: (Start)
%F A019427 a(n) = (-1+3*(-1)^n-4*(-1+(-1)^n)*n)/2 for n>1.
%F A019427 a(n) = 2*a(n-2)-a(n-4) for n>5.
%F A019427 G.f.: x*(x^4-x^3+4*x^2+x+3) / ((x-1)^2*(x+1)^2). (End)
%e A019427 0.25534192122103626650448223... = 0 + 1/(3 + 1/(1 + 1/(10 + 1/(1 + ...)))). - _Harry J. Smith_, Jun 13 2009
%t A019427  Join[{0, 3}, LinearRecurrence[{0, 2, 0, -1}, {1, 10, 1, 18}, 100]] (* _Vincenzo Librandi_, Jan 03 2016 *)
%o A019427 (PARI) { allocatemem(932245000); default(realprecision, 91000); x=contfrac(tan(1/4)); for (n=0, 20000, write("b019427.txt", n, " ", x[n+1])); } \\ _Harry J. Smith_, Jun 13 2009
%o A019427 (PARI) Vec(x*(x^4-x^3+4*x^2+x+3)/((x-1)^2*(x+1)^2) + O(x^100)) \\ _Colin Barker_, Sep 08 2013
%o A019427 (Magma) [0,3] cat [(-1+3*(-1)^n-4*(-1+(-1)^n)*n)/2: n in [2..80]]; // _Vincenzo Librandi_, Jan 03 2016
%Y A019427 Cf. A161013 (decimal expansion). - _Harry J. Smith_, Jun 13 2009
%K A019427 nonn,cofr,easy
%O A019427 0,2
%A A019427 _David W. Wilson_