| Author |
Comment/Response |
Bill Simpson
|
 06/01/12 12:26pm
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;
Show[
Table[
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"];
Plot[{(Evaluate[\[Beta][T]/.s])/bb},{T,a,b},PlotRange->Automatic],
{M,0.,.8,.1}
]
]
URL: , |
|
