Re: FindRoots?
- To: mathgroup at smc.vnet.net
- Subject: [mg112137] Re: FindRoots?
- From: Andrzej Kozlowski <akozlowski at gmail.com>
- Date: Wed, 1 Sep 2010 06:28:52 -0400 (EDT)
On 30 Aug 2010, at 12:19, David Park wrote: > I don't know why Wolfram doesn't acquire the rights to this package and > incorporate it as part of regular Mathematica. Finding all the roots of a > real function on a 1-dimensional domain is the most common case of root > finding, and it is the one thing regular Mathematica is very poor at. Ted's > package is quite robust and returns all the roots in order. Mathematica's Reduce can do this (in the case of complex analytic functions, which covers most important practical cases) since version 7. I have posted examples of this several times, but somehow nobody seems to notice and we keep reading strange remarks like the above. Well, I don't want it to sound like I am trying to diminish Ted's achievement, but quite frankly, from the mathematical point of view, the methods used in the RootSearch package are rather primitive, relatively to the current state of knowledge in this field, of course. Not only does Reduce solve most transcendental equations much faster but it also does not require setting the values of any options by hand, and, what I think is most important, the results thus obtained are provably correct, which is not the case with Ted's package. In fact, much more general methods of solving transcendental equations exist (in general the functions need not be analytic and can involve systems of n-real equations in n-variab l! es with non-vanishing Jacobian). An example of such a method can be found in one of my demonstrations on the demonstrations site (the algorithm can easily be converted to a practical method for real life probl ems. In view of this, how could there be any justification for the suggestion that Wolfram incorporates what is actually an impressive but amateurish package into what is supposed to be a sophisticated mathematical program? Sometimes certain capabilities are missing from Mathematica simply because implementing "bleeding edge" algorithms in a way that is satisfactory for a program that has the kind of aspirations that Mathematica obviously does takes time. These sort of "helpful suggesting", particularly when coming from non-experts, have no chance whatever of being accepted and even can produce the impression of an ulterior motive being involved.