%Program Dist.m
%
function d=Dist(xP,yP,xQ,yQ);
d=sqrt((xP-xQ)^2+(yP-yQ)^2);
end