Difference between revisions of "File:Sphericalbessely,n=0glyph.png"

Latest revision as of 23:14, 10 June 2016

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

def sphericalbessely(n,z):
        return ((pi/(2*z))**(1/2))*bessely(n+1/2,z)

x1=np.arange(0.01,20,0.01)
x2=np.arange(-25,25,0.1)

f=np.vectorize(lambda x: sphericalbessely(0,x))
y1=f(x1)

def g(x):
        return 0
h=np.vectorize(g)
y2=h(x2)
y3=f(x3)

fig, ax = plt.subplots()

ax.spines['top'].set_color('none')
ax.spines['left'].set_color('none')
ax.spines['right'].set_color('none')
ax.spines['bottom'].set_color('none')

plt.tick_params(axis='x', which='both', bottom='off', top='off', labelbottom='off')
plt.tick_params(axis='y', which='both', left='off', right='off', labelleft='off')

plt.plot(x1,y1,linewidth=5,color='Black')
plt.plot(x2,y2,linewidth=1,color='Black')
plt.plot(y2,x2,linewidth=1,color='Black')

plt.ylim([-0.6,0.6])
plt.xlim([0,20.3])

plt.savefig('sphericalbessely,n=0glyph.png',bbox_inches='tight',pad_inches=0.0)