function d=perpDist(x,y)

d=abs((x(2)-x(1))*(y(1)-y(3))-(x(1)-x(3))*(y(2)-y(1)))/sqrt((x(2)-x(1))^2+(y(2)-y(1))^2);