On Sat, Mar 20, 2010 at 2:47 PM, Ondrej Certik <[email protected]> wrote: > Thanks Aaron! > > It may be, that the new polynomials can't handle floating point > numbers (for some algorithmic reasons), in which case imho we should > provide the user with the following options: > > a) convert floating point numbers to Rational and run solve(), and at > the end, convert all roots to floating point numbers again, it has the > advantage of finding all the roots (I mean as long as they can be > found analytically) > b) use some numerical algorithms > > Maybe for polynomial equations, there are numerical algorithms that > return (find) all roots automatically? nsolve() currently needs an > initial guess.
Wolfram alpha seems to be able to find all 4 roots numerically (and automatically): http://www.wolframalpha.com/input/?i=solve(-12*(2*sqrt(3)+%2B+9)+*+12++%2B+(78+%2B+24.3*sqrt(3)+)*x**2+-+3*sqrt(2)+*+x**3+-+x**4,+x) so we should do it too. Ondrej -- You received this message because you are subscribed to the Google Groups "sympy" group. To post to this group, send email to [email protected]. To unsubscribe from this group, send email to [email protected]. For more options, visit this group at http://groups.google.com/group/sympy?hl=en.
