;$Id: axel_check.pro,v 1.2 2016/12/16 03:25:38 xiangyu Exp $
;
;  This script is included as "@axel_check" after the first plot in
;  moments_shima_Nx32_a24.pro (current copy moments_shima_Nx32_a24_axel.pro)
;
;  mv idl.ps ~/tex/xiangyu/SwarmSmolu_numerics/gamma_distribution/fig/axel_check.ps
;
s=10./1e10
r=grange(1.,4000.,2000)
debug=0
;
if debug then wait,.4
n=1.0e+10 & mu=23.42 & lam=2.394
f=s*n*r^mu*exp(-lam*r)*lam^(mu+1)/gamma(mu+1)
oplot,r,f,li=2,thick=3
;
if debug then wait,.4
n=2.0e+9 & mu=3.50 & lam=0.300
f=s*n*r^mu*exp(-lam*r)*lam^(mu+1)/gamma(mu+1)
oplot,r,f,li=2,thick=3
;
;if debug then wait,.4
;n=3.4e+8 & mu=1.25 & lam=0.098
;f=s*n*r^mu*exp(-lam*r)*lam^(mu+1)/gamma(mu+1)
;oplot,r,f,li=1
;
if debug then wait,.4
n=6.7e+07 & mu=-0.23 & lam=0.034
f=s*n*r^mu*exp(-lam*r)*lam^(mu+1)/gamma(mu+1)
oplot,r,f,li=2,thick=3