docs: def i4_sobol_generate_std_normal(dim_num, n, skip=1):
      """
      Generates multivariate standard normal quasi-random variables.
  
      Parameters:
        Input, integer dim_num, the spatial dimension.
        Input, integer n, the number of points to generate.
        Input, integer SKIP, the number of initial points to skip.
  
        Output, real np array of shape (n, dim_num).
      """
  
      sobols = i4_sobol_generate(dim_num, n, skip)
  
      normals = norm.ppf(sobols)
  
      return normals