MathGroup Archive 2002

[Date Index] [Thread Index] [Author Index]

Search the Archive

Bug in NSolve?

  • To: mathgroup at smc.vnet.net
  • Subject: [mg33788] Bug in NSolve?
  • From: rodgerroSPAMNOT.siteconnect.com at library2.airnews.net
  • Date: Tue, 16 Apr 2002 03:50:44 -0400 (EDT)
  • Sender: owner-wri-mathgroup at wolfram.com

The Mathematica book says in section 2.9.5:

"NSolve will always give you the complete set of numerical solutions to any
polynomial equation in one variable."

However,  Mathematica version 4.1.2  on a Pentium III gives this result:

NSolve[4877361379 x^2 -9754525226 x + 4877163849==0,x,20]

{{ x---> 0.9999797501 }}

On the other hand, 

NSolve[4877361379 x^2 -9754525226 x + 4877163849==0,x,10]

gives:

{{ x ---> 0.99998 } , { x---> 0.99998 } }

Machine precision works, but arbitrary precision doesn't?

At least now we get two roots; the correct solution is given by:

N[Solve[4877361379 x^2 -9754525226 x + 4877163849==0,x ] ,10]

which is:

{{ x ---> 0.99998 - 2.89955 x 10^-10 i },{ x---> 0.99998 + 2.89955 x 10^-10i }}

Does anybody know what's going on here?




  • Prev by Date: Is it possible?
  • Next by Date: RE: Subscript[x, y] and the Symbol x
  • Previous by thread: Re: Is it possible?
  • Next by thread: Re: Bug in NSolve?