In[]:=
ClearAll[M,m,Mi,fd,l,F,a,Mcr,Υ,k,gd,g,f,z,t,Gy,Z,ul,unl]
In[]:=
G=6.67*
-11
10
;​​c=3*
8
10
;​​k=1.38*
-23
10
;​​hs=6.63*
-34
10
/(2*Pi);​​σ=5.67*
-8
10
;​​σT=66.52*
-30
10
;​​mp=1.67*
-27
10
;​​Gy=3.16*
16
10
;​​tu=13.77*Gy;
In[]:=
fd=(hs*
4
c
)/
2
G
*f;​​gd=(hs*
4
c
)/
2
G
*g;​​Z=
2/3
(14*Gy/t)
;
In[]:=
λ=1.12;​​mH=1.67*
-27
10
;​​cs=6*
3
10
;​​n=0.25*
0
10
*
3
Z
;​​ρ=n*mH;​​unl=.25
1.5
(620*
3
10
)
*
0
10
*Pi*λ*mH*
2
G
*
-3/2
cs
;
In[]:=
s1=NDSolveM'[t]-1*fd
2
M[t]
+1*unl*
4.3
(14*Gy/t)
*1+
1/2
(1-
2
a
)
*
2
M[t]
/.{a0.1,f4.684*
-6
10
},M[
-20
10
*Gy]1*
11
10
,M[t],{t,1,2*tu}
NDSolve
:At t == 0.000316​, step size is effectively zero; singularity or stiff system suspected.
NDSolve
:At t == 0.000316​, step size is effectively zero; singularity or stiff system suspected.
Out[]=
M[t]InterpolatingFunction
Domain: {{1.,0.000316}}
Output: scalar
[t]
In[]:=
p4=LogLogPlot[Evaluate[M[t]/.s1],{t,1,2*tu}]
InterpolatingFunction
:Input value {1.00084} lies outside the range of data in the interpolating function. Extrapolation will be used.
Out[]=