I get several results and want to check, which one is a real number. From trial and error, I can see that maxSol[2][x] is the only plausable solution.
Is there something like isReal()?
my attempt (don't laugh) ;)
for i in maxSol: print RealField(50)(i)
where I obviously get: TypeError: Unable to convert x (='.....
I think "a in RR" is what you want to use. Your first attempt could work too, using a try/except block.
For testing integrality, or rationality, "in ZZ" and "in QQ" work similarly.
And if you know the numbers you're checking are complex, you could also test to see whether their imaginary parts are zero.
It is not so easy: arccos(2) in RR is evaluated as true, while arccos(2).n() is NaN. The best we found with my student is a.n() != NaN and a.n().is_real() -- is this too complicated? See also http://bit.ly/LndAoI
again, this works perfect :)
for a in maxSol:
if a[x] in RR:
print RealField(50)(a[x])
Asked: 2012-05-04 12:39:24 +0100
Seen: 3,244 times
Last updated: May 04 '12
Copyright Sage, 2010. Some rights reserved under creative commons license. Content on this site is licensed under a Creative Commons Attribution Share Alike 3.0 license.