import matplotlib.pyplot as plt
import numpy as np

from mpl_toolkits.axes_grid1.inset_locator import zoomed_inset_axes, mark_inset

blue="#007ec1"

xs1 = np.linspace(-3, 3, 10001)
xs2=1/(np.arange(1.,1001,1.)*np.pi/2)

xs = np.union1d(np.union1d(xs1,xs2),-xs2)
ys = np.sin(1/xs)

f, ax = plt.subplots()

ax.plot(xs,ys, color=blue)
ax.set_xlim(-2.5,2.5)
ax.set_ylim(-2,2)

axins = zoomed_inset_axes(ax,1.7,loc=4)
axins.plot(xs,ys)
axins.set_xlim(-0.3,0.3)
axins.set_ylim(-1.1,1.1)

mark_inset(ax, axins, loc1=3, loc2=2, fc="0.95", ec="0.7")

plt.xticks(visible=False)
plt.yticks(visible=False)

plt.savefig("The function sin(1 over x).svg")