from spectrum import *
p = pburg(marple_data, 15, NFFT=4096)
p.plot(sides='centerdc')