fig = plt.figure(figsize=(10, 7))
x = np.linspace(0, 1.75, 1000)
bins = plt.hist(stats[:, 0], bins=20, color='grey', alpha=0.5, ec='black')
plt.plot(x, gaus(x, *fit1)+gaus(x, *fit2), color='black', linestyle='dashed')
plt.xlabel('Frequencey of Max Amplitude', fontsize=20)

plt.annotate(r'$\mu_{{1}}={:0.2f}\pm{:0.2f}$ Hz'.format(fit1[0], err1[0]), xy=(0.5, 6), fontsize=15)
plt.annotate(r'$\sigma_{{1}}={:0.2f}\pm{:0.2f}$ Hz'.format(fit1[1], err1[1]), xy=(0.5, 5.5), fontsize=15)

plt.annotate(r'$\mu_{{2}}={:0.2f}\pm{:0.2f}$ Hz'.format(fit2[0], err2[0]), xy=(1, 6), fontsize=15)
plt.annotate(r'$\sigma_{{2}}={:0.2f}\pm{:0.2f}$ Hz'.format(fit2[1], err2[1]), xy=(1, 5.5), fontsize=15)

plt.savefig('Figures/FrequencyDist.pdf', bbox_inches='tight')