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