login
Continued fraction for sqrt(39).
2

%I #26 Dec 27 2023 00:40:16

%S 6,4,12,4,12,4,12,4,12,4,12,4,12,4,12,4,12,4,12,4,12,4,12,4,12,4,12,4,

%T 12,4,12,4,12,4,12,4,12,4,12,4,12,4,12,4,12,4,12,4,12,4,12,4,12,4,12,

%U 4,12,4,12,4,12,4,12,4,12,4

%N Continued fraction for sqrt(39).

%H Harry J. Smith, <a href="/A040032/b040032.txt">Table of n, a(n) for n = 0..20000</a>

%H G. Xiao, <a href="http://wims.unice.fr/~wims/en_tool~number~contfrac.en.html">Contfrac</a>

%H <a href="/index/Con#confC">Index entries for continued fractions for constants</a>

%H <a href="/index/Rec#order_02">Index entries for linear recurrences with constant coefficients</a>, signature (0,1).

%e 6.2449979983983982058468931... = 6 + 1/(4 + 1/(12 + 1/(4 + 1/(12 + ...)))). - _Harry J. Smith_, Jun 05 2009

%p Digits := 100: convert(evalf(sqrt(N)),confrac,90,'cvgts'):

%t ContinuedFraction[Sqrt[39],300] (* _Vladimir Joseph Stephan Orlovsky_, Mar 06 2011 *)

%t PadRight[{6},100,{12,4}] (* or *) Join[{6},LinearRecurrence[{0,1},{4,12},100]] (* _Harvey P. Dale_, Feb 09 2015 *)

%o (PARI) { allocatemem(932245000); default(realprecision, 35000); x=contfrac(sqrt(39)); for (n=0, 20000, write("b040032.txt", n, " ", x[n+1])); } \\ _Harry J. Smith_, Jun 05 2009

%Y Cf. A010493 Decimal expansion. - _Harry J. Smith_, Jun 05 2009

%K nonn,cofr,easy

%O 0,1

%A _N. J. A. Sloane_