from spectrum import pev, marple_data
p = pev(marple_data, 15, NFFT=4096)
p.plot()