In [11]:
import numpy as np
%matplotlib nbagg
import matplotlib.pyplot as plt
# the plot says a an effective kM = 0.8 is reasonable
k = np.arange(0, 2.**10)
kM = k[-1]
fig = plt.figure(figsize=(4,3))
a = fig.add_subplot(111)
a.plot(k/kM, np.exp(-36*(k / kM)**32))
fig.tight_layout()
In [ ]: