MathGroup Archive 2003

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

Search the Archive

Re: bug???

  • To: mathgroup at smc.vnet.net
  • Subject: [mg45238] Re: bug???
  • From: drbob at bigfoot.com (Bobby R. Treat)
  • Date: Tue, 23 Dec 2003 05:13:31 -0500 (EST)
  • References: <bs3nhd$io4$1@smc.vnet.net>
  • Sender: owner-wri-mathgroup at wolfram.com

NSolve tries (I believe) to solve the problem ALGEBRAICALLY. There's
no closed form solution for most polynomials of order greater than 5,
so... good luck with order 802. Try this instead:

Timing[FindRoot[
   x^802 + 4*x^403 + 2*x^2 - 
     3 == 0, {x, 1}, 
   WorkingPrecision -> 40]]
{0.009999999999999995*Second, 
  {x -> 0.9964552050290507443\
02715516924973612445054031242\
178101144`39.99999999999999}}

Fast enough for you?

Bobby


"ANDY" <wangzhengyao at hotmail.com> wrote in message news:<bs3nhd$io4$1 at smc.vnet.net>...
> hello,
> 
> please running the next Expression:
> 
> Timing[NSolve[x^802 + 4x^403 + 2x^2 - 3 == 0, x, 4]]
> 
> in Mathematica4.1 and Mathematica 5.
> 
> Why???


  • Prev by Date: Re: Heart Plot
  • Next by Date: polynomialsorting.html + tensor
  • Previous by thread: Re: bug???
  • Next by thread: Re: bug???