from numpy import *
from matplotlib.pyplot import *
from mpl_toolkits.axes_grid.axislines import SubplotZero
 
pdf = lambda c,mu,x: sqrt(c/(2.*pi)) * exp(-0.5*c / (x-mu)) / (x-mu)**1.5

fig = figure(figsize=(6,5))
ax = SubplotZero(fig,111)
fig.add_subplot(ax)
ax.grid(True)
ax.minorticks_on()
 
x = logspace(-1.6,1.6)
ax.loglog(x,pdf(.5,0.,x),c='r',lw=1.5,label=r'$c=1/2$')#'$c={}^1\!/_2$')
ax.loglog(x,pdf(1.,0.,x),c='k',lw=1.5,label=r'$c=1$')
ax.loglog(x,pdf(2.,0.,x),c='c',lw=1.5,label=r'$c=2$')
ax.loglog(x,pdf(4.,0.,x),c='b',lw=1.5,label=r'$c=4$')
ax.loglog(x,pdf(8.,0.,x),c='g',lw=1.5,label=r'$c=8$')

ax.set_ylim(1e-4,1)
ax.set_xlim(0.05,15)
ax.legend(loc=4,frameon=False)
ax.set_xticks((0.1,1.0,10.))
ax.set_xticklabels((r'0.1',r'1.0',r'10.0'))
ax.set_xlabel(r'x $-\!\rightarrow$')

fig.savefig("Levy0_LdistributionPDF.svg",bbox_inches="tight",pad_inches=.15)