Comment #3 on issue 2813 by smi...@gmail.com: solve(x**2 - x - 0.1) doesn't give both solutions
http://code.google.com/p/sympy/issues/detail?id=2813
In pull 656 this gives [1/2 + sqrt(35)/10, -sqrt(35)/10 + 1/2] It would probably be better to return Floats, but for now, unless rational is False, any floats will be made into Rationals. If the rational flag is used then the same problem occurs but it appears to be related to the precision requested in n(): [Dbg]>>> for i in range(1,21): [Dbg]... print i, val.evalf(36).n(i,chop=True) [Dbg]... 1 0 2 0 3 0 4 0 5 0 6 0 7 0 8 0 9 0 10 0 11 0 12 0 13 0 14 0 15 0 16 0 17 -5.2041704279304213e-18 18 -5.20417042793042128e-18 19 -5.204170427930421283e-18 20 -5.2041704279304212832e-18 b/c of this an immediate return should be deferred to see if other simplification can prove the solution or better numerical zero detections implemented. -- You received this message because you are subscribed to the Google Groups "sympy-issues" group. To post to this group, send email to sympy-issues@googlegroups.com. To unsubscribe from this group, send email to sympy-issues+unsubscr...@googlegroups.com. For more options, visit this group at http://groups.google.com/group/sympy-issues?hl=en.