File:Gammaplot.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

x0=np.arange(0.03,5.2,0.001); #right piece
x1=np.arange(-0.97,-0.03,0.001); #left 1 piece
x2=np.arange(-1.99,-1.001,0.001); #left 2 piece
x3=np.arange(-2.99999,-2.001,0.001); #left 3 piece
x4=np.arange(-3.9999,-3.001,0.001); #left 4 piece
x5=np.arange(-4.999999,-4.001,0.001); #left 5 piece
x6=np.arange(-33,33,0.1); #vertical

f=np.vectorize(gamma)
y0=f(x0)
y1=f(x1)
y2=f(x2)
y3=f(x3)
y4=f(x4)
y5=f(x5)
def g(x):
        return 0
def h(x):
        return -1
def k(x):
        return -2
def l(x):
        return -3
def m(x):
        return -4
def n(x):
        return -5
o=np.vectorize(g)
p=np.vectorize(h)
q=np.vectorize(k)
r=np.vectorize(l)
s=np.vectorize(m)
t=np.vectorize(n)

y6=o(x6)
y7=p(x6)
y8=q(x6)
y9=r(x6)
y10=s(x6)
y11=t(x6)

fig, ax = plt.subplots()

xlabels = [item.get_text() for item in ax.get_xticklabels()]
xlabels=[r'$-5$',r'$-4$',r'$-3$',r'$-2$',r'$-1$',r'$0$',r'$1$',r'$2$',r'$3$',r'$4$',r'$5$']
ax.set_xticks([-5,-4,-3,-2,-1,0,1,2,3,4,5])
ax.set_xticklabels(xlabels,horizontalalignment='center')

plt.xlabel(r'$x$')
plt.ylabel(r'$\Gamma(x)$')

plt.plot(x0,y0,linewidth=2,color='Black')
plt.plot(x1,y1,linewidth=2,color='Black')
plt.plot(x2,y2,linewidth=2,color='Black')
plt.plot(x3,y3,linewidth=2,color='Black')
plt.plot(x4,y4,linewidth=2,color='Black')
plt.plot(x5,y5,linewidth=2,color='Black')
plt.plot(y6,x6,'--',linewidth=2,color='Black')
plt.plot(y7,x6,'--',linewidth=2,color='Black')
plt.plot(y8,x6,'--',linewidth=2,color='Black')
plt.plot(y9,x6,'--',linewidth=2,color='Black')
plt.plot(y10,x6,'--',linewidth=2,color='Black')
plt.plot(y11,x6,'--',linewidth=2,color='Black')


plt.xlim([-5.3,5.3])
plt.ylim([-20,30])

plt.savefig('gammaplot.png',bbox_inches='tight',pad_inches=0.15)