code: def angle(u, v, w=None):
      
      if w is not None:
          u = u - v
          v = w - v
      c = np.dot(u, v) / norm(u) / norm(v)
      return np.arccos(np.clip(c, -1, 1))
docs: Returns the angle in radians between two vectors.