[R] Solving 100th order equation

Duncan Murdoch murdoch at stats.uwo.ca
Sat May 24 13:56:18 CEST 2008


Shubha Vishwanath Karanth wrote:
> To apply uniroot I don't even know the interval values... Does numerical methods help me? Or any other method?
>   
How about plotting the function?  If you write it as f(x) = x^100 - ... 
- 6*x - 4000, then clearly f(0) is negative, and large enough x will be 
positive.  Just find
out what "large enough" means (it depends on the ...) and use that to 
start uniroot().

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?
>>
>>  
>>
>> For example my equation is:
>>
>>  
>>
>> (x^100 )- (2*x^99) +(10*x^50)+.............. +(6*x ) = 4000
>>
>>  
>>
>> I have only one unknown value and that is x. How do I solve for this?
>>
>>   
>>     
> uniroot() will find one root.  If you want all of them, I don't know 
> what is available.
>
> Duncan Murdoch
>   
>>  
>>
>>  
>>
>>  
>>
>> BR, Shubha
>>
>> Shubha Karanth | Amba Research
>>
>> Ph +91 80 3980 8031 | Mob +91 94 4886 4510 
>>
>> Bangalore * Colombo * London * New York * San José * Singapore * www.ambaresearch.com
>>
>>  
>>
>> This e-mail may contain confidential and/or privileged i...{{dropped:13}}
>>
>>   
>> ------------------------------------------------------------------------
>>
>> ______________________________________________
>> 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.
>>   
>>     
>
> This e-mail may contain confidential and/or privileged i...{{dropped:10}}
>
> ______________________________________________
> 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.
>



More information about the R-help mailing list