% ##################################################################
% # Function to calculate Euclidean distance from all locations specified
% # in vectors: lat & lon
% # to a single location specified with: x & y
% ##################################################################

function dist=dist_vec_ED(x,y,lat,lon,ii)
dist = zeros(length(lat),1);
xx=x.*pi./180;
yy=y.*pi./180;
latr=lat.*pi./180;
lonr=lon.*pi./180;
R = 6371;                           % Earth's radius in km
for i=1:ii
    delta_lat = latr(i) - yy;        % difference in latitude
    delta_lon = lonr(i) - xx;        % difference in longitude
    b = sqrt(delta_lat^2+delta_lon^2);
    dist(i)=R * b * 1000;  % distance in metres
end