Divergence: one function for all of these three probability distributions

We should have one function distd() for the divergence (Rényi, Bhattacharyya or Hellinger) for MGGD, MCD or MTD. Similarly: one function kld() for the Kullback-Leibler divergence.