       Re: Plotting a series of Roots

• To: mathgroup at smc.vnet.net
• Subject: [mg128800] Re: Plotting a series of Roots
• From: William Duhe <wjduhe at loyno.edu>
• Date: Tue, 27 Nov 2012 03:31:49 -0500 (EST)
• Delivered-to: l-mathgroup@mail-archive0.wolfram.com
• Delivered-to: l-mathgroup@wolfram.com
• Delivered-to: mathgroup-newout@smc.vnet.net
• Delivered-to: mathgroup-newsend@smc.vnet.net

```So here I have an algorithm that gives me how my roots change as a function of lambda or t1[lambda,1]:

s1[lambda_, t_] =
x[t] /. DSolve[{x'[t] == lambda, x == 0}, x[t], t][];

lambda t;

t1[lambda_, value_] = t /. Solve[s1[lambda, t] == value, t][];

value/lambda;

s1[lambda, t1[lambda, value]];

value;

Plot[t1[lambda, 1], {lambda, .1, 1}, PlotRange -> Automatic]

What I would like to do is in a similar fashion solve the differential equation bellow and plot how a particular root of that equation changes with \[Alpha]. It is important to note here that I have to use NDSolve rather than DSolve for this type of equation.What I would like is analogous to the above code modified to be t1[\[Alpha],1]with the modified equations.

M = 10000;
g = 1;
\[Alpha] = 1;
A = 10;
a = .01;(*initial plot temp*)
b = 1.2; (*final plot temp*)
c = 0.00000001; (*initial temp*)
d = 0; (*initial \[Beta]*)
m = 1;
bb = (g/(2*\[Pi])^(3/2)*(m/T)^(3/2)*E^(-m/T));

s = NDSolve[{\[Beta]'[T] == (3*\[Alpha]^2*M*T^(5/2))/(
Sqrt[2*g]*m^(9/2))*(bb^2 - \[Beta][T]^2), \[Beta][c] ==
d}, \[Beta], {T, a, b}, Method -> "BDF"];

```

• Prev by Date: Re: Numerical expression
• Next by Date: Re: inserting an equation and plotting it
• Previous by thread: Re: Plotting a series of Roots
• Next by thread: Re: Plotting a series of Roots