function nei = RD_neighbors(locs)
%

[num_pts,dummy] = size(locs);
nei = zeros(num_pts,num_pts);
for p1 = 1:num_pts-1
    pt1 = locs(p1,:);
    for p2 = p1+1:num_pts
        pt2 = locs(p2,:);
        nei(p1,p2) = norm(pt1-pt2);
        nei(p2,p1) = nei(p1,p2);
    end
end
