def pathl(cop_dat): # to calculate COP path length delt = np.diff(cop_dat[:,(0,1)], axis = 0) sqs = np.square(delt) sum_s = np.add(sqs[:,0],sqs[:,1]) lgths = np.sqrt(sum_s) return np.sum(lgths)