code: def get_unf(N=100):
      
  #
  # get uniform directions  [dec,inc]
      z = random.uniform(-1., 1., size=N)
      t = random.uniform(0., 360., size=N)  # decs
      i = np.arcsin(z) * 180. / np.pi  # incs
      return np.array([t, i]).transpose()
  
  # def get_unf(N): #Jeff's way
      """
       subroutine to retrieve N uniformly distributed directions
      """
docs: Generates N uniformly distributed directions
      using the way described in Fisher et al. (1987).
      Parameters
      __________
      N : number of directions, default is 100
  
      Returns
      ______
      array of nested dec,inc pairs