Comment #21 on issue 694 by smichr: Add tests from "Review of CAS  
mathematical capabilities", by Michael Wester
http://code.google.com/p/sympy/issues/detail?id=694

What is it that is not being worked on...the tests or this particular test?  
As far as
this solution "filling a screen", that is no longer the case. Solutions  
don't
simplify down to zero, however:

>>> l=[];s=solve(x**4+x**3+x**2+x+1,x)
>>> for si in s:
...     l.append(simplify(si))
...     
>>> l
[-1/4 - 5**(1/2)/4 + (-40 + 8*5**(1/2))**(1/2)/8, -1/4 + 5**(1/2)/4 + (-40 -
8*5**(1/2))**(1/2)/8, -1/4 - (-40 - 8*5**(1/2))**(1/2)/8 + 5**(1/2)/4, -1/4  
-
5**(1/2)/4 - (-40 + 8*5**(1/2))**(1/2)/8]
>>> for si in s:
...     print simplify((x**4+x**3+x**2+x+1).subs(x,si))
...     
5*(-40 + 8*5**(1/2))**(1/2)/64 - 5*5**(1/2)*(-40 + 8*5**(1/2))**(1/2)/64 -
5**(1/2)*(-40 + 8*5**(1/2))**(3/2)/512
5*(-40 - 8*5**(1/2))**(1/2)/64 + 5*5**(1/2)*(-40 - 8*5**(1/2))**(1/2)/64 +
5**(1/2)*(-40 - 8*5**(1/2))**(3/2)/512
-5*(-40 - 8*5**(1/2))**(1/2)/64 - 5*5**(1/2)*(-40 - 8*5**(1/2))**(1/2)/64 -
5**(1/2)*(-40 - 8*5**(1/2))**(3/2)/512
-5*(-40 + 8*5**(1/2))**(1/2)/64 + 5*5**(1/2)*(-40 + 8*5**(1/2))**(1/2)/64 +
5**(1/2)*(-40 + 8*5**(1/2))**(3/2)/512

But the solutions are numerically correct:

>>> for si in s:
...     print N(((x**4+x**3+x**2+x+1).subs(x,si)))
...     
.0e-126 + .0e-125*I
.0e-125 + .0e-126*I
.0e-125 - .0e-126*I
.0e-126 - .0e-125*I

--
You received this message because you are listed in the owner
or CC fields of this issue, or because you starred this issue.
You may adjust your issue notification preferences at:
http://code.google.com/hosting/settings

--

You received this message because you are subscribed to the Google Groups 
"sympy-issues" group.
To post to this group, send email to sympy-iss...@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.


Reply via email to