%I #8 Mar 09 2014 13:40:27
%S 1,4,91,240,840,989,5396,117719,312120,1089720,1284121,7003804,
%T 152799571,405130920,1414456320,1666787669,9090932396,198333725039,
%U 525859622640,1835963213040,2163489110641,11800023246004,257437022301451
%N Denominators a(n) of Pythagorean approximations b(n)/a(n) to 2/3.
%C See A195500 for a discussion and references.
%t r = 2/3; z = 25;
%t p[{f_, n_}] := (#1[[2]]/#1[[
%t 1]] &)[({2 #1[[1]] #1[[2]], #1[[1]]^2 - #1[[
%t 2]]^2} &)[({Numerator[#1], Denominator[#1]} &)[
%t Array[FromContinuedFraction[
%t ContinuedFraction[(#1 + Sqrt[1 + #1^2] &)[f], #1]] &, {n}]]]];
%t {a, b} = ({Denominator[#1], Numerator[#1]} &)[
%t p[{r, z}]] (* A195559, A195560 *)
%t Sqrt[a^2 + b^2] (* A195561 *)
%t (* _Peter J. C. Moses_, Sep 02 2011 *)
%Y Cf. A195500, A195560, A195561.
%K nonn
%O 1,2
%A _Clark Kimberling_, Sep 21 2011
|