Author 
Comment/Response 
William Duhe

06/06/12 11:24am
g = 1;
\[Alpha] = 1;
A = 10;
a = .2;(*initial plot temp*)
b = 2.5;(*final plot temp*)
c = 0.00000001;(*initial temp*)
d = 0;(*initial \[Beta]*)
m = 1;
M = 1;
Show[Plot[(
E^(m/T)*(T/m)^2)/(((2 \[Pi])^(3/2)*
Sqrt[g]*\[Alpha]^2)/((6*g*\[Alpha]^2)*(m/M))), {T, a, b},
PlotStyle > {Red}],
Table[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], {M, 1000, 1000, 100}]]
I want to plot both of these functions shown in the 2D plot above together in a 3D plot. I want to vary M simultaneously in both functions from 110,000. I am not sure the best way to go about this. Any thoughts? I want the two grids laid down on top of one another in order to analyze intersections in the data. Thanks a lot for any help! Attached is a copy of the file along with some other 3D plots I produced for single functions.
Attachment: Cosmology_Plots_3_4.nb, URL: , 
