Solving a Complex Diff Eq.
|
 05/30/12 08:03am
M = 100;
g = 1;
\[Alpha] = 1;
A = 10;
a = 1;(*initial plot temp*)
b = 10; (*final plot temp*)
c = .0001; (*initial temp*)
d = 100; (*initial \[Beta]*)
m = 1;
s = NDSolve[{\[Beta]'[T] == (3*\[Alpha]^2*M*T^(5/2))/(
Sqrt[2*g]*m^(
9/2))*((g/(2*pi)^(3/2)*(m/T)^(3/2)*E^(-m/T))^2 - \[Beta][
T]^2), \[Beta][c] == d}, \[Beta], {T, a, b}, Method -> "BDF"];
Plot[{Evaluate[\[Beta][T] /. s]}, {T, a, b}, PlotRange -> Automatic]
Here is a copy of my code, its generation errors and I'm not really sure how to solve the issue. I tried switching methods of solving but it doesn't seem to help. I have also attached a copy of the file with the errors up. Thanks for any help!
Attachment: Cosmology_Plots_3_Unfinished.nb, URL: , |
