dist = numpy.linalg.norm(a-b, ord=2) #ord=2 is default and means Euclidean distance, but I'm showing here that you can specify it