%I #29 Nov 12 2023 05:58:32
%S 3,1,6,1,6,1,6,1,6,1,6,1,6,1,6,1,6,1,6,1,6,1,6,1,6,1,6,1,6,1,6,1,6,1,
%T 6,1,6,1,6,1,6,1,6,1,6,1,6,1,6,1,6,1,6,1,6,1,6,1,6,1,6,1,6,1,6,1,6,1,
%U 6,1,6,1,6,1,6,1,6,1,6,1,6
%N Continued fraction for sqrt(15).
%H Harry J. Smith, <a href="/A040011/b040011.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).
%F From _Amiram Eldar_, Nov 12 2023: (Start)
%F Multiplicative with a(2^e) = 6, and a(p^e) = 1 for an odd prime p.
%F Dirichlet g.f.: zeta(s) * (1 + 5/2^s). (End)
%e 3.872983346207416885179265399... = 3 + 1/(1 + 1/(6 + 1/(1 + 1/(6 + ...)))). - _Harry J. Smith_, Jun 03 2009
%p Digits := 100: convert(evalf(sqrt(N)),confrac,90,'cvgts'):
%t ContinuedFraction[Sqrt[15],300] (* _Vladimir Joseph Stephan Orlovsky_, Mar 05 2011 *)
%t PadRight[{3},120,{6,1}] (* _Harvey P. Dale_, Apr 14 2020 *)
%o (PARI) { allocatemem(932245000); default(realprecision, 19000); x=contfrac(sqrt(15)); for (n=0, 20000, write("b040011.txt", n, " ", x[n+1])); } \\ _Harry J. Smith_, Jun 03 2009
%Y Cf. A010472 (decimal expansion).
%K nonn,cofr,easy,mult
%O 0,1
%A _N. J. A. Sloane_