from spectrum import *
p =  pminvar(data_cosine(), 15)
p()
p.plot()