;$Id: axel_check2.pro,v 1.2 2016/12/16 03:25:39 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=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=2,thick=3
;