#!/usr/bin/env python
 
"""This generates the graph depicting Re sinc(z)."""
# Author: Ika. 2015-08-13

from matplotlib.pyplot import *
import matplotlib.pyplot as plt
import numpy as np

low = -7.0
high = 7.0
tic = 0.1

X = np.arange(low, high, tic)

Z = np.outer(X,X)

for i in range(len(X)):
	for j in range(len(X)):
		x = low + tic*j
		y = low + tic*i
		if np.abs(x)<tic and np.abs(y)<tic:
			Z[i][j] = 0.0
		else:
			Z[i][j] = ((np.exp(1j*(x+1j*y))-np.exp(-1j*(x+1j*y)))/(2j*(x+1j*y))).imag

plt.figure(figsize=(4,3.2))
plt.axes([0.18,0.08,0.76,0.9])

plt.imshow(Z, origin='lower', \
			cmap=plt.cm.gist_rainbow,\
			extent=(low, high, low, high))
plt.colorbar(shrink=0.84)
for a in plt.gca().get_xticklabels():
    a.set_size(10)
for a in plt.gca().get_yticklabels():
    a.set_size(10)
plt.savefig("sinc_re.svg")