function d = pointdist(x1,y1,x2,y2)

d = sqrt((x1-x2).^2 + (y1-y2).^2);