Dist.m first

finish Dist.m. Now we can compute the called dist(P,C(P))