Create HR.py

Phase plot and time series, no power spectra or coupled time series yet.