File:Besseljn=0plot.png

#!/usr/bin/python
import numpy as np
import matplotlib.pyplot as plt
from mpmath import *
from pylab import rcParams
rcParams['figure.figsize'] = 10, 10

x=np.arange(-20,20,0.01);
def f(x):
        return besselj(0,x)
g=np.vectorize(f)
y=g(x)
def h(x):
        return 0
k=np.vectorize(h)
y2=k(x)

fig, ax = plt.subplots()

xlabels = [item.get_text() for item in ax.get_xticklabels()]
xlabels=[r'$-20$',r'$-15$',r'$-10$',r'$-5$',r'$0$',r'$j_{0,1}$',r'$j_{0,2}$',r'$j_{0,3}$',r'$j_{0,4}$',r'$j_{0,5}$',r'$j_{0,6}$',r'$20$']
ax.set_xticks([-20,-15,-10,-5,0,2.40483,5.52008,8.65373,11.7915,14.9309,18.0711,20])
ax.set_xticklabels(xlabels,horizontalalignment='center')

plt.xlabel(r'$x$')
plt.ylabel(r'$\mathrm{J}_0(x)$')

plt.plot(x,y,linewidth=2,color='Black')
plt.plot(x,y2,'--',linewidth=2,color='Black')

plt.savefig('besseljn=0plot.png',bbox_inches='tight',pad_inches=0.15)