spectrum.py: add a simple script for plotting power spectra