Below is the outcome I get when I use myEqConstants.n(100): third variant: [-0.192261149873352 + .0e-21⋅ⅈ, -29.4436689512432 + .0e-19⋅ⅈ, 27.1235615040264 - .0e-19⋅ⅈ, 0, 0.192368597090154 - .0e-21 ⋅ⅈ]
the imaginary parts are zeros, but they are still in the solution. Should they not disappear? Thanks Regars, anartz On Wed, Aug 12, 2009 at 1:24 PM, Fredrik Johansson < fredrik.johans...@gmail.com> wrote: > > 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 -~----------~----~----~----~------~----~------~--~---
