In[]:=
f[r_]:=1-
2M
r
+
2
Q
2
r
+
8
3
πP
2
r
​​mass=
3
2
Q
+3
2
r
+8Pπ
4
r
6r
/.rrh;​​T[r_]:=
-
2
Q
+
2
r
+8Pπ
4
r
4π
3
r
​​veff[r_]:=
f[r]
2
r
In[]:=
rp=(r/.Last[NSolve[D[veff[r],r]==0,r]])/.Mmass;
In[]:=
rs=
rp
f[rp]
;
In[]:=
rs/.{M1.2,rh0.4,Q0.5,P0.1}
Out[]=
1.9333
In[]:=
min=First@FindMinimum[T[rh]/.{Q0.5,P0.1}&&0<rh<1,{rh,0.6}]​​max=First@FindMaximum[T[rh]/.{Q0.5,P0.1}&&0<rh<1,{rh,0.6}]
Out[]=
0.160526
Out[]=
0.224907
In[]:=
colorf1=Blend[{{min,Yellow},{max,Red}},#1]&;​​bl=BarLegend[{colorf1,{min,max}},10,LegendMarkerSize300,LegendLabelStyle["legend",20,"Panel"]];
In[]:=
pl=ParametricPlot[{rsCos[θ],rsSin[θ]}/.{M1.2,Q0.5,P0.1},{rh,0,1},{θ,0,2π},PlotStylecolorf1];​​Legended[pl,Placed[bl,Right]]
Out[]=
legend