Time-domain simulation for fractional relaxation of Havriliak-Negami type