docs: def rand(n, **kwargs):
      """Random vector from the n-Sphere
  
      This function will return a random vector which is an element of the n-Sphere.
      The n-Sphere is defined as a vector in R^n+1 with a norm of one. 
  
      Basically, we'll find a random vector in R^n+1 and normalize it. 
      This uses the method of Marsaglia 1972.
  
      Parameters
      ----------
      None
  
      Returns
      -------
      rvec 
          Random (n+1,) numpy vector with a norm of 1
  
      """
      rvec = np.random.randn(3)
      rvec = rvec / np.linalg.norm(rvec)
      return rvec