Re: plotting
- To: mathgroup at smc.vnet.net
- Subject: [mg124870] Re: plotting
- From: David Bailey <dave at removedbailey.co.uk>
- Date: Fri, 10 Feb 2012 05:53:15 -0500 (EST)
- Delivered-to: l-mathgroup@mail-archive0.wolfram.com
- References: <jh08a4$i5l$1@smc.vnet.net>
On 09/02/2012 10:48, David wrote:
> Dear All,
>
> I'm trying to plot the function r[x] as a function of x (x=Ï?/Ï?0) but when I put values for x couldn't get number. rewrite
> This is kindly ask you for your comments on how to plot r[x] as a function of x.
>
> Thank you in advance for your help. :)
>
> \[Gamma]m = 10^11;
> \[Gamma]e = 2.73*10^13;
> \[Omega]e = 1.37*10^16;
> \[Omega]m = 0.5 \[Omega]e;
> \[Omega]0 = 0.4092 \[Omega]e;
> \[Kappa]0 = 10^4;
> x[\[Omega]_] := \[Omega]/\[Omega]e
> y1 = 1.971*10^-3;
> y2 = 0.23*10^-5;
> a := 1 - 1/(x[\[Omega]]^2 + y1^2)
>
>
> b := \[Gamma]e/(\[Omega]e * x[\[Omega]] (x[\[Omega]]^2 + y1^2))
>
>
> c := 1 - (\[Omega]m^2/\[Omega]e^2)/(x[\[Omega]]^2 + y2^2)
>
> d := (\[Omega]m^2/\[Omega]e^2*\[Gamma]m)/(\[Omega]e*
> x[\[Omega]] (x[\[Omega]]^2 + y2^2))
>
> r[x_] := \[Sqrt]((\[Omega]e^2*x[\[Omega]]^2)/(
> 9*10^16 (a^2 - b^2 - 1.69) +
> 4 a^2*b^2 ) (1.3 (a^2 - b^2) - 1.69 (a*c - b*d))*(a^2 - b^2 -
> 1.69) + ((2.6 a*b - 1.69 (a*d + b*c)) 2 a*b))
>
> Plot[r[x]/\[Kappa]0, {x, 0, 0.6}]
>
You have defined a function x[\[Omega]_], but you didn't give a value
for \[Omega] itself (i.e. it doesn't have a global value) - so the
function r[x] isn't a number. The easiest way to see that is to evaluate
something like r[0] .
David Bailey
http://www.dbaileyconsultancy.co.uk