Re: [R] Solving 100th order equation

From: Spencer Graves <spencer.graves_at_pdf.com>
Date: Sat, 24 May 2008 08:32:59 -0700

Robert A LaBudde wrote:
> At 07:58 AM 5/24/2008, Duncan Murdoch wrote:

>> Shubha Vishwanath Karanth wrote:
>>> To apply uniroot I don't even know the interval values... Does 
>>> numerical methods help me? Or any other method?
>>>
>>
>> I forgot:  we also have polyroot().
>> Duncan Murdoch
>>> Thanks and Regards,
>>> Shubha
>>>
>>> -----Original Message-----
>>> From: Duncan Murdoch [mailto:murdoch_at_stats.uwo.ca] Sent: Saturday, 
>>> May 24, 2008 5:08 PM
>>> To: Shubha Vishwanath Karanth
>>> Subject: Re: [R] Solving 100th order equation
>>>
>>> Shubha Vishwanath Karanth wrote:
>>>
>>>> Hi R,
>>>>
>>>>
>>>>
>>>> I have a 100th order equation for which I need to solve the value 
>>>> for x. Is there a package to do this?
>>>> <snip>
>

> Finding all of the roots of a 100-th order polynomial on a computer is
> not an easy task. It would take move than 100 digits of precision to
> do so. A similar problem to finding all 100 eigenvalues of a 100x100
> ill-conditioned matrix.
>

> The suggestion to use graphics to find small intervals localizing the
> roots of interest first would make the most sense, at least for the
> real roots.
>

> It is highly unlikely that all 100 roots are needed. Usually the ones
> of interest are the one with the largest real part, or the ones inside
> the unit circle, or the real roots, etc.
And uniroot, polyroot or solve.polynom should provide marvelous starting values.
>

> Is this a homework problem?
>

> ================================================================
> Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: ral_at_lcfltd.com
> Least Cost Formulations, Ltd. URL: http://lcfltd.com/
> 824 Timberlake Drive Tel: 757-467-0954
> Virginia Beach, VA 23464-3239 Fax: 757-467-2947
>

> "Vere scire est per causas scire"
>

> ______________________________________________
> R-help_at_r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide
> http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.


R-help_at_r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. Received on Sat 24 May 2008 - 21:03:22 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Sat 24 May 2008 - 21:30:58 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive