x_min, x_max, y_min, y_max = -2 * np.pi, 2 * np.pi, -2 * np.pi, 2 * np.pi
f = [[f1(x, y) for x in np.linspace(x_min, x_max, num=500)]
for y in np.linspace(y_min, y_max, num=500)]
plt.figure(figsize=(10, 7))
plt.imshow(f,
interpolation="bicubic",
origin="lower",
extent=[x_min, x_max, y_min, y_max])
plt.colorbar()
plt.title(r"Title here (remove for papers)")
plt.xlabel(r"Description of $x$ coordinate (units)")
plt.ylabel(r"Description of $y$ coordinate (units)")
plt.tight_layout()
plt.show()