A Grunwald–Letnikov scheme for fractional operators of Havriliak–Negami type