function d=distance2D(x1,x2,y1,y2)
    d=sqrt((x2-x1).^2+(y2-y1).^2)
end


% d=distance2D(0,3,0,4)