code: def mnl_utility_transform(sys_utility_array, *args, **kwargs):
      
      # Return a 2D array of systematic utility values
      if len(sys_utility_array.shape) == 1:
          systematic_utilities = sys_utility_array[:, np.newaxis]
      else:
          systematic_utilities = sys_utility_array
  
      return systematic_utilities
docs: Parameters
      ----------
      sys_utility_array : ndarray.
          Should have 1D or 2D. Should have been created by the dot product of a
          design matrix and an array of index coefficients.
  
      Returns
      -------
          systematic_utilities : 2D ndarray.
              The input systematic utilities. If `sys_utility_array` is 2D, then
              `sys_utility_array` is returned. Else, returns
              `sys_utility_array[:, None]`.