;$Id: axel_check.pro,v 1.1 2016/12/12 14:53:23 brandenb 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./1e11
r=grange(1.,4000.,2000)
debug=0
;
if debug then wait,.4
n=1.0e+11 & mu=23.61 & lam=2.412
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=2.0e+10 & mu=2.16 & lam=0.242
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=3.8e+8 & mu=-.69 & lam=0.017
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=1.1e+07 & mu=-0.92 & lam=0.003
f=s*n*r^mu*exp(-lam*r)*lam^(mu+1)/gamma(mu+1)
oplot,r,f,li=1