code: def vecangle(v1, v2, deg=True):
      
      if np.array_equal(v1, v2):
          return 0.0
      dm = np.dot(v1, v2)
      cm = np.linalg.norm(v1) * np.linalg.norm(v2)
      angle = np.arccos(round(dm / cm, 10))  # Round here to prevent floating point errors
      return np.degrees([angle, ])[0] if deg else angle
docs: Calculate the angle between two vectors
      :param v1: coordinates of vector v1
      :param v2: coordinates of vector v2
      :returns : angle in degree or rad