import matplotlib.pyplot as plt
import numpy as np

s, intensity = np.load('data/powder.npy')
fig = plt.figure()
ax = fig.add_subplot(111)
ax.plot(s, intensity, 'k')
ax.set_xlabel('s ($4 \pi / \AA$)')
ax.set_ylabel('Diffracted intensity (counts)')
ax.set_title('Background-subtracted diffraction pattern of rutile VO$_2$')
plt.show()