Mathematica 9 is now available
Services & Resources / Wolfram Forums / MathGroup Archive
-----

MathGroup Archive 2012

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

Search the Archive

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




  • Prev by Date: Re: ListInterpolate and missing values
  • Next by Date: Re: ListInterpolate and missing values
  • Previous by thread: Re: plotting
  • Next by thread: Re: plotting