Author 
Comment/Response 
William Duhe

05/31/12 11:06am
I want to vary parameters in this differential equation.
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"];
Plot[{(Evaluate[\[Beta][T] /. s])/bb}, {T, a, b},
PlotRange > Automatic]
Notice that what I am actually evaluating is the ratio of the differential equation to bb. I would like to track the changes in this ratio as I vary M from 110000. How would I do this? I will attach a copy of the work I have tried so far using the Table syntax.
Attachment: Cosmology_Plots_3.nb, URL: , 
