function dist = pd(d1, d2)
dist = sqrt(sum((d1-d2).^2));