[R] How to solve difficult equations?

Robert A LaBudde ral at lcfltd.com
Wed Apr 25 16:47:50 CEST 2007


At 03:15 AM 4/25/2007, francogrex wrote:

>This below is not solvable with uniroot to find "a":
>fn=function(a){
>b=(0.7/a)-a
>(1/(a+b+1))-0.0025
>}
>uniroot(fn,c(-500,500))  gives
>"Error in uniroot(fn, c(-500, 500)) : f() values at end points not of
>opposite sign"
>
>I read R-help posts and someone wrote a function:
>http://finzi.psych.upenn.edu/R/Rhelp02a/archive/92407.html
>but it is not very precise. Is there any '"standard" function in R that can
>solve this? thanks.

Actually, if you're solving fn(a)==0, then some trivial algebra leads 
to a linear equation with a=0.001754.

Why are you trying to solve this numerically? Is it a single instance 
of a larger, more general 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"



More information about the R-help mailing list