Updates:
Cc: matt...@gmail.com
Comment #5 on issue 1857 by smi...@gmail.com: evalf/simplify accuracy
http://code.google.com/p/sympy/issues/detail?id=1857
Here's another example:
z
(-27*12**(1/3)*sqrt(31)*I +
27*2**(2/3)*3**(1/3)*sqrt(31)*I)/(-2511*2**(2/3)*3**
(1/3) + (29*18**(1/3) + 9*2**(1/3)*3**(2/3)*sqrt(31)*I +
87*2**(1/3)*3**(1/6)*I)
**2)
pprint(z)
3 ____ ____ 2/3 3 ___ ____
- 27*\/ 12 *\/ 31 *I + 27*2 *\/ 3 *\/ 31 *I
---------------------------------------------------------------------------
2
2/3 3 ___ / 3 ____ 3 ___ 2/3 ____ 3 ___ 6 ___ \
- 2511*2 *\/ 3 + \29*\/ 18 + 9*\/ 2 *3 *\/ 31 *I + 87*\/ 2 *\/ 3 *I/
N(z)
-.0e+3 + .0e+3*I #..........those are such large zeros
numer(z).n()
.0e-134*I #.........so the numer is clearly a small zero
denom(z).n()
-69122.3265868904 + 39971.2403652584*I
factor(z).n()
-.0e+0 + .0e+0*I #............after factoring, the zeros are "smaller"
--
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.