code: def angular_errors(hyp_axes):
      
      # Not quite sure why this is sqrt but it is empirically correct
      ax = N.sqrt(hyp_axes)
      return tuple(N.arctan2(ax[-1],ax[:-1]))
docs: Minimum and maximum angular errors
      corresponding to 1st and 2nd axes
      of PCA distribution.
  
      Ordered as [minimum, maximum] angular error.