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 

      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 


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) ==


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

John Jowett
Electronic mail: World Wide Web: 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