Status: New
Owner: ----
Labels: Type-Defect Priority-Medium

New issue 1439 by smichr: polynomial solution missing real roots
http://code.google.com/p/sympy/issues/detail?id=1439

The equation -1936 - 5056*x - 7592*x**2 + 2704*x**3 - 49*x**4 has two real
roots at about 3.6 and 52 but only imaginary roots (and approx 0 at that)
are being found.

>>> solve(_,x)
.
.etc...
.
3176523*256871224631296524288**(1/2))**(5/12)*(1086740688896/3176523 +
256871224631296524288**(1/2)/117649)**(1/6))]
>>> [tmp.evalf() for tmp in _]
[.0e-16 + .0e-16*I, .0e-16 + .0e-16*I, -.0e-18 - .0e-18*I, .0e-17  
+ .0e-17*I]

Is the stiffness of this eqn giving the solver problems?

/c

--
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-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
-~----------~----~----~----~------~----~------~--~---

Reply via email to