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

x1=np.arange(-21,0.9999,0.001);
x2=np.arange(1.0001,10,0.001);
xvert=np.arange(-33,33,0.1);

def g(x):
        return 1;
h=np.vectorize(g)
yvert=h(xvert)

def k(x):
        return 0;
l=np.vectorize(k)
zeroline=l(xvert)

f=np.vectorize(lambda x: zeta(x))
y1=f(x1)
y2=f(x2)

fig, ax = plt.subplots()

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

plt.xlim(-19,10)
plt.ylim(-4,4)

plt.plot(x1,y1,linewidth=2,color='Black')
plt.plot(x2,y2,linewidth=2,color='Black')
plt.plot(yvert,xvert,'--',linewidth=2,color='Black')
plt.plot(xvert,yvert,'--',linewidth=2,color='Black')
plt.plot(xvert,zeroline,'--',linewidth=0.5,color='Black')

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