code: def anglesep(lon0: float, lat0: float,
               lon1: float, lat1: float, deg: bool = True) -> float:
      
      if angular_separation is None:
          raise ImportError('angledist requires AstroPy. Try angledis_meeus')
  
      if deg:
          lon0 = radians(lon0)
          lat0 = radians(lat0)
          lon1 = radians(lon1)
          lat1 = radians(lat1)
  
      sep_rad = angular_separation(lon0, lat0, lon1, lat1)
  
      if deg:
          return degrees(sep_rad)
      else:
          return sep_rad
docs: Parameters
      ----------
  
      lon0 : float
          longitude of first point
      lat0 : float
          latitude of first point
      lon1 : float
          longitude of second point
      lat1 : float
          latitude of second point
      deg : bool, optional
            degrees input/output  (False: radians in/out)
  
      Returns
      -------
  
      sep_rad : float or numpy.ndarray of float
          angular separation
  
      For reference, this is from astropy astropy/coordinates/angle_utilities.py
      Angular separation between two points on a sphere.