| Author |
Comment/Response |
Bill Simpson
|
 06/01/12 12:10pm
One method, trying to make the smallest changes to what you have and accomplish what I think you are asking for. There should be a cleaner simpler way of doing this, but perhaps this will do.
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;
zed = {};
For[M = 0., M <= .8, M += .1,
bb = (g/(2*\[Pi])^(3/2)*(m/T)^(3/2)*\[ExponentialE]^(-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"];
zed = Append[zed, Plot[{(Evaluate[\[Beta][T] /. s])/bb}, {T, a, b}, PlotRange->Automatic]]
];
Show[zed]
URL: , |
|
