get_hrf function

Generate a double gamma model of the HRF as used in SPM.