import matplotlib.pyplot as plt
import numpy as np
from crystals import Crystal
graphite = Crystal.from_database('C')
from skued import powdersim
q = np.linspace(1, 10, 1024)
diff = powdersim(graphite, q)
plt.figure()
plt.plot(q, diff/diff.max())
plt.xlim([q.min(), q.max()])
plt.xlabel('$q (1/\AA)$')
plt.ylabel('Diffracted intensity (A.u.)')
plt.title('Polycrystalline graphite diffraction')