function eucDist = eucDist(x1,y1,x2,y2)

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

