Just for the record, I checked these with the recent Sage beta (8.1.beta9),
and everything works, no errors.
On Friday, February 18, 2011 at 4:07:38 PM UTC, D. S. McNeil wrote:
>
> A somewhat simpler test case, which I think preserves the qualitative
> issue:
>
> sage: from sage.rings.polynomial
A somewhat simpler test case, which I think preserves the qualitative issue:
sage: from sage.rings.polynomial.real_roots import real_roots
sage:
sage: x = polygen(QQ)
sage: f = 2503841067*x^13 - 15465014877*x^12 + 37514382885*x^11 -
44333754994*x^10 + 24138665092*x^9 - 2059014842*x^8 - 3197810701*
Just for the record. The problem seems to be related to RIF. For the
inexact ring RR, it works:
len(e.roots(ring=RR))
13
len(e.real_roots())
13
numeric approximations of the two missing roots are:
0.953956769342757, 0.957223630414975
This pair of roots is exactly the pair of most close roots am
I was not sufficiently careful in posting my polynomial e(x) and
apparently some bad line breaks and spaces were introduced. This reply
is to post a properly-wrapped copy of the polynomial.
Zach Teitler
594813163982683595176185452130468519129644968918377817372539082858201843
58576518918331892948