Derivation of the relaxation time distribution function for the Havriliak-Negami model