def continued_fraction(r):
while True:
i = int(r)
yield i
if i == r:
break
r = 1 / (r - i)
Thilina Rathnayake於 2013年6月28日星期五UTC+8上午4時16分34秒寫道：
>
> Hi all,
>
> Is there a function in sympy which gives the continued fraction
> representation
> of quadratic irrationalities like, (sqrt(17) + 5 ) / 4?
>
> I found an implementation of a continued fraction function in Shor's
> algorithm:
> http://docs.sympy.org/dev/modules/physics/quantum/shor.html
>
> But it doesn't seem to give correct results when irrationalities are
> involved.
>
> >>> from sympy.physics.quantum.shor import continued_fraction
>> >>> continued_fraction(sqrt(17)+5, 4)
>> [2, 4]
>>
>
> Actual answer should be 2, [3, 1, 1] with the numbers inside [] repeating
> forever.
> Any help would be appreciated.
>
> Regards,
> Thilina.
> .
>
--
You received this message because you are subscribed to the Google Groups
"sympy" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to [email protected].
To post to this group, send email to [email protected].
Visit this group at http://groups.google.com/group/sympy.
For more options, visit https://groups.google.com/groups/opt_out.
