def METHOD_NAME(XYZ: ArrayLike) -> NDArrayFloat:
    """
    Convert from *whitepoint* *CIE XYZ* tristimulus values to
    *Hunter L,a,b* :math:`K_{a}` and :math:`K_{b}` chromaticity
    coefficients.
    Parameters
    ----------
    XYZ
        *Whitepoint* *CIE XYZ* tristimulus values.
    Returns
    -------
    :class:`numpy.ndarray`
        *Hunter L,a,b* :math:`K_{a}` and :math:`K_{b}` chromaticity
        coefficients.
    References
    ----------
    :cite:`HunterLab2008c`
    Examples
    --------
    >>> XYZ = np.array([109.850, 100.000, 35.585])
    >>> XYZ_to_K_ab_HunterLab1966(XYZ)  # doctest: +ELLIPSIS
    array([ 185.2378721...,   38.4219142...])
    """
    X, _Y, Z = tsplit(XYZ)
    K_a = 175 * np.sqrt(X / 98.043)
    K_b = 70 * np.sqrt(Z / 118.115)
    K_ab = tstack([K_a, K_b])
    return K_ab