MathGroup Archive 1998

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

Search the Archive

Re: calculation



Saeed Esmaily Rashid wrote:
> 
> Hello!
> 
> My name is Saeed and i'm studying physics. I have a problem which ihope
> someone can help me with it. I have an equation
> 
> .001422409738*R^(73/80)*(R^(1/2)-1)^(-13/40) == .04
> 
> i'm using the Solve function in Mathematica 3.0 to solve it for R, but
> it calculats endlessly and takes very long time. the question is that
> is there any way to optimize this equation or using another function 


Hello, 
      This is a numerical problem so use FindRoot, not Solve, 

FindRoot[.001422409738*R^(73/80)*(R^(1/2)-1)^(-13/40) == .04,{R,2}]

and you'll get 

{R->81.268}

Plotting the function also helps.  Try

Plot[.001422409738*R^(73/80)*(R^(1/2)-1)^(-13/40) - .04,{R,0,100}]

and you'll find another root close to the singularity

FindRoot[.001422409738*R^(73/80)*(R^(1/2)-1)^(-13/40) ==
.04,{R,1.00001}]

{R->1.00007}

You can also try using the option WorkingPrecision (if you know  that it
is justified).

John Jowett
-- 
----------------------------------------------------------------------
Electronic mail: John.Jowett@cern.ch World Wide Web: 
http://wwwslap.cern.ch/~jowett/ Telephone:       +41 22 76 76643     
<-----[you can leave voice mail] Fax:             +41 22 767 8480
Postal address:  SL Division,   CERN,   CH-1211 Geneva 23, Switzerland
Location:        Room 1-D18, Building 865, CERN Prevessin Site, France
Secretary:       +41 (0)22 76 75272 or 73076  <-----[changed 12/12/97]
----------------------------------------------------------------------



  • Prev by Date: Re: Graphics Output
  • Next by Date: Re: Numbered Equations?
  • Prev by thread: Re: calculation
  • Next by thread: Re: calculation