OFFSET
1,1
COMMENTS
See A195500 for a discussion and references.
MATHEMATICA
r = 3/2; z = 21;
p[{f_, n_}] := (#1[[2]]/#1[[
1]] &)[({2 #1[[1]] #1[[2]], #1[[1]]^2 - #1[[
2]]^2} &)[({Numerator[#1], Denominator[#1]} &)[
Array[FromContinuedFraction[
ContinuedFraction[(#1 + Sqrt[1 + #1^2] &)[f], #1]] &, {n}]]]];
{a, b} = ({Denominator[#1], Numerator[#1]} &)[
Sqrt[a^2 + b^2] (* A195552 *)
(* Peter J. C. Moses, Sep 02 2011 *)
CROSSREFS
KEYWORD
nonn,frac
AUTHOR
Clark Kimberling, Sep 21 2011
STATUS
approved