On Tue, Aug 11, 2009 at 10:58 PM, Vinzent Steinberg<vinzent.steinb...@googlemail.com> wrote: > > # mpmath has a solver for polynomials, but we have to convert it to a > list of > # coefficients (please not that the results are not very accurate, you > can refine > # them using an iterative solver)
polyroots should give full accuracy unless there are repeated roots. > Indeed it's strange that the imaginary parts don't vanish, even if you > use higher precision for evaluating. This smells like a bug (assuming > the roots are really real). Fredrik, what do you think? > Please note that the second variant is somewhat inaccurate, you can > however use it as starting points for the first variant (see > documentation). It should work if you pass myEqConstants.n(50) instead of myEqConstants as input. I think it assumes that the input is only 15-digit accurate. Fredrik --~--~---------~--~----~------------~-------~--~----~ You received this message because you are subscribed to the Google Groups "sympy" group. To post to this group, send email to sympy@googlegroups.com To unsubscribe from this group, send email to sympy+unsubscr...@googlegroups.com For more options, visit this group at http://groups.google.com/group/sympy?hl=en -~----------~----~----~----~------~----~------~--~---