MathGroup Archive 2004

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

Search the Archive

Re: Plotting an implicit function

  • To: mathgroup at smc.vnet.net
  • Subject: [mg50598] Re: [mg50559] Plotting an implicit function
  • From: Bob Hanlon <hanlonr at cox.net>
  • Date: Sat, 11 Sep 2004 06:45:25 -0400 (EDT)
  • Reply-to: hanlonr at cox.net
  • Sender: owner-wri-mathgroup at wolfram.com

r[M_]:=r /. FindRoot[
        1+(-1-0.03*M^2+1.2*M)*r-
            (0.05*E^(r*M)/r-0.05/r+E^(r*M)*(1-0.05*M)) == 0,
        {r,0.01}];

data=Table[{M,r[M]}, {M,5,20,.1}];

ListPlot[data,
    Frame->True,
    Axes->False, 
    PlotJoined->True, 
    PlotStyle->RGBColor[1, 0, 0], 
    PlotRange->All];

data2=Table[{M,r[M]}, {M,5.855,5.86,.0002}];

ListPlot[data2,
    Frame->True,
    Axes->False, 
    PlotJoined->True, 
    PlotStyle->RGBColor[1, 0, 0], 
    PlotRange->All];


Bob Hanlon

> 
> From: david.friskin at upe.ac.za (David Friskin)
To: mathgroup at smc.vnet.net
> Date: 2004/09/10 Fri AM 04:05:59 EDT
> To: mathgroup at smc.vnet.net
> Subject: [mg50598] [mg50559] Plotting an implicit function
> 
> I have the following implicit function:
> f[r, M]=1 + (-1 - 0.03*M^2 + 
>     1.2*M)*r - (0.05*E^(r*M)/r - 0.05/r + E^(r*M)*(1 - 0.05*M))==0
> 
> where r is a function of M i.e. r(M).
> 
> I would like to plot r(M) vs M, for {M,0,20}.
> 
> I tried using
>  ImplicitPlot[{f[r,M] == 0},{M,0,20}], but Solve is having problems
> finding the roots to my function.
> 
> Could someone help me with the code to do this, perhaps using FindRoot
> somehow?
> 
> Thanks
> David
> 
> P.S. r > 0
> 
> 


  • Prev by Date: Re: ColorFunctions again (making z=0 be different from z=1)
  • Next by Date: RE: Plotting an implicit function
  • Previous by thread: Re: Plotting an implicit function
  • Next by thread: RE: Plotting an implicit function