Add spectrum() function to obsutils and create make_spectrum.py example script to demonstrate its usage.