import numpy as np
import matplotlib.pyplot as plt

b = 2.5
x = np.linspace(0, 20, 200)

f = np.exp(-(x - b) ** 2 / 2)
fig, ax = plt.subplots(1, 1, figsize=(9, 3))
ax.plot(x, f, "k")
ax.set_ylim([-0.5, 1.5])
ax.set_xticks(np.linspace(0, 20, 21))
ax.fill_between(x, f, np.zeros_like(f), color="k")
ax.set_title("Seismic Wavelet", fontsize=13)
ax.set_xlabel("Time (ms)")

# fig.savefig('SeisWavelet.png',dpi=300)
