cosmo pdist hdr

function d = cosmo_pdist(x, distance)