distance_matrix function

Distance matrix calculation