code: def major_axis(points):
      
      U, S, V = np.linalg.svd(points)
      axis = util.unitize(np.dot(S, V))
      return axis
docs: Returns an approximate vector representing the major axis of points
  
      Parameters
      -------------
      points: (n, dimension) float, points in space
  
      Returns
      -------------
      axis: (dimension,) float, vector along approximate major axis