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¶