Method¶

refer to workplanJul29.org. Updated at Dec 3.

y(r)=y0(1+r^2/θ^2_c)^{−(3β−1)/2}

For β=1 case, we have: Y_{int}(r)=πθ^2_cy_0\ln[1+(r/θ_c)^2].

For β=2/3 case, we have: Y_{int}(r)=πθ^2_cy_0\times2[\sqrt{1+(r/θ_c)^2}−1].

For R500, it is normally 3 times of Rcore.

Source code¶

Result¶